Re: [Paraview] Cataylst pipeline not updating

2017-02-14 Thread Andy Bauer
Hi Tim,

I'm glad this fixed the issue. You'll want to always pass in the time to
UpdatePipeline() for Catalyst. The issue was that the pipeline thought that
everything it had was current since the time didn't change and none if the
filters' parameters changed and therefore it skipped recomputing what it
thought would be the same result. I do wonder if the start of the pipeline
should mark that it is modified. I'll have to think about that some more.
Even if that were true though you'd still want to pass in the time value to
UpdatePipeline(). This is true as well for ParaView Python pipelines in
general for time dependent sources.

No problem about the delay in reporting -- I was more than happy to let
your computers do the proof rather than verifying it myself :)

Cheers,
Andy

On Tue, Feb 14, 2017 at 11:00 AM, Gallagher, Timothy P <
tim.gallag...@gatech.edu> wrote:

> Hi Andy,
>
>
> Looks like passing the current time step into the UpdatePipeline calls
> made it work as expected. Sorry for the delay in reporting, it took a long
> time for the job to move through the queue!
>
>
> Anyway, looks like it's working properly and I am getting the data as
> expected.
>
>
> For future reference, should I always just pass the time in? Or are the
> possibly other, unintended side effects from that as well to watch out for?
>
>
> Thanks again,
>
>
> Tim
>
>
> --
> *From:* Andy Bauer <andy.ba...@kitware.com>
> *Sent:* Monday, February 13, 2017 11:30 AM
>
> *To:* Gallagher, Timothy P
> *Cc:* paraview@paraview.org
> *Subject:* Re: [Paraview] Cataylst pipeline not updating
>
> Hi Tim,
>
> It's possible that there's something funky going on when the data set
> being passed through parts of the pipeline are empty. I'm hoping for the
> easy fix with passing in the time since I've know that has caused
> unexpected behavior before. If that doesn't solve your issue then I'll
> probably have to do something similar and run it locally to see what's
> going on.
>
> Cheers,
> Andy
>
> On Mon, Feb 13, 2017 at 11:06 AM, Gallagher, Timothy P <
> tim.gallag...@gatech.edu> wrote:
>
>> The RequestDataDescription is the same as the default output from the
>> script generator, so I didn't include it. The DoCoProcessing gets called
>> every time step and the volumeIntegrals data is updated and written
>> correctly every delta_t, so I'm confident in that part.
>>
>>
>> I will add the calls you suggested and re-run. It takes some time for it
>> to evolve to the point something will happen -- I don't have it saving data
>> frequently enough to run a quick test, I was expecting the pipeline to save
>> it for me!
>>
>>
>> The other thing that occurred to me -- is there something about the
>> Contour filter that turns it off/disables it when the pipeline is built and
>> there is no surface at that contour level?
>>
>>
>> Thanks,
>>
>>
>> Tim
>>
>>
>> --
>> *From:* Andy Bauer <andy.ba...@kitware.com>
>> *Sent:* Monday, February 13, 2017 10:56 AM
>> *To:* Gallagher, Timothy P
>> *Cc:* paraview@paraview.org
>> *Subject:* Re: [Paraview] Cataylst pipeline not updating
>>
>> Hi Tim,
>>
>> Try adding in the current time to the UpdatePipeline() calls, e.g.
>> *.UpdatePipeline(datadescription.GetCurrentTime())
>>
>> It's tough to tell with RequestDataDescription() missing but I assume
>> that DoCoProcessing() is getting called when you want, correct?
>>
>> Let me know if that doesn't solve the problem.
>>
>> Best,
>> Andy
>>
>> On Mon, Feb 13, 2017 at 5:45 AM, Gallagher, Timothy P <
>> tim.gallag...@gatech.edu> wrote:
>>
>>> Back with another question!
>>>
>>>
>>> I have a pipeline that has 1 input and 4 outputs. The simplest is a
>>> IntegrateVariables filter hooked up to the input. Another branch of the
>>> pipeline goes through several PythonCalculator filters and a
>>> CellDataToPointData filter and then arrives at the three outputs:
>>>
>>>
>>> 1. The CD2PD filter is connected to a Contour filter with a single
>>> level, X. This is then connected to an IntegrateVariables filter (used to
>>> compute the area of the Contour) and a Writer.
>>>
>>> 2. The CD2PD filter is also connected to another Contour filter with a
>>> single level Y. This is connected to an IntegrateVariables and Writer just
>>> as in 1.
>>>
>>> 3. The CD2PD filter is sliced two times with each slice in a different
>>> ren

Re: [Paraview] Cataylst pipeline not updating

2017-02-14 Thread Gallagher, Timothy P
Hi Andy,


Looks like passing the current time step into the UpdatePipeline calls made it 
work as expected. Sorry for the delay in reporting, it took a long time for the 
job to move through the queue!


Anyway, looks like it's working properly and I am getting the data as expected.


For future reference, should I always just pass the time in? Or are the 
possibly other, unintended side effects from that as well to watch out for?


Thanks again,


Tim



From: Andy Bauer <andy.ba...@kitware.com>
Sent: Monday, February 13, 2017 11:30 AM
To: Gallagher, Timothy P
Cc: paraview@paraview.org
Subject: Re: [Paraview] Cataylst pipeline not updating

Hi Tim,

It's possible that there's something funky going on when the data set being 
passed through parts of the pipeline are empty. I'm hoping for the easy fix 
with passing in the time since I've know that has caused unexpected behavior 
before. If that doesn't solve your issue then I'll probably have to do 
something similar and run it locally to see what's going on.

Cheers,
Andy

On Mon, Feb 13, 2017 at 11:06 AM, Gallagher, Timothy P 
<tim.gallag...@gatech.edu<mailto:tim.gallag...@gatech.edu>> wrote:

The RequestDataDescription is the same as the default output from the script 
generator, so I didn't include it. The DoCoProcessing gets called every time 
step and the volumeIntegrals data is updated and written correctly every 
delta_t, so I'm confident in that part.


I will add the calls you suggested and re-run. It takes some time for it to 
evolve to the point something will happen -- I don't have it saving data 
frequently enough to run a quick test, I was expecting the pipeline to save it 
for me!


The other thing that occurred to me -- is there something about the Contour 
filter that turns it off/disables it when the pipeline is built and there is no 
surface at that contour level?


Thanks,


Tim



From: Andy Bauer <andy.ba...@kitware.com<mailto:andy.ba...@kitware.com>>
Sent: Monday, February 13, 2017 10:56 AM
To: Gallagher, Timothy P
Cc: paraview@paraview.org<mailto:paraview@paraview.org>
Subject: Re: [Paraview] Cataylst pipeline not updating

Hi Tim,

Try adding in the current time to the UpdatePipeline() calls, e.g. 
*.UpdatePipeline(datadescription.GetCurrentTime())

It's tough to tell with RequestDataDescription() missing but I assume that 
DoCoProcessing() is getting called when you want, correct?

Let me know if that doesn't solve the problem.

Best,
Andy

On Mon, Feb 13, 2017 at 5:45 AM, Gallagher, Timothy P 
<tim.gallag...@gatech.edu<mailto:tim.gallag...@gatech.edu>> wrote:

Back with another question!


I have a pipeline that has 1 input and 4 outputs. The simplest is a 
IntegrateVariables filter hooked up to the input. Another branch of the 
pipeline goes through several PythonCalculator filters and a 
CellDataToPointData filter and then arrives at the three outputs:


1. The CD2PD filter is connected to a Contour filter with a single level, X. 
This is then connected to an IntegrateVariables filter (used to compute the 
area of the Contour) and a Writer.

2. The CD2PD filter is also connected to another Contour filter with a single 
level Y. This is connected to an IntegrateVariables and Writer just as in 1.

3. The CD2PD filter is sliced two times with each slice in a different render 
view to save images.


I have the coprocessor running every delta_t seconds. The initial 
IntegrateVariables filter is updated with UpdatePipeline and then the values 
are saved to disk. This works every step, no issues.


I then call UpdatePipeline on both of the Contour and both of the 
IntegrateVariables filters in 1 and 2. I check if the 'Area' field in the 
IntegrateVariables exceeds a threshold, and if it does, I save the Contour's to 
disk as a VTK file. This is not working correctly.


If I start the simulation from a point where the Contour's exist in the flow 
already, then everything works great and the Area field changes in time and it 
saves the Contour's.


However, if the simulation starts without the Contours already in the flow, the 
Area fields stay zero forever. It's as if the pipeline isn't updating or 
something. I have attached the DoCoProcessing function used. I don't see 
anything immediately obvious that would prevent it from updating or getting the 
correct values for the Area.


Am I missing something about how the pipeline or contour filters update?


Thanks,


Tim

___
Powered by www.kitware.com<http://www.kitware.com>

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscrib

Re: [Paraview] Cataylst pipeline not updating

2017-02-13 Thread Andy Bauer
Hi Tim,

It's possible that there's something funky going on when the data set being
passed through parts of the pipeline are empty. I'm hoping for the easy fix
with passing in the time since I've know that has caused unexpected
behavior before. If that doesn't solve your issue then I'll probably have
to do something similar and run it locally to see what's going on.

Cheers,
Andy

On Mon, Feb 13, 2017 at 11:06 AM, Gallagher, Timothy P <
tim.gallag...@gatech.edu> wrote:

> The RequestDataDescription is the same as the default output from the
> script generator, so I didn't include it. The DoCoProcessing gets called
> every time step and the volumeIntegrals data is updated and written
> correctly every delta_t, so I'm confident in that part.
>
>
> I will add the calls you suggested and re-run. It takes some time for it
> to evolve to the point something will happen -- I don't have it saving data
> frequently enough to run a quick test, I was expecting the pipeline to save
> it for me!
>
>
> The other thing that occurred to me -- is there something about the
> Contour filter that turns it off/disables it when the pipeline is built and
> there is no surface at that contour level?
>
>
> Thanks,
>
>
> Tim
>
>
> --
> *From:* Andy Bauer <andy.ba...@kitware.com>
> *Sent:* Monday, February 13, 2017 10:56 AM
> *To:* Gallagher, Timothy P
> *Cc:* paraview@paraview.org
> *Subject:* Re: [Paraview] Cataylst pipeline not updating
>
> Hi Tim,
>
> Try adding in the current time to the UpdatePipeline() calls, e.g.
> *.UpdatePipeline(datadescription.GetCurrentTime())
>
> It's tough to tell with RequestDataDescription() missing but I assume that
> DoCoProcessing() is getting called when you want, correct?
>
> Let me know if that doesn't solve the problem.
>
> Best,
> Andy
>
> On Mon, Feb 13, 2017 at 5:45 AM, Gallagher, Timothy P <
> tim.gallag...@gatech.edu> wrote:
>
>> Back with another question!
>>
>>
>> I have a pipeline that has 1 input and 4 outputs. The simplest is a
>> IntegrateVariables filter hooked up to the input. Another branch of the
>> pipeline goes through several PythonCalculator filters and a
>> CellDataToPointData filter and then arrives at the three outputs:
>>
>>
>> 1. The CD2PD filter is connected to a Contour filter with a single level,
>> X. This is then connected to an IntegrateVariables filter (used to compute
>> the area of the Contour) and a Writer.
>>
>> 2. The CD2PD filter is also connected to another Contour filter with a
>> single level Y. This is connected to an IntegrateVariables and Writer just
>> as in 1.
>>
>> 3. The CD2PD filter is sliced two times with each slice in a different
>> render view to save images.
>>
>>
>> I have the coprocessor running every delta_t seconds. The initial
>> IntegrateVariables filter is updated with UpdatePipeline and then the
>> values are saved to disk. This works every step, no issues.
>>
>>
>> I then call UpdatePipeline on both of the Contour and both of the
>> IntegrateVariables filters in 1 and 2. I check if the 'Area' field in the
>> IntegrateVariables exceeds a threshold, and if it does, I save the
>> Contour's to disk as a VTK file. This is not working correctly.
>>
>>
>> If I start the simulation from a point where the Contour's exist in the
>> flow already, then everything works great and the Area field changes in
>> time and it saves the Contour's.
>>
>>
>> However, if the simulation starts without the Contours already in the
>> flow, the Area fields stay zero forever. It's as if the pipeline isn't
>> updating or something. I have attached the DoCoProcessing function used. I
>> don't see anything immediately obvious that would prevent it from updating
>> or getting the correct values for the Area.
>>
>>
>> Am I missing something about how the pipeline or contour filters update?
>>
>>
>> Thanks,
>>
>>
>> Tim
>>
>> ___
>> Powered by www.kitware.com
>>
>> Visit other Kitware open-source projects at
>> http://www.kitware.com/opensource/opensource.html
>>
>> Please keep messages on-topic and check the ParaView Wiki at:
>> http://paraview.org/Wiki/ParaView
>>
>> Search the list archives at: http://markmail.org/search/?q=ParaView
>>
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/mailman/listinfo/paraview
>>
>>
>
___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Cataylst pipeline not updating

2017-02-13 Thread Gallagher, Timothy P
The RequestDataDescription is the same as the default output from the script 
generator, so I didn't include it. The DoCoProcessing gets called every time 
step and the volumeIntegrals data is updated and written correctly every 
delta_t, so I'm confident in that part.


I will add the calls you suggested and re-run. It takes some time for it to 
evolve to the point something will happen -- I don't have it saving data 
frequently enough to run a quick test, I was expecting the pipeline to save it 
for me!


The other thing that occurred to me -- is there something about the Contour 
filter that turns it off/disables it when the pipeline is built and there is no 
surface at that contour level?


Thanks,


Tim



From: Andy Bauer <andy.ba...@kitware.com>
Sent: Monday, February 13, 2017 10:56 AM
To: Gallagher, Timothy P
Cc: paraview@paraview.org
Subject: Re: [Paraview] Cataylst pipeline not updating

Hi Tim,

Try adding in the current time to the UpdatePipeline() calls, e.g. 
*.UpdatePipeline(datadescription.GetCurrentTime())

It's tough to tell with RequestDataDescription() missing but I assume that 
DoCoProcessing() is getting called when you want, correct?

Let me know if that doesn't solve the problem.

Best,
Andy

On Mon, Feb 13, 2017 at 5:45 AM, Gallagher, Timothy P 
<tim.gallag...@gatech.edu<mailto:tim.gallag...@gatech.edu>> wrote:

Back with another question!


I have a pipeline that has 1 input and 4 outputs. The simplest is a 
IntegrateVariables filter hooked up to the input. Another branch of the 
pipeline goes through several PythonCalculator filters and a 
CellDataToPointData filter and then arrives at the three outputs:


1. The CD2PD filter is connected to a Contour filter with a single level, X. 
This is then connected to an IntegrateVariables filter (used to compute the 
area of the Contour) and a Writer.

2. The CD2PD filter is also connected to another Contour filter with a single 
level Y. This is connected to an IntegrateVariables and Writer just as in 1.

3. The CD2PD filter is sliced two times with each slice in a different render 
view to save images.


I have the coprocessor running every delta_t seconds. The initial 
IntegrateVariables filter is updated with UpdatePipeline and then the values 
are saved to disk. This works every step, no issues.


I then call UpdatePipeline on both of the Contour and both of the 
IntegrateVariables filters in 1 and 2. I check if the 'Area' field in the 
IntegrateVariables exceeds a threshold, and if it does, I save the Contour's to 
disk as a VTK file. This is not working correctly.


If I start the simulation from a point where the Contour's exist in the flow 
already, then everything works great and the Area field changes in time and it 
saves the Contour's.


However, if the simulation starts without the Contours already in the flow, the 
Area fields stay zero forever. It's as if the pipeline isn't updating or 
something. I have attached the DoCoProcessing function used. I don't see 
anything immediately obvious that would prevent it from updating or getting the 
correct values for the Area.


Am I missing something about how the pipeline or contour filters update?


Thanks,


Tim

___
Powered by www.kitware.com<http://www.kitware.com>

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Cataylst pipeline not updating

2017-02-13 Thread Andy Bauer
Oops, datadescription.GetTime() instead of datadescription.GetCurrentTime().

On Mon, Feb 13, 2017 at 10:56 AM, Andy Bauer  wrote:

> Hi Tim,
>
> Try adding in the current time to the UpdatePipeline() calls, e.g.
> *.UpdatePipeline(datadescription.GetCurrentTime())
>
> It's tough to tell with RequestDataDescription() missing but I assume that
> DoCoProcessing() is getting called when you want, correct?
>
> Let me know if that doesn't solve the problem.
>
> Best,
> Andy
>
> On Mon, Feb 13, 2017 at 5:45 AM, Gallagher, Timothy P <
> tim.gallag...@gatech.edu> wrote:
>
>> Back with another question!
>>
>>
>> I have a pipeline that has 1 input and 4 outputs. The simplest is a
>> IntegrateVariables filter hooked up to the input. Another branch of the
>> pipeline goes through several PythonCalculator filters and a
>> CellDataToPointData filter and then arrives at the three outputs:
>>
>>
>> 1. The CD2PD filter is connected to a Contour filter with a single level,
>> X. This is then connected to an IntegrateVariables filter (used to compute
>> the area of the Contour) and a Writer.
>>
>> 2. The CD2PD filter is also connected to another Contour filter with a
>> single level Y. This is connected to an IntegrateVariables and Writer just
>> as in 1.
>>
>> 3. The CD2PD filter is sliced two times with each slice in a different
>> render view to save images.
>>
>>
>> I have the coprocessor running every delta_t seconds. The initial
>> IntegrateVariables filter is updated with UpdatePipeline and then the
>> values are saved to disk. This works every step, no issues.
>>
>>
>> I then call UpdatePipeline on both of the Contour and both of the
>> IntegrateVariables filters in 1 and 2. I check if the 'Area' field in the
>> IntegrateVariables exceeds a threshold, and if it does, I save the
>> Contour's to disk as a VTK file. This is not working correctly.
>>
>>
>> If I start the simulation from a point where the Contour's exist in the
>> flow already, then everything works great and the Area field changes in
>> time and it saves the Contour's.
>>
>>
>> However, if the simulation starts without the Contours already in the
>> flow, the Area fields stay zero forever. It's as if the pipeline isn't
>> updating or something. I have attached the DoCoProcessing function used. I
>> don't see anything immediately obvious that would prevent it from updating
>> or getting the correct values for the Area.
>>
>>
>> Am I missing something about how the pipeline or contour filters update?
>>
>>
>> Thanks,
>>
>>
>> Tim
>>
>> ___
>> Powered by www.kitware.com
>>
>> Visit other Kitware open-source projects at
>> http://www.kitware.com/opensource/opensource.html
>>
>> Please keep messages on-topic and check the ParaView Wiki at:
>> http://paraview.org/Wiki/ParaView
>>
>> Search the list archives at: http://markmail.org/search/?q=ParaView
>>
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/mailman/listinfo/paraview
>>
>>
>
___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Cataylst pipeline not updating

2017-02-13 Thread Andy Bauer
Hi Tim,

Try adding in the current time to the UpdatePipeline() calls, e.g.
*.UpdatePipeline(datadescription.GetCurrentTime())

It's tough to tell with RequestDataDescription() missing but I assume that
DoCoProcessing() is getting called when you want, correct?

Let me know if that doesn't solve the problem.

Best,
Andy

On Mon, Feb 13, 2017 at 5:45 AM, Gallagher, Timothy P <
tim.gallag...@gatech.edu> wrote:

> Back with another question!
>
>
> I have a pipeline that has 1 input and 4 outputs. The simplest is a
> IntegrateVariables filter hooked up to the input. Another branch of the
> pipeline goes through several PythonCalculator filters and a
> CellDataToPointData filter and then arrives at the three outputs:
>
>
> 1. The CD2PD filter is connected to a Contour filter with a single level,
> X. This is then connected to an IntegrateVariables filter (used to compute
> the area of the Contour) and a Writer.
>
> 2. The CD2PD filter is also connected to another Contour filter with a
> single level Y. This is connected to an IntegrateVariables and Writer just
> as in 1.
>
> 3. The CD2PD filter is sliced two times with each slice in a different
> render view to save images.
>
>
> I have the coprocessor running every delta_t seconds. The initial
> IntegrateVariables filter is updated with UpdatePipeline and then the
> values are saved to disk. This works every step, no issues.
>
>
> I then call UpdatePipeline on both of the Contour and both of the
> IntegrateVariables filters in 1 and 2. I check if the 'Area' field in the
> IntegrateVariables exceeds a threshold, and if it does, I save the
> Contour's to disk as a VTK file. This is not working correctly.
>
>
> If I start the simulation from a point where the Contour's exist in the
> flow already, then everything works great and the Area field changes in
> time and it saves the Contour's.
>
>
> However, if the simulation starts without the Contours already in the
> flow, the Area fields stay zero forever. It's as if the pipeline isn't
> updating or something. I have attached the DoCoProcessing function used. I
> don't see anything immediately obvious that would prevent it from updating
> or getting the correct values for the Area.
>
>
> Am I missing something about how the pipeline or contour filters update?
>
>
> Thanks,
>
>
> Tim
>
> ___
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at http://www.kitware.com/
> opensource/opensource.html
>
> Please keep messages on-topic and check the ParaView Wiki at:
> http://paraview.org/Wiki/ParaView
>
> Search the list archives at: http://markmail.org/search/?q=ParaView
>
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/mailman/listinfo/paraview
>
>
___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


[Paraview] Cataylst pipeline not updating

2017-02-13 Thread Gallagher, Timothy P
Back with another question!


I have a pipeline that has 1 input and 4 outputs. The simplest is a 
IntegrateVariables filter hooked up to the input. Another branch of the 
pipeline goes through several PythonCalculator filters and a 
CellDataToPointData filter and then arrives at the three outputs:


1. The CD2PD filter is connected to a Contour filter with a single level, X. 
This is then connected to an IntegrateVariables filter (used to compute the 
area of the Contour) and a Writer.

2. The CD2PD filter is also connected to another Contour filter with a single 
level Y. This is connected to an IntegrateVariables and Writer just as in 1.

3. The CD2PD filter is sliced two times with each slice in a different render 
view to save images.


I have the coprocessor running every delta_t seconds. The initial 
IntegrateVariables filter is updated with UpdatePipeline and then the values 
are saved to disk. This works every step, no issues.


I then call UpdatePipeline on both of the Contour and both of the 
IntegrateVariables filters in 1 and 2. I check if the 'Area' field in the 
IntegrateVariables exceeds a threshold, and if it does, I save the Contour's to 
disk as a VTK file. This is not working correctly.


If I start the simulation from a point where the Contour's exist in the flow 
already, then everything works great and the Area field changes in time and it 
saves the Contour's.


However, if the simulation starts without the Contours already in the flow, the 
Area fields stay zero forever. It's as if the pipeline isn't updating or 
something. I have attached the DoCoProcessing function used. I don't see 
anything immediately obvious that would prevent it from updating or getting the 
correct values for the Area.


Am I missing something about how the pipeline or contour filters update?


Thanks,


Tim

#  Processing method 
def WriteDataFile(comm, object, filename, curTime, CellData=True):
  object.UpdatePipeline()

  if CellData:
data = object.CellData
  else:
data = object.PointData

  if comm.Get_rank() == 0:
with  open(filename, 'a+') as fid:
  # Rewind the file to see if it is empty or not
  fid.seek(0)
  if len(fid.read(1)) == 0:
# The file is empty, we need to create the header
outputStr = '# Time '
for var in data.keys():
  outputStr += '%s  ' % var
outputStr += '\n'
fid.write(outputStr)

  outputStr = '%15.7e ' % curTime
  for var in data.keys():
value = data[var].GetRange()[0]
outputStr += '%15.7e' % value
  outputStr += '\n'
  fid.write(outputStr)

def DoCoProcessing(datadescription):
"Callback to do co-processing for current timestep"
global coprocessor
global lastTime
global deltaT
global comm
global paths

# Update the coprocessor by providing it the newly generated simulation data.
# If the pipeline hasn't been setup yet, this will setup the pipeline.
coprocessor.UpdateProducers(datadescription)

curTime = datadescription.GetTime()
timestep = datadescription.GetTimeStep()

if curTime >= lastTime + deltaT:
  lastTime = curTime

  coprocessor.Pipeline.volumeIntegrals.UpdatePipeline()

  WriteDataFile(comm, coprocessor.Pipeline.volumeIntegrals,
paths['outdir']+'/post/volumeIntegrals.dat',
curTime,
CellData = True)

  # Check if we need to write images and the first flameArea datasets to disk
  coprocessor.Pipeline.flameArea.UpdatePipeline()
  area = np.zeros(1)

  if comm.Get_rank() == 0:
try:
  area[0] = coprocessor.Pipeline.flameArea.CellData['Area'].GetRange()[0]
except AttributeError:
  area[0] = 0.0
  else:
area[0] = 0.0

  comm.Bcast(area, root=0)
  if area[0] > 1e-9:
# Write image capture (Last arg: rescale lookup table), if appropriate.
coprocessor.WriteImages(datadescription, rescale_lookuptable=False)

# Write out the surface and the area
coprocessor.Pipeline.flameSurfaceWriter.FileName = paths['outdir']+'/post/flameSurface_%07i.vtm' % timestep
coprocessor.Pipeline.flameSurfaceWriter.UpdatePipeline(curTime)

WriteDataFile(comm, coprocessor.Pipeline.flameArea,
  paths['outdir']+'/post/flameArea.dat',
  curTime,
  CellData = False)


  # Check for the second flame surface
  coprocessor.Pipeline.flameArea2.UpdatePipeline()
  area = np.zeros(1)

  if comm.Get_rank() == 0:
try:
  area[0] = coprocessor.Pipeline.flameArea2.CellData['Area'].GetRange()[0]
except AttributeError:
  area[0] = 0.0
  else:
area[0] = 0.0

  comm.Bcast(area, root=0)
  if area[0] > 1e-9:
# Write image capture (Last arg: rescale lookup table), if appropriate.