<div dir="ltr"><div><div><div>Hi Alexandre,<br><br></div>I'm glad it's working for you and thanks for letting other people know what you did.<br><br></div>In my opinion I think that just having the adaptor produce a vtkUnstructuredGrid is the way to go. Certain filters will only take in vtkUnstructuredGrids and others will only take in vtkPolyData. You may be forced to used different scripts for 3D simulations than 1D and 2D.<br><br></div><div>If you get some nice results with Catalyst, feel free to share them with me. I'm always interested in what people are using Catalyst for.<br></div><div><br></div>Regards,<br>Andy<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Nov 27, 2014 at 5:40 AM, Alexandre Ancel <span dir="ltr"><<a href="mailto:alexandre.ancel@cemosis.fr" target="_blank">alexandre.ancel@cemosis.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div><div><div>Hello,<br><br></div>Regarding to the vtkout_type, it is defined in the corresponding header:<br><a href="https://github.com/feelpp/feelpp/blob/develop/feel/feelfilters/exporterVTK.hpp" target="_blank">https://github.com/feelpp/feelpp/blob/develop/feel/feelfilters/exporterVTK.hpp</a><br></div>I guess it is not really an issue as it is a typedef to vtkUnstructuredGrid. Initially I intended to use vtkPolyData for 1D and 2D meshes; and vtkUnstructuredGrid for 3D meshes, but I guess that I will totally switch to vtkUnstucturedGrid as it would be able to handle all the cases.<br><br></div>Thanks for the advice regarding to ForceOutput, I'll change it.<br><br></div>Now back to my issue, I managed to make it work. In the end, it was only me misusing Catalyst.<br></div>Just how I make it work now: (for people that might have a similar issue to mine)<br></div>- Launch Paraview and allow connections for Catalyst<br></div>- Pause the simulation (It will thus pause at the first timestep)<br></div>- Launch the simulation and visualize the data in Paraview<br><br></div><div>Thanks again for your time and help,<br><br>Best regards,<br></div><div>Alexandre Ancel<br></div><div><div><div><div> <br></div></div></div></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Nov 26, 2014 at 6:06 PM, Andy Bauer <span dir="ltr"><<a href="mailto:andy.bauer@kitware.com" target="_blank">andy.bauer@kitware.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Hi,<br><br></div>My guess is that the templating isn't working the way you think it is. I couldn't see where vtkout_type was being set but in general you probably don't want to template on the VTK data sets since the API where the work is done is different for each of them. Also, input->GetClassName() should be returning the name of the actual class you want to use, e.g. vtkPolyData, vtkUnstructuredGrid, etc.<br><br></div>Another thing I noticed is that you probably want to get rid of the SetForceOutput(true) in:<span><br><table><tbody><tr><td><span>if</span>(inSituProcessor-><span>RequestDataDescription</span>(dataDescription.<span>GetPointer</span>()) != <span>0</span>)</td>
</tr>
<tr>
<td> {</td>
</tr>
<tr>
<td> dataDescription-><span>GetInputDescriptionByName</span>(<span><span>"</span>input<span>"</span></span>)-><span>SetGrid</span>(out);</td>
</tr>
<tr>
<td> dataDescription-><span>SetForceOutput</span>(<span>true</span>);</td>
</tr>
<tr>
<td> std::cout << <span><span>"</span>CoProcess <span>"</span></span> << inSituProcessor-><span>CoProcess</span>(dataDescription.<span>GetPointer</span>())<< std::endl;</td>
</tr>
<tr>
<td> }</td></tr></tbody></table><br></span></div>This will ensure that all output from Catalyst is generated when possibly only a single output was requested in the pipelines.<br><br></div>Regards,<br>Andy<br></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Nov 26, 2014 at 11:33 AM, Alexandre Ancel <span dir="ltr"><<a href="mailto:alexandre.ancel@cemosis.fr" target="_blank">alexandre.ancel@cemosis.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div>Hello,<br><br></div>An update on this issue that I'm still having:<br></div>I
tried to pinpoint the source of the problem. I went through Paraview
and VTK code, and found that there is a check that fails to pass
apparently, resulting in the errors previously mentionned. This check fails on the Paraview side.<br></div><div>In
VTK/Common/ExecutionModel/vtkDemandDrivenPipeline.cxx (from which vtkPVPostFilterExecutive inherits), the
InputTypeIsValid (l.763) function seems to check for data type, and fails to find a
matching type when I try to visualize the data extracted from the
PVTrivialProducer.<br></div><div>The recorded type (with
input->GetClassName()) is vtkDataObject and the IsA function is
unable to match it with any of the type in the list iterated over.<br><br></div><div>Did I went to far in the code ? Could the information about the Unstructured grid nature be lost when I extract the data in Paraview ?<br><br></div><div>By
the way, this list is probably not the right place to talk about it, since it is related to VTK,
but isn't the error printed when no match is found a bit confusing for
the user ?<br></div><div>I mean if the function does not find a valid
type, it will by default choose the first type of the list, resulting in
the kind of errors I had (vtkPolyData in my case):<span><br>vtkPVPostFilterExecutive (0x3efc890): Input for connection index 0 on
input port index 0 for algorithm vtkPVPostFilter(0x3efbe70) is of type
vtkDataObject, but a vtkPolyData is required.<br><br></span></div><div>Whereas the type checked was not matched with any type, thus not a vtkPolyData.<br><br></div>By the way, if you want to have a look at the C++ code used in our library, it is available in:<br><a href="https://github.com/feelpp/feelpp/blob/develop/feel/feelfilters/exporterVTK_impl.hpp" target="_blank">https://github.com/feelpp/feelpp/blob/develop/feel/feelfilters/exporterVTK_impl.hpp</a><br></div>With:<br>- the initialization part starting at line 115<br>- the coprocessing part starting at line 406<br><br></div>Thanks again for the help,<br></div>Best regards,<br>Alexandre Ancel<br><div><div><div><div><br></div></div></div></div></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Nov 24, 2014 at 4:56 PM, Alexandre Ancel <span dir="ltr"><<a href="mailto:alexandre.ancel@cemosis.fr" target="_blank">alexandre.ancel@cemosis.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Hello,<br><br></div>I tested writing data from the python script, and the data is apparently written correctly<br></div> when opening the files with Paraview. The issue is only present I try to extract the data in Paraview.<br></div><div>I'll continue to investigate both in my code and Paraview. If you have other pieces of advice as to where I could be looking, I would be interested.<br></div><div><br></div>I will also test if I can reproduce the issue with one of the example.<br><br></div>Regards,<br>Alexandre Ancel<br></div><div class="gmail_extra"><div><div><br><div class="gmail_quote">On Fri, Nov 21, 2014 at 10:57 PM, Andy Bauer <span dir="ltr"><<a href="mailto:andy.bauer@kitware.com" target="_blank">andy.bauer@kitware.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>Hi Alexandre,<br><br></div>If you could reproduce the errors that you're getting in one of the examples, we could try to track that down. I tried using your script and didn't have any issues with the CxxFullExample in the Catalyst example code and PV 4.2.<br><br></div><div><br></div>Regards,<br>Andy<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Fri, Nov 21, 2014 at 12:14 PM, Dan Lipsa <span dir="ltr"><<a href="mailto:dan.lipsa@kitware.com" target="_blank">dan.lipsa@kitware.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><br></div><div>When I click on the icon located near the producer to visualize the data, I get an Extract component added to the pipeline. Trying to visualize this, then throws me the following errors:<br><br>---<br>ERROR: In /home/ancel/Downloads/ParaView-v4.2.0-source/VTK/Common/ExecutionModel/vtkDemandDrivenPipeline.cxx, line 810<br>vtkPVPostFilterExecutive (0x3efc890): Input for connection index 0 on input port index 0 for algorithm vtkPVPostFilter(0x3efbe70) is of type vtkDataObject, but a vtkPolyData is required.<br><br><br>ERROR: In /home/ancel/Downloads/ParaView-v4.2.0-source/VTK/Common/ExecutionModel/vtkDemandDrivenPipeline.cxx, line 810<br>vtkPVPostFilterExecutive (0x3efc890): Input for connection index 0 on input port index 0 for algorithm vtkPVPostFilter(0x3efbe70) is of type vtkDataObject, but a vtkPolyData is required.<br><br><br>ERROR: In /home/ancel/Downloads/ParaView-v4.2.0-source/ParaViewCore/ServerManager/Rendering/vtkSMParaViewPipelineControllerWithRendering.cxx, line 504<br>vtkSMParaViewPipelineControllerWithRendering (0x38bbc00): Data cannot be shown in the defaulted render view!!<br></div><div>---<br></div></div></blockquote><div><br></div></span><div>I have not run an example for an unstructured grid. It should work similarly though.</div><div>I recently wrote a Catalyst adapter for an unstructured grid. I will run Catalyst Live with that to see how that works.</div><div><br></div><div>You could try however to save images or entire data from the pipeline script as it is done in the Catalyst examples.</div><span><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>A few more questions come in my mind now:<br></div><div>- If I understand correctly and if the user wants to pause the simulation, for example when he has only 1 timestep, he would have to do the following steps:<br></div><div>* Connect with catalyst in Paraview<br></div><div>* Set a breakpoint for time step 0<br></div><div>* launch the simulation<br></div><div>Is that alright ?<br></div><div>Would it be possible to do this with the play/next frame paraview buttons ?<br></div></div></blockquote><div><br></div></span><div>The communication between catalyst and the simulation happens only at time steps. Because of that, Catalyst could only stop the simulation at time step 2. </div><span><div> </div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>- Looking at the last lines of the blog post concerning pvserver. Is pvserver catalyst enabled by default ?<br></div></div></blockquote></span><div>Yes.</div><span><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div><div>Does it mean that I can build the following pipeline ?<br><span>Paraview (local computer) <-> pvserver (cluster frontal node) <-> catalyst enabled simulation (cluster nodes)</span> <br></div></div></blockquote><div><br></div></span><div> You need to enable MPI to be able to run this. pvserver is the paraview server that does all the rendering, only images are sent to ParaView (GUI).</div><span><font color="#888888"><div><br></div><div>Dan</div><div><br></div></font></span></div></div></div>
<br></div></div><span>_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/paraview" target="_blank">http://public.kitware.com/mailman/listinfo/paraview</a><br>
<br></span></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><br></div></div><span>-- <br><div><div dir="ltr"><div><div dir="ltr"><span><font color="#888888"><div dir="ltr"><font face="courier new, monospace">Alexandre Ancel <br>Docteur, Ingénieur de recherche / Phd, Research Engineer<br></font><div>
<font face="courier new, monospace"><a href="http://www.cemosis.fr" target="_blank">Ce<font color="#888888">mosis</font></a> - <a href="mailto:alexandre.ancel@cemosis.fr" target="_blank">alexandre.ancel@cemosis.fr</a></font></div><div><font face="courier new, monospace">Tel: +33 (0)3 68 8</font><span><font color="#888888"><font face="courier new, monospace"><b>5 02 06</b></font></font></span></div>
<div><font face="courier new, monospace">IRMA - 7, rue René Descartes</font></div><div><font face="courier new, monospace">67 000 Strasbourg, France<br></font></div></div></font></span></div></div></div></div>
</span></div>
</blockquote></div><br><br clear="all"><br>-- <br><div><div dir="ltr"><div><div dir="ltr"><span><font color="#888888"><div dir="ltr"><font face="courier new, monospace">Alexandre Ancel <br>Docteur, Ingénieur de recherche / Phd, Research Engineer<br></font><div>
<font face="courier new, monospace"><a href="http://www.cemosis.fr" target="_blank">Ce<font color="#888888">mosis</font></a> - <a href="mailto:alexandre.ancel@cemosis.fr" target="_blank">alexandre.ancel@cemosis.fr</a></font></div><div><font face="courier new, monospace">Tel: +33 (0)3 68 8</font><span><font color="#888888"><font face="courier new, monospace"><b>5 02 06</b></font></font></span></div>
<div><font face="courier new, monospace">IRMA - 7, rue René Descartes</font></div><div><font face="courier new, monospace">67 000 Strasbourg, France<br></font></div></div></font></span></div></div></div></div>
</div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div><div dir="ltr"><div><div dir="ltr"><span><font color="#888888"><div dir="ltr"><font face="courier new, monospace">Alexandre Ancel <br>Docteur, Ingénieur de recherche / Phd, Research Engineer<br></font><div>
<font face="courier new, monospace"><a href="http://www.cemosis.fr" target="_blank">Ce<font color="#888888">mosis</font></a> - <a href="mailto:alexandre.ancel@cemosis.fr" target="_blank">alexandre.ancel@cemosis.fr</a></font></div><div><font face="courier new, monospace">Tel: +33 (0)3 68 8</font><span><font color="#888888"><font face="courier new, monospace"><b>5 02 06</b></font></font></span></div>
<div><font face="courier new, monospace">IRMA - 7, rue René Descartes</font></div><div><font face="courier new, monospace">67 000 Strasbourg, France<br></font></div></div></font></span></div></div></div></div>
</div>
</div></div></blockquote></div><br></div>