Did you also include the w.UpdatePipeline() part after the writer was created? It won't create the file unless that's in there.<br><br>Andy<br><br><div class="gmail_quote">On Thu, Mar 21, 2013 at 3:17 PM, Pettey . Lucas <span dir="ltr"><<a href="mailto:LPettey@drc.com" target="_blank">LPettey@drc.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Andy,<br>
<br>
Changing to Scalars=['POINTS','psi01'] did not change anything. I still get the same error and blank image.<br>
<br>
I tried adding the XMLPImageDataWriter command, but no file 'fullgrid.pvti' was produced so I can't tell what is being populated into the grid.<br>
<br>
I hardwired the C++ adaptor to 'psi01' rather than the character array name in case their was a problem with the way the string was passed. I had checked that name contained the correct characters, but I figured hardwiring would be safer at this stage.<br>
<br>
Anything else that you think I should try?<br>
<br>
Thanks,<br>
Lucas<br>
<br>
________________________________<br>
From: <a href="mailto:paraview-bounces@paraview.org">paraview-bounces@paraview.org</a> [<a href="mailto:paraview-bounces@paraview.org">paraview-bounces@paraview.org</a>] on behalf of Andy Bauer [<a href="mailto:andy.bauer@kitware.com">andy.bauer@kitware.com</a>]<br>
Sent: Thursday, March 21, 2013 12:49 PM<br>
To: Pettey . Lucas<br>
<div class="im">Cc: paraview [<a href="mailto:paraview@paraview.org">paraview@paraview.org</a>]<br>
Subject: Re: [Paraview] CoProcessing with Fortran and ParaView 3.98<br>
<br>
</div><div class="im">I think the problem is that:<br>
Clip1 = Clip( guiName="Clip1", InsideOut=0, Scalars=['POINTS', ''], Value=95.0, UseValueAsOffset=0, Crinkleclip=0, ClipType="Scalar" )<br>
<br>
Should probably have the name of the point data field set in it. Do this with (assuming psi01 is the name of the field you want to clip with respect to:<br>
Scalars=['POINTS', 'psi01']<br>
<br>
<br>
My guess is that the scalar field isn't getting set properly. Try adding in:<br>
w = XMLPImageDataWriter(Input=psi01_, FileName = 'fullgrid.pvti')<br>
w.UpdatePipeline()<br>
<br>
right after the CreateProducer() method to see what's in your grid.<br>
<br>
By the way, does adding in the "//char(0)" work for converting Fortran strings to C/C++ strings (e.g. in "call addfield(psi,"psi01"//char(0))") ? If it does, will it also work for Fortran77? That's a neat little trick!<br>
<br>
Andy<br>
<br>
<br>
</div><div class="im">On Wed, Mar 20, 2013 at 5:29 PM, Pettey . Lucas <<a href="mailto:LPettey@drc.com">LPettey@drc.com</a><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a>>> wrote:<br>
I have solved the C++ to Fortran array issues, but I am still not getting a meaningful visualization. I have attached the latest files, including a screenshot of the exported state and the resulting image that is produced by the coprocessor.<br>
<br>
The coprocessor is generating an error:<br>
<br>
ERROR: In /Users/lucaspettey/ParaView-git/ParaView/VTK/Filters/General/vtkTableBasedClipDataSet.cxx, line 1956<br>
vtkPVClipDataSet (0x7feaa4998550): no input scalars.<br>
<br>
This error is not appearing at the time=0 co-process, but it is appearing each time step the co-processor is being asked to write an image. All of the images are the same blank axes as in the attached image0.png.<br>
<br>
I'm not sure what to try next so any help is appreciated.<br>
<br>
Thanks,<br>
Lucas<br>
<br>
________________________________________<br>
</div>From: <a href="mailto:paraview-bounces@paraview.org">paraview-bounces@paraview.org</a><mailto:<a href="mailto:paraview-bounces@paraview.org">paraview-bounces@paraview.org</a>> [<a href="mailto:paraview-bounces@paraview.org">paraview-bounces@paraview.org</a><mailto:<a href="mailto:paraview-bounces@paraview.org">paraview-bounces@paraview.org</a>>] on behalf of Pettey . Lucas [<a href="mailto:LPettey@drc.com">LPettey@drc.com</a><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a>>]<br>
<div class="im">Sent: Wednesday, March 20, 2013 11:45 AM<br>
To: Andy Bauer<br>
</div><div class="im">Cc: paraview [<a href="mailto:paraview@paraview.org">paraview@paraview.org</a><mailto:<a href="mailto:paraview@paraview.org">paraview@paraview.org</a>>]<br>
Subject: Re: [Paraview] CoProcessing with Fortran and ParaView 3.98<br>
<br>
</div><div class="im">So after some troubleshooting it seems that C++ doesn't like the complex data type the way Fortran is sending it.<br>
<br>
The code runs when I just send a real value. Images appear, but of course they are meaningless since it isn't the correct data. I may try passing the real and imaginary separate. I may also look into row/column order issues although this grid is 100x100x100 so it shouldn't affect anything.<br>
<br>
Lucas<br>
________________________________________<br>
</div>From: <a href="mailto:paraview-bounces@paraview.org">paraview-bounces@paraview.org</a><mailto:<a href="mailto:paraview-bounces@paraview.org">paraview-bounces@paraview.org</a>> [<a href="mailto:paraview-bounces@paraview.org">paraview-bounces@paraview.org</a><mailto:<a href="mailto:paraview-bounces@paraview.org">paraview-bounces@paraview.org</a>>] on behalf of Pettey . Lucas [<a href="mailto:LPettey@drc.com">LPettey@drc.com</a><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a>>]<br>
<div class="im">Sent: Wednesday, March 20, 2013 11:12 AM<br>
To: Andy Bauer<br>
</div><div class="im">Cc: paraview [<a href="mailto:paraview@paraview.org">paraview@paraview.org</a><mailto:<a href="mailto:paraview@paraview.org">paraview@paraview.org</a>>]<br>
Subject: Re: [Paraview] CoProcessing with Fortran and ParaView 3.98<br>
<br>
</div><div class="im">Hi Andy and Berk,<br>
<br>
Here is the error that I am receiving.<br>
<br>
I know that it is initializing ok and failing at the first call to the co-processor. I have attached my python script output as well as my test code and adaptors.<br>
<br>
CoProcessingExample(19247,0x7fff77470180) malloc: *** mmap(size=3739605888335872) failed (error code=12)<br>
*** error: can't allocate region<br>
*** set a breakpoint in malloc_error_break to debug<br>
Operating system error: Cannot allocate memory<br>
Memory allocation failed<br>
<br>
Thanks for the help,<br>
Lucas<br>
<br>
</div>From: Andy Bauer [<a href="mailto:andy.bauer@kitware.com">andy.bauer@kitware.com</a><mailto:<a href="mailto:andy.bauer@kitware.com">andy.bauer@kitware.com</a>>]<br>
<div class="im">Sent: Wednesday, March 20, 2013 9:08 AM<br>
To: Pettey . Lucas<br>
</div>Cc: Berk Geveci; paraview [<a href="mailto:paraview@paraview.org">paraview@paraview.org</a><mailto:<a href="mailto:paraview@paraview.org">paraview@paraview.org</a>>]<br>
<div class="im">Subject: Re: [Paraview] CoProcessing with Fortran and ParaView 3.98<br>
<br>
</div><div class="im">Hi Lucas,<br>
<br>
Can you share your co-processing script and any run output? There may be warnings in the output or info in the script that will help figure that out.<br>
<br>
Thanks,<br>
Andy<br>
<br>
</div><div class="im">On Wed, Mar 20, 2013 at 10:04 AM, Pettey . Lucas <<a href="mailto:LPettey@drc.com">LPettey@drc.com</a><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a>><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a>>>> wrote:<br>
Hi Berk,<br>
<br>
Thanks for the help.<br>
<br>
I found that this example is written based on an earlier version of VTK so some of the function calls have changed. Based on your advice and some of the examples contained in the ParaView wiki, I have gotten it to compile and run.<br>
<br>
Now I just have to figure out why it isn't generating any VTK objects....<br>
<br>
Probably something wrong with my python co-processing script.<br>
<br>
Thanks<br>
Lucas<br>
________________________________<br>
</div>From: Berk Geveci [<a href="mailto:berk.geveci@kitware.com">berk.geveci@kitware.com</a><mailto:<a href="mailto:berk.geveci@kitware.com">berk.geveci@kitware.com</a>><mailto:<a href="mailto:berk.geveci@kitware.com">berk.geveci@kitware.com</a><mailto:<a href="mailto:berk.geveci@kitware.com">berk.geveci@kitware.com</a>>>]<br>
<div class="im">Sent: Tuesday, March 19, 2013 6:48 PM<br>
To: Pettey . Lucas<br>
</div>Cc: paraview [<a href="mailto:paraview@paraview.org">paraview@paraview.org</a><mailto:<a href="mailto:paraview@paraview.org">paraview@paraview.org</a>><mailto:<a href="mailto:paraview@paraview.org">paraview@paraview.org</a><mailto:<a href="mailto:paraview@paraview.org">paraview@paraview.org</a>>>]<br>
<div class="im">Subject: Re: [Paraview] CoProcessing with Fortran and ParaView 3.98<br>
<br>
You can simply delete this line:<br>
<br>
Grid->SetNumberOfScalarComponents(2);<br>
<br>
It is redundant given the following later:<br>
<br>
field->SetNumberOfComponents(2);<br>
<br>
You also need to #include <vtkSmartPointer.h>.<br>
<br>
Best,<br>
-berk<br>
<br>
<br>
<br>
</div><div><div class="h5">On Tue, Mar 19, 2013 at 2:31 PM, Pettey . Lucas <<a href="mailto:LPettey@drc.com">LPettey@drc.com</a><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a>><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a>>><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a>><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a><mailto:<a href="mailto:LPettey@drc.com">LPettey@drc.com</a>>>>> wrote:<br>
Hello everyone,<br>
<br>
I am trying to run a simple Fortran coprocessing example with a complex data type. I am following the example here:<br>
<br>
<a href="http://en.wikibooks.org/wiki/Parallel_Spectral_Numerical_Methods/Visualization_with_ParaView_CoProcessing#ParaView_CoProcessing_Resources" target="_blank">http://en.wikibooks.org/wiki/Parallel_Spectral_Numerical_Methods/Visualization_with_ParaView_CoProcessing#ParaView_CoProcessing_Resources</a><br>
<br>
I have written a basic Fortran code that outputs the data in binary, opened that in ParaView, output the coprocessing script and then written the adaptor following the example in the above link. My adaptor looks like:<br>
<br>
extern "C" void createcpimagedata_(int* nx, int* ny, int* nz)<br>
{<br>
if (!ParaViewCoProcessing::GetCoProcessorData()) {<br>
vtkGenericWarningMacro("Unable to access CoProcessorData.");<br>
return;<br>
}<br>
<br>
// The simulation grid is a 2-dimensional topologically and geometrically<br>
// regular grid. In VTK/ParaView, this is considered an image data set.<br>
vtkImageData* Grid = vtkImageData::New();<br>
<br>
// assuming dimZ == 1 for now<br>
Grid->SetDimensions(*nx, *ny, *nz);<br>
Grid->SetNumberOfScalarComponents(2);<br>
<br>
// Setting the Origin and Spacing are also options.<br>
<br>
// Name should be consistent between here, Fortran and Python client script.<br>
ParaViewCoProcessing::GetCoProcessorData()->GetInputDescriptionByName("input")->SetGrid(Grid);<br>
}<br>
<br>
// Add field(s) to the data container.<br>
// Separate from above because this will be dynamic, grid is static.<br>
// Might be an issue, VTK probably assumes row major, but<br>
// omeg probably passed column major...<br>
// by hand name mangling for fortran<br>
extern "C" void addfield_(double* scalars, char* name)<br>
{<br>
vtkCPInputDataDescription *idd = ParaViewCoProcessing::GetCoProcessorData()->GetInputDescriptionByName("input");<br>
<br>
vtkImageData* Image = vtkImageData::SafeDownCast(idd->GetGrid());<br>
<br>
if (!Image) {<br>
vtkGenericWarningMacro("No adaptor grid to attach field data to.");<br>
return;<br>
}<br>
<br>
<br>
// field name must match that in the fortran code.<br>
if (idd->IsFieldNeeded(name)) {<br>
vtkSmartPointer<vtkDoubleArray> field = vtkSmartPointer<vtkDoubleArray>::New();<br>
field->SetNumberOfComponents(2);<br>
field->SetName(name);<br>
field->SetArray(scalars, 2* Image->GetNumberOfPoints(), 1);<br>
Image->GetPointData()->AddArray(field);<br>
<br>
}<br>
}<br>
<br>
I am getting some errors from the Makefile generated by Cmake:<br>
<br>
/Users/lucaspettey/codes/BEC/coprocVTK.cxx:40:38: error: too few arguments to function call, expected 2, have 1<br>
Grid->SetNumberOfScalarComponents(2);<br>
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ^<br>
/Users/lucaspettey/ParaView-git/ParaView/VTK/Common/DataModel/vtkImageData.h:298:3: note: 'SetNumberOfScalarComponents' declared here<br>
static void SetNumberOfScalarComponents( int n, vtkInformation* meta_data);<br>
^<br>
/Users/lucaspettey/codes/BEC/coprocVTK.cxx:67:5: error: unknown type name 'vtkSmartPointer'; did you mean 'vtkSmartPointerBase'?<br>
vtkSmartPointer<vtkDoubleArray> field = vtkSmartPointer<vtkDoubleArray>::New();<br>
^~~~~~~~~~~~~~~<br>
vtkSmartPointerBase<br>
/Users/lucaspettey/ParaView-git/ParaView/VTK/Common/Core/vtkOStreamWrapper.h:36:7: note: 'vtkSmartPointerBase' declared here<br>
class vtkSmartPointerBase;<br>
^<br>
/Users/lucaspettey/codes/BEC/coprocVTK.cxx:67:20: error: expected unqualified-id<br>
vtkSmartPointer<vtkDoubleArray> field = vtkSmartPointer<vtkDoubleArray>::New();<br>
^<br>
/Users/lucaspettey/codes/BEC/coprocVTK.cxx:68:5: error: use of undeclared identifier 'field'<br>
field->SetNumberOfComponents(2);<br>
^<br>
/Users/lucaspettey/codes/BEC/coprocVTK.cxx:69:5: error: use of undeclared identifier 'field'<br>
field->SetName(name);<br>
^<br>
/Users/lucaspettey/codes/BEC/coprocVTK.cxx:70:5: error: use of undeclared identifier 'field'<br>
field->SetArray(scalars, 2* Image->GetNumberOfPoints(), 1);<br>
^<br>
/Users/lucaspettey/codes/BEC/coprocVTK.cxx:71:37: error: use of undeclared identifier 'field'<br>
Image->GetPointData()->AddArray(field);<br>
^<br>
<br>
Any help is appreciated.<br>
<br>
Thanks,<br>
Lucas<br>
<br>
________________________________<br>
This electronic message transmission and any attachments that accompany it contain information from DRCŪ (Dynamics Research Corporation) or its subsidiaries, or the intended recipient, which is privileged, proprietary, business confidential, or otherwise protected from disclosure and is the exclusive property of DRC and/or the intended recipient. The information in this email is solely intended for the use of the individual or entity that is the intended recipient. If you are not the intended recipient, any use, dissemination, distribution, retention, or copying of this communication, attachments, or substance is prohibited. If you have received this electronic transmission in error, please immediately reply to the author via email that you received the message by mistake and also promptly and permanently delete this message and all copies of this email and any attachments. We thank you for your assistance and apologize for any inconvenience.<br>
_______________________________________________<br>
</div></div>Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><<a href="http://www.kitware.com" target="_blank">http://www.kitware.com</a>><<a href="http://www.kitware.com" target="_blank">http://www.kitware.com</a>><<a href="http://www.kitware.com" target="_blank">http://www.kitware.com</a>><br>
<div class="HOEnZb"><div class="h5"><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://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><br>
<br>
_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><<a href="http://www.kitware.com" target="_blank">http://www.kitware.com</a>><<a href="http://www.kitware.com" target="_blank">http://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://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><br>
<br>
_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><<a href="http://www.kitware.com" target="_blank">http://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://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><br>
<br>
</div></div></blockquote></div><br>