[Paraview] CoProcessing

Andy Bauer andy.bauer at kitware.com
Tue Aug 17 15:49:12 EDT 2010


For 1:  Can you send me your
Utilities/VTKPythonWrapping/paraview/servermanager.py in the PV source as
well as the
Utilities/VTKPythonWrapping/site-packages/paraview/servermanager.py file in
your PV build?

For 2: I'll take a look at that since I have to do the same thing on a
remote machine I'll be working on.

For 3: No idea why.  If you can set up a simple problem I can take a look at
it.

Andy

On Tue, Aug 17, 2010 at 10:06 AM, Jacques Papper
<jacques.papper at gmail.com>wrote:

> Hi All,
>
> Although I've got the coprocessing to work somehow, I still have a couple
> of problems :
>
> -> I have pulled from git yesterday.
>
> 1. None inside the python scripts generated by the plugin throws up some
> issues :
>
> Traceback (most recent call last):
>   File "<string>", line 2, in <module>
>   File "/users/boreas01/jacques/TESTING/A6_saloon_50k/pvCopro.py", line
> 138, in DoCoProcessing
>     DataRepresentation1.SelectOrientationVectors = [None, '']
>   File
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/Utilities/VTKPythonWrapping/site-packages/paraview/servermanager.py",
> line 207, in __setattr__
>     setter(self, value)
>   File
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/Utilities/VTKPythonWrapping/site-packages/paraview/servermanager.py",
> line 2367, in setProperty
>     return self.GetPropertyValue(propName)
>   File
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/Utilities/VTKPythonWrapping/site-packages/paraview/servermanager.py",
> line 256, in SetPropertyWithName
>     def SetPropertyWithName(self, pname, arg):
>   File
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/Utilities/VTKPythonWrapping/site-packages/paraview/servermanager.py",
> line 799, in SetData
>     if len(values) == 1:
> UnboundLocalError: local variable 'val' referenced before assignment
>
>
> 2. When Compiling with MESA, and switching offscreen rendering ON I get
> problems at runtime when coprocessing (I have no problems saving screenshots
> directly from within paraview) :
>
>
> --- Some of my CMAKE settings :
> //Path to a file.
> OPENGL_INCLUDE_DIR:PATH=/usr/include
>
> //Path to a library.
> OPENGL_gl_LIBRARY:FILEPATH=/usr/lib64/libGL.so
>
> //Path to a library.
> OPENGL_glu_LIBRARY:FILEPATH=/usr/lib64/libGLU.so
>
> //Path to a file.
> OPENGL_xmesa_INCLUDE_DIR:PATH=/usr/include
>
> //Path to a file.
> OSMESA_INCLUDE_DIR:PATH=/usr/include
>
> //Path to a library.
> OSMESA_LIBRARY:FILEPATH=/usr/lib64/libOSMesa.so
>
> ----- RunTime error message ----
>
>
> #0  Foam::error::printStack(Foam::Ostream&) in
> "/users/boreas01/jacques/FOAMpro-1.6-dev/FOAM-dev/lib/linux64GccDPOpt/libFOAM.so"
> #1  Foam::sigSegv::sigSegvHandler(int) in
> "/users/boreas01/jacques/FOAMpro-1.6-dev/FOAM-dev/lib/linux64GccDPOpt/libFOAM.so"
> #2  ?? in "/lib64/libc.so.6"
> #3  __rawmemchr in "/lib64/libc.so.6"
> #4  ?? in "/lib64/libc.so.6"
> #5  vsscanf in "/lib64/libc.so.6"
> #6  _IO_sscanf in "/lib64/libc.so.6"
> #7  vtkOpenGLExtensionManager::ReadOpenGLExtensions() in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkRendering.so.pv3.9"
> #8  vtkOpenGLExtensionManager::Update() in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkRendering.so.pv3.9"
> #9  vtkOpenGLExtensionManager::ExtensionSupported(char const*) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkRendering.so.pv3.9"
> #10  vtkOpenGLRenderWindow::OpenGLInit() in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkRendering.so.pv3.9"
> #11  vtkXOpenGLRenderWindow::Start() in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkRendering.so.pv3.9"
> #12  vtkRenderWindow::DoStereoRender() in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkRendering.so.pv3.9"
> #13  vtkRenderWindow::DoFDRender() in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkRendering.so.pv3.9"
> #14  vtkRenderWindow::DoAARender() in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkRendering.so.pv3.9"
> #15  vtkRenderWindow::Render() in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkRendering.so.pv3.9"
> #16  vtkXOpenGLRenderWindow::Render() in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkRendering.so.pv3.9"
> #17  vtkXOpenGLRenderWindowCommand(vtkClientServerInterpreter*,
> vtkObjectBase*, char const*, vtkClientServerStream const&,
> vtkClientServerStream&) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkRenderingCS.so"
> #18  vtkClientServerInterpreter::ProcessCommandInvoke(vtkClientServerStream
> const&, int) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkClientServer.so"
> #19  vtkClientServerInterpreter::ProcessOneMessage(vtkClientServerStream
> const&, int) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkClientServer.so"
> #20  vtkClientServerInterpreter::ProcessStream(vtkClientServerStream
> const&) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkClientServer.so"
> #21  vtkSelfConnection::ProcessStreamLocally(vtkClientServerStream&) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
> #22
> vtkMPISelfConnection::SendStreamToDataServerRoot(vtkClientServerStream&) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
> #23  vtkProcessModuleConnection::SendStream(unsigned int,
> vtkClientServerStream&) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
> #24  vtkProcessModuleConnectionManager::SendStream(long long, unsigned int,
> vtkClientServerStream&, int) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
> #25  vtkProcessModule::SendStream(long long, unsigned int,
> vtkClientServerStream&, int) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
> #26  vtkSMRenderViewProxy::PerformRender() in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerManager.so"
> #27  vtkSMViewProxy::StillRender() in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerManager.so"
> #28  vtkSMRenderViewProxy::CaptureWindow(int) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerManager.so"
> #29  vtkSMRenderViewProxy::WriteImage(char const*, char const*, int) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerManager.so"
> #30  ?? at vtkSMRenderViewProxyPython.cxx:0
> #31  vtkPythonUtil::CallOverloadedMethod(PyMethodDef*, _object*, _object*)
> in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPythonCore.so.pv3.9"
> #32  PyObject_Call in "/usr/lib64/libpython2.6.so.1.0"
> #33  PyEval_EvalFrameEx in "/usr/lib64/libpython2.6.so.1.0"
> #34  PyEval_EvalCodeEx in "/usr/lib64/libpython2.6.so.1.0"
> #35  PyEval_EvalFrameEx in "/usr/lib64/libpython2.6.so.1.0"
> #36  PyEval_EvalCodeEx in "/usr/lib64/libpython2.6.so.1.0"
> #37  PyEval_EvalFrameEx in "/usr/lib64/libpython2.6.so.1.0"
> #38  PyEval_EvalFrameEx in "/usr/lib64/libpython2.6.so.1.0"
> #39  PyEval_EvalCodeEx in "/usr/lib64/libpython2.6.so.1.0"
> #40  PyEval_EvalCode in "/usr/lib64/libpython2.6.so.1.0"
> #41  ?? in "/usr/lib64/libpython2.6.so.1.0"
> #42  PyRun_StringFlags in "/usr/lib64/libpython2.6.so.1.0"
> #43  PyRun_SimpleStringFlags in "/usr/lib64/libpython2.6.so.1.0"
> #44  vtkPVPythonInterpretor::RunSimpleString(char const*) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVPythonInterpretor.so"
> #45  vtkCPPythonScriptPipeline::CoProcess(vtkCPDataDescription*) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkCoProcessor.so"
> #46  vtkCPProcessor::CoProcess(vtkCPDataDescription*) in
> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkCoProcessor.so"
>
>
> When compiling without MESA and OFFSCREEN then it works fine but I get an
> annoying window that pops up all the time...
>
> 3. I am able to use the Parallel Multi Block Data Set Writer, on a very
> simple hex mesh. But when I use it on a more complex mesh (including
> POLYHEDRA cells) I get the following error :
>
> When opening the file in paraview :
>
>   ERROR: In
> /users/boreas01/jacques/PARAVIEW/ParaView/VTK/IO/vtkXMLParser.cxx, line 483
>
> vtkXMLDataParser (0x18aa700): Error parsing XML in stream at line 1, column
> 7, byte index 7: not well-formed (invalid token)
>
>  ERROR: In
> /users/boreas01/jacques/PARAVIEW/ParaView/VTK/IO/vtkXMLReader.cxx, line 362
>
> vtkXMLUnstructuredGridReader (0x18cbbc0): Error parsing input file.
> ReadXMLInformation aborting.
>
>  ERROR: In
> /users/boreas01/jacques/PARAVIEW/ParaView/VTK/Filtering/vtkExecutive.cxx,
> line 756
>
> vtkStreamingDemandDrivenPipeline (0xdadf60): Algorithm
> vtkXMLUnstructuredGridReader(0x18cbbc0) returned failure for request:
> vtkInformation (0x127b170)
>
> Debug: Off
>
> Modified Time: 80976
>
> Reference Count: 1
>
> Registered Events: (none)
>
> Request: REQUEST_INFORMATION
>
> ALGORITHM_AFTER_FORWARD: 1
>
> FORWARD_DIRECTION: 0
>
>
>
> When running the coprocessor  I do not get any error message. And I am
> still able to ouput png screenshots (when not using OFFSCREEN)... Is there
> an incompatibility between the writers and the polyhedra cell type ? Or is
> it possible that somewhere I am not filling in the cells / points properly ?
>
>
> Best,
> Jacques
>
>
>
>
>
> 2010/8/16 pat marion <pat.marion at kitware.com>
>
>> cc'ing the list.
>>
>> Pat
>>
>>
>> ---------- Forwarded message ----------
>> From: Jacques Papper <jacques.papper at gmail.com>
>> Date: Mon, Aug 16, 2010 at 1:55 PM
>> Subject: Re: [Paraview] CoProcessing
>> To: pat marion <pat.marion at kitware.com>
>>
>>
>> Once again... My mistake .. I'm not used to submodules !
>>
>> 2010/8/16 Jacques Papper <jacques.papper at gmail.com>
>>
>> Hi, I just pulled my git..
>>>
>>> CMake Error in Utilities/VTKClientServer/Wrapping/CMakeLists.txt:
>>>   Cannot find source file "vtkParsePreprocess.c".  Tried extensions .c .C
>>>   .c++ .cc .cpp .cxx .m .M .mm .h .hh .h++ .hm .hpp .hxx .in .txx
>>>
>>> I tried to configure and generate cmake but I get the same error when
>>> generating...
>>>
>>> Should I start from a clean cmake cache ?
>>>
>>> Jacques
>>>
>>> 2010/8/16 pat marion <pat.marion at kitware.com>
>>>
>>> If you can't upgrade your paraview version you could try hand applying
>>>> the patch to /users/boreas01/jacques/
>>>>
>>>>>
>>>>> PARAVIEW/ParaView-bin/Utilities/VTKPythonWrapping/site-packages/paraview/servermanager.py
>>>>
>>>>
>>>> The patch:
>>>>
>>>>
>>>> http://paraview.org/gitweb?p=ParaView.git;a=commitdiff;h=23027770320168c3f7877c6e4150d5ca834f3b9b
>>>>
>>>> Pat
>>>>
>>>>
>>>> On Mon, Aug 16, 2010 at 11:01 AM, Jacques Papper <
>>>> jacques.papper at gmail.com> wrote:
>>>>
>>>>> Hi,
>>>>>
>>>>> I'm sorry I was just rerunning to reproduce the error, and did not
>>>>> realize I missed changing my input back to input.foam instead of toto.foam
>>>>> ....
>>>>>
>>>>> Here is the actual error message :
>>>>>
>>>>> Traceback (most recent call last):
>>>>>   File "<string>", line 2, in <module>
>>>>>   File "/users/boreas01/jacques/TESTING/A6_saloon_50k/pvCopro.py", line
>>>>> 138, in DoCoProcessing
>>>>>     DataRepresentation1.SelectOrientationVectors = [None, '']
>>>>>   File
>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/Utilities/VTKPythonWrapping/site-packages/paraview/servermanager.py",
>>>>> line 207, in __setattr__
>>>>>     setter(self, value)
>>>>>   File
>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/Utilities/VTKPythonWrapping/site-packages/paraview/servermanager.py",
>>>>> line 2367, in setProperty
>>>>>     return self.SetPropertyWithName(propName, value)
>>>>>   File
>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/Utilities/VTKPythonWrapping/site-packages/paraview/servermanager.py",
>>>>> line 256, in SetPropertyWithName
>>>>>     prop.SetData(arg)
>>>>>   File
>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/Utilities/VTKPythonWrapping/site-packages/paraview/servermanager.py",
>>>>> line 799, in SetData
>>>>>     self.SMProperty.SetElement(3,  str(val))
>>>>> UnboundLocalError: local variable 'val' referenced before assignment
>>>>>
>>>>>
>>>>> Jacques
>>>>>
>>>>> 2010/8/16 Andy Bauer <andy.bauer at kitware.com>
>>>>>
>>>>> The other thing is that there was a problem when outputting images from
>>>>>> coprocessing.  Try updating to commit 2302777 to see if that's it.
>>>>>>
>>>>>> Andy
>>>>>>
>>>>>>
>>>>>> On Mon, Aug 16, 2010 at 9:07 AM, pat marion <pat.marion at kitware.com>wrote:
>>>>>>
>>>>>>> Hi Jacques,
>>>>>>>
>>>>>>> That error message seems to be unrelated to the lines containing
>>>>>>> 'None' that you copied in a previous email.  It looks to me like your code
>>>>>>> that calls the coprocessor is not inserting the vtk data set under the name
>>>>>>> 'test.foam'.  In the calling code you need something like:
>>>>>>>
>>>>>>> dataDesc->AddInput("test.foam")
>>>>>>>
>>>>>>> dataDesc->GetInputDescriptionByName("test.foam")->SetGrid(dataObject);
>>>>>>>
>>>>>>> Pat
>>>>>>>
>>>>>>>
>>>>>>> On Mon, Aug 16, 2010 at 8:35 AM, Jacques Papper <
>>>>>>> jacques.papper at gmail.com> wrote:
>>>>>>>
>>>>>>>> Hi Pat,
>>>>>>>>
>>>>>>>> This is what shows up for me :
>>>>>>>>
>>>>>>>>
>>>>>>>> datadescription.GetInputDescriptionByName(input_name).AllFieldsOff()
>>>>>>>> AttributeError: 'NoneType' object has no attribute 'AllFieldsOff'
>>>>>>>>
>>>>>>>>
>>>>>>>> Inside the Python script :
>>>>>>>>
>>>>>>>> def RequestDataDescription(datadescription):
>>>>>>>>     "Callback to populate the request for current timestep"
>>>>>>>>     timestep = datadescription.GetTimeStep()
>>>>>>>>
>>>>>>>>
>>>>>>>>     input_name = 'test.foam'
>>>>>>>>     if (timestep % 0.0004 == 0) :
>>>>>>>>
>>>>>>>> datadescription.GetInputDescriptionByName(input_name).AllFieldsOn()
>>>>>>>>
>>>>>>>> datadescription.GetInputDescriptionByName(input_name).GenerateMeshOn()
>>>>>>>>     else:
>>>>>>>>
>>>>>>>> datadescription.GetInputDescriptionByName(input_name).AllFieldsOff()
>>>>>>>>
>>>>>>>> datadescription.GetInputDescriptionByName(input_name).GenerateMeshOff()
>>>>>>>>
>>>>>>>>
>>>>>>>> Best,
>>>>>>>>
>>>>>>>>
>>>>>>>> Jacques
>>>>>>>>
>>>>>>>>
>>>>>>>> 2010/8/13 pat marion <pat.marion at kitware.com>
>>>>>>>>
>>>>>>>> Hi Jacques,
>>>>>>>>>
>>>>>>>>> My python trace files have similar lines, and I am able to run them
>>>>>>>>> without error.  I've only just tested with latest git master, not 3.8.  Can
>>>>>>>>> you copy the specific error message you're getting so I can take a look?
>>>>>>>>>
>>>>>>>>> Pat
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Fri, Aug 13, 2010 at 6:20 PM, Jacques Papper <
>>>>>>>>> jacques.papper at gmail.com> wrote:
>>>>>>>>>
>>>>>>>>>>    Sorry I meant None... not Null :
>>>>>>>>>>
>>>>>>>>>> For example :
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>     DataRepresentation1.SelectOrientationVectors = [None, '']
>>>>>>>>>>     DataRepresentation1.SelectedMapperIndex = 'Projected tetra'
>>>>>>>>>>     DataRepresentation1.CubeAxesTickLocation = 'Inside'
>>>>>>>>>>     DataRepresentation1.CubeAxesXAxisMinorTickVisibility = 1
>>>>>>>>>>     DataRepresentation1.CubeAxesYAxisVisibility = 1
>>>>>>>>>>     DataRepresentation1.SelectionPointLabelFontFamily = 'Arial'
>>>>>>>>>>     DataRepresentation1.Source.ShaftResolution = 6
>>>>>>>>>>     DataRepresentation1.CubeAxesFlyMode = 'Closest Triad'
>>>>>>>>>>     DataRepresentation1.SelectScaleArray = [None, '']
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> 2010/8/13 Andy Bauer <andy.bauer at kitware.com>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> On Fri, Aug 13, 2010 at 3:39 AM, Jacques Papper <
>>>>>>>>>>> jacques.papper at gmail.com> wrote:
>>>>>>>>>>>
>>>>>>>>>>>> I found the problem... I was inserting PointData instead of
>>>>>>>>>>>> CellData :) ,,,
>>>>>>>>>>>> Anyway it seems to be starting to work now.. 2 things though :
>>>>>>>>>>>>
>>>>>>>>>>>> -> Whenever I try to run directly a python script exported from
>>>>>>>>>>>> paraview (Through the Coprocessing export state) it won't run until I remove
>>>>>>>>>>>> all references to LIC and all lines containing Null...
>>>>>>>>>>>>
>>>>>>>>>>>> Now that you mention it, I didn't put in anything special for
>>>>>>>>>>> plugins.  Setting the plugins to auto-load may work but my guess is that it
>>>>>>>>>>> won't.  I'll look into the plugin thing and the best way to do it in the
>>>>>>>>>>> script generator.
>>>>>>>>>>>
>>>>>>>>>>> What do you mean by Null?  I don't know what that is so you'll
>>>>>>>>>>> have to explain it to me.
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>> For the LIC, is it because the python script doesn't
>>>>>>>>>>>> automatically load the LIC plugin ?
>>>>>>>>>>>> For the Null, do you have this problem as well ?
>>>>>>>>>>>>
>>>>>>>>>>>> -> Is the Live Source Data code already available in git ? I'd
>>>>>>>>>>>> love to test it out even if it is not robust yet.
>>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> We don't yet have permission to release the Live Source Data code
>>>>>>>>>>> yet but I think it should be available in the next couple of months.
>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> Best,
>>>>>>>>>>>> Jacques
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> 2010/8/12 Andy Bauer <andy.bauer at kitware.com>
>>>>>>>>>>>>
>>>>>>>>>>>>> You probably want to add the field data to the specific grids
>>>>>>>>>>>>> in the multiblock data set and not just the block itself.  As an example:
>>>>>>>>>>>>>
>>>>>>>>>>>>> create multiblock
>>>>>>>>>>>>> for each grid:
>>>>>>>>>>>>>   create grid points and cells
>>>>>>>>>>>>>   add point and cell data to grid
>>>>>>>>>>>>>   add grid to block
>>>>>>>>>>>>>
>>>>>>>>>>>>> Note that in parallel for coprocessing that the block must have
>>>>>>>>>>>>> the same structure on each process but whether a leaf node is empty or not
>>>>>>>>>>>>> can change.  For example:
>>>>>>>>>>>>>
>>>>>>>>>>>>> proc 0 mutliblock:
>>>>>>>>>>>>>   block 0 -> unstructured grid
>>>>>>>>>>>>>   block 1 -> uniform grid
>>>>>>>>>>>>>   block 2 -> NULL
>>>>>>>>>>>>> proc 1 multiblock:
>>>>>>>>>>>>>   block 0 -> NULL
>>>>>>>>>>>>>   block 1 -> NULL
>>>>>>>>>>>>>   block 2 -> unstructured grid
>>>>>>>>>>>>>
>>>>>>>>>>>>> I'm not sure of all of the details but what I do is that I set
>>>>>>>>>>>>> a block to only be non-NULL on a single process.
>>>>>>>>>>>>>
>>>>>>>>>>>>> Hope this helps,
>>>>>>>>>>>>> Andy
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>> On Thu, Aug 12, 2010 at 2:38 PM, Jacques Papper <
>>>>>>>>>>>>> jacques.papper at gmail.com> wrote:
>>>>>>>>>>>>>
>>>>>>>>>>>>>> Hi All,
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> I am making some more progress on this. Everything was working
>>>>>>>>>>>>>> until I tried adding field data... I now get an explosion :
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> #2  ?? in "/lib64/libc.so.6"
>>>>>>>>>>>>>> #3  vtkDataArrayTemplate<double>::ComputeScalarRange(int) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkCommon.so.pv3.9"
>>>>>>>>>>>>>> #4  vtkDataArray::ComputeRange(int) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkCommon.so.pv3.9"
>>>>>>>>>>>>>> #5  vtkPVArrayInformation::CopyFromObject(vtkObject*) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
>>>>>>>>>>>>>> #6
>>>>>>>>>>>>>> vtkPVDataSetAttributesInformation::CopyFromDataSetAttributes(vtkDataSetAttributes*)
>>>>>>>>>>>>>> in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
>>>>>>>>>>>>>> #7  vtkPVDataInformation::CopyFromDataSet(vtkDataSet*) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
>>>>>>>>>>>>>> #8  vtkPVDataInformation::CopyFromObject(vtkObject*) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
>>>>>>>>>>>>>> #9  vtkPVCompositeDataInformation::CopyFromObject(vtkObject*)
>>>>>>>>>>>>>> in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
>>>>>>>>>>>>>> #10
>>>>>>>>>>>>>> vtkPVDataInformation::CopyFromCompositeDataSet(vtkCompositeDataSet*) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
>>>>>>>>>>>>>> #11  vtkPVDataInformation::CopyFromObject(vtkObject*) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
>>>>>>>>>>>>>> #12  vtkSelfConnection::GatherInformation(unsigned int,
>>>>>>>>>>>>>> vtkPVInformation*, vtkClientServerID) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
>>>>>>>>>>>>>> #13
>>>>>>>>>>>>>> vtkSynchronousMPISelfConnection::GatherInformation(unsigned int,
>>>>>>>>>>>>>> vtkPVInformation*, vtkClientServerID) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerCommon.so"
>>>>>>>>>>>>>> #14  vtkSMOutputPort::GatherDataInformation() in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerManager.so"
>>>>>>>>>>>>>> #15  vtkSMOutputPort::GetDataInformation() in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerManager.so"
>>>>>>>>>>>>>> #16  vtkSMDataTypeDomain::IsInDomain(vtkSMSourceProxy*, int)
>>>>>>>>>>>>>> in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerManager.so"
>>>>>>>>>>>>>> #17  vtkSMDataTypeDomain::IsInDomain(vtkSMProperty*) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerManager.so"
>>>>>>>>>>>>>> #18  vtkSMProperty::IsInDomains() in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerManager.so"
>>>>>>>>>>>>>> #19
>>>>>>>>>>>>>> vtkSMRenderViewProxy::CreateDefaultRepresentation(vtkSMProxy*, int) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVServerManager.so"
>>>>>>>>>>>>>> #20  ?? at vtkSMRenderViewProxyPython.cxx:0
>>>>>>>>>>>>>> #21  PyEval_EvalFrameEx in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #22  PyEval_EvalCodeEx in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #23  PyEval_EvalFrameEx in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #24  PyEval_EvalCodeEx in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #25  PyEval_EvalFrameEx in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #26  PyEval_EvalCodeEx in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #27  PyEval_EvalFrameEx in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #28  PyEval_EvalCodeEx in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #29  PyEval_EvalFrameEx in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #30  PyEval_EvalFrameEx in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #31  PyEval_EvalCodeEx in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #32  PyEval_EvalCode in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #33  ?? in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #34  PyRun_StringFlags in "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #35  PyRun_SimpleStringFlags in
>>>>>>>>>>>>>> "/usr/lib64/libpython2.6.so.1.0"
>>>>>>>>>>>>>> #36  vtkPVPythonInterpretor::RunSimpleString(char const*) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkPVPythonInterpretor.so"
>>>>>>>>>>>>>> #37
>>>>>>>>>>>>>> vtkCPPythonScriptPipeline::CoProcess(vtkCPDataDescription*) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkCoProcessor.so"
>>>>>>>>>>>>>> #38  vtkCPProcessor::CoProcess(vtkCPDataDescription*) in
>>>>>>>>>>>>>> "/users/boreas01/jacques/PARAVIEW/ParaView-bin/bin/libvtkCoProcessor.so"
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Any ideas where this can be coming from ?
>>>>>>>>>>>>>> When I do not add the field data it just complains about
>>>>>>>>>>>>>> missing fields so I know that the code is not completely wrong...
>>>>>>>>>>>>>> I am adding the data to the first block of my
>>>>>>>>>>>>>> multiblockdataset... Is that a problem ?
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Thanks !
>>>>>>>>>>>>>> Jacques
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> 2010/8/12 Jacques Papper <jacques.papper at gmail.com>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Hi Andy, Pat,
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> I have the skeleton of the coprocessing working (for the
>>>>>>>>>>>>>>> moment I just tried by populating the point field - no cells, no data).
>>>>>>>>>>>>>>> I am trying to figure out the most efficient way of
>>>>>>>>>>>>>>> transferring my data structure to VTK with minimal copying and looping...
>>>>>>>>>>>>>>> I would like to use the vtkUnstructuredGrid::SetCells()
>>>>>>>>>>>>>>> method, which means I probably need to use the vtkCellArray::SetCells()
>>>>>>>>>>>>>>> method, and therefore the vtkIdType::SetArray() method...
>>>>>>>>>>>>>>> The problem is I am not sure what the format of the array
>>>>>>>>>>>>>>> should be ?
>>>>>>>>>>>>>>> I am dealing with unstructured polyhedra cells (most of which
>>>>>>>>>>>>>>> are known types like hex prism etc..).
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Of course, I could also use the InsertNextCell() methods, but
>>>>>>>>>>>>>>> I was wondering if there wasn't a more efficient way...
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Best,
>>>>>>>>>>>>>>> Jacques
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> 2010/8/9 Andy Bauer <andy.bauer at kitware.com>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Hi Jacques,
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> I put an example on the wiki based on Pat's example for
>>>>>>>>>>>>>>>> running a coprocessing.  It's at:
>>>>>>>>>>>>>>>> http://paraview.org/Wiki/Coprocessing_example#Python_Scripts
>>>>>>>>>>>>>>>> Let me know if you have any problems with it so that I can
>>>>>>>>>>>>>>>> fix any errors that I have.
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Also, I just made changes so that the coprocessing library
>>>>>>>>>>>>>>>> gets installed properly.
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>>>> Andy
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> On Wed, Aug 4, 2010 at 7:21 PM, Jacques Papper <
>>>>>>>>>>>>>>>> jacques.papper at gmail.com> wrote:
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Thanks Pat,
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> I am trying to compile my first attempt.
>>>>>>>>>>>>>>>>> Although I do find the headers I need, I am having some
>>>>>>>>>>>>>>>>> trouble identifying which libraries I need to link to.
>>>>>>>>>>>>>>>>> For example, I assume that I should be linking to
>>>>>>>>>>>>>>>>> libvtkCoProcessor.so ... Although I find this library in the compilation
>>>>>>>>>>>>>>>>> folder (under bin), I do not find it in the make installed version under
>>>>>>>>>>>>>>>>> lib...
>>>>>>>>>>>>>>>>> Is this normal ?
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Best,
>>>>>>>>>>>>>>>>> Jacques
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> 2010/8/4 pat marion <pat.marion at kitware.com>
>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>> Below is a simple cxx program I have used for testing.  It
>>>>>>>>>>>>>>>>>> creates one sphere per process, positions the sphere as a function of the
>>>>>>>>>>>>>>>>>> process id, and the spheres grow/shrink over time.
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> Pat
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> #include "vtkCPProcessor.h"
>>>>>>>>>>>>>>>>>> #include "vtkCPPythonScriptPipeline.h"
>>>>>>>>>>>>>>>>>> #include "vtkMultiProcessController.h"
>>>>>>>>>>>>>>>>>> #include "vtkXMLUnstructuredGridReader.h"
>>>>>>>>>>>>>>>>>> #include "vtkUnstructuredGrid.h"
>>>>>>>>>>>>>>>>>> #include "vtkCPDataDescription.h"
>>>>>>>>>>>>>>>>>> #include "vtkCPInputDataDescription.h"
>>>>>>>>>>>>>>>>>> #include "vtkSmartPointer.h"
>>>>>>>>>>>>>>>>>> #include "vtkPolyData.h"
>>>>>>>>>>>>>>>>>> #include "vtkSphereSource.h"
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> #include <stdio.h>
>>>>>>>>>>>>>>>>>> #include <string>
>>>>>>>>>>>>>>>>>> #include <sstream>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> static unsigned int procId;
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> void myprint(const std::string& str)
>>>>>>>>>>>>>>>>>> {
>>>>>>>>>>>>>>>>>>   printf("driver (%u): %s\n", procId, str.c_str());
>>>>>>>>>>>>>>>>>> }
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> class DataGenerator {
>>>>>>>>>>>>>>>>>> public:
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   DataGenerator()
>>>>>>>>>>>>>>>>>>     {
>>>>>>>>>>>>>>>>>>     this->Sphere =
>>>>>>>>>>>>>>>>>> vtkSmartPointer<vtkSphereSource>::New();
>>>>>>>>>>>>>>>>>>     this->Sphere->SetThetaResolution(30);
>>>>>>>>>>>>>>>>>>     this->Sphere->SetPhiResolution(30);
>>>>>>>>>>>>>>>>>>     this->Sphere->SetCenter(procId*4.0, 0, 0);
>>>>>>>>>>>>>>>>>>     this->Index = 0;
>>>>>>>>>>>>>>>>>>     }
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   vtkSmartPointer<vtkPolyData> GetNext()
>>>>>>>>>>>>>>>>>>     {
>>>>>>>>>>>>>>>>>>     double radius = fabs(sin(0.1 * this->Index));
>>>>>>>>>>>>>>>>>>     this->Index++;
>>>>>>>>>>>>>>>>>>     this->Sphere->SetRadius(1.0 + radius);
>>>>>>>>>>>>>>>>>>     this->Sphere->Update();
>>>>>>>>>>>>>>>>>>     vtkSmartPointer<vtkPolyData> ret =
>>>>>>>>>>>>>>>>>> vtkSmartPointer<vtkPolyData>::New();
>>>>>>>>>>>>>>>>>>     ret->DeepCopy(this->Sphere->GetOutput());
>>>>>>>>>>>>>>>>>>     return ret;
>>>>>>>>>>>>>>>>>>     }
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> protected:
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   int Index;
>>>>>>>>>>>>>>>>>>   vtkSmartPointer<vtkSphereSource> Sphere;
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> };
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> int main(int argc, char* argv[])
>>>>>>>>>>>>>>>>>> {
>>>>>>>>>>>>>>>>>>   if (argc < 3)
>>>>>>>>>>>>>>>>>>     {
>>>>>>>>>>>>>>>>>>     printf("Usage: %s <cp python file> <number of
>>>>>>>>>>>>>>>>>> steps>\n", argv[0]);
>>>>>>>>>>>>>>>>>>     return 1;
>>>>>>>>>>>>>>>>>>     }
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   std::string cpPythonFile = argv[1];
>>>>>>>>>>>>>>>>>>   int nSteps = atoi(argv[2]);
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   myprint("starting coprocessor");
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   vtkCPProcessor* processor = vtkCPProcessor::New();
>>>>>>>>>>>>>>>>>>   processor->Initialize();
>>>>>>>>>>>>>>>>>>   vtkCPPythonScriptPipeline* pipeline =
>>>>>>>>>>>>>>>>>> vtkCPPythonScriptPipeline::New();
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   // mpi was initialized when we called
>>>>>>>>>>>>>>>>>> vtkCPPythonScriptPipeline::New()
>>>>>>>>>>>>>>>>>>   procId =
>>>>>>>>>>>>>>>>>> vtkMultiProcessController::GetGlobalController()->GetLocalProcessId();
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   // read the coprocessing python file
>>>>>>>>>>>>>>>>>>   myprint("loading pipeline python file: " +
>>>>>>>>>>>>>>>>>> cpPythonFile);
>>>>>>>>>>>>>>>>>>   int success =
>>>>>>>>>>>>>>>>>> pipeline->Initialize(cpPythonFile.c_str());
>>>>>>>>>>>>>>>>>>   if (!success)
>>>>>>>>>>>>>>>>>>     {
>>>>>>>>>>>>>>>>>>     myprint("aborting");
>>>>>>>>>>>>>>>>>>     return 1;
>>>>>>>>>>>>>>>>>>     }
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   processor->AddPipeline(pipeline);
>>>>>>>>>>>>>>>>>>   pipeline->Delete();
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   if (nSteps == 0)
>>>>>>>>>>>>>>>>>>     {
>>>>>>>>>>>>>>>>>>     return 0;
>>>>>>>>>>>>>>>>>>     }
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   // create a data source
>>>>>>>>>>>>>>>>>>   DataGenerator generator;
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   // do coprocessing
>>>>>>>>>>>>>>>>>>   double tStart = 0.0;
>>>>>>>>>>>>>>>>>>   double tEnd = 1.0;
>>>>>>>>>>>>>>>>>>   double stepSize = (tEnd - tStart)/nSteps;
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   vtkCPDataDescription* dataDesc =
>>>>>>>>>>>>>>>>>> vtkCPDataDescription::New();
>>>>>>>>>>>>>>>>>>   dataDesc->AddInput("input");
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   for (int i = 0; i < nSteps; ++i)
>>>>>>>>>>>>>>>>>>     {
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>     double currentTime = tStart + stepSize*i;
>>>>>>>>>>>>>>>>>>     std::stringstream timeStr;
>>>>>>>>>>>>>>>>>>     timeStr << "time(" << i << ", " << currentTime << ")";
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>     dataDesc->SetTimeData(currentTime, i);
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>     myprint("call RequestDataDescription, " +
>>>>>>>>>>>>>>>>>> timeStr.str());
>>>>>>>>>>>>>>>>>>     int do_coprocessing =
>>>>>>>>>>>>>>>>>> processor->RequestDataDescription(dataDesc);
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>     if (do_coprocessing)
>>>>>>>>>>>>>>>>>>       {
>>>>>>>>>>>>>>>>>>       myprint("calling CoProcess, " + timeStr.str());
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>       vtkSmartPointer<vtkDataObject> dataObject =
>>>>>>>>>>>>>>>>>>         generator.GetNext();
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> dataDesc->GetInputDescriptionByName("input")->SetGrid(dataObject);
>>>>>>>>>>>>>>>>>>       processor->CoProcess(dataDesc);
>>>>>>>>>>>>>>>>>>       }
>>>>>>>>>>>>>>>>>>     }
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   myprint("finalizing");
>>>>>>>>>>>>>>>>>>   dataDesc->Delete();
>>>>>>>>>>>>>>>>>>   processor->Finalize();
>>>>>>>>>>>>>>>>>>   processor->Delete();
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>   return 0;
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> }
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>> On Wed, Aug 4, 2010 at 12:14 PM, Andy Bauer <
>>>>>>>>>>>>>>>>>> andy.bauer at kitware.com> wrote:
>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Hi Jacques,
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> There is a polyhedra cell type in VTK now --
>>>>>>>>>>>>>>>>>>> http://www.vtk.org/doc/nightly/html/classvtkPolyhedron.html
>>>>>>>>>>>>>>>>>>> As far as I know it works with all of the proper filters
>>>>>>>>>>>>>>>>>>> but since I haven't tried it yet I won't promise that.  The good news is
>>>>>>>>>>>>>>>>>>> that Will Schroeder had a high interest in it and probably worked on some of
>>>>>>>>>>>>>>>>>>> it so I'd assume that it's working quite well right now.
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> As for Phasta, it does run in parallel (as props to their
>>>>>>>>>>>>>>>>>>> developers it was a finalist for the 2009 Gordon Bell prize).  The grid is
>>>>>>>>>>>>>>>>>>> already partitioned and each process runs the adaptor and creates an
>>>>>>>>>>>>>>>>>>> unstructured grid from its portion of the partitioned mesh.  Thus, there
>>>>>>>>>>>>>>>>>>> isn't any need for mpi calls in the adaptor code.  If you had ghost cell
>>>>>>>>>>>>>>>>>>> information in your partitioned mesh and wanted to get fancy you should be
>>>>>>>>>>>>>>>>>>> able to add that to your partitioned grid to make some of the filters faster
>>>>>>>>>>>>>>>>>>> but I haven't tried that.
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Andy
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> On Wed, Aug 4, 2010 at 11:58 AM, Jacques Papper <
>>>>>>>>>>>>>>>>>>> jacques.papper at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Thanks a lot Andy, Takuya,
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> I'm using the PhastaAdaptor, and the FortranAdaptorAPI
>>>>>>>>>>>>>>>>>>>> as a guide for the moment.
>>>>>>>>>>>>>>>>>>>> I know there were talks of getting POLYHEDRAL cell
>>>>>>>>>>>>>>>>>>>> support in VTK. Do you know if this is there yet?
>>>>>>>>>>>>>>>>>>>> My dataset is multiregion unstructured polyhedral mesh
>>>>>>>>>>>>>>>>>>>> domain decomposed amongst each processor.
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Is the Phasta code parallelized ? If so, I do not see
>>>>>>>>>>>>>>>>>>>> any MPI statements in the adaptor code ?
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Jacques
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> 2010/8/4 Andy Bauer <andy.bauer at kitware.com>
>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>> Hi Jacques,
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> What type of data set do you have?  Even though the
>>>>>>>>>>>>>>>>>>>>> PHASTA adaptor (
>>>>>>>>>>>>>>>>>>>>> ParaView/CoProcessing/Adaptors/FortranAdaptors/PhastaAdaptor) is for fortran
>>>>>>>>>>>>>>>>>>>>> code it may give you an idea.  Also stepping through the example in
>>>>>>>>>>>>>>>>>>>>> ParaView/CoProcessing/CoProcessor/Testing/Cxx/PythonScriptCoProcessingExample.cxx
>>>>>>>>>>>>>>>>>>>>> may help as well.
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> I'll spend some time this week putting up a skeleton of
>>>>>>>>>>>>>>>>>>>>> a simulation code on the coprocessing wiki that should hopefully be easier
>>>>>>>>>>>>>>>>>>>>> to follow.  I'll let you know when it's done.
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> Andy
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>> On Wed, Aug 4, 2010 at 8:02 AM, Jacques Papper <
>>>>>>>>>>>>>>>>>>>>> jacques.papper at gmail.com> wrote:
>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>  Hi All,
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> Sorry for my last post, I figured out that I had
>>>>>>>>>>>>>>>>>>>>>> wrongly set my PYTHONPATH..
>>>>>>>>>>>>>>>>>>>>>> All the tests work ok now. Still interested in
>>>>>>>>>>>>>>>>>>>>>> CoProcessing adaptors examples though :)
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> Thanks
>>>>>>>>>>>>>>>>>>>>>> Jacques
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> 2010/8/4 Jacques Papper <jacques.papper at gmail.com>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> Hi All,
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> I'm starting to look into the ParaView CoProcessing
>>>>>>>>>>>>>>>>>>>>>>> libraries.
>>>>>>>>>>>>>>>>>>>>>>> I just pulled from git today, and compiled it all up
>>>>>>>>>>>>>>>>>>>>>>> following the guidelines in :
>>>>>>>>>>>>>>>>>>>>>>> http://www.paraview.org/Wiki/CoProcessing
>>>>>>>>>>>>>>>>>>>>>>> I didn't find :
>>>>>>>>>>>>>>>>>>>>>>> *BUILD_PYTHON_COPROCESSING_ADAPTOR
>>>>>>>>>>>>>>>>>>>>>>> *but instead :
>>>>>>>>>>>>>>>>>>>>>>> PARAVIEW_BUILD_PLUGIN_CoProcessingScriptGenerator
>>>>>>>>>>>>>>>>>>>>>>> anyway the compilation went through without any
>>>>>>>>>>>>>>>>>>>>>>> issues.
>>>>>>>>>>>>>>>>>>>>>>> I then tried :
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> ctest -R CoProcessing
>>>>>>>>>>>>>>>>>>>>>>> Test project
>>>>>>>>>>>>>>>>>>>>>>> /users/boreas01/jacques/PARAVIEW/ParaView-bin
>>>>>>>>>>>>>>>>>>>>>>>     Start 491: CoProcessingTestPythonScript
>>>>>>>>>>>>>>>>>>>>>>> 1/3 Test #491: CoProcessingTestPythonScript
>>>>>>>>>>>>>>>>>>>>>>> ...........   Passed    0.45 sec
>>>>>>>>>>>>>>>>>>>>>>>     Start 492: CoProcessingPythonScriptGridPlot
>>>>>>>>>>>>>>>>>>>>>>> 2/3 Test #492: CoProcessingPythonScriptGridPlot
>>>>>>>>>>>>>>>>>>>>>>> .......***Failed    0.09 sec
>>>>>>>>>>>>>>>>>>>>>>>     Start 493: CoProcessingPythonScriptPressurePlot
>>>>>>>>>>>>>>>>>>>>>>> 3/3 Test #493: CoProcessingPythonScriptPressurePlot
>>>>>>>>>>>>>>>>>>>>>>> ...***Failed    0.09 sec
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> 33% tests passed, 2 tests failed out of 3
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> Total Test time (real) =   0.68 sec
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> The following tests FAILED:
>>>>>>>>>>>>>>>>>>>>>>>         492 - CoProcessingPythonScriptGridPlot
>>>>>>>>>>>>>>>>>>>>>>> (Failed)
>>>>>>>>>>>>>>>>>>>>>>>         493 - CoProcessingPythonScriptPressurePlot
>>>>>>>>>>>>>>>>>>>>>>> (Failed)
>>>>>>>>>>>>>>>>>>>>>>> Errors while running CTest
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> Is this a problem in my current installation or on
>>>>>>>>>>>>>>>>>>>>>>> the master branch ?
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> Finally, I would like to start writing an adaptor to
>>>>>>>>>>>>>>>>>>>>>>> a C++ parallelised code. Can you tell me what is the closest code I can
>>>>>>>>>>>>>>>>>>>>>>> inspire myself from ?
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> Thank you,
>>>>>>>>>>>>>>>>>>>>>>> Jacques
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>> PS sorry mixed up the subjects..
>>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> _______________________________________________
>>>>>>>>>>>>>>>>>>>>>> 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
>>>>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>>>>> Follow this link to subscribe/unsubscribe:
>>>>>>>>>>>>>>>>>>>>>> http://www.paraview.org/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
>>>>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>>>> Follow this link to subscribe/unsubscribe:
>>>>>>>>>>>>>>>>>>> http://www.paraview.org/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
>>
>> Follow this link to subscribe/unsubscribe:
>> http://www.paraview.org/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
>
> Follow this link to subscribe/unsubscribe:
> http://www.paraview.org/mailman/listinfo/paraview
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.paraview.org/pipermail/paraview/attachments/20100817/776c3e57/attachment-0001.htm>


More information about the ParaView mailing list