[Paraview] CoProcessing
Andy Bauer
andy.bauer at kitware.com
Mon Aug 16 10:51:41 EDT 2010
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
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.paraview.org/pipermail/paraview/attachments/20100816/1a72695f/attachment-0001.htm>
More information about the ParaView
mailing list