<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Andy,
<div class=""><br class="">
</div>
<div class="">Thank you for your quick answer.</div>
<div class=""><br class="">
</div>
<div class="">You probably already read my previous answer about this issue.</div>
<div class="">The problem was located in NeedToCreateGrid(...), which also clears any field data associated with an existing grid.</div>
<div class="">The solution consists in calling simply this function only once for every new time step.</div>
<div class=""><br class="">
</div>
<div class="">That said, you raised a good point about IsFieldNeeded(…), which was also on my radar.</div>
<div class=""><br class="">
</div>
<div class="">I already observed that RequestDataDescription(datadescription) in the python script sets AllFields in the data description object to On. As a consequence, all the fields specified in the adaptor are passed to Catalyst since IsFieldNeeded() always
 returns true, whether the corresponding field is used in the Catalyst pipeline or not. Since we have improved our adaptor (we should commit it back now) and increased the number of fields we are potentially interest in for coprocesing purpose, this can indeed
 leads to additional memory usage and cpu time.</div>
<div class=""><br class="">
</div>
<div class="">I definitely agree it would therefore be quite useful to have the possibility to request only the desired field variables through the IsFieldNeeded() function in the adaptor.</div>
<div class="">If you have any advice regarding this feature, I would be very interested in trying that out.</div>
<div class=""><br class="">
</div>
<div class="">Thank you for your help!</div>
<div class=""><br class="">
</div>
<div class="">Cheers,</div>
<div class=""><br class="">
</div>
<div class="">Michel</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div>
<blockquote type="cite" class="">
<div class="">On Oct 26, 2015, at 5:57 PM, Andy Bauer <<a href="mailto:andy.bauer@kitware.com" class="">andy.bauer@kitware.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">
<div class="">
<div class="">
<div class="">
<div class="">
<div class="">Hi Michel,<br class="">
<br class="">
</div>
You should be able to pass a single field at a time to Catalyst. I'm not sure where the problem is but my first guess is that maybe you're giving the same name to all of the fields. What does the code that's calling addfield() look like?
<br class="">
<br class="">
</div>
Note that the Catalyst API uses things like idd->IsFieldNeeded("pressure") to check if a field is needed by the pipeline. This has been in the API since nearly the beginning but we've never had a chance to generate Python scripts which can take advantage of
 loading only desired fields. This can potentially save on both execution time and memory usage. This is on my radar again but I'm not sure when it will get done. You can modify the Python scripts though to just request the desired field variables in the RequestDataDescription()
 method and everything should work as desired. Let us know if you want to try that out and need help with it.<br class="">
<br class="">
</div>
Cheers,<br class="">
</div>
Andy<br class="">
</div>
<div class="">
<div class="">
<div class="">
<div class="">
<div class="">
<div class="">
<div class="gmail_extra"><br class="">
<div class="gmail_quote">On Mon, Oct 26, 2015 at 11:22 AM, Michel Rasquin <span dir="ltr" class="">
<<a href="mailto:michel.rasquin@colorado.edu" target="_blank" class="">michel.rasquin@colorado.edu</a>></span> wrote:<br class="">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Hi everyone,<br class="">
<br class="">
I am trying to add some fields to a vtkCPAdaptorAPI object for coprocessing with Catalyst.<br class="">
I rely for that purpose on the successful implementation of the Phasta adaptor provided along with ParaView.<br class="">
See ParaView-v4.4.0-source/CoProcessing/Adaptors/PhastaAdaptor/PhastaAdaptor.cxx.<br class="">
After the initialization of the coprocessing objects and the generation of the grid, the current implementation to add fields in the phasta adaptor relies on the following function:<br class="">
<br class="">
void addfields(… double* dofArray, double* vortArray, double * otherFieldOfInterest … )<br class="">
{<br class="">
  vtkCPInputDataDescription* idd = vtkCPAdaptorAPI::GetCoProcessorData()->GetInputDescriptionByName("input”);<br class="">
  vtkUnstructuredGrid* UnstructuredGrid = vtkUnstructuredGrid::SafeDownCast(idd->GetGrid());<br class="">
  if(!UnstructuredGrid) {<br class="">
    vtkGenericWarningMacro("No unstructured grid to attach field data to.");<br class="">
    return;<br class="">
  }<br class="">
<br class="">
  // now add numerical field data<br class="">
  //velocity<br class="">
  vtkIdType NumberOfNodes = UnstructuredGrid->GetNumberOfPoints();<br class="">
  if(idd->IsFieldNeeded("velocity"))<br class="">
  {<br class="">
    vtkDoubleArray* velocity = vtkDoubleArray::New();<br class="">
    velocity->SetName("velocity");<br class="">
    velocity->SetNumberOfComponents(3);<br class="">
    velocity->SetNumberOfTuples(NumberOfNodes);<br class="">
    for (vtkIdType idx=0; idx<NumberOfNodes; idx++) {<br class="">
      velocity->SetTuple3(idx, dofArray[idx],<br class="">
                          dofArray[idx+ *nshg],<br class="">
                          dofArray[idx+ *nshg*2]);<br class="">
    }<br class="">
    UnstructuredGrid->GetPointData()->AddArray(velocity);<br class="">
    velocity->Delete();<br class="">
  }<br class="">
<br class="">
  if(idd->IsFieldNeeded(“vorticity"))<br class="">
  {<br class="">
    vtkDoubleArray* vorticity = vtkDoubleArray::New();<br class="">
    velocity->SetName(“vorticity");<br class="">
    velocity->SetNumberOfComponents(3);<br class="">
    velocity->SetNumberOfTuples(NumberOfNodes);<br class="">
    for (vtkIdType idx=0; idx<NumberOfNodes; idx++) {<br class="">
      velocity->SetTuple3(idx, vortArray[idx],<br class="">
                          vortArray[idx+ *nshg],<br class="">
                          vortArray[idx+ *nshg*2]);<br class="">
    }<br class="">
    UnstructuredGrid->GetPointData()->AddArray(vorticity);<br class="">
    vorticity->Delete();<br class="">
    }<br class="">
<br class="">
  // etc for any the other fields of interest for Catalyst<br class="">
}<br class="">
<br class="">
Currently, all the fields requested for coprocessing needs to be attached in this function at the same time, using the same pointer to vtkUnstructuredGrid resulting from the SafeDownCast mentioned above. However, I need a more flexible implementation so that
 I can call addfield (with no “s”) as many times as needed and attach a single field to the vtkCPAdaptorAPI object each time this function is called.<br class="">
<br class="">
Concretely, my first implementation is simply the following:<br class="">
<br class="">
void addfield(std::string fieldName, int* NumberOfComp, double* fieldArray)<br class="">
{<br class="">
  vtkCPInputDataDescription* idd = vtkCPAdaptorAPI::GetCoProcessorData()->GetInputDescriptionByName("input");<br class="">
  vtkUnstructuredGrid* UnstructuredGrid = vtkUnstructuredGrid::SafeDownCast(idd->GetGrid());<br class="">
  if(!UnstructuredGrid) {<br class="">
    vtkGenericWarningMacro("No unstructured grid to attach field data to.");<br class="">
    return;<br class="">
  }<br class="">
<br class="">
  // Get number of nodes<br class="">
  vtkIdType NumberOfNodes = UnstructuredGrid->GetNumberOfPoints();<br class="">
<br class="">
  // Add field<br class="">
  if(idd->IsFieldNeeded(fieldName.c_str())) {<br class="">
    vtkDoubleArray* dataArray = vtkDoubleArray::New();<br class="">
    dataArray->SetName(fieldName.c_str());<br class="">
    dataArray->SetNumberOfComponents(*NumberOfComp);<br class="">
    dataArray->SetNumberOfTuples(NumberOfNodes);<br class="">
    // fill in dataArray from fieldArray, NumberOfNodes and NumberOfComp<br class="">
    …<br class="">
    UnstructuredGrid->GetPointData()->AddArray(dataArray);<br class="">
    dataArray->Delete();<br class="">
  }<br class="">
}<br class="">
<br class="">
The problem is that only the last field passed to this new addfield() function can be actually used by Catalyst for coprocessing.<br class="">
Indeed, it appears that all other fields previously passed to addfield() cannot be retrieved from the vtkCPAdaptorAPI object.<br class="">
Consequently, any filter in the Catalyst pipeline that relies on the N-1 first fields (out of N in total) passed to addfields() will be ignored because relevant data is missing.<br class="">
<br class="">
I suspect the issue is in one of the first two lines of the addfield() function, namely<br class="">
<br class="">
  vtkCPInputDataDescription* idd = vtkCPAdaptorAPI::GetCoProcessorData()->GetInputDescriptionByName("input");<br class="">
  vtkUnstructuredGrid* UnstructuredGrid = vtkUnstructuredGrid::SafeDownCast(idd->GetGrid());<br class="">
<br class="">
Could you please let me know if it is possible to pass one single field at a time to the Catalyst adaptor from different locations of the code, or if all the fields must be passed in one shot?<br class="">
<br class="">
Thank you for your help.<br class="">
<br class="">
Best regards,<br class="">
<br class="">
Michel<br class="">
<br class="">
<br class="">
<br class="">
_______________________________________________<br class="">
Powered by <a href="http://www.kitware.com/" rel="noreferrer" target="_blank" class="">
www.kitware.com</a><br class="">
<br class="">
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank" class="">
http://www.kitware.com/opensource/opensource.html</a><br class="">
<br class="">
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" rel="noreferrer" target="_blank" class="">
http://paraview.org/Wiki/ParaView</a><br class="">
<br class="">
Search the list archives at: <a href="http://markmail.org/search/?q=ParaView" rel="noreferrer" target="_blank" class="">
http://markmail.org/search/?q=ParaView</a><br class="">
<br class="">
Follow this link to subscribe/unsubscribe:<br class="">
<a href="http://public.kitware.com/mailman/listinfo/paraview" rel="noreferrer" target="_blank" class="">http://public.kitware.com/mailman/listinfo/paraview</a><br class="">
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>