[Paraview] Program to write .VTK
Hello everybody, I'm searching for a program to write some .DAT files in .VTK format. (for * Paraview* visualization) These .DATs contains grid structure(bathymetry), lat, long, scalars and vectors. I'm writing them *manually* in .VTK, like these: # vtk DataFile Version 3.0 vtk output ASCII DATASET STRUCTURED_POINTS *DIMENSIONS 387 312 1 * SPACING 1 1 1 ORIGIN 0 0 0 *POINT_DATA 120744 * SCALARS batimetria float *LOOKUP_TABLE default 0.000e+000 0.000e+000 0.000e+000 0.000e+000 0.000e+000... 0.000e+000 0.000e+000 0.000e+000 0.000e+000 0.000e+000...* ... The DAT files holds the grid (bathymetry), then I spent lots of time copying and pasting the grid from .DAT to *LOOKUP_TABLE default*, and writing the *dimensions and point data* of each file... Since I'm working with large time data-steps files, I would like to known if anybody use or recommend a program to speed up this process. Thanks a lot, Renan. -- Renan Leser de Medeiros Laboratory of Computational Methods in Engineering Federal University of Rio de Janeiro - COPPE P.O.Box 68552 - ZIP 21949-900 Rio de Janeiro, RJ, Brasil ___ 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
[Paraview] [newbie] Animating time-series MFIX data
Hi all, Please excuse my newbie question - I'm a new ParaVIEW user this week. I have a time-series dataset generated from the MFIX CFD code that I want to visualise. I have no problem loading and displaying the data. I can generate an isosurface, for example, and animate such that the isosurface value changes over time. But I cannot seem to animate over the time series data (e.g. for a constant isosurface value). For example if I change the Time Step slider in the Object Inspector and apply, the view updates to reflect the data at that time step. However if I add an animation track (Cell Arrays) the display does not update during the anim. Hopefully this is a simple newbie question :) Thanks in advance for your help, Regards, Mark ___ 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
Re: [Paraview] [newbie] Animating time-series MFIX data
Actually, this seems to be a bug in the reader. Creating the contour and then hitting play should just work. I tested and it does not. We will look into it. -berk On Thu, Feb 12, 2009 at 5:24 AM, Mark S. Bentley m...@lunartech.org wrote: Hi all, Please excuse my newbie question - I'm a new ParaVIEW user this week. I have a time-series dataset generated from the MFIX CFD code that I want to visualise. I have no problem loading and displaying the data. I can generate an isosurface, for example, and animate such that the isosurface value changes over time. But I cannot seem to animate over the time series data (e.g. for a constant isosurface value). For example if I change the Time Step slider in the Object Inspector and apply, the view updates to reflect the data at that time step. However if I add an animation track (Cell Arrays) the display does not update during the anim. Hopefully this is a simple newbie question :) Thanks in advance for your help, Regards, Mark ___ 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
Re: [Paraview] particle tracking
Berk Oops. I didn't notice that it was a "steady state" flow field. The TemporalStremTRacer expects unsteady flows, however it would be quite simple to modify it to use a single step. (I think easier than taking streamlines and animating particles along them - though thinking about it, the streamtracer code does have all that conversion of cell units to time which would be useful). One could subclass the temporaltreamtracer, suppress the time requests and interpolate over a single step, then output time dependent data where the particles animate along the path (- how many steps would one want?) - clicking the play button would just advect them along. (Just thinking out loud). JB I don't think this answers Pei's question. To my knowledge the temporal stream tracer (i.e. particle tracker) works only for time-dependent data. Am I wrong? To animate particles in a steady-state flow field, I'd think that you would generate streamlines and then somehow animate particles along those. I can't think of anything that would do this out-of-box but writing a filter that does this would be pretty easy. Is this something more than a few users would want? -berk On Thu, Feb 12, 2009 at 12:22 AM, David E DeMarle dave.dema...@kitware.com wrote: Hi Pei-Ying, You may want to try ParaView Meshless (https://twiki.cscs.ch/twiki/bin/view/ParaViewMeshless). It is a version of ParaView with the cutting edge of John Biddiscombe's particle tracking work in it. Some of those features will be/have been integrated into the main paraview code, but I can not vouch for when that will happen. cheers, Dave DeMarle 2008/12/25 Pei-Ying Hsieh phsieh2...@yahoo.com: Dear PVers: I have a steady state flow field. I am wondering if PV can do particle tracking animation. what I am looking for is similar to pathlines, but, instead of showing the lines, I would like to do an animation of seed particles "flowing" through the flow field. Is this possible? Thanks! and wish everyone here a Happy New Year! Pei ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview -- David E DeMarle Kitware, Inc. RD Engineer 28 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-371-3971 x109 ___ 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 -- John Biddiscombe,email:biddisco @ cscs.ch http://www.cscs.ch/ CSCS, Swiss National Supercomputing Centre | Tel: +41 (91) 610.82.07 Via Cantonale, 6928 Manno, Switzerland | Fax: +41 (91) 610.82.82 ___ 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
Re: [Paraview] particle tracking
Berk Geveci wrote: To animate particles in a steady-state flow field, I'd think that you would generate streamlines and then somehow animate particles along those. this is exactly the technique I use. Generate streamlines. Then iso-contour the streamline object with the scalar field IntegrationTime. Use a single threshold. Then animate the threshold value. Use Mode=Sequence, get a ramp between minimum time and maximun time. Use Glyphs (small spheres or arrows) attached to the iso-valued contours and they will animate their position along the streamlines. Jean Swiss national Supercomputing Center ___ 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
Re: [Paraview] particle tracking
Interesting. I guess the biggest challenge of doing this in our pipeline would be finding the time range. Ideally, the integration should continue until all particles leave the domain or get stuck in a stagnant region. That would require integrating them first and then looking at the largest integration time. While I was writing this, Jean's response showed up. It is pretty ingenious :-) This is why I love ParaView. It allows you to do things in ways developers did not foresee. -berk On Thu, Feb 12, 2009 at 8:42 AM, John Biddiscombe biddi...@cscs.ch wrote: Berk Oops. I didn't notice that it was a steady state flow field. The TemporalStremTRacer expects unsteady flows, however it would be quite simple to modify it to use a single step. (I think easier than taking streamlines and animating particles along them - though thinking about it, the streamtracer code does have all that conversion of cell units to time which would be useful). One could subclass the temporaltreamtracer, suppress the time requests and interpolate over a single step, then output time dependent data where the particles animate along the path (- how many steps would one want?) - clicking the play button would just advect them along. (Just thinking out loud). JB I don't think this answers Pei's question. To my knowledge the temporal stream tracer (i.e. particle tracker) works only for time-dependent data. Am I wrong? To animate particles in a steady-state flow field, I'd think that you would generate streamlines and then somehow animate particles along those. I can't think of anything that would do this out-of-box but writing a filter that does this would be pretty easy. Is this something more than a few users would want? -berk On Thu, Feb 12, 2009 at 12:22 AM, David E DeMarle dave.dema...@kitware.com wrote: Hi Pei-Ying, You may want to try ParaView Meshless (https://twiki.cscs.ch/twiki/bin/view/ParaViewMeshless). It is a version of ParaView with the cutting edge of John Biddiscombe's particle tracking work in it. Some of those features will be/have been integrated into the main paraview code, but I can not vouch for when that will happen. cheers, Dave DeMarle 2008/12/25 Pei-Ying Hsieh phsieh2...@yahoo.com: Dear PVers: I have a steady state flow field. I am wondering if PV can do particle tracking animation. what I am looking for is similar to pathlines, but, instead of showing the lines, I would like to do an animation of seed particles flowing through the flow field. Is this possible? Thanks! and wish everyone here a Happy New Year! Pei ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview -- David E DeMarle Kitware, Inc. RD Engineer 28 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-371-3971 x109 ___ 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 -- John Biddiscombe,email:biddisco @ cscs.ch http://www.cscs.ch/ CSCS, Swiss National Supercomputing Centre | Tel: +41 (91) 610.82.07 Via Cantonale, 6928 Manno, Switzerland | Fax: +41 (91) 610.82.82 ___ 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
[Paraview] Why is PKdTree used for IMAGE Volume Rendering?
I just tried to use ParaView from Tuesday's CVS (The network was down last night.) on both CPU cores to volume render a large (1/4 the size of available RAM) image. As soon as I selected volume render for the representation I saw the status bar show PKdTree and moments later all of the system memory was consumed and the system started swapping. I had to kill paraview and the pvserver processes from the console to end the swap-death. I reran paraview with only one CPU core and proceeded as normal. My question is: why is PKdTree needed at all for image data? For unstructured data it of course makes sense but for image data it just wastes a lot of RAM. I haven't tried parallel image volume rendering in a while so I don't know when this showed up. signature.asc Description: This is a digitally signed message part ___ 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
Re: [Paraview] Why is PKdTree used for IMAGE Volume Rendering?
Do you have any unstructured data being rendering in the scene? KdTree is used to distribute that unstructured data in accordance with the partitions provided by the structured data. Utkarsh On Thu, Feb 12, 2009 at 2:09 PM, Kevin H. Hobbs hob...@ohiou.edu wrote: I just tried to use ParaView from Tuesday's CVS (The network was down last night.) on both CPU cores to volume render a large (1/4 the size of available RAM) image. As soon as I selected volume render for the representation I saw the status bar show PKdTree and moments later all of the system memory was consumed and the system started swapping. I had to kill paraview and the pvserver processes from the console to end the swap-death. I reran paraview with only one CPU core and proceeded as normal. My question is: why is PKdTree needed at all for image data? For unstructured data it of course makes sense but for image data it just wastes a lot of RAM. I haven't tried parallel image volume rendering in a while so I don't know when this showed up. ___ 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
Re: [Paraview] Why is PKdTree used for IMAGE Volume Rendering?
On Thu, 2009-02-12 at 14:26 -0500, Utkarsh Ayachit wrote: Do you have any unstructured data being rendering in the scene? No, just one big .vti image. signature.asc Description: This is a digitally signed message part ___ 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
Re: [Paraview] Why is PKdTree used for IMAGE Volume Rendering?
Kevin, The KdTree is still built (using the partitions from the structured dataset set itself -- hence isn't a very compute intensive task) just in case it's needed to partition any unstructured data, hence you're seeing the PKdTree progress. But I don't think that's where it is hanging. Is the progress bar saying PKtree:100 ? It may actually have frozen in the volume rendering code. Is it possible to attach a debugger to any of the server processes and see where it's stuck (provide the stack trace if possible)? Utkarsh On Thu, Feb 12, 2009 at 2:37 PM, Kevin H. Hobbs hob...@ohiou.edu wrote: On Thu, 2009-02-12 at 14:26 -0500, Utkarsh Ayachit wrote: Do you have any unstructured data being rendering in the scene? No, just one big .vti image. ___ 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
Re: [Paraview] Why is PKdTree used for IMAGE Volume Rendering?
Ah and there's the bad news. In case of image data for the if condition at vtkPKdTree.cxx:435, I'd expect the condition to succeed and go to this-ProcessUserDefinedCuts(..) however, it seems like that condition is failing and it's going to this-MultiProcessBuildLocator(..) which is used for unstructured datasets. Are you able to reproduce the behavior with a simpler dataset (say wavelet) or something that;s already in ParaViewData, or some dataset that I can try it with? Also give me all the steps to reproduce the problem (since I am not able to reproduce it with 5 procs using Wavelet source). Utkarsh On Thu, Feb 12, 2009 at 3:24 PM, Kevin H. Hobbs hob...@ohiou.edu wrote: On Thu, 2009-02-12 at 14:56 -0500, Utkarsh Ayachit wrote: Kevin, The KdTree is still built (using the partitions from the structured dataset set itself -- hence isn't a very compute intensive task) just in case it's needed to partition any unstructured data, hence you're seeing the PKdTree progress. But I don't think that's where it is hanging. Is the progress bar saying PKtree:100 ? No, it only made it to 58 It may actually have frozen in the volume rendering code. Is it possible to attach a debugger to any of the server processes and see where it's stuck (provide the stack trace if possible)? Paraview stops completely as soon as I attach but here is the first bt before it really hung up: (gdb) bt #0 vtkDataArrayTemplatedouble::DataChanged (this=0x2e89fa0) at /home/kevin/kitware/ParaView3/VTK/Common/vtkDataArrayTemplate.txx:1202 #1 0x7f035a9cc221 in vtkPoints::SetPoint () at /home/kevin/kitware/ParaView3/VTK/Common/vtkPoints.h:126 #2 vtkImageData::GetCell (this=0x2e897d0, cellId=48799904) at /home/kevin/kitware/ParaView3/VTK/Filtering/vtkImageData.cxx:392 #3 0x7f035aa00bed in vtkKdTree::ComputeCellCenters (this=0x28853d0, set=0x0) at /home/kevin/kitware/ParaView3/VTK/Filtering/vtkKdTree.cxx:677 #4 0x7f035eee71eb in vtkPKdTree::MultiProcessBuildLocator (this=0x28853d0, volBounds=0x2e91140) at /home/kevin/kitware/ParaView3/VTK/Parallel/vtkPKdTree.cxx:498 #5 0x7f035eee7683 in vtkPKdTree::BuildLocator (this=0x28853d0) at /home/kevin/kitware/ParaView3/VTK/Parallel/vtkPKdTree.cxx:441 #6 0x00769725 in vtkKdTreeManager::Update (this=0x2877c50) at /home/kevin/kitware/ParaView3/Servers/Filters/vtkKdTreeManager.cxx:143 #7 0x078fc403 in vtkKdTreeManagerCommand (arlu=0x24cfac0, ob=0x2877c50, method=value optimized out, m...@0x7fff6aeabc10, resultstre...@0x24ceb40) at /home/kevin/kitware/ParaView3_Build/Servers/Filters/vtkKdTreeManagerClientServer.cxx:123 #8 0x7f035f40e36f in vtkClientServerInterpreter::ProcessCommandInvoke (this=0x24cfac0, c...@0x7fff6aeac080, midx=0) at /home/kevin/kitware/ParaView3/Utilities/VTKClientServer/vtkClientServerInterpreter.cxx:384 #9 0x7f035f41058e in vtkClientServerInterpreter::ProcessOneMessage (this=0x24cfac0, c...@0x7fff6aeac080, message=0) at /home/kevin/kitware/ParaView3/Utilities/VTKClientServer/vtkClientServerInterpreter.cxx:214 #10 0x7f035f4107fd in vtkClientServerInterpreter::ProcessStream (this=0x24cfac0, c...@0x7fff6aeac080) at /home/kevin/kitware/ParaView3/Utilities/VTKClientServer/vtkClientServerInterpreter.cxx:183 #11 0x001ae220 in vtkSelfConnection::ProcessStreamLocally (this=value optimized out, stre...@0x7fff6aeac080) at /home/kevin/kitware/ParaView3/Servers/Common/vtkSelfConnection.cxx:120 #12 0x0015e7ca in vtkMPISelfConnection::SendStreamToServerNodeInternal (this=0x2985a40, remoteId=-1, stre...@0x7fff6aeac080) at /home/kevin/kitware/ParaView3/Servers/Common/vtkMPISelfConnection.cxx:246 #13 0x0015dc25 in vtkMPISelfConnection::SendStreamToDataServer (this=0x2e89fa0, stream=value optimized out) at /home/kevin/kitware/ParaView3/Servers/Common/vtkMPISelfConnection.cxx:215 ---Type return to continue, or q return to quit--- #14 0x0015fb87 in vtkProcessModuleConnection::SendStream (this=0x2985a40, servers=value optimized out, stre...@0x7fff6aeac080) at /home/kevin/kitware/ParaView3/Servers/Common/vtkProcessModuleConnection.cxx:225 #15 0x001629a9 in vtkProcessModuleConnectionManager::SendStream (this=value optimized out, connectionID=value optimized out, serverFlags=1, stre...@0x7fff6aeac080, reset=1) at /home/kevin/kitware/ParaView3/Servers/Common/vtkProcessModuleConnectionManager.cxx:536 #16 0x0016a3ca in vtkProcessModule::SendStream (this=0x24d14c0, connectionID=1, server=1, stre...@0x7fff6aeac080, resetStream=1) at /home/kevin/kitware/ParaView3/Servers/Common/vtkProcessModule.cxx:870 #17 0x0015ae85 in vtkClientConnectionRMI (remoteArg=0x2f1c3c0, remoteArgLength=value optimized out) at /home/kevin/kitware/ParaView3/Servers/Common/vtkClientConnection.cxx:57 #18 0x7f035eebef2d in vtkMultiProcessController::ProcessRMI
Re: [Paraview] particle tracking
That's a really cool solution. Bravo Jean. -Ken On 2/12/09 6:58 AM, Jean Favre jfa...@cscs.ch wrote: Berk Geveci wrote: To animate particles in a steady-state flow field, I'd think that you would generate streamlines and then somehow animate particles along those. this is exactly the technique I use. Generate streamlines. Then iso-contour the streamline object with the scalar field IntegrationTime. Use a single threshold. Then animate the threshold value. Use Mode=Sequence, get a ramp between minimum time and maximun time. Use Glyphs (small spheres or arrows) attached to the iso-valued contours and they will animate their position along the streamlines. Jean Swiss national Supercomputing Center ___ 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 Kenneth Moreland *** Sandia National Laboratories *** *** *** *** email: kmo...@sandia.gov ** *** ** phone: (505) 844-8919 *** web: http://www.cs.unm.edu/~kmorel ___ 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
Re: [Paraview] particle tracking
This conversation has basically become academic, but it would also be pretty easy to create a filter that reported a bunch of time steps and just passed the same data every time. The particle tracer would thing it was a time varying data even though it was not. Jean's solution is still easier (and more clever), though. -Ken On 2/12/09 7:00 AM, Berk Geveci berk.gev...@kitware.com wrote: Interesting. I guess the biggest challenge of doing this in our pipeline would be finding the time range. Ideally, the integration should continue until all particles leave the domain or get stuck in a stagnant region. That would require integrating them first and then looking at the largest integration time. While I was writing this, Jean's response showed up. It is pretty ingenious :-) This is why I love ParaView. It allows you to do things in ways developers did not foresee. -berk On Thu, Feb 12, 2009 at 8:42 AM, John Biddiscombe biddi...@cscs.ch wrote: Berk Oops. I didn't notice that it was a steady state flow field. The TemporalStremTRacer expects unsteady flows, however it would be quite simple to modify it to use a single step. (I think easier than taking streamlines and animating particles along them - though thinking about it, the streamtracer code does have all that conversion of cell units to time which would be useful). One could subclass the temporaltreamtracer, suppress the time requests and interpolate over a single step, then output time dependent data where the particles animate along the path (- how many steps would one want?) - clicking the play button would just advect them along. (Just thinking out loud). JB I don't think this answers Pei's question. To my knowledge the temporal stream tracer (i.e. particle tracker) works only for time-dependent data. Am I wrong? To animate particles in a steady-state flow field, I'd think that you would generate streamlines and then somehow animate particles along those. I can't think of anything that would do this out-of-box but writing a filter that does this would be pretty easy. Is this something more than a few users would want? -berk On Thu, Feb 12, 2009 at 12:22 AM, David E DeMarle dave.dema...@kitware.com wrote: Hi Pei-Ying, You may want to try ParaView Meshless (https://twiki.cscs.ch/twiki/bin/view/ParaViewMeshless). It is a version of ParaView with the cutting edge of John Biddiscombe's particle tracking work in it. Some of those features will be/have been integrated into the main paraview code, but I can not vouch for when that will happen. cheers, Dave DeMarle 2008/12/25 Pei-Ying Hsieh phsieh2...@yahoo.com: Dear PVers: I have a steady state flow field. I am wondering if PV can do particle tracking animation. what I am looking for is similar to pathlines, but, instead of showing the lines, I would like to do an animation of seed particles flowing through the flow field. Is this possible? Thanks! and wish everyone here a Happy New Year! Pei ___ ParaView mailing list ParaView@paraview.org http://www.paraview.org/mailman/listinfo/paraview -- David E DeMarle Kitware, Inc. RD Engineer 28 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-371-3971 x109 ___ 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 -- John Biddiscombe,email:biddisco @ cscs.ch http://www.cscs.ch/ CSCS, Swiss National Supercomputing Centre | Tel: +41 (91) 610.82.07 Via Cantonale, 6928 Manno, Switzerland | Fax: +41 (91) 610.82.82 ___ 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 Kenneth Moreland *** Sandia National Laboratories *** *** *** *** email: kmo...@sandia.gov ** *** ** phone: (505) 844-8919 *** web: http://www.cs.unm.edu/~kmorel ___ Powered by www.kitware.com Visit other Kitware open-source projects at
Re: [Paraview] Why is PKdTree used for IMAGE Volume Rendering?
Kevin, Attached is a patch. Can you verify that it works after applying the patch? If so, I'll commit it to CVS. Thanks Utkarsh On Thu, Feb 12, 2009 at 4:45 PM, Kevin H. Hobbs hob...@ohiou.edu wrote: On Thu, 2009-02-12 at 15:34 -0500, Utkarsh Ayachit wrote: Ah and there's the bad news. In case of image data for the if condition at vtkPKdTree.cxx:435, I'd expect the condition to succeed and go to this-ProcessUserDefinedCuts(..) however, it seems like that condition is failing and it's going to this-MultiProcessBuildLocator(..) which is used for unstructured datasets. Are you able to reproduce the behavior with a simpler dataset (say wavelet) or something that;s already in ParaViewData, or some dataset that I can try it with? Yes, I was able to reproduce the problem with the wavelet source data. and I got a new datum as well. It seems I have to load the large image from a .vti file for the error to occur. The error does not occur when I create a large image in parallel with the wavelet source. Also give me all the steps to reproduce the problem (since I am not able to reproduce it with 5 procs using Wavelet source). Utkarsh In paraview connected to the builtin server I created a wavelet source with extents 0 1958 0 1498 0 77 to match my file. I then saved that to a vti file. I then restarted paraview and started a 2 core server that the dialog runs with mpirun -np 2 \ /home/kevin/kitware/ParaView3_Build/bin/pvserver \ --use-offscreen-rendering I loaded the image then selected volume rendering for the representation. Then BLAMO swap-death. Index: Servers/Filters/vtkKdTreeManager.cxx === RCS file: /cvsroot/ParaView3/ParaView3/Servers/Filters/vtkKdTreeManager.cxx,v retrieving revision 1.7 diff -u -3 -p -r1.7 vtkKdTreeManager.cxx --- Servers/Filters/vtkKdTreeManager.cxx 24 Nov 2008 14:32:41 - 1.7 +++ Servers/Filters/vtkKdTreeManager.cxx 12 Feb 2009 22:24:40 - @@ -24,6 +24,7 @@ #include vtkPoints.h #include vtkPVUpdateSuppressor.h #include vtkSmartPointer.h +#include vtkSphereSource.h #include vtkUnstructuredGrid.h #define VTK_CREATE(type, name) \ @@ -38,7 +39,6 @@ class vtkKdTreeManager::vtkAlgorithmSet vtkStandardNewMacro(vtkKdTreeManager); vtkCxxRevisionMacro(vtkKdTreeManager, $Revision: 1.7 $); vtkCxxSetObjectMacro(vtkKdTreeManager, StructuredProducer, vtkAlgorithm); -vtkCxxSetObjectMacro(vtkKdTreeManager, KdTree, vtkPKdTree); // vtkKdTreeManager::vtkKdTreeManager() { @@ -46,6 +46,7 @@ vtkKdTreeManager::vtkKdTreeManager() this-StructuredProducer = 0; this-KdTree = 0; this-NumberOfPieces = 1; + this-KdTreeInitialized = false; } // @@ -95,6 +96,16 @@ void vtkKdTreeManager::RemoveAllProducer } // +void vtkKdTreeManager::SetKdTree(vtkPKdTree* tree) +{ + if (this-KdTree != tree) +{ +vtkSetObjectBodyMacro(KdTree, vtkPKdTree, tree); +this-KdTreeInitialized = false; +} +} + +// void vtkKdTreeManager::Update() { vtkAlgorithmSet::iterator iter; @@ -132,17 +143,26 @@ void vtkKdTreeManager::Update() } this-KdTree-RemoveAllDataSets(); - for (dsIter = outputs.begin(); dsIter != outputs.end(); ++dsIter) -{ -this-AddDataSetToKdTree(*dsIter); -} - - static bool initialized = false; - if (!initialized) + if (!this-KdTreeInitialized) { +// HACK: This hack fixes the following issue: +// * create wavelet (num procs = 4) +// * volume render -- broken!!! +// * change some wavelet parameter (force the KdTree to rebuild) and all's fine! +// Seems like something doesn't get initialized correctly, I have no idea +// what. This seems to overcome the issue. +vtkSphereSource* sphere = vtkSphereSource::New(); +sphere-Update(); +this-KdTree-AddDataSet(sphere-GetOutput()); +sphere-Delete(); this-KdTree-BuildLocator(); -initialized = true; +this-KdTree-RemoveAllDataSets(); +this-KdTreeInitialized = true; } + for (dsIter = outputs.begin(); dsIter != outputs.end(); ++dsIter) +{ +this-AddDataSetToKdTree(*dsIter); +} if (this-StructuredProducer) { Index: Servers/Filters/vtkKdTreeManager.h === RCS file: /cvsroot/ParaView3/ParaView3/Servers/Filters/vtkKdTreeManager.h,v retrieving revision 1.2 diff -u -3 -p -r1.2 vtkKdTreeManager.h --- Servers/Filters/vtkKdTreeManager.h 7 Sep 2007 13:34:15 - 1.2 +++ Servers/Filters/vtkKdTreeManager.h 12 Feb 2009 22:24:40 - @@ -65,6 +65,7 @@ protected: void AddDataSetToKdTree(vtkDataSet *data); + bool
Re: [Paraview] Why is PKdTree used for IMAGE Volume Rendering?
On Thu, 2009-02-12 at 17:27 -0500, Utkarsh Ayachit wrote: Kevin, Attached is a patch. Can you verify that it works after applying the patch? If so, I'll commit it to CVS. Thanks Utkarsh Yup that fixed it. The PKdTree completes immediately and then I get my volume rendering. Thank you! signature.asc Description: This is a digitally signed message part ___ 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
Re: [Paraview] particletracer
Sorry, but this doesnt make any sense to me. I have a 3D dataset with magnitude and velocity (X,Y,Z) for 13 time frames. And as you wrote, paraview internally converts multiple vktDataSets into a TemporalDataSet. But particle tracer requires multiple time steps, why does it then give error that a vktDataSet is required? A vktDataSet, according to what you wrote, is just a single timestep. Martin Martin A single timestep is a vtkDataSet, StremTracer generates streamline through a vector field for this data. The ParticleTracer requires multiple time steps to advect particles through the field over time. The particle tracer requests 2 time steps per iteration from it's input, internally paraview converts multiple vtkDataSets from the input into a TemporalDataSet. JB Thanks you, I can now conenct the two sources. I tried on Windows and found this alternative if I want to use input or source in the particle tracer. When choosing source I could select both dataset and seed source. Problem on my Mac was that the dialog box was to small so I didnt find this option. Btw what is the difference between vtkTemporalDataSet and vktDataSet? ParticleTracer requires vktDataSet which I dont have, is it possible to convert the dataset somehow. I can use filter StreamTracer on my dataset and because StreamTracer and ParticleTracer is based on the same physics, I should be able to visualize particles. ERROR: In /Users/kitware/Berk/ParaView3/VTK/Filtering/vtkDemandDrivenPipeline.cxx, line 822 vtkCompositeDataPipeline (0x225e2020): Input for connection index 0 on input port index 1 for algorithm vtkTemporalStreamTracer(0x225e20b0) is of type vtkMultiBlockDataSet, but a vtkDataSet is required. Cheers, Martin Uppman Can you tell us exactly (detailed, step-by-step) what you are doing? That dialog works fine for me (running OS X 10.5). -berk On Tue, Feb 10, 2009 at 4:10 AM, Martin Uppman wrote: Hi again, thank you for fast answer John! I'm aware of the problem that my particle tracer doesnt have any seed source, but I have been working for quite awhile sloving this problem. I did as you said, created temporalcache and seed (point seed) but I cant connect the 2 with my particle tracer. I found your (John's) great powerpoint presentation on the web explaining everything there is to explain. Though I cant make it work!!! In one of the pictures in the presentation there are some boxes where you can choose which sources to use in the particle tracer (in object inspector - properties - source). I dont have that option in my program. Does someone know why (using v3.4.0 on a Mac)? About v3.4.0, isnt it the lastest version? Cant fint anything newer on paraview.org. Cheers, Martin Uppman Martin ERROR: In /Users/kitware/Berk/ParaView3/VTK/Filtering/vtkDemandDrivenPipeline.cxx, line 722 vtkCompositeDataPipeline (0x217c7fb0): Input port 1 of algorithm vtkTemporalStreamTracer(0x1e345e00) has 0 connections but is not optional. Input 0 should be a time dependent dataset. It is best to use a TemporalCache before the particle tracer to stop multiple updates for the same time step. Input 1 should be some seed points. Ensight - TemporalCache - ParticleTracer Some Other data eg.line/plane - the ignore pipeline time option is provided to get around some animation issues that used to plague paraview. You should not need it, but the correct way to use it, is set IgnorePipelineTime to true, then set an animation keyframe to animate TimeStep from 0-N - the particle tracer will then use the timestep to generate update times for the Ensight input, instead of usigin the GUI generated time. If you are using the particle tracer from prior to 3.4 then you should switch to a newer one as it is contiunually being improved. I will soon (months, not days) be checking in a lot of changes which enable it to handle dynamic meshes better. JB I've been looking for information on this and a found a thread which I think concerns the same problem. In that thread Berk made a patch to fix this problem. I also tried to ignore pipeline time, but this doent work. I'm working with .case files and I want to animate the blood flow in a human heart using particle tracers. I'm using a Mac with ParaView version 3.4.0. Thanks for looking into this. Cheers, Martin Uppman ___ 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 -- John Biddiscombe, email:biddisco @ cscs.ch http://www.cscs.ch/ CSCS, Swiss National Supercomputing Centre | Tel: +41 (91)