Re: [Paraview] Numpy import error
*works since I had to change PYTHON_LIBRARY to point to the lib subdirectory of /Library/Frameworks/Python.framework/Versions/2.7/, similar to the executable and include_dir versions On Thu, Mar 21, 2013 at 10:59 PM, Sohail Shafii sssha...@ucdavis.eduwrote: Importing numpy into python works. Same with import sys.path. What's interesting is that ParaView uses python stored in /usr/local/ while the version of Python that I use (i.e. default python from the terminal) resides in /Library/Frameworks. The one in /usr/local does not have numpy. I didn't think that ParaView was using the one in /urs/local, since it was using the python binaries (PYTHON_EXECUTABLE and PYTHON_INCLUDE_DIR) from the Frameworks directory according to ccmake. However, it was using the /usr/local version for the PYTHON_LIBRARY, which didn't catch my attention when I was browsing the ccmake configuration. So it was effectively using both for different tasks. Not surprising since my library path pointed to /usr/lib and my path had the Frameworks folder as the first item. I recompiled and it works now. Thanks. On Thu, Mar 21, 2013 at 5:59 PM, Utkarsh Ayachit utkarsh.ayac...@kitware.com wrote: Run python from terminal and import numpy. Does that work? If so, import sys and print sys.path. Now try the same in Python shell from ParaView. Does that work? If not, compare the sys.path from python and ParaView's Python shell. Any differences? On Thu, Mar 21, 2013 at 8:53 PM, Sohail Shafii sssha...@ucdavis.edu wrote: Hi, I have numpy installed installed on Python 2.7.3. for Mac OS 10.8.3. I compiled against the same python while building ParaView 3.98.1 64-bit, but when I create a Python programmable filter in ParaView, I cannot import numpy -- I get a message that says ImportError: No module named numpy. Python was enabled during compilation as well. Any advice? Thanks, Sohail ___ 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
[Paraview] Numpy import error
Hi, I have numpy installed installed on Python 2.7.3. for Mac OS 10.8.3. I compiled against the same python while building ParaView 3.98.1 64-bit, but when I create a Python programmable filter in ParaView, I cannot import numpy -- I get a message that says ImportError: No module named numpy. Python was enabled during compilation as well. Any advice? Thanks, Sohail ___ 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] Numpy import error
Importing numpy into python works. Same with import sys.path. What's interesting is that ParaView uses python stored in /usr/local/ while the version of Python that I use (i.e. default python from the terminal) resides in /Library/Frameworks. The one in /usr/local does not have numpy. I didn't think that ParaView was using the one in /urs/local, since it was using the python binaries (PYTHON_EXECUTABLE and PYTHON_INCLUDE_DIR) from the Frameworks directory according to ccmake. However, it was using the /usr/local version for the PYTHON_LIBRARY, which didn't catch my attention when I was browsing the ccmake configuration. So it was effectively using both for different tasks. Not surprising since my library path pointed to /usr/lib and my path had the Frameworks folder as the first item. I recompiled and it works now. Thanks. On Thu, Mar 21, 2013 at 5:59 PM, Utkarsh Ayachit utkarsh.ayac...@kitware.com wrote: Run python from terminal and import numpy. Does that work? If so, import sys and print sys.path. Now try the same in Python shell from ParaView. Does that work? If not, compare the sys.path from python and ParaView's Python shell. Any differences? On Thu, Mar 21, 2013 at 8:53 PM, Sohail Shafii sssha...@ucdavis.edu wrote: Hi, I have numpy installed installed on Python 2.7.3. for Mac OS 10.8.3. I compiled against the same python while building ParaView 3.98.1 64-bit, but when I create a Python programmable filter in ParaView, I cannot import numpy -- I get a message that says ImportError: No module named numpy. Python was enabled during compilation as well. Any advice? Thanks, Sohail ___ 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] ParaView thinks that it does not have a wind blade reader
Thanks for the help. On Sunday, January 20, 2013, Andy Bauer wrote: FYI: this finally has made its way into ParaView such that it will be available in the GUI for single process runs. Andy On Thu, Dec 20, 2012 at 8:08 PM, Sohail Shafii sssha...@ucdavis.eduwrote: Thanks. Please take your time. If there any other problems with this reader then feel free to consult me. 1) typically I make a separate build directory for VTK or ParaView. The test I ran looked for ../../../../vtkData or something like that, so it seems like the test executable has to be four directories deep. So I guess there is a common parent to the data folders and paraview, not sure what this parent would be or if there is a reference for this hierarchy. 2) I guess in this case I would have to experiment when I do make the switch to the newer VTK On Thu, Dec 20, 2012 at 9:45 AM, Andy Bauer andy.ba...@kitware.comwrote: Hi Sohail, Changing the number of fields from 23 to 16 fixed the warnings. I'll get these changes into VTK and then change ParaView so that it can read in WindBlade files with a single process. That may take a week or so though. As for the other questions: 1) If the data directories are called VTKData, VTKLargeData, and ParaViewData and they share the same parent directory as VTK or ParaView, when you build VTK or ParaView it will automatically pick up the path to each of them. 2) For the library names, after modularization they've been changed and relocated from the bin subdirectory to the lib subdirectory. Many of the old libraries have been subdivided into better organized libraries. I don't know any specific document that specifies what to do for this though. Andy On Sun, Dec 16, 2012 at 9:39 PM, Sohail Shafii sssha...@ucdavis.eduwrote: I have some time now, so I addressed the issue. I don't the code has a problem; the problem is related to the .wind file. Near the bottom, it describes how many fields exist in the comp.out files (DATA_VARIABLES). For your test data, each file is about 64.5 MB. The .wind says that each file should have 23 fields. For n scalar fields, 1 is reserved for UVW; (n-1) for scalar data. So for this one, 1 for UWV while 22 for scalar data. The UVW block is about 10 MB (160x80x70x(4 bytes)x(3 components)); each scalar field is about 3.6 MB (160x80x70x(4 bytes)). So 10 MB + 22 scalar fields*3.6 gives you about 90 MB. That's way too big compared to the actual sizes of the comp.out files that you guys have for testing. * The test works when you reduce the DATA_VARIABLES item from 23 to 16*(which will give you close to 64.5 MB). I'm not sure why this .wind file is bad. Maybe older wind blade readers just ignored errors like these. Anyway, the test seems to launch a window and close right away. I didn't see any other errors after I fixed the wind file that I am using. A couple of questions: 1. How do you organize the VTK build directory with the testdata? The VTK tests seem to expect the data to be four directories UP the hierarchy, so I had to place my vtkData directory in my home folder. 2. I tried the new versions of VTK with QT and it appears that the -l commands that I put in the .pro file for linking (i.e. -lvtkIO) don't seem to work...is there a document somewhere that helps in this regard? I reverted to an older version of VTK for now but I want to eventually upgrade. Sohail On Sat, Nov 3, 2012 at 1:16 PM, Sohail Shafii sssha...@ucdavis.eduwrote: Also I'm removing Jimmy from this list since he may not be working with wind files at this point. Sohail On Sat, Nov ___ 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] ParaView thinks that it does not have a wind blade reader
I have some time now, so I addressed the issue. I don't the code has a problem; the problem is related to the .wind file. Near the bottom, it describes how many fields exist in the comp.out files (DATA_VARIABLES). For your test data, each file is about 64.5 MB. The .wind says that each file should have 23 fields. For n scalar fields, 1 is reserved for UVW; (n-1) for scalar data. So for this one, 1 for UWV while 22 for scalar data. The UVW block is about 10 MB (160x80x70x(4 bytes)x(3 components)); each scalar field is about 3.6 MB (160x80x70x(4 bytes)). So 10 MB + 22 scalar fields*3.6 gives you about 90 MB. That's way too big compared to the actual sizes of the comp.out files that you guys have for testing. * The test works when you reduce the DATA_VARIABLES item from 23 to 16* (which will give you close to 64.5 MB). I'm not sure why this .wind file is bad. Maybe older wind blade readers just ignored errors like these. Anyway, the test seems to launch a window and close right away. I didn't see any other errors after I fixed the wind file that I am using. A couple of questions: 1. How do you organize the VTK build directory with the testdata? The VTK tests seem to expect the data to be four directories UP the hierarchy, so I had to place my vtkData directory in my home folder. 2. I tried the new versions of VTK with QT and it appears that the -l commands that I put in the .pro file for linking (i.e. -lvtkIO) don't seem to work...is there a document somewhere that helps in this regard? I reverted to an older version of VTK for now but I want to eventually upgrade. Sohail On Sat, Nov 3, 2012 at 1:16 PM, Sohail Shafii sssha...@ucdavis.edu wrote: Also I'm removing Jimmy from this list since he may not be working with wind files at this point. Sohail On Sat, Nov 3, 2012 at 1:15 PM, Sohail Shafii sssha...@ucdavis.eduwrote: At this point the changes will have to wait, unfortunately. For the past two weeks I've been involved with a TopoInVis submission so I couldn't get around to it. I'll look at it this weekend. No need to wait for me. On Thu, Nov 1, 2012 at 10:29 AM, Andy Bauer andy.ba...@kitware.comwrote: We're getting ready for a pre-release of PV 4.0 now so we'll need this soon. If we need to update to a newer WindBlade file for the test, we can do that too. You'll have to think though on how to handle backwards compatibility. Andy On Wed, Oct 31, 2012 at 12:37 PM, Sohail Shafii sssha...@ucdavis.eduwrote: Well you probably know already, but I won't be able to make the changes right now. I have a submission due this Friday. The reader still works for newer data sets but fails for the outdated one that you guys use for testing. So the fixes will have to come later. On Fri, Oct 26, 2012 at 2:39 PM, Andy Bauer andy.ba...@kitware.comwrote: This change is only in VTK in gerrit so you need to follow the directions there to check out the code and run the test. The directions are at http://review.source.kitware.com/#/t/1396/ I'm not sure why paraview is crashing though. I'll check on that. andy On Fri, Oct 26, 2012 at 5:06 PM, Sohail Shafii sssha...@ucdavis.eduwrote: Current master version can't launch (MacOS Mountain Lion): paraview(575,0x7fff72e33180) malloc: *** error for object 0x111058450: pointer being freed was not allocated *** set a breakpoint in malloc_error_break to debug Program received signal SIGABRT, Aborted. 0x7fff82432212 in __pthread_kill () (gdb) bt 14 #0 0x7fff82432212 in __pthread_kill () #1 0x7fff87964af4 in pthread_kill () #2 0x7fff879a8dce in abort () #3 0x7fff8797c959 in free () #4 0x7fff5fc0762a in __dyld__ZN4dyldL12registerDOFsERKNSt3__16vectorIN11ImageLoader7DOFInfoENS0_9allocatorIS3_ () #5 0x7fff5fc0f06b in __dyld__ZN11ImageLoader4linkERKNS_11LinkContextEbbRKNS_10RPathChainE () #6 0x7fff5fc04ce4 in __dyld__ZN4dyld4linkEP11ImageLoaderbRKNS0_10RPathChainE () #7 0x7fff5fc0a81a in __dyld_dlopen () #8 0x7fff8ac96dd8 in dlopen () #9 0x7fff83dcfa10 in gfxInitializeLibrary () #10 0x00011e6014fe in gliInitializeLibrary () #11 0x7fff8a4f63e6 in glcPluginCount () #12 0x7fff8a4f62a8 in glcPluginCount () #13 0x7fff8a4f6091 in glcPluginCount () (More stack frames follow...) Which version of the master works, so that I can debug? Sohail On Wed, Oct 24, 2012 at 10:04 AM, Andy Bauer andy.ba...@kitware.comwrote: A beta version will be released by mid-November. In order for the desired changes to make it in to that, the VTK part that's waiting your changes needs to go in early next week so that I can make the changes to ParaView to enable the WindBlade reader when running in either built-in mode or connecting to a serial server. Andy On Wed, Oct 24, 2012 at 12:13 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Haven't had a chance yet due to a paper. I'll look at it this week. How soon will it be released? Also
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
At this point the changes will have to wait, unfortunately. For the past two weeks I've been involved with a TopoInVis submission so I couldn't get around to it. I'll look at it this weekend. No need to wait for me. On Thu, Nov 1, 2012 at 10:29 AM, Andy Bauer andy.ba...@kitware.com wrote: We're getting ready for a pre-release of PV 4.0 now so we'll need this soon. If we need to update to a newer WindBlade file for the test, we can do that too. You'll have to think though on how to handle backwards compatibility. Andy On Wed, Oct 31, 2012 at 12:37 PM, Sohail Shafii sssha...@ucdavis.eduwrote: Well you probably know already, but I won't be able to make the changes right now. I have a submission due this Friday. The reader still works for newer data sets but fails for the outdated one that you guys use for testing. So the fixes will have to come later. On Fri, Oct 26, 2012 at 2:39 PM, Andy Bauer andy.ba...@kitware.comwrote: This change is only in VTK in gerrit so you need to follow the directions there to check out the code and run the test. The directions are at http://review.source.kitware.com/#/t/1396/ I'm not sure why paraview is crashing though. I'll check on that. andy On Fri, Oct 26, 2012 at 5:06 PM, Sohail Shafii sssha...@ucdavis.eduwrote: Current master version can't launch (MacOS Mountain Lion): paraview(575,0x7fff72e33180) malloc: *** error for object 0x111058450: pointer being freed was not allocated *** set a breakpoint in malloc_error_break to debug Program received signal SIGABRT, Aborted. 0x7fff82432212 in __pthread_kill () (gdb) bt 14 #0 0x7fff82432212 in __pthread_kill () #1 0x7fff87964af4 in pthread_kill () #2 0x7fff879a8dce in abort () #3 0x7fff8797c959 in free () #4 0x7fff5fc0762a in __dyld__ZN4dyldL12registerDOFsERKNSt3__16vectorIN11ImageLoader7DOFInfoENS0_9allocatorIS3_ () #5 0x7fff5fc0f06b in __dyld__ZN11ImageLoader4linkERKNS_11LinkContextEbbRKNS_10RPathChainE () #6 0x7fff5fc04ce4 in __dyld__ZN4dyld4linkEP11ImageLoaderbRKNS0_10RPathChainE () #7 0x7fff5fc0a81a in __dyld_dlopen () #8 0x7fff8ac96dd8 in dlopen () #9 0x7fff83dcfa10 in gfxInitializeLibrary () #10 0x00011e6014fe in gliInitializeLibrary () #11 0x7fff8a4f63e6 in glcPluginCount () #12 0x7fff8a4f62a8 in glcPluginCount () #13 0x7fff8a4f6091 in glcPluginCount () (More stack frames follow...) Which version of the master works, so that I can debug? Sohail On Wed, Oct 24, 2012 at 10:04 AM, Andy Bauer andy.ba...@kitware.comwrote: A beta version will be released by mid-November. In order for the desired changes to make it in to that, the VTK part that's waiting your changes needs to go in early next week so that I can make the changes to ParaView to enable the WindBlade reader when running in either built-in mode or connecting to a serial server. Andy On Wed, Oct 24, 2012 at 12:13 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Haven't had a chance yet due to a paper. I'll look at it this week. How soon will it be released? Also please respond to my sssha...@ucdavis.edu as this address is not on the paraview list anymore. Sohail -- *From:* Andy Bauer andy.ba...@kitware.com *To:* Sohail Shafii sohailsha...@yahoo.com *Cc:* David E DeMarle dave.dema...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov *Sent:* Wednesday, October 24, 2012 8:09 AM *Subject:* Re: [Paraview] ParaView thinks that it does not have a wind blade reader Hi Sohail, Did you get a chance to look at the warnings for the vtkIOGeometryCxx-TestWindBladeReader test for the http://review.source.kitware.com/#/t/1396/ gerrit topic? It's going to need to be fixed pretty soon if it's going to make it into the ParaView 4.0 release. Thanks, Andy On Wed, Oct 3, 2012 at 10:22 PM, Andy Bauer andy.ba...@kitware.comwrote: It's available with git clone git://vtk.org/VTKLargeData.gitVTKLargeData. It's in the Data/WindBladeReader subdirectory there. On Wed, Oct 3, 2012 at 6:27 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Which data set are you using for the test? If it's the single turbine data set with the curvy terrain, I don't have it (I use other data sets) so I can't do a similar test. Let me know where I can download it. Sohail -- *From:* Andy Bauer andy.ba...@kitware.com *To:* David E DeMarle dave.dema...@kitware.com *Cc:* Sohail Shafii sohailsha...@yahoo.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov *Sent:* Wednesday, October 3, 2012 3:16 PM *Subject:* Re: [Paraview] ParaView thinks that it does not have a wind blade reader I finally got around to doing the VTK part of the WindBladeReader work. It's on gerrit at http
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
Also I'm removing Jimmy from this list since he may not be working with wind files at this point. Sohail On Sat, Nov 3, 2012 at 1:15 PM, Sohail Shafii sssha...@ucdavis.edu wrote: At this point the changes will have to wait, unfortunately. For the past two weeks I've been involved with a TopoInVis submission so I couldn't get around to it. I'll look at it this weekend. No need to wait for me. On Thu, Nov 1, 2012 at 10:29 AM, Andy Bauer andy.ba...@kitware.comwrote: We're getting ready for a pre-release of PV 4.0 now so we'll need this soon. If we need to update to a newer WindBlade file for the test, we can do that too. You'll have to think though on how to handle backwards compatibility. Andy On Wed, Oct 31, 2012 at 12:37 PM, Sohail Shafii sssha...@ucdavis.eduwrote: Well you probably know already, but I won't be able to make the changes right now. I have a submission due this Friday. The reader still works for newer data sets but fails for the outdated one that you guys use for testing. So the fixes will have to come later. On Fri, Oct 26, 2012 at 2:39 PM, Andy Bauer andy.ba...@kitware.comwrote: This change is only in VTK in gerrit so you need to follow the directions there to check out the code and run the test. The directions are at http://review.source.kitware.com/#/t/1396/ I'm not sure why paraview is crashing though. I'll check on that. andy On Fri, Oct 26, 2012 at 5:06 PM, Sohail Shafii sssha...@ucdavis.eduwrote: Current master version can't launch (MacOS Mountain Lion): paraview(575,0x7fff72e33180) malloc: *** error for object 0x111058450: pointer being freed was not allocated *** set a breakpoint in malloc_error_break to debug Program received signal SIGABRT, Aborted. 0x7fff82432212 in __pthread_kill () (gdb) bt 14 #0 0x7fff82432212 in __pthread_kill () #1 0x7fff87964af4 in pthread_kill () #2 0x7fff879a8dce in abort () #3 0x7fff8797c959 in free () #4 0x7fff5fc0762a in __dyld__ZN4dyldL12registerDOFsERKNSt3__16vectorIN11ImageLoader7DOFInfoENS0_9allocatorIS3_ () #5 0x7fff5fc0f06b in __dyld__ZN11ImageLoader4linkERKNS_11LinkContextEbbRKNS_10RPathChainE () #6 0x7fff5fc04ce4 in __dyld__ZN4dyld4linkEP11ImageLoaderbRKNS0_10RPathChainE () #7 0x7fff5fc0a81a in __dyld_dlopen () #8 0x7fff8ac96dd8 in dlopen () #9 0x7fff83dcfa10 in gfxInitializeLibrary () #10 0x00011e6014fe in gliInitializeLibrary () #11 0x7fff8a4f63e6 in glcPluginCount () #12 0x7fff8a4f62a8 in glcPluginCount () #13 0x7fff8a4f6091 in glcPluginCount () (More stack frames follow...) Which version of the master works, so that I can debug? Sohail On Wed, Oct 24, 2012 at 10:04 AM, Andy Bauer andy.ba...@kitware.comwrote: A beta version will be released by mid-November. In order for the desired changes to make it in to that, the VTK part that's waiting your changes needs to go in early next week so that I can make the changes to ParaView to enable the WindBlade reader when running in either built-in mode or connecting to a serial server. Andy On Wed, Oct 24, 2012 at 12:13 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Haven't had a chance yet due to a paper. I'll look at it this week. How soon will it be released? Also please respond to my sssha...@ucdavis.edu as this address is not on the paraview list anymore. Sohail -- *From:* Andy Bauer andy.ba...@kitware.com *To:* Sohail Shafii sohailsha...@yahoo.com *Cc:* David E DeMarle dave.dema...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov *Sent:* Wednesday, October 24, 2012 8:09 AM *Subject:* Re: [Paraview] ParaView thinks that it does not have a wind blade reader Hi Sohail, Did you get a chance to look at the warnings for the vtkIOGeometryCxx-TestWindBladeReader test for the http://review.source.kitware.com/#/t/1396/ gerrit topic? It's going to need to be fixed pretty soon if it's going to make it into the ParaView 4.0 release. Thanks, Andy On Wed, Oct 3, 2012 at 10:22 PM, Andy Bauer andy.ba...@kitware.comwrote: It's available with git clone git://vtk.org/VTKLargeData.gitVTKLargeData. It's in the Data/WindBladeReader subdirectory there. On Wed, Oct 3, 2012 at 6:27 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Which data set are you using for the test? If it's the single turbine data set with the curvy terrain, I don't have it (I use other data sets) so I can't do a similar test. Let me know where I can download it. Sohail -- *From:* Andy Bauer andy.ba...@kitware.com *To:* David E DeMarle dave.dema...@kitware.com *Cc:* Sohail Shafii sohailsha...@yahoo.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov *Sent:* Wednesday, October 3, 2012 3:16 PM
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
Well you probably know already, but I won't be able to make the changes right now. I have a submission due this Friday. The reader still works for newer data sets but fails for the outdated one that you guys use for testing. So the fixes will have to come later. On Fri, Oct 26, 2012 at 2:39 PM, Andy Bauer andy.ba...@kitware.com wrote: This change is only in VTK in gerrit so you need to follow the directions there to check out the code and run the test. The directions are at http://review.source.kitware.com/#/t/1396/ I'm not sure why paraview is crashing though. I'll check on that. andy On Fri, Oct 26, 2012 at 5:06 PM, Sohail Shafii sssha...@ucdavis.eduwrote: Current master version can't launch (MacOS Mountain Lion): paraview(575,0x7fff72e33180) malloc: *** error for object 0x111058450: pointer being freed was not allocated *** set a breakpoint in malloc_error_break to debug Program received signal SIGABRT, Aborted. 0x7fff82432212 in __pthread_kill () (gdb) bt 14 #0 0x7fff82432212 in __pthread_kill () #1 0x7fff87964af4 in pthread_kill () #2 0x7fff879a8dce in abort () #3 0x7fff8797c959 in free () #4 0x7fff5fc0762a in __dyld__ZN4dyldL12registerDOFsERKNSt3__16vectorIN11ImageLoader7DOFInfoENS0_9allocatorIS3_ () #5 0x7fff5fc0f06b in __dyld__ZN11ImageLoader4linkERKNS_11LinkContextEbbRKNS_10RPathChainE () #6 0x7fff5fc04ce4 in __dyld__ZN4dyld4linkEP11ImageLoaderbRKNS0_10RPathChainE () #7 0x7fff5fc0a81a in __dyld_dlopen () #8 0x7fff8ac96dd8 in dlopen () #9 0x7fff83dcfa10 in gfxInitializeLibrary () #10 0x00011e6014fe in gliInitializeLibrary () #11 0x7fff8a4f63e6 in glcPluginCount () #12 0x7fff8a4f62a8 in glcPluginCount () #13 0x7fff8a4f6091 in glcPluginCount () (More stack frames follow...) Which version of the master works, so that I can debug? Sohail On Wed, Oct 24, 2012 at 10:04 AM, Andy Bauer andy.ba...@kitware.comwrote: A beta version will be released by mid-November. In order for the desired changes to make it in to that, the VTK part that's waiting your changes needs to go in early next week so that I can make the changes to ParaView to enable the WindBlade reader when running in either built-in mode or connecting to a serial server. Andy On Wed, Oct 24, 2012 at 12:13 PM, Sohail Shafii sohailsha...@yahoo.comwrote: Haven't had a chance yet due to a paper. I'll look at it this week. How soon will it be released? Also please respond to my sssha...@ucdavis.edu as this address is not on the paraview list anymore. Sohail -- *From:* Andy Bauer andy.ba...@kitware.com *To:* Sohail Shafii sohailsha...@yahoo.com *Cc:* David E DeMarle dave.dema...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov *Sent:* Wednesday, October 24, 2012 8:09 AM *Subject:* Re: [Paraview] ParaView thinks that it does not have a wind blade reader Hi Sohail, Did you get a chance to look at the warnings for the vtkIOGeometryCxx-TestWindBladeReader test for the http://review.source.kitware.com/#/t/1396/ gerrit topic? It's going to need to be fixed pretty soon if it's going to make it into the ParaView 4.0 release. Thanks, Andy On Wed, Oct 3, 2012 at 10:22 PM, Andy Bauer andy.ba...@kitware.comwrote: It's available with git clone git://vtk.org/VTKLargeData.gitVTKLargeData. It's in the Data/WindBladeReader subdirectory there. On Wed, Oct 3, 2012 at 6:27 PM, Sohail Shafii sohailsha...@yahoo.comwrote: Which data set are you using for the test? If it's the single turbine data set with the curvy terrain, I don't have it (I use other data sets) so I can't do a similar test. Let me know where I can download it. Sohail -- *From:* Andy Bauer andy.ba...@kitware.com *To:* David E DeMarle dave.dema...@kitware.com *Cc:* Sohail Shafii sohailsha...@yahoo.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov *Sent:* Wednesday, October 3, 2012 3:16 PM *Subject:* Re: [Paraview] ParaView thinks that it does not have a wind blade reader I finally got around to doing the VTK part of the WindBladeReader work. It's on gerrit at http://review.source.kitware.com/#/c/7734/. It still needs some cleanup but at least both the serial test and the parallel test with a single process pass. For the serial test, vtkIOGeometryCxx-TestWindBladeReader, I'm still getting the following warning which needs to be fixed as well: Warning: In /home/acbauer/CODE/VTK/VTK/IO/Geometry/vtkWindBladeReader.cxx, line 568 vtkWindBladeReader (0x1373840): WindBladeReader error reading file: /home/acbauer/DATA/VTKLargeData/Data/WindBladeReader/test1_topo.wind Premature EOF while reading block of data. Expected 896000 but got 0 Sohail, can you check this out and fix it? Thanks
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
Current master version can't launch (MacOS Mountain Lion): paraview(575,0x7fff72e33180) malloc: *** error for object 0x111058450: pointer being freed was not allocated *** set a breakpoint in malloc_error_break to debug Program received signal SIGABRT, Aborted. 0x7fff82432212 in __pthread_kill () (gdb) bt 14 #0 0x7fff82432212 in __pthread_kill () #1 0x7fff87964af4 in pthread_kill () #2 0x7fff879a8dce in abort () #3 0x7fff8797c959 in free () #4 0x7fff5fc0762a in __dyld__ZN4dyldL12registerDOFsERKNSt3__16vectorIN11ImageLoader7DOFInfoENS0_9allocatorIS3_ () #5 0x7fff5fc0f06b in __dyld__ZN11ImageLoader4linkERKNS_11LinkContextEbbRKNS_10RPathChainE () #6 0x7fff5fc04ce4 in __dyld__ZN4dyld4linkEP11ImageLoaderbRKNS0_10RPathChainE () #7 0x7fff5fc0a81a in __dyld_dlopen () #8 0x7fff8ac96dd8 in dlopen () #9 0x7fff83dcfa10 in gfxInitializeLibrary () #10 0x00011e6014fe in gliInitializeLibrary () #11 0x7fff8a4f63e6 in glcPluginCount () #12 0x7fff8a4f62a8 in glcPluginCount () #13 0x7fff8a4f6091 in glcPluginCount () (More stack frames follow...) Which version of the master works, so that I can debug? Sohail On Wed, Oct 24, 2012 at 10:04 AM, Andy Bauer andy.ba...@kitware.com wrote: A beta version will be released by mid-November. In order for the desired changes to make it in to that, the VTK part that's waiting your changes needs to go in early next week so that I can make the changes to ParaView to enable the WindBlade reader when running in either built-in mode or connecting to a serial server. Andy On Wed, Oct 24, 2012 at 12:13 PM, Sohail Shafii sohailsha...@yahoo.comwrote: Haven't had a chance yet due to a paper. I'll look at it this week. How soon will it be released? Also please respond to my sssha...@ucdavis.edu as this address is not on the paraview list anymore. Sohail -- *From:* Andy Bauer andy.ba...@kitware.com *To:* Sohail Shafii sohailsha...@yahoo.com *Cc:* David E DeMarle dave.dema...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov *Sent:* Wednesday, October 24, 2012 8:09 AM *Subject:* Re: [Paraview] ParaView thinks that it does not have a wind blade reader Hi Sohail, Did you get a chance to look at the warnings for the vtkIOGeometryCxx-TestWindBladeReader test for the http://review.source.kitware.com/#/t/1396/ gerrit topic? It's going to need to be fixed pretty soon if it's going to make it into the ParaView 4.0 release. Thanks, Andy On Wed, Oct 3, 2012 at 10:22 PM, Andy Bauer andy.ba...@kitware.comwrote: It's available with git clone git://vtk.org/VTKLargeData.gitVTKLargeData. It's in the Data/WindBladeReader subdirectory there. On Wed, Oct 3, 2012 at 6:27 PM, Sohail Shafii sohailsha...@yahoo.comwrote: Which data set are you using for the test? If it's the single turbine data set with the curvy terrain, I don't have it (I use other data sets) so I can't do a similar test. Let me know where I can download it. Sohail -- *From:* Andy Bauer andy.ba...@kitware.com *To:* David E DeMarle dave.dema...@kitware.com *Cc:* Sohail Shafii sohailsha...@yahoo.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov *Sent:* Wednesday, October 3, 2012 3:16 PM *Subject:* Re: [Paraview] ParaView thinks that it does not have a wind blade reader I finally got around to doing the VTK part of the WindBladeReader work. It's on gerrit at http://review.source.kitware.com/#/c/7734/. It still needs some cleanup but at least both the serial test and the parallel test with a single process pass. For the serial test, vtkIOGeometryCxx-TestWindBladeReader, I'm still getting the following warning which needs to be fixed as well: Warning: In /home/acbauer/CODE/VTK/VTK/IO/Geometry/vtkWindBladeReader.cxx, line 568 vtkWindBladeReader (0x1373840): WindBladeReader error reading file: /home/acbauer/DATA/VTKLargeData/Data/WindBladeReader/test1_topo.wind Premature EOF while reading block of data. Expected 896000 but got 0 Sohail, can you check this out and fix it? Thanks, Andy On Fri, Jul 20, 2012 at 12:02 PM, David E DeMarle dave.dema...@kitware.com wrote: Thanks. In today's git master it now sits in VTK/IO/Parallel/vtkWindbladeReader. The serial version should go in VTK/IO/Geometry. The MPI version should go in a new module VTK/IO/MPIParallel. (VTK/IO/Parallel depends only on vtkSocketController now and should not depend on vtkMPIController). Finally when that is done we can remove the flag in paraview's xml that prevents it from being used in serial. David E DeMarle Kitware, Inc. RD Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909 On Thu, Jul 19, 2012 at 7:08 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Well
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
MPI is only used if it's available. It tries to see if there is a MPIController available. If not, then the number of ranks is 1, which means that the parallel class calls the super class's requestdata (super class is the serial wind blade reader). However, the parallel class overrides certain functions to make use of MPI. So even the parallel instantiation will call the superclass's requestdata, the overriden functions that are called from that point on assume that MPI is available...this is a problem when the MPI controller is not available. On another note, it appears that the rectilinear grid that I am reading (VTK file, structured points I think) is read by the vtkFileSeriesReader, which makes use of the vtkPDataSetReader class. Strange that if I use vtkPDataSetReader myself in C++ code that uses VTK and not paraview, bounds/cell information is not generated. So I think that ParaView does something else besides using that reader...hmmm. On Mon, Aug 27, 2012 at 7:32 PM, Andy Bauer andy.ba...@kitware.com wrote: On Mon, Aug 27, 2012 at 6:53 PM, Sohail Shafii sssha...@ucdavis.eduwrote: I see...so if MPI is running then it will create a parallel version of the class. Otherwise, it will create the normal reader. Technically, if ParaView is built with MPI enabled it will create the parallel version of the class every time whether or not the client, server, or python scripts are actually run with mpiexec, mpirun, or whatever is used to start the executable. So you can't rely on MPI being initialized and should not initialize it yourself in your class as that can cause problems as well. I think that's what you meant but wanted to be as clear as possible for the implications for the classes. As far as the check goes for a multiprocess controller; while it does call the parent request data in the serial case, the inherited parallel class overrides some of the methods. Which is a problem because it refers to its own file pointer in those methods (it uses MPI_File instead of a FILE*). Hmm. Sohail On Thu, Aug 23, 2012 at 3:20 PM, Andy Bauer andy.ba...@kitware.comwrote: Hi Sohail, Sorry for the slow reply but I wanted to make sure how things should work first and that took a bit to figure out. Anyways, this is closer to what's needed but the way it is supposed to work in VTK is that we use vtkObjectFactory to determine which one should be created when vtkWindBladeReader::New() is called. This is done at build time based on whether or not VTK is built with MPI. So if VTK and ParaView are built in parallel then there is code added that creates a vtkPWindBladeReader when vtkWindBladeReader::New() is called. This is regardless of whether or not the executable is run with mpi initialized or not. Because of this, the parallel version of the reader is responsible for checking whether mpi is initialized, ideally by doing vtkMultiProcessController::GetGlobalController()-IsA(vtkMPIController) is true, and if it isn't it should just call the parent class's methods for reading in the files. You may want to look at vtkPNrrdReader and the CMakeLists.txt file in that directory to see how they did things. Note though that if it's a single process running that it doesn't do any MPI calls even in the parallel version so your code will be slightly different than theirs. As for other parts of the email below: 1) yes, please keep cc'ing the list as there are some details in here that will be helpful to other people trying to implement their own parallel readers or filters 2) I'm not sure what's going on with the gradient filter but am guessing that it has to do with the multiblock data set and that there probably isn't that same array in each block. But that's just a guess. Did you try extracting the air block and do the vorticity and/or q criterion computation on that? 3) vtkSMReaderFactory may be the class you're looking for to determine which reader is getting used. ParaViewReaders.xml has a list of readers that it will try to read in a file. Andy On Tue, Aug 21, 2012 at 7:11 PM, Sohail Shafii sssha...@ucdavis.eduwrote: sent via another e-mail address; yahoo prevents me from uploading stuff because it's stupid Here's the first version. Do you want me to keep ccing the list? Not sure if that's necessary at this point because this is related to the development of ParaView and not help. I've run it with an older version of the master, and it seems to read in a test .wind file just fine. For now I had to store the vtkPWindBladeReader* files into VTK/Parallel/MPI since I don't have VTK/IO/MPIParallel module yet. Now when I run the gradient filter after the field, the gradient filter thinks that the input array does not have the appropriate number of components for Q criterion or vorticity...but it should (UVW does have three components per tuple). Might be my outdated copy of the git repo. Haven't debugged that gradient code yet. *I
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
Well I guess it's not a problem in the derived functions if I check to see if MPI is available; if not, call the superclass's version of the virtual function. On Tue, Aug 28, 2012 at 11:32 AM, Sohail Shafii sssha...@ucdavis.eduwrote: MPI is only used if it's available. It tries to see if there is a MPIController available. If not, then the number of ranks is 1, which means that the parallel class calls the super class's requestdata (super class is the serial wind blade reader). However, the parallel class overrides certain functions to make use of MPI. So even the parallel instantiation will call the superclass's requestdata, the overriden functions that are called from that point on assume that MPI is available...this is a problem when the MPI controller is not available. On another note, it appears that the rectilinear grid that I am reading (VTK file, structured points I think) is read by the vtkFileSeriesReader, which makes use of the vtkPDataSetReader class. Strange that if I use vtkPDataSetReader myself in C++ code that uses VTK and not paraview, bounds/cell information is not generated. So I think that ParaView does something else besides using that reader...hmmm. On Mon, Aug 27, 2012 at 7:32 PM, Andy Bauer andy.ba...@kitware.comwrote: On Mon, Aug 27, 2012 at 6:53 PM, Sohail Shafii sssha...@ucdavis.eduwrote: I see...so if MPI is running then it will create a parallel version of the class. Otherwise, it will create the normal reader. Technically, if ParaView is built with MPI enabled it will create the parallel version of the class every time whether or not the client, server, or python scripts are actually run with mpiexec, mpirun, or whatever is used to start the executable. So you can't rely on MPI being initialized and should not initialize it yourself in your class as that can cause problems as well. I think that's what you meant but wanted to be as clear as possible for the implications for the classes. As far as the check goes for a multiprocess controller; while it does call the parent request data in the serial case, the inherited parallel class overrides some of the methods. Which is a problem because it refers to its own file pointer in those methods (it uses MPI_File instead of a FILE*). Hmm. Sohail On Thu, Aug 23, 2012 at 3:20 PM, Andy Bauer andy.ba...@kitware.comwrote: Hi Sohail, Sorry for the slow reply but I wanted to make sure how things should work first and that took a bit to figure out. Anyways, this is closer to what's needed but the way it is supposed to work in VTK is that we use vtkObjectFactory to determine which one should be created when vtkWindBladeReader::New() is called. This is done at build time based on whether or not VTK is built with MPI. So if VTK and ParaView are built in parallel then there is code added that creates a vtkPWindBladeReader when vtkWindBladeReader::New() is called. This is regardless of whether or not the executable is run with mpi initialized or not. Because of this, the parallel version of the reader is responsible for checking whether mpi is initialized, ideally by doing vtkMultiProcessController::GetGlobalController()-IsA(vtkMPIController) is true, and if it isn't it should just call the parent class's methods for reading in the files. You may want to look at vtkPNrrdReader and the CMakeLists.txt file in that directory to see how they did things. Note though that if it's a single process running that it doesn't do any MPI calls even in the parallel version so your code will be slightly different than theirs. As for other parts of the email below: 1) yes, please keep cc'ing the list as there are some details in here that will be helpful to other people trying to implement their own parallel readers or filters 2) I'm not sure what's going on with the gradient filter but am guessing that it has to do with the multiblock data set and that there probably isn't that same array in each block. But that's just a guess. Did you try extracting the air block and do the vorticity and/or q criterion computation on that? 3) vtkSMReaderFactory may be the class you're looking for to determine which reader is getting used. ParaViewReaders.xml has a list of readers that it will try to read in a file. Andy On Tue, Aug 21, 2012 at 7:11 PM, Sohail Shafii sssha...@ucdavis.eduwrote: sent via another e-mail address; yahoo prevents me from uploading stuff because it's stupid Here's the first version. Do you want me to keep ccing the list? Not sure if that's necessary at this point because this is related to the development of ParaView and not help. I've run it with an older version of the master, and it seems to read in a test .wind file just fine. For now I had to store the vtkPWindBladeReader* files into VTK/Parallel/MPI since I don't have VTK/IO/MPIParallel module yet. Now when I run the gradient filter after the field, the gradient filter thinks
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
As far as the gradient filter is concerned: I pick the UVW array in the GUI. It does compute a gradient, but only for one array -- A-Scale turbulence (all variables are loaded). So picking a different array in the GUI doesn't seem to make a difference as it keeps computing gradient for that one input (point data) array. That's why it can't compute vorticity or q-criterion because a-scale turbulence is a 1-component array. Might be because I'm using an older version of the repo from late July. Sohail On Tue, Aug 28, 2012 at 11:51 AM, Sohail Shafii sssha...@ucdavis.eduwrote: Well I guess it's not a problem in the derived functions if I check to see if MPI is available; if not, call the superclass's version of the virtual function. On Tue, Aug 28, 2012 at 11:32 AM, Sohail Shafii sssha...@ucdavis.eduwrote: MPI is only used if it's available. It tries to see if there is a MPIController available. If not, then the number of ranks is 1, which means that the parallel class calls the super class's requestdata (super class is the serial wind blade reader). However, the parallel class overrides certain functions to make use of MPI. So even the parallel instantiation will call the superclass's requestdata, the overriden functions that are called from that point on assume that MPI is available...this is a problem when the MPI controller is not available. On another note, it appears that the rectilinear grid that I am reading (VTK file, structured points I think) is read by the vtkFileSeriesReader, which makes use of the vtkPDataSetReader class. Strange that if I use vtkPDataSetReader myself in C++ code that uses VTK and not paraview, bounds/cell information is not generated. So I think that ParaView does something else besides using that reader...hmmm. On Mon, Aug 27, 2012 at 7:32 PM, Andy Bauer andy.ba...@kitware.comwrote: On Mon, Aug 27, 2012 at 6:53 PM, Sohail Shafii sssha...@ucdavis.eduwrote: I see...so if MPI is running then it will create a parallel version of the class. Otherwise, it will create the normal reader. Technically, if ParaView is built with MPI enabled it will create the parallel version of the class every time whether or not the client, server, or python scripts are actually run with mpiexec, mpirun, or whatever is used to start the executable. So you can't rely on MPI being initialized and should not initialize it yourself in your class as that can cause problems as well. I think that's what you meant but wanted to be as clear as possible for the implications for the classes. As far as the check goes for a multiprocess controller; while it does call the parent request data in the serial case, the inherited parallel class overrides some of the methods. Which is a problem because it refers to its own file pointer in those methods (it uses MPI_File instead of a FILE*). Hmm. Sohail On Thu, Aug 23, 2012 at 3:20 PM, Andy Bauer andy.ba...@kitware.comwrote: Hi Sohail, Sorry for the slow reply but I wanted to make sure how things should work first and that took a bit to figure out. Anyways, this is closer to what's needed but the way it is supposed to work in VTK is that we use vtkObjectFactory to determine which one should be created when vtkWindBladeReader::New() is called. This is done at build time based on whether or not VTK is built with MPI. So if VTK and ParaView are built in parallel then there is code added that creates a vtkPWindBladeReader when vtkWindBladeReader::New() is called. This is regardless of whether or not the executable is run with mpi initialized or not. Because of this, the parallel version of the reader is responsible for checking whether mpi is initialized, ideally by doing vtkMultiProcessController::GetGlobalController()-IsA(vtkMPIController) is true, and if it isn't it should just call the parent class's methods for reading in the files. You may want to look at vtkPNrrdReader and the CMakeLists.txt file in that directory to see how they did things. Note though that if it's a single process running that it doesn't do any MPI calls even in the parallel version so your code will be slightly different than theirs. As for other parts of the email below: 1) yes, please keep cc'ing the list as there are some details in here that will be helpful to other people trying to implement their own parallel readers or filters 2) I'm not sure what's going on with the gradient filter but am guessing that it has to do with the multiblock data set and that there probably isn't that same array in each block. But that's just a guess. Did you try extracting the air block and do the vorticity and/or q criterion computation on that? 3) vtkSMReaderFactory may be the class you're looking for to determine which reader is getting used. ParaViewReaders.xml has a list of readers that it will try to read in a file. Andy On Tue, Aug 21, 2012 at 7:11 PM, Sohail Shafii sssha...@ucdavis.eduwrote
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
I see...so if MPI is running then it will create a parallel version of the class. Otherwise, it will create the normal reader. As far as the check goes for a multiprocess controller; while it does call the parent request data in the serial case, the inherited parallel class overrides some of the methods. Which is a problem because it refers to its own file pointer in those methods (it uses MPI_File instead of a FILE*). Hmm. Sohail On Thu, Aug 23, 2012 at 3:20 PM, Andy Bauer andy.ba...@kitware.com wrote: Hi Sohail, Sorry for the slow reply but I wanted to make sure how things should work first and that took a bit to figure out. Anyways, this is closer to what's needed but the way it is supposed to work in VTK is that we use vtkObjectFactory to determine which one should be created when vtkWindBladeReader::New() is called. This is done at build time based on whether or not VTK is built with MPI. So if VTK and ParaView are built in parallel then there is code added that creates a vtkPWindBladeReader when vtkWindBladeReader::New() is called. This is regardless of whether or not the executable is run with mpi initialized or not. Because of this, the parallel version of the reader is responsible for checking whether mpi is initialized, ideally by doing vtkMultiProcessController::GetGlobalController()-IsA(vtkMPIController) is true, and if it isn't it should just call the parent class's methods for reading in the files. You may want to look at vtkPNrrdReader and the CMakeLists.txt file in that directory to see how they did things. Note though that if it's a single process running that it doesn't do any MPI calls even in the parallel version so your code will be slightly different than theirs. As for other parts of the email below: 1) yes, please keep cc'ing the list as there are some details in here that will be helpful to other people trying to implement their own parallel readers or filters 2) I'm not sure what's going on with the gradient filter but am guessing that it has to do with the multiblock data set and that there probably isn't that same array in each block. But that's just a guess. Did you try extracting the air block and do the vorticity and/or q criterion computation on that? 3) vtkSMReaderFactory may be the class you're looking for to determine which reader is getting used. ParaViewReaders.xml has a list of readers that it will try to read in a file. Andy On Tue, Aug 21, 2012 at 7:11 PM, Sohail Shafii sssha...@ucdavis.eduwrote: sent via another e-mail address; yahoo prevents me from uploading stuff because it's stupid Here's the first version. Do you want me to keep ccing the list? Not sure if that's necessary at this point because this is related to the development of ParaView and not help. I've run it with an older version of the master, and it seems to read in a test .wind file just fine. For now I had to store the vtkPWindBladeReader* files into VTK/Parallel/MPI since I don't have VTK/IO/MPIParallel module yet. Now when I run the gradient filter after the field, the gradient filter thinks that the input array does not have the appropriate number of components for Q criterion or vorticity...but it should (UVW does have three components per tuple). Might be my outdated copy of the git repo. Haven't debugged that gradient code yet. *I was wondering if this is the type of code that you want? Or should the classes be designed differently?* Also, I have only run it parallel because if I try to instantiate it inside of the python shell in serial it tries to run RequestInformation (which assumes that the file path exists but does not). It's a work-in-progress but it's a start. A random question: if ParaView reads a file, where would I put a breakpoint that tells me which reader it ends up picking?? I have some separate VTK code that reads in a rectilinear grid, but does not generate cell/bounds information. ParaView does generate information about the bounds though. Sohail -- *From:* David E DeMarle dave.dema...@kitware.com *To:* Sohail Shafii sohailsha...@yahoo.com *Cc:* Andy Bauer andy.ba...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov *Sent:* Friday, July 20, 2012 10:02 AM *Subject:* Re: [Paraview] ParaView thinks that it does not have a wind blade reader Thanks. In today's git master it now sits in VTK/IO/Parallel/vtkWindbladeReader. The serial version should go in VTK/IO/Geometry. The MPI version should go in a new module VTK/IO/MPIParallel. (VTK/IO/Parallel depends only on vtkSocketController now and should not depend on vtkMPIController). Finally when that is done we can remove the flag in paraview's xml that prevents it from being used in serial. David E DeMarle Kitware, Inc. RD Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909
Re: [Paraview] Numpy masking (via programm filter) not quite working in parallel
Thanks. I'll look into it and let you know if I have any questions. Sohail From: George Zagaris george.zaga...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: David Thompson david.thomp...@kitware.com; paraview@paraview.org paraview@paraview.org Sent: Monday, August 20, 2012 7:36 PM Subject: Re: [Paraview] Numpy masking (via programm filter) not quite working in parallel Hello, I have done some work in VTK towards that. Currently, this works only with structured data. If your data is structured, you can use vtkPStructuredGridConnectivity to generate a bitmask (an unsigned char associated with each mesh node). The bitmask can then be queried using the vtkGhostArray (See http://www.vtk.org/doc/nightly/html/classvtkGhostArray.html) to see which properties are turned on for the given mesh node. If the IGNORE property is ON, that means that the node is a duplicate with respect to the given process and should not be taken into account. The implementation assigns ownership to the lowest rank of a mesh node -- on all other processes, the IGNORE property is set to ON. TestPStructuredGridConnectivity illustrates how this can be used. The test ensures that: (1) Ensuring that the total number of nodes computed is invariant to the number of processes/partitions (2) Ensuring that computing the average for a given field is also invariant to the numer of processes/partitions If this could be useful in your work, I can elaborate more if you have any questions on how to use it. Best, George On Fri, Aug 17, 2012 at 6:18 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Very well. I think that should be fine for now. Sohail From: David Thompson david.thomp...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: David Thompson david.thomp...@kitware.com; paraview@paraview.org paraview@paraview.org Sent: Friday, August 17, 2012 4:15 PM Subject: Re: [Paraview] Numpy masking (via programm filter) not quite working in parallel Hi Sohail, ... Will this point sharing information become available in later versions of ParaView? i.e. would it ever be easy to identify and count duplicates? I don't know of any plans to implement this (but that doesn't mean there aren't). David ___ 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
[Paraview] Numpy masking (via programm filter) not quite working in parallel
Hi, I need to use Numpy in a lot of the programmable filters that I write, and I've run into differences in how its masking feature works in serial and parallel. Masking allows one to filter out portions of an array that do not pass some condition. As an example, I've created a stock paraview wavelet, and saved it as a pvd file. I then load it in, and run this inside of a programmable filter: --- import numpy data = inputs[0].PointData['RTData'] # create a mask that tells us which points are equal to one mask = numpy.ma.masked_equal(data, 1) # filter data array by the mask conditions (so that other points are excluded) maskedPnts = numpy.extract(mask, data) print len(maskedPnts) --- In serial mode, I get 9261 points. With two processes, I get 2 x 4851 or 9702. So masking always produces more points. Any ideas to why that is? Is there anything I can do/print out to see why masking doesn't quite work in parallel? Thanks, Sohail ___ 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] Numpy masking (via programm filter) not quite working in parallel
Hi, I guess I forgot about the data duplication issue. For most statistics that I calculate in parallel, the basic results are not affected that much (i.e. mean, standard dev). But for the Pearson correlation coefficient, which is a bit more complicated, there is a noticeable difference (0.833 for parallel and 0.99 for serial). Will this point sharing information become available in later versions of ParaView? i.e. would it ever be easy to identify and count duplicates? Thanks, Sohail From: David Thompson david.thomp...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org paraview@paraview.org Sent: Friday, August 17, 2012 3:40 PM Subject: Re: [Paraview] Numpy masking (via programm filter) not quite working in parallel Hi Sohail, This is likely caused by points shared on several processes. While ParaView splits the cells of a mesh across processes, cells on the boundary between processes share vertices. Thus if a vertex bounds cells split across 3 processes, that vertex will appear in 3 different lists of local vertices and thus be counted 3 times instead of once. There is nothing in ParaView for determining which vertices are shared (and by how many processes). This also affects the statistics filters when run in parallel on point-centered data. For most large data, the number of vertices on inter-process boundaries is small compared to the total so we documented the behavior for the statistics filters but did not implement a solution (because the expected skew is small). David On Aug 17, 2012, at 5:05 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I need to use Numpy in a lot of the programmable filters that I write, and I've run into differences in how its masking feature works in serial and parallel. Masking allows one to filter out portions of an array that do not pass some condition. As an example, I've created a stock paraview wavelet, and saved it as a pvd file. I then load it in, and run this inside of a programmable filter: --- import numpy data = inputs[0].PointData['RTData'] # create a mask that tells us which points are equal to one mask = numpy.ma.masked_equal(data, 1) # filter data array by the mask conditions (so that other points are excluded) maskedPnts = numpy.extract(mask, data) print len(maskedPnts) --- In serial mode, I get 9261 points. With two processes, I get 2 x 4851 or 9702. So masking always produces more points. Any ideas to why that is? Is there anything I can do/print out to see why masking doesn't quite work in parallel? Thanks, Sohail ___ 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] Numpy masking (via programm filter) not quite working in parallel
Very well. I think that should be fine for now. Sohail From: David Thompson david.thomp...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: David Thompson david.thomp...@kitware.com; paraview@paraview.org paraview@paraview.org Sent: Friday, August 17, 2012 4:15 PM Subject: Re: [Paraview] Numpy masking (via programm filter) not quite working in parallel Hi Sohail, ... Will this point sharing information become available in later versions of ParaView? i.e. would it ever be easy to identify and count duplicates? I don't know of any plans to implement this (but that doesn't mean there aren't). David___ 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] ParaView thinks that it does not have a wind blade reader
When running serially, vtkSMReaderFactory::CanReadFile sees that the session has a single process. As such, on line 660 it realizes that the source proxy's process support (in this case the source is a vtkWindBladeReader proxy) is multiple processes, which causes a return of false from that function. This makes ParaView think that the wind blade reader can't read a .wind file in serial. According to Andy, the reader works find in parallel. I assume that a fix for this is pretty trivial? It never did this before in previous versions of ParaView and I think that the wind blade reader code has remained mostly the same. Andy mentioned some details in person when he was talking with Jon but at that time I didn't know as much as I do now. There is a mention below about enabling the filter in serial and parallel -- is this done in a xml file? Sohail From: David E DeMarle dave.dema...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: Andy Bauer andy.ba...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov Sent: Tuesday, July 17, 2012 7:00 PM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader Each file that has an association with that file extension, made via the XML, has its CanReadFile method called. Put a breakpoint there to determine if ParaView even knows about that reader. On Tuesday, July 17, 2012, Sohail Shafii wrote: Sorry, but I meant to ask a different question. Here's what Andy said before, which relates to my concern: I went through the code and it jogged my memory about previous problems with it. It had an explicit MPI call in the constructor that would cause problems when paraview wasn't getting built with MPI. Any time a user tried to open an unknown file it would try to create a windblade reader to see if the reader could handle the file and this would cause a crash. If you launch a separate pvserver the wind blade reader will be available in the gui. What I meant was; which ParaView class creates or tests a reader to see if it can handle a file? The test code you specified has a predefined data set that is loaded in, and I think is more of an example to your average user. Correct me if I'm wrong. So basically ParaView will recognize the file extension, tries to load it, crashes somewhere, then asks you specify the data format from a dialog. I'm trying to figure out where it crashes so that I can see how that relates to the wind blade reader class and any MPI calls that are in there. Sohail Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader kargad:ParaView demarle$ find . -name *WindBlade* -print --- ./VTK/IO/Parallel/Testing/Cxx/TestWindBladeReader.cxx ./VTK/IO/Parallel/vtkWindBladeReader.cxx ./VTK/IO/Parallel/vtkWindBladeReader.h David E DeMarle Kitware, Inc. RD Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909 On Tue, Jul 17, 2012 at 2:09 PM, Sohail Shafii sohailsha...@yahoo.com wrote: What is the class that is responsible for testing to reader to see if it can open the .wind file? Sohail From: Andy Bauer andy.ba...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: Woodring, Jonathan L woodr...@lanl.gov; Sebastien Jourdain sebastien.jourd...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org Sent: Wednesday, July 11, 2012 10:53 AM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader I was reluctant to enable the reader in serial (i.e. not using MPI IO which I think is getting done by default with the parallel build) because of the following warnings that I got: Warning: In /home/acbauer/CODE/ParaView/ ParaView/VTK/IO/Parallel/vtkWindBladeReader.cxx, line 942 vtkWindBladeReader (0x3a0c800): WindBladeReader error reading file: /home/acbauer/DATA/VTKLargeData/Data/WindBladeReader/test1_topo.wind Premature EOF while reading block of data. Expected 896000 but got 0 I just assumed that it was getting bad results so I didn't check the serial output with the parallel output. In either case the code should be changed before enabling it in serial so that it reads the files correctly and doesn't produce any warning messages. Andy On Wed, Jul 11, 2012 at 12:04 PM, Sohail Shafii sohailsha...@yahoo.com wrote: I'm not aware of any option to build/link against mpi-io. This is just a standard build with mpi support that suffers from this issue. The ParaView website's builds also have the same bug. From: Woodring, Jonathan L woodr...@lanl.gov To: Sohail Shafii sohailsha...@yahoo.com; Andy Bauer andy.ba...@kitware.com; Sebastien Jourdain sebastien.jourd...@kitware.com Cc: migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
Well it seems like #1 and #2 can be solved by writing a new parallel wind blade reader class. I had similar experience by creating distributed versions of class that I created in other vtk directories. The normal wind blade reader would have no parallel calls, which would address #1. What I find interesting is that in your directories, the wind blade reader is located in VTK/IO/Parallel. For the ParaView snapshot that I downloaded, as well as the 3.14 source, it's still in VTK/Parallel. No Parallel folder in VTK/IO. Is this a recent change? So maybe the ideal folder layout would be: VTK/IO/WindBladeReader.* VTK/IO/Parallel/DistributedWindBladeReader.* I can make changes to the WindBladeReader to remove the MPI-IO calls and make a distributed wind blade reader class in case you want to add them to the repository. The one template I have for the parallel/serial split is the vtkStreamTracer and its distributed counterparts so I can use those as a guide. If you have better recommendations please let me know. Sohail From: David E DeMarle dave.dema...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: Andy Bauer andy.ba...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov Sent: Thursday, July 19, 2012 3:20 PM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader Two issues need to be fixed. 1) that behavior in 3.14 is there because of this fix which we had to rush in just before the 3.14 release. Committer: Utkarsh Ayachit utkarsh.ayac...@kitware.com 2012-02-08 17:12:18 Parent: e401d569f6e54c49c3497852463f6b13a01da8b6 (Update version after 3.14.0-RC1 tag.) Child: 99fdb1aeb5997c29b43393f45a478c4b53a23109 (Merge topic 'fix_windblade_segfault') Branches: many (31) Follows: v3.14.0-RC1 Precedes: v3.14.0-RC2 Fix segfault with ParaView built with MPI when opening *.vtr files. The problem was that vtkWindBlade reader was getting created to test CanReadFile and that would cause segfault unless MPI was initialized. Solved the problem by added support to vtkSMReaderFactory to respect the ProcessSupport flag on SourceProxy and change windblade reader XML to say that it only works in parallel. In other words, Windblade reader's Constructor should not take it upon itself to initialize MPI. Once that is fixed, we get to issue 2 before paraview 4.0 2) VTK's new modularization has the rule that no module can behave differently with parallel (or most other cmake config options) turned configured on or off. This is so because it makes it very hard to package VTK when the libs behave differently in different situations. Ie, we can't have any more #ifdef VTK_USE_MPI do something #else do something different #endif in VTK classes. Now if we want a class to behave differently when MPI is turned on or off we have to make a serial only version in a non MPI dependent module and a subclass of that in an MPI dependent module. David E DeMarle Kitware, Inc. RD Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909 On Thu, Jul 19, 2012 at 3:52 PM, Sohail Shafii sohailsha...@yahoo.com wrote: When running serially, vtkSMReaderFactory::CanReadFile sees that the session has a single process. As such, on line 660 it realizes that the source proxy's process support (in this case the source is a vtkWindBladeReader proxy) is multiple processes, which causes a return of false from that function. This makes ParaView think that the wind blade reader can't read a .wind file in serial. According to Andy, the reader works find in parallel. I assume that a fix for this is pretty trivial? It never did this before in previous versions of ParaView and I think that the wind blade reader code has remained mostly the same. Andy mentioned some details in person when he was talking with Jon but at that time I didn't know as much as I do now. There is a mention below about enabling the filter in serial and parallel -- is this done in a xml file? Sohail From: David E DeMarle dave.dema...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: Andy Bauer andy.ba...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov Sent: Tuesday, July 17, 2012 7:00 PM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader Each file that has an association with that file extension, made via the XML, has its CanReadFile method called. Put a breakpoint there to determine if ParaView even knows about that reader. On Tuesday, July 17, 2012, Sohail Shafii wrote: Sorry, but I meant to ask a different question. Here's what Andy said before, which relates to my concern: I went through the code and it jogged my memory about previous problems with it. It had an explicit
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
What is the class that is responsible for testing to reader to see if it can open the .wind file? Sohail From: Andy Bauer andy.ba...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: Woodring, Jonathan L woodr...@lanl.gov; Sebastien Jourdain sebastien.jourd...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org Sent: Wednesday, July 11, 2012 10:53 AM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader I was reluctant to enable the reader in serial (i.e. not using MPI IO which I think is getting done by default with the parallel build) because of the following warnings that I got: Warning: In /home/acbauer/CODE/ParaView/ ParaView/VTK/IO/Parallel/vtkWindBladeReader.cxx, line 942 vtkWindBladeReader (0x3a0c800): WindBladeReader error reading file: /home/acbauer/DATA/VTKLargeData/Data/WindBladeReader/test1_topo.wind Premature EOF while reading block of data. Expected 896000 but got 0 I just assumed that it was getting bad results so I didn't check the serial output with the parallel output. In either case the code should be changed before enabling it in serial so that it reads the files correctly and doesn't produce any warning messages. Andy On Wed, Jul 11, 2012 at 12:04 PM, Sohail Shafii sohailsha...@yahoo.com wrote: I'm not aware of any option to build/link against mpi-io. This is just a standard build with mpi support that suffers from this issue. The ParaView website's builds also have the same bug. From: Woodring, Jonathan L woodr...@lanl.gov To: Sohail Shafii sohailsha...@yahoo.com; Andy Bauer andy.ba...@kitware.com; Sebastien Jourdain sebastien.jourd...@kitware.com Cc: migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org Sent: Tuesday, July 10, 2012 10:06 AM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader I added the MPI calls to be in sync with what was being added in the simulation. In the SuperComputing demo last year, we updated both the HiGrad FireTec to use MPI-IO and likewise ParaView to do the same thing. It was working when I last checked it in November. It could be that you didn't build or link against MPI-IO? It might require that the build needs to be guarded if you don't have MPI-IO. I thought I did that but it might have gotten changed during the VTK overhaul… That's the only thing that I can think of at the moment, because I'm pretty sure I had it working in serial with MPI-IO (it shouldn't matter.) From: Sohail Shafii sohailsha...@yahoo.com Reply-To: Sohail Shafii sohailsha...@yahoo.com Date: Tue, 10 Jul 2012 08:44:01 -0700 To: Andy Bauer andy.ba...@kitware.com, Sebastien Jourdain sebastien.jourd...@kitware.com Cc: migic...@gmail.com migic...@gmail.com, woodring woodr...@lanl.gov, paraview@paraview.org paraview@paraview.org Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader I'm afraid my knowledge here is a bit limited. All of my edits have been in relation to loading the blade data. The MPI calls were put in by someone else. By the way, is the wind blade reader updated in VTK as well? (i.e. the VTK package) I remember downloading VTK a while ago and there was a mismatch between that windbladereader and the one in ParaView. Sohail From: Andy Bauer andy.ba...@kitware.com To: Sebastien Jourdain sebastien.jourd...@kitware.com Cc: Sohail Shafii sohailsha...@yahoo.com; migic...@gmail.com migic...@gmail.com; woodr...@lanl.gov woodr...@lanl.gov; paraview@paraview.org paraview@paraview.org Sent: Tuesday, July 3, 2012 11:54 AM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader Hi Sohail, I went through the code and it jogged my memory about previous problems with it. It had an explicit MPI call in the constructor that would cause problems when paraview wasn't getting built with MPI. Any time a user tried to open an unknown file it would try to create a windblade reader to see if the reader could handle the file and this would cause a crash. If you launch a separate pvserver the wind blade reader will be available in the gui. If someone fixes this problem in the reader we can enable it properly through the gui. Right now though I'm getting a bunch of warnings when running in serial like: Warning: In /home/acbauer/CODE/ParaView/ParaView/VTK/IO/Parallel/vtkWindBladeReader.cxx, line 942 vtkWindBladeReader (0x3a0c800): WindBladeReader error reading file: /home/acbauer/DATA/VTKLargeData/Data/WindBladeReader/test1_topo.wind Premature EOF while reading block of data. Expected 896000 but got 0 Andy On Tue, Jul 3, 2012 at 12:09 PM, Andy Bauer andy.ba...@kitware.com wrote: I get the same problem in both 3.14.1 and master. The reader is still available but just not conveniently through
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
Sorry, but I meant to ask a different question. Here's what Andy said before, which relates to my concern: I went through the code and it jogged my memory about previous problems with it. It had an explicit MPI call in the constructor that would cause problems when paraview wasn't getting built with MPI. Any time a user tried to open an unknown file it would try to create a windblade reader to see if the reader could handle the file and this would cause a crash. If you launch a separate pvserver the wind blade reader will be available in the gui. What I meant was; which ParaView class creates or tests a reader to see if it can handle a file? The test code you specified has a predefined data set that is loaded in, and I think is more of an example to your average user. Correct me if I'm wrong. So basically ParaView will recognize the file extension, tries to load it, crashes somewhere, then asks you specify the data format from a dialog. I'm trying to figure out where it crashes so that I can see how that relates to the wind blade reader class and any MPI calls that are in there. Sohail From: David E DeMarle dave.dema...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: Andy Bauer andy.ba...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org; Woodring, Jonathan L woodr...@lanl.gov Sent: Tuesday, July 17, 2012 12:43 PM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader kargad:ParaView demarle$ find . -name *WindBlade* -print --- ./VTK/IO/Parallel/Testing/Cxx/TestWindBladeReader.cxx ./VTK/IO/Parallel/vtkWindBladeReader.cxx ./VTK/IO/Parallel/vtkWindBladeReader.h David E DeMarle Kitware, Inc. RD Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909 On Tue, Jul 17, 2012 at 2:09 PM, Sohail Shafii sohailsha...@yahoo.com wrote: What is the class that is responsible for testing to reader to see if it can open the .wind file? Sohail From: Andy Bauer andy.ba...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: Woodring, Jonathan L woodr...@lanl.gov; Sebastien Jourdain sebastien.jourd...@kitware.com; migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org Sent: Wednesday, July 11, 2012 10:53 AM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader I was reluctant to enable the reader in serial (i.e. not using MPI IO which I think is getting done by default with the parallel build) because of the following warnings that I got: Warning: In /home/acbauer/CODE/ParaView/ ParaView/VTK/IO/Parallel/vtkWindBladeReader.cxx, line 942 vtkWindBladeReader (0x3a0c800): WindBladeReader error reading file: /home/acbauer/DATA/VTKLargeData/Data/WindBladeReader/test1_topo.wind Premature EOF while reading block of data. Expected 896000 but got 0 I just assumed that it was getting bad results so I didn't check the serial output with the parallel output. In either case the code should be changed before enabling it in serial so that it reads the files correctly and doesn't produce any warning messages. Andy On Wed, Jul 11, 2012 at 12:04 PM, Sohail Shafii sohailsha...@yahoo.com wrote: I'm not aware of any option to build/link against mpi-io. This is just a standard build with mpi support that suffers from this issue. The ParaView website's builds also have the same bug. From: Woodring, Jonathan L woodr...@lanl.gov To: Sohail Shafii sohailsha...@yahoo.com; Andy Bauer andy.ba...@kitware.com; Sebastien Jourdain sebastien.jourd...@kitware.com Cc: migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org Sent: Tuesday, July 10, 2012 10:06 AM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader I added the MPI calls to be in sync with what was being added in the simulation. In the SuperComputing demo last year, we updated both the HiGrad FireTec to use MPI-IO and likewise ParaView to do the same thing. It was working when I last checked it in November. It could be that you didn't build or link against MPI-IO? It might require that the build needs to be guarded if you don't have MPI-IO. I thought I did that but it might have gotten changed during the VTK overhaul… That's the only thing that I can think of at the moment, because I'm pretty sure I had it working in serial with MPI-IO (it shouldn't matter.) From: Sohail Shafii sohailsha...@yahoo.com Reply-To: Sohail Shafii sohailsha...@yahoo.com Date: Tue, 10 Jul 2012 08:44:01 -0700 To: Andy Bauer andy.ba...@kitware.com, Sebastien Jourdain sebastien.jourd...@kitware.com Cc: migic...@gmail.com migic...@gmail.com, woodring woodr...@lanl.gov, paraview@paraview.org paraview@paraview.org Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader I'm afraid
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
I'm not aware of any option to build/link against mpi-io. This is just a standard build with mpi support that suffers from this issue. The ParaView website's builds also have the same bug. From: Woodring, Jonathan L woodr...@lanl.gov To: Sohail Shafii sohailsha...@yahoo.com; Andy Bauer andy.ba...@kitware.com; Sebastien Jourdain sebastien.jourd...@kitware.com Cc: migic...@gmail.com migic...@gmail.com; paraview@paraview.org paraview@paraview.org Sent: Tuesday, July 10, 2012 10:06 AM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader I added the MPI calls to be in sync with what was being added in the simulation. In the SuperComputing demo last year, we updated both the HiGrad FireTec to use MPI-IO and likewise ParaView to do the same thing. It was working when I last checked it in November. It could be that you didn't build or link against MPI-IO? It might require that the build needs to be guarded if you don't have MPI-IO. I thought I did that but it might have gotten changed during the VTK overhaul… That's the only thing that I can think of at the moment, because I'm pretty sure I had it working in serial with MPI-IO (it shouldn't matter.) From: Sohail Shafii sohailsha...@yahoo.com Reply-To: Sohail Shafii sohailsha...@yahoo.com Date: Tue, 10 Jul 2012 08:44:01 -0700 To: Andy Bauer andy.ba...@kitware.com, Sebastien Jourdain sebastien.jourd...@kitware.com Cc: migic...@gmail.com migic...@gmail.com, woodring woodr...@lanl.gov, paraview@paraview.org paraview@paraview.org Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader I'm afraid my knowledge here is a bit limited. All of my edits have been in relation to loading the blade data. The MPI calls were put in by someone else. By the way, is the wind blade reader updated in VTK as well? (i.e. the VTK package) I remember downloading VTK a while ago and there was a mismatch between that windbladereader and the one in ParaView. Sohail From: Andy Bauer andy.ba...@kitware.com To: Sebastien Jourdain sebastien.jourd...@kitware.com Cc: Sohail Shafii sohailsha...@yahoo.com; migic...@gmail.com migic...@gmail.com; woodr...@lanl.gov woodr...@lanl.gov; paraview@paraview.org paraview@paraview.org Sent: Tuesday, July 3, 2012 11:54 AM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader Hi Sohail, I went through the code and it jogged my memory about previous problems with it. It had an explicit MPI call in the constructor that would cause problems when paraview wasn't getting built with MPI. Any time a user tried to open an unknown file it would try to create a windblade reader to see if the reader could handle the file and this would cause a crash. If you launch a separate pvserver the wind blade reader will be available in the gui. If someone fixes this problem in the reader we can enable it properly through the gui. Right now though I'm getting a bunch of warnings when running in serial like: Warning: In /home/acbauer/CODE/ParaView/ParaView/VTK/IO/Parallel/vtkWindBladeReader.cxx, line 942 vtkWindBladeReader (0x3a0c800): WindBladeReader error reading file: /home/acbauer/DATA/VTKLargeData/Data/WindBladeReader/test1_topo.wind Premature EOF while reading block of data. Expected 896000 but got 0 Andy On Tue, Jul 3, 2012 at 12:09 PM, Andy Bauer andy.ba...@kitware.com wrote: I get the same problem in both 3.14.1 and master. The reader is still available but just not conveniently through the GUI. To access it though you can bring up the python shell and do the following: w = WindBladereader() w.Filename = path/test1_topo.wind Show() Render() This shouldn't be too hard of a fix for ParaView. I'll let you know when it's done. Andy On Mon, Jul 2, 2012 at 6:03 PM, Sebastien Jourdain sebastien.jourd...@kitware.com wrote: Hi Sohail, a similar bug was fixed in master not so long ago. Could you try master to see if the issue remains ? Thanks, Seb On Mon, Jul 2, 2012 at 5:36 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, We are trying to open a .wind data set. In the open file menu, ParaView lists the .wind format as a possible option. Once we click open, however, it complains that a reader for that data format cannot be found. This is true for the compiled-from-source and pre-compiled versions of ParaView 3.14.1. The wind blade source still exists in the VTK/Parallel folder and it is referenced in the corresponding CMakeLists.txt file. Any thoughts? Thanks, Sohail ___ 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
Re: [Paraview] ParaView thinks that it does not have a wind blade reader
Jimmy said that it works. Thanks. Sohail From: Andy Bauer andy.ba...@kitware.com To: Sebastien Jourdain sebastien.jourd...@kitware.com Cc: Sohail Shafii sohailsha...@yahoo.com; migic...@gmail.com migic...@gmail.com; woodr...@lanl.gov woodr...@lanl.gov; paraview@paraview.org paraview@paraview.org Sent: Tuesday, July 3, 2012 10:09 AM Subject: Re: [Paraview] ParaView thinks that it does not have a wind blade reader I get the same problem in both 3.14.1 and master. The reader is still available but just not conveniently through the GUI. To access it though you can bring up the python shell and do the following: w = WindBladereader() w.Filename = path/test1_topo.wind Show() Render() This shouldn't be too hard of a fix for ParaView. I'll let you know when it's done. Andy On Mon, Jul 2, 2012 at 6:03 PM, Sebastien Jourdain sebastien.jourd...@kitware.com wrote: Hi Sohail, a similar bug was fixed in master not so long ago. Could you try master to see if the issue remains ? Thanks, Seb On Mon, Jul 2, 2012 at 5:36 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, We are trying to open a .wind data set. In the open file menu, ParaView lists the .wind format as a possible option. Once we click open, however, it complains that a reader for that data format cannot be found. This is true for the compiled-from-source and pre-compiled versions of ParaView 3.14.1. The wind blade source still exists in the VTK/Parallel folder and it is referenced in the corresponding CMakeLists.txt file. Any thoughts? Thanks, Sohail ___ 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
[Paraview] ParaView thinks that it does not have a wind blade reader
Hi, We are trying to open a .wind data set. In the open file menu, ParaView lists the .wind format as a possible option. Once we click open, however, it complains that a reader for that data format cannot be found. This is true for the compiled-from-source and pre-compiled versions of ParaView 3.14.1. The wind blade source still exists in the VTK/Parallel folder and it is referenced in the corresponding CMakeLists.txt file. Any thoughts? Thanks, Sohail ___ 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] Compile error for vtkWrapPython (locally installed python)
Hi, I am currently compiling paraview 3.14.1 and using a locally installed version of python2.7. My version of cmake is 2.8.8. The OS is Red Hat Enterprise Linux Client release 5.8. I get the following error during compilation: [ 0%] Built target vtkWrapPythonInit make[2]: *** No rule to make target `/home/***/local/lib/libpython2.7.a', needed by `bin/vtkWrapPython'. Stop. make[1]: *** [VTK/Wrapping/CMakeFiles/vtkWrapPython.dir/all] Error 2 make[1]: *** Waiting for unfinished jobs [ 0%] Built target vtkWrapHierarchy [ 0%] Built target vtksys [ 1%] Built target vtkzlib make: *** [all] Error 2 With verbose one, it fails after this line: make -f VTK/Wrapping/CMakeFiles/vtkWrapPython.dir/build.make VTK/Wrapping/CMakeFiles/vtkWrapPython.dir/build In the build.make file, it has the following under vtkWrapPython_EXTERNAL_OBJECTS = bin/vtkWrapPython: VTK/Wrapping/CMakeFiles/vtkWrapPython.dir/build.make bin/vtkWrapPython: /home/***/local/lib/libpython2.7.a ... Any ideas? Thanks, Sohail ___ 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] vtkDataSet algorithm that takes composite data set as input
This is my custom print self: Debug: Off Modified Time: 122500 Reference Count: 1 Registered Events: (none) Request: REQUEST_DATA_OBJECT FORWARD_DIRECTION: 0 ALGORITHM_AFTER_FORWARD: 1 FROM_OUTPUT_PORT: 0 For some odd reason gdb does not allow me to step through vtkDataSetAlgorithm's ProcessRequest. But it does return an error, as the executive prints the following: vtkPVCompositeDataPipeline (0x1dd7640): Algorithm vtkFastMarcher(0x1e2e9d0) returned failure for request: vtkInformation (0x2d0d590) Debug: Off Modified Time: 276195 Reference Count: 1 Registered Events: (none) Request: REQUEST_DATA_OBJECT FORWARD_DIRECTION: 0 FROM_OUTPUT_PORT: 0 ALGORITHM_AFTER_FORWARD: 1 Debug: Off Modified Time: 276205 Reference Count: 1 Registered Events: (none) Request: REQUEST_DATA_OBJECT FORWARD_DIRECTION: 0 FROM_OUTPUT_PORT: 0 ALGORITHM_AFTER_FORWARD: 1 It appears to fail in the superclass's ProcessRequest here: // create the output if(request-Has(vtkDemandDrivenPipeline::REQUEST_DATA_OBJECT())) { return this-RequestDataObject(request, inputVector, outputVector); } Basically RequestDataObject can't cast the inInfo into a vtkDataSet, because it's a composite. I guess I have to override it? Sohail From: David E DeMarle dave.dema...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org paraview@paraview.org Sent: Monday, January 23, 2012 7:25 AM Subject: Re: [Paraview] vtkDataSet algorithm that takes composite data set as input I'ld sit down with a debugger and see which request exactly is failing and why. Try overriding vtkAlgorithm::ProcessRequest() in your filter so that you can easily get a breakpoint in an instance of your class. In that method call PrintSelf(cerr, vtkIndent(0)) on the request vtkInformation. Put the breakpoint right after that. After the breakpoint call this-Superclass::ProcessRequest(). When the PrintSelf tells you which request fails, rerun and step in until you see why the handling method returns. David E DeMarle Kitware, Inc. RD Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909 On Fri, Jan 20, 2012 at 3:34 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Thanks. I get these sorts of errors: vtkPVCompositeDataPipeline (0x2a831a0): Algorithm vtkFastMarcher(0x2a80860) returned failure for request: vtkInformation (0x339b0f0 Before RequestData is reached... Sohail From: David E DeMarle dave.dema...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org paraview@paraview.org Sent: Friday, January 20, 2012 11:42 AM Subject: Re: [Paraview] vtkDataSet algorithm that takes composite data set as input Try something like this. In this case the mapper takes in either type, which you may or may not want. // int vtkCompositePolyDataMapper2::FillInputPortInformation( int vtkNotUsed(port), vtkInformation* info) { info-Set(vtkAlgorithm::INPUT_REQUIRED_DATA_TYPE(), vtkPolyData); info-Append(vtkAlgorithm::INPUT_REQUIRED_DATA_TYPE(), vtkCompositeDataSet); return 1; } David E DeMarle Kitware, Inc. RD Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909 On Fri, Jan 20, 2012 at 2:36 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I am developing an algorithm that requires a composite data set as input while creating a vtkDataSet as output -- basically it inherit from vtkDataSetAlgorithm. I know that one must provide a new definition to the FillInputPortInformation(int port, vtkInformation *info) function, but I'm not sure how. Please advise, thank you. Sohail ___ 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
[Paraview] vtkDataSet algorithm that takes composite data set as input
Hi, I am developing an algorithm that requires a composite data set as input while creating a vtkDataSet as output -- basically it inherit from vtkDataSetAlgorithm. I know that one must provide a new definition to the FillInputPortInformation(int port, vtkInformation *info) function, but I'm not sure how. Please advise, thank you. Sohail___ 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] vtkDataSet algorithm that takes composite data set as input
Thanks. I get these sorts of errors: vtkPVCompositeDataPipeline (0x2a831a0): Algorithm vtkFastMarcher(0x2a80860) returned failure for request: vtkInformation (0x339b0f0 Before RequestData is reached... Sohail From: David E DeMarle dave.dema...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org paraview@paraview.org Sent: Friday, January 20, 2012 11:42 AM Subject: Re: [Paraview] vtkDataSet algorithm that takes composite data set as input Try something like this. In this case the mapper takes in either type, which you may or may not want. // int vtkCompositePolyDataMapper2::FillInputPortInformation( int vtkNotUsed(port), vtkInformation* info) { info-Set(vtkAlgorithm::INPUT_REQUIRED_DATA_TYPE(), vtkPolyData); info-Append(vtkAlgorithm::INPUT_REQUIRED_DATA_TYPE(), vtkCompositeDataSet); return 1; } David E DeMarle Kitware, Inc. RD Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909 On Fri, Jan 20, 2012 at 2:36 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I am developing an algorithm that requires a composite data set as input while creating a vtkDataSet as output -- basically it inherit from vtkDataSetAlgorithm. I know that one must provide a new definition to the FillInputPortInformation(int port, vtkInformation *info) function, but I'm not sure how. Please advise, thank you. Sohail ___ 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] Data does not load off of a partition from saved (.wind) state
The 555 MB data set would probably suffice. I could debug this problem with my data set. In the state file, I see a few places where UVW is mentioned: Property name=SelectionPointFieldDataArrayName id=227.SelectionPointFieldDataArrayName number_of_elements=1 Element index=0 value=vtkOriginalPointIds/ Domain name=array_list id=227.SelectionPointFieldDataArrayName.array_list String text=Density/ String text=UVW/ String text=Vorticity/ /Domain /Property and: Property name=PointArrayInfo id=190.PointArrayInfo number_of_elements=52 Element index=48 value=UVW/ Element index=49 value=1/ and: Property name=PointArrayStatus id=190.PointArrayStatus number_of_elements=52 ... Element index=48 value=UVW/ Element index=49 value=1/ and: Domain name=array_list id=190.PointArrayStatus.array_list ... String text=UVW/ So I would think that it would get loaded in. I guess I could debug the part of the ParaView could that reads the state file...because once the wind blade reader starts doing processing it doesn't even read in UVW. Sohail From: David E DeMarle dave.dema...@kitware.com To: Utkarsh Ayachit utkarsh.ayac...@kitware.com Cc: Sohail Shafii sohailsha...@yahoo.com Sent: Monday, November 21, 2011 7:36 AM Subject: Re: [Paraview] Data does not load off of a partition from saved (.wind) state Probably not. I have the 555MB kitware.wind dataset from Pat Fasel back in 3/2010 and the 124MB test1_topo.wind data from Sohail (which is in VTKLargeData) here. David E DeMarle Kitware, Inc. RD Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909 On Mon, Nov 21, 2011 at 10:28 AM, Utkarsh Ayachit utkarsh.ayac...@kitware.com wrote: Dave, Do you have this dataset? Utkarsh On Sun, Nov 20, 2011 at 8:21 PM, Sohail Shafii sohailsha...@yahoo.com wrote: It appears that this problem also happens when I put the data on the root partition. For some reason when ParaView restores the state, wind blade reader thinks it's being told to load in the Density array, and not the UVW array as I specified. I can't seem to load in the UVW array after paraview starts either...this only happens when trying to restore a state for this data set (i.e. not when I open the data set without loading a state). This doesn't happen with other .wind data sets, only this new fairly large one (which is about 600 MB per time step). Sohail From: Sohail Shafii sohailsha...@yahoo.com To: paraview@paraview.org paraview@paraview.org Sent: Thursday, November 17, 2011 2:52 PM Subject: Data does not load off of a partition from saved (.wind) state Hi, I've noticed that if I open a data set from a separate partition (after mounting it in Linux), and if I save a state file related to said data set, the data does not load at all when I load the statefile if after I start ParaView from a clean state. Please note that it doesn't matter if the state file is saved on my root partition or on the separate partition. Also, this is specific to the wind blade reader. I've debugged it and it appears to be reading the data when loading the saved state...for some reason that information does not get passed to the ParaView pipeline. The end result is that I have a Field, Blade and Ground ports except all three have zero points and cells. This is for the release branch of ParaView. Any thoughts? Thanks. Sohail ___ 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] Data does not load off of a partition from saved (.wind) state
Addendum: What might help is to add another filter after the Field portion, like do a Gradient of Unstructured Data Sets and select UVW as the input. When you reload the state file you might notice that something goes wrong (i.e. UVW is not read in and an error message pops up). If that doesn't happen then there might be some problem with my paraview installation or my new data set. Sohail From: Sohail Shafii sohailsha...@yahoo.com To: David E DeMarle dave.dema...@kitware.com; Utkarsh Ayachit utkarsh.ayac...@kitware.com Cc: paraview@paraview.org paraview@paraview.org Sent: Monday, November 21, 2011 10:04 AM Subject: Re: [Paraview] Data does not load off of a partition from saved (.wind) state The 555 MB data set would probably suffice. I could debug this problem with my data set. In the state file, I see a few places where UVW is mentioned: Property name=SelectionPointFieldDataArrayName id=227.SelectionPointFieldDataArrayName number_of_elements=1 Element index=0 value=vtkOriginalPointIds/ Domain name=array_list id=227.SelectionPointFieldDataArrayName.array_list String text=Density/ String text=UVW/ String text=Vorticity/ /Domain /Property and: Property name=PointArrayInfo id=190.PointArrayInfo number_of_elements=52 Element index=48 value=UVW/ Element index=49 value=1/ and: Property name=PointArrayStatus id=190.PointArrayStatus number_of_elements=52 ... Element index=48 value=UVW/ Element index=49 value=1/ and: Domain name=array_list id=190.PointArrayStatus.array_list ... String text=UVW/ So I would think that it would get loaded in. I guess I could debug the part of the ParaView could that reads the state file...because once the wind blade reader starts doing processing it doesn't even read in UVW. Sohail From: David E DeMarle dave.dema...@kitware.com To: Utkarsh Ayachit utkarsh.ayac...@kitware.com Cc: Sohail Shafii sohailsha...@yahoo.com Sent: Monday, November 21, 2011 7:36 AM Subject: Re: [Paraview] Data does not load off of a partition from saved (.wind) state Probably not. I have the 555MB kitware.wind dataset from Pat Fasel back in 3/2010 and the 124MB test1_topo.wind data from Sohail (which is in VTKLargeData) here. David E DeMarle Kitware, Inc. RD Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909 On Mon, Nov 21, 2011 at 10:28 AM, Utkarsh Ayachit utkarsh.ayac...@kitware.com wrote: Dave, Do you have this dataset? Utkarsh On Sun, Nov 20, 2011 at 8:21 PM, Sohail Shafii sohailsha...@yahoo.com wrote: It appears that this problem also happens when I put the data on the root partition. For some reason when ParaView restores the state, wind blade reader thinks it's being told to load in the Density array, and not the UVW array as I specified. I can't seem to load in the UVW array after paraview starts either...this only happens when trying to restore a state for this data set (i.e. not when I open the data set without loading a state). This doesn't happen with other .wind data sets, only this new fairly large one (which is about 600 MB per time step). Sohail From: Sohail Shafii sohailsha...@yahoo.com To: paraview@paraview.org paraview@paraview.org Sent: Thursday, November 17, 2011 2:52 PM Subject: Data does not load off of a partition from saved (.wind) state Hi, I've noticed that if I open a data set from a separate partition (after mounting it in Linux), and if I save a state file related to said data set, the data does not load at all when I load the statefile if after I start ParaView from a clean state. Please note that it doesn't matter if the state file is saved on my root partition or on the separate partition. Also, this is specific to the wind blade reader. I've debugged it and it appears to be reading the data when loading the saved state...for some reason that information does not get passed to the ParaView pipeline. The end result is that I have a Field, Blade and Ground ports except all three have zero points and cells. This is for the release branch of ParaView. Any thoughts? Thanks. Sohail ___ 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
Re: [Paraview] Data does not load off of a partition from saved (.wind) state
It appears that this problem also happens when I put the data on the root partition. For some reason when ParaView restores the state, wind blade reader thinks it's being told to load in the Density array, and not the UVW array as I specified. I can't seem to load in the UVW array after paraview starts either...this only happens when trying to restore a state for this data set (i.e. not when I open the data set without loading a state). This doesn't happen with other .wind data sets, only this new fairly large one (which is about 600 MB per time step). Sohail From: Sohail Shafii sohailsha...@yahoo.com To: paraview@paraview.org paraview@paraview.org Sent: Thursday, November 17, 2011 2:52 PM Subject: Data does not load off of a partition from saved (.wind) state Hi, I've noticed that if I open a data set from a separate partition (after mounting it in Linux), and if I save a state file related to said data set, the data does not load at all when I load the statefile if after I start ParaView from a clean state. Please note that it doesn't matter if the state file is saved on my root partition or on the separate partition. Also, this is specific to the wind blade reader. I've debugged it and it appears to be reading the data when loading the saved state...for some reason that information does not get passed to the ParaView pipeline. The end result is that I have a Field, Blade and Ground ports except all three have zero points and cells. This is for the release branch of ParaView. Any thoughts? Thanks. Sohail___ 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] Data does not load off of a partition from saved (.wind) state
Hi, I've noticed that if I open a data set from a separate partition (after mounting it in Linux), and if I save a state file related to said data set, the data does not load at all when I load the statefile if after I start ParaView from a clean state. Please note that it doesn't matter if the state file is saved on my root partition or on the separate partition. Also, this is specific to the wind blade reader. I've debugged it and it appears to be reading the data when loading the saved state...for some reason that information does not get passed to the ParaView pipeline. The end result is that I have a Field, Blade and Ground ports except all three have zero points and cells. This is for the release branch of ParaView. Any thoughts? Thanks. Sohail ___ 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] Linking error when trying to use vtkGLSL/vtkXMLShader
Thanks. That seems to do the trick. The only other issue is the rendering context. I've used the line integral convolution/lic filters as an example, and it seems like I need a SetContext function, which is passed a pointer to an (new) allocated vtkRenderWindow object. For reference, I am looking at vtkImageDataLIC2D.cxx. I will then need to downcast this object to a vtkOpenGLRenderWindow and pass it to vtkShaderProgram2, which requires a proper context...well I've done all of that, and so far I've noticed that a new (small, blacked out) window with the title Visualization Toolkit pops up when I start my filter. Now I have to do something with this fragment shader (which currently does nothing). Am I on the right track here? Or I am issuing an improper context? Sohail From: Utkarsh Ayachit utkarsh.ayac...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org paraview@paraview.org Sent: Wednesday, October 26, 2011 8:43 AM Subject: Re: [Paraview] Linking error when trying to use vtkGLSL/vtkXMLShader Classes in graphics cannot depend on classes in rendering. You can move your filter to hybrid. Utkarsh On Oct 26, 2011, at 12:11 AM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I am currently trying to do some shader operations in a standard C++ filter. This filter (inherits from vtkMultiBlockDataSetAlgorithm, resides in VTK/Graphics) will be used to voxelize a polygonal object (basically the result of the contour filter) and write the results to a 3D texture. I tried to do the following: vtkXMLShader *xmlShader = vtkXMLShader::New(); vtkGLSLShader *glslShader = vtkGLSLShader::New(); Which results in: Undefined symbols for architecture x86_64: vtkXMLShader::New(), referenced from: vtkRepulsiveForceFieldGenerator::RequestData(vtkInformation*, vtkInformationVector**, vtkInformationVector*)in vtkRepulsiveForceFieldGenerator.cxx.o vtkGLSLShader::New(), referenced from: vtkRepulsiveForceFieldGenerator::RequestData(vtkInformation*, vtkInformationVector**, vtkInformationVector*)in vtkRepulsiveForceFieldGenerator.cxx.o ld: symbol(s) not found for architecture x86_64 collect2: ld returned 1 exit status I'm going to guess that standard filters can't use code like this, and that my filter should reside inside of the VTK/Rendering folder and follow the examples in that folder. Thanks, Sohail ___ 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
[Paraview] Linking error when trying to use vtkGLSL/vtkXMLShader
Hi, I am currently trying to do some shader operations in a standard C++ filter. This filter (inherits from vtkMultiBlockDataSetAlgorithm, resides in VTK/Graphics) will be used to voxelize a polygonal object (basically the result of the contour filter) and write the results to a 3D texture. I tried to do the following: vtkXMLShader *xmlShader = vtkXMLShader::New(); vtkGLSLShader *glslShader = vtkGLSLShader::New(); Which results in: Undefined symbols for architecture x86_64: vtkXMLShader::New(), referenced from: vtkRepulsiveForceFieldGenerator::RequestData(vtkInformation*, vtkInformationVector**, vtkInformationVector*)in vtkRepulsiveForceFieldGenerator.cxx.o vtkGLSLShader::New(), referenced from: vtkRepulsiveForceFieldGenerator::RequestData(vtkInformation*, vtkInformationVector**, vtkInformationVector*)in vtkRepulsiveForceFieldGenerator.cxx.o ld: symbol(s) not found for architecture x86_64 collect2: ld returned 1 exit status I'm going to guess that standard filters can't use code like this, and that my filter should reside inside of the VTK/Rendering folder and follow the examples in that folder. Thanks, Sohail___ 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] Multiple polydata in same pipeline -- don't always seem centered around same origin
Thanks. Changing the setting Coincident Topology Resolution to Do Nothing made it work. Sohail From: David E DeMarle dave.dema...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org paraview@paraview.org Sent: Wednesday, October 12, 2011 1:02 PM Subject: Re: [Paraview] Multiple polydata in same pipeline -- don't always seem centered around same origin ParaView draws points and lines slightly offset so that you see wireframes on top of surfaces. Try playing with Edit-Settings-Render View-General-Coincident Topology Resolution to change this behavior. David E DeMarle Kitware, Inc. RD Engineer 21 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-881-4909 On Wed, Oct 12, 2011 at 3:34 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I have been writing some code that basically computes points on certain portions of line segments (kind of line intersection points), and I can't help but notice that when I zoom in, these points appear disconnected from lines the should lie on top of. To illustrate my point, I have created a simple example using stock paraview filters (see attachments). In centroid_zoomed_out, I basically create a line using the line source filter, then using a programmable filter I compute its centroid and create a vertex. When I render both pieces of polydata, the result looks correct when one is zoomed out. But zoomed-in (centroid_zoomed_in), it appears to be disconnected from the line. Is there any reason for this? I have never rendered points and lines together in paraview, but I would expect them to exist around the same origin. Thanks, Sohail ___ 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] Accessing OpenGL from standard VTK filter
This is personal code, by maybe one day I will contribute it to VTK, as like a plugin or something There are pieces of code that special things that do access OpenGL directly however, like the LIC code. But I could do things another way...I so see classes related to shaders (vtkShader) and setting up the view volume (vtkPersp*), so those classes may come in handy Sohail From: Andy Bauer andy.ba...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org paraview@paraview.org Sent: Wednesday, October 5, 2011 8:56 AM Subject: Re: [Paraview] Accessing OpenGL from standard VTK filter Are you hoping that this will make its way into the main VTK repo or is this personal code that you plan on maintaining yourself and plan on linking with VTK? There are people that know more about this than me but this sounds like a bad idea for the main VTK repo. A lot of the OpenGL specific calls are hidden by factories which create them as needed. For example, you'd create a vtkPolyDataMapper and the factory would create the vtkOpenGLPolyDataMapper, vtkMesaPolyDataMapper, etc as needed by the build. Andy On Tue, Oct 4, 2011 at 6:34 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, If I were to write a typical VTK filter C++ (i.e. the type with VTKRequestData, not existing in the VTK/Rendering folder), would it be fine to have OpenGL calls in there? I'm guessing based on the way VTK is designed, most classes that access OpenGL directly might exist in the Rendering folder. This is more of a design question anyway and relates to how to VTK classes are designed. My class will add data structures to the pipeline but it may need a fragment shader for computations. Sohail ___ 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
[Paraview] Accessing OpenGL from standard VTK filter
Hi, If I were to write a typical VTK filter C++ (i.e. the type with VTKRequestData, not existing in the VTK/Rendering folder), would it be fine to have OpenGL calls in there? I'm guessing based on the way VTK is designed, most classes that access OpenGL directly might exist in the Rendering folder. This is more of a design question anyway and relates to how to VTK classes are designed. My class will add data structures to the pipeline but it may need a fragment shader for computations. Sohail___ 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] Comparative view and wind blade reader
Thank you. I can wait for this fix to be updated in the release, I guess. Sohail From: Utkarsh Ayachit utkarsh.ayac...@kitware.com To: Andy Bauer andy.ba...@kitware.com Cc: Sohail Shafii sohailsha...@yahoo.com; paraview@paraview.org paraview@paraview.org Sent: Tuesday, August 30, 2011 8:38 AM Subject: Re: [Paraview] Comparative view and wind blade reader Sohail, I've tracked this issue down. Attached is the patch. I'll get this fix in 3.12. Utkarsh On Tue, Aug 9, 2011 at 2:30 PM, Andy Bauer andy.ba...@kitware.com wrote: I forgot to mention that that one only has 2 time steps. I have a sample one with 4 time steps too if you need it. On Tue, Aug 9, 2011 at 2:28 PM, Andy Bauer andy.ba...@kitware.com wrote: Yes, it's in VTKLargeData/Data/WindBladeReader. Select the test1_topo.wind file. Andy On Tue, Aug 9, 2011 at 2:06 PM, Sohail Shafii sohailsha...@yahoo.com wrote: I believe Andy should have at least one lying around. Sohail From: Utkarsh Ayachit utkarsh.ayac...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org paraview@paraview.org Sent: Tuesday, August 9, 2011 11:07 AM Subject: Re: [Paraview] Comparative view and wind blade reader Ah, it's possible that the view has issues with multiple output ports, I don't think I ever verified that. I forget, is there any small dataset I can test this with already in ParaView (or where can I get it) ? Utkarsh On Tue, Aug 9, 2011 at 11:39 AM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I have recently become interested in comparing different time steps using the windblade reader. One approach involves a 3D comparative view. After loading my .wind data set, I created a 2X1 Comparative View. In the Comparative View Inspector, I clicked the + symbol when Time was selected. I can view the Field portion without problems in either view (with different timesteps for each view), but when I select Blade, the first view renders it correctly while the second view renders a surface around the data set. When I select Ground, the first view renders the item correctly while the second renders the data set's outline. I wonder if Comparative view has issues with multiple ports. Any ideas? Thank you, Sohail ___ 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] Unstructured to structured...
Thanks. I had a hunch that this was the answer. In my case, the surface would represent values that smaller than the requested iso-value, so inner cells would be below, while bordering cells would have points that are equal to and above. I needed to do the conversion because the paper specifically makes use of the iso-surface boundary in 3D space. I suppose I could have a multi-block data set with two data sets (using Group data sets): one unstructured (contour) and another structured (original). I could find out where the points of the iso-surface exist in the structured grid, although that part will take some thought. Sohail From: Moreland, Kenneth kmo...@sandia.gov To: Sohail Shafii sohailsha...@yahoo.com; paraview@paraview.org paraview@paraview.org Sent: Thursday, August 11, 2011 4:32 PM Subject: Re: [Paraview] Unstructured to structured... In general, no. It is not possible to convert an unstructured grid to a structured grid because there are structures representable with unstructured grids that are non-representable with structured grids. For your question, it would be much easier to get the information you want by just looking at the scalar value on the original structured grid. You don't need to look at the isosurface. For example, all cells inside the isosurface will have all scalar values either above or below the isovalue (depending on the orientation of the surface). Bordering cells have scalar values both above and below the isovalue. -Ken Kenneth Moreland *** Sandia National Laboratories *** *** *** *** email: kmo...@sandia.gov ** *** ** phone: (505) 844-8919 *** web: http://www.cs.unm.edu/~kmorel From: Sohail Shafii sohailsha...@yahoo.com Reply-To: Sohail Shafii sohailsha...@yahoo.com Date: Thu, 11 Aug 2011 14:54:18 -0700 To: paraview@paraview.org paraview@paraview.org Subject: [Paraview] Unstructured to structured... Hi, I've been wondering: would it be possible to take an unstructured grid using existing filters and make it into structured? Here's why: I compute an iso-surface for a data set (resulting surface is unstructured subset of original structured data set), which needs to be compared against the cells of the original structured grid. Specifically, I need to find out which cells of the original data set border the iso-surface, and which are inside. All filters that do appending/combination of the contour and original (structured) data set create unstructured grids. Thanks, Sohail ___ 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
[Paraview] Unstructured to structured...
Hi, I've been wondering: would it be possible to take an unstructured grid using existing filters and make it into structured? Here's why: I compute an iso-surface for a data set (resulting surface is unstructured subset of original structured data set), which needs to be compared against the cells of the original structured grid. Specifically, I need to find out which cells of the original data set border the iso-surface, and which are inside. All filters that do appending/combination of the contour and original (structured) data set create unstructured grids. Thanks, Sohail___ 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] Comparative view and wind blade reader
Hi, I have recently become interested in comparing different time steps using the windblade reader. One approach involves a 3D comparative view. After loading my .wind data set, I created a 2X1 Comparative View. In the Comparative View Inspector, I clicked the + symbol when Time was selected. I can view the Field portion without problems in either view (with different timesteps for each view), but when I select Blade, the first view renders it correctly while the second view renders a surface around the data set. When I select Ground, the first view renders the item correctly while the second renders the data set's outline. I wonder if Comparative view has issues with multiple ports. Any ideas? Thank you, Sohail___ 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] Comparative view and wind blade reader
I believe Andy should have at least one lying around. Sohail From: Utkarsh Ayachit utkarsh.ayac...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org paraview@paraview.org Sent: Tuesday, August 9, 2011 11:07 AM Subject: Re: [Paraview] Comparative view and wind blade reader Ah, it's possible that the view has issues with multiple output ports, I don't think I ever verified that. I forget, is there any small dataset I can test this with already in ParaView (or where can I get it) ? Utkarsh On Tue, Aug 9, 2011 at 11:39 AM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I have recently become interested in comparing different time steps using the windblade reader. One approach involves a 3D comparative view. After loading my .wind data set, I created a 2X1 Comparative View. In the Comparative View Inspector, I clicked the + symbol when Time was selected. I can view the Field portion without problems in either view (with different timesteps for each view), but when I select Blade, the first view renders it correctly while the second view renders a surface around the data set. When I select Ground, the first view renders the item correctly while the second renders the data set's outline. I wonder if Comparative view has issues with multiple ports. Any ideas? Thank you, Sohail ___ 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] Accessing field and wind blade data (together) in python filter
Thank you for the suggestions. Since the data sets were on separate ports, here's what I did: I did a re-sample with data set to map the input (field) to the source (blade). But the blade cell data items were missing in the result, so I then did a Group data set with the result of the mapping and the blade. This means that I do not have to make any changes in the vtkWindBladeReader code. Sohail From: David E DeMarle dave.dema...@kitware.com To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org paraview@paraview.org Sent: Thursday, July 28, 2011 9:25 AM Subject: Re: [Paraview] Accessing field and wind blade data (together) in python filter One thought that comes to mind is the Resample with DataSet filter, which is a filter which has the effect of merging pipeline branches together. Set the Input to be the field and the Source to be the blade and the result will be a copy of the blade with the data values of the field sampled onto the points of the blade. David E DeMarle Kitware, Inc. RD Engineer 28 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-371-3971 x109 On Tue, Jul 26, 2011 at 6:29 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, Right now the field and blade data are separate for the wind blade reader (as expected). I was wondering if, after creating a python filter below the blade object in the pipeline, I could somehow access the field data? I'm guessing not based on some experimenting. I need to do this because I have to correlate points on the blade to the same areas in the field to calculate certain quantities (which requires interpolated field values). If the python route is not possible then I can just make changes in vtkWindBladeReader instead. Sohail ___ 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
[Paraview] Accessing field and wind blade data (together) in python filter
Hi, Right now the field and blade data are separate for the wind blade reader (as expected). I was wondering if, after creating a python filter below the blade object in the pipeline, I could somehow access the field data? I'm guessing not based on some experimenting. I need to do this because I have to correlate points on the blade to the same areas in the field to calculate certain quantities (which requires interpolated field values). If the python route is not possible then I can just make changes in vtkWindBladeReader instead. Sohail___ 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] Possible issue with velocity interpolator class
Hi, I've been developer a vortex core tracer, which is basically a customized version of vtkStreamTracer. I use the interpolation code quite heavily in the class, in this test case I use the point interpolator (i.e. vtkInterpolatedVelocityField class). Now in one part of the code, I've noticed that doing two calls to the interpolator's FunctionValues for a bad boundary point produces strange results. After the first call to FunctionValues I obtain a valid cell (using GetCell) and after the second call I get a NULL cell. The second call is correct because the point in question's y coordinate is -0.0012773338068285955, which outside of the data set's 0-boundary. Anyway, I will look into this further to see what's up; I was wondering if anybody has ideas. Thank you, Sohail Shafii___ 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] Interesting bug when using windbladereader and multiple time steps
By current development version are you referring to the master branch? Sohail --- On Wed, 7/13/11, Andy Bauer andy.ba...@kitware.com wrote: From: Andy Bauer andy.ba...@kitware.com Subject: Re: [Paraview] Interesting bug when using windbladereader and multiple time steps To: Sohail Shafii sohailsha...@yahoo.com Cc: Patricia KFasel p...@lanl.gov, paraview@paraview.org Date: Wednesday, July 13, 2011, 8:41 AM I'm not sure why this isn't working for you. I just tried it on my machine with 2 processes and was able to use the gradient of unstructured grid filter on the field grid. Can you try it out with the current development version of paraview to see if you can reproduce it there? Andy On Mon, Jul 11, 2011 at 9:38 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, There is actually another bug (which may or may not be related) which did not happen before with ParaView 3.8.0 but appears with the release branch of the ParaView that I am using. If I load a wind data set after connecting to 2-process localhost (mpirun -np 2 pvserver), the gradient of unstructured grid appears greyed out in the filters menu so I can't use that (which is necessary for me). If I run mpirun -np 1 pvserver, this is not an issue. It's only an issue with 2 processes or above. This doesn't happen if I load in a *foam file (i.e. openFOAM file); only with .wind and the most up-to-date release version of ParaView as opposed to 3.8.0. Sohail --- On Fri, 7/8/11, Andy Bauer andy.ba...@kitware.com wrote: From: Andy Bauer andy.ba...@kitware.com Subject: Re: [Paraview] Interesting bug when using windbladereader and multiple time steps To: Sohail Shafii sohailsha...@yahoo.com Cc: Patricia KFasel p...@lanl.gov, paraview@paraview.org Date: Friday, July 8, 2011, 10:22 AM 3 or 4 time steps should be fine. Thanks, Andy On Fri, Jul 8, 2011 at 1:18 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Thank you for looking into this. I only modified the windbladereader to handle a slight format change, so my knowledge here is a bit limited (in regards to its original design) but I'll help out when necessary. From what I see there are two outputs (one field, one blade). I guess the nice thing about having a separate output for the blade is that I can render it separately from the field. I can send you a data set with more timesteps (similar to what you are using now, actually), problem is my DSL connection has a poor upload. Would three or four timesteps do? Sohail --- On Fri, 7/8/11, Andy Bauer andy.ba...@kitware.com wrote: From: Andy Bauer andy.ba...@kitware.com Subject: Re: [Paraview] Interesting bug when using windbladereader and multiple time steps To: Sohail Shafii sohailsha...@yahoo.com, Fasel, Patricia K p...@lanl.gov Cc: paraview@paraview.org Date: Friday, July 8, 2011, 9:59 AM Ok, I'm able to replicate the bug. The test file I have (http://vtk.org/gitweb?p=VTKLargeData.git;a=tree;f=Data/WindBladeReader;h=3459e8ded59eb05796e922d90200bd4e83231df0;hb=HEAD) only has 2 time steps so I may need another set of files from you with more time steps. Playing around with the reader in ParaView, I think it may be worse than you realize. Only the time information for the field is getting updated properly and if you go back in time steps, the blade time step either still increases or stays the same if it has hit the last time step. I need to talk to our temporal pipeline expert to figure out the proper behavior. One question though -- is there a reason that the reader doesn't return a multiblock with the 3 grids in it instead of having 3 separate output ports? I'm thinking that if the output was changed to a multiblock then these time issues would go away. Andy On Thu, Jul 7, 2011 at 10:28 PM, Sohail Shafii sohailsha...@yahoo.com wrote: I believe that data set went you (related to the updates I made for the wind blade reader)...the one with ten timesteps, should do the trick. Do you still have it? Sohail --- On Thu, 7/7/11, Andy Bauer andy.ba...@kitware.com wrote: From: Andy Bauer andy.ba...@kitware.com Subject: Re: [Paraview] Interesting bug when using windbladereader and multiple time steps To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org Date: Thursday, July 7, 2011, 7:09 PM I'll take a look at it. I may need a time dependent data set if I can't figure it out though. Andy On Thu, Jul 7, 2011 at 9:13 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I've noticed a bug when using the windbladereader; the portion of the code that reads the blade file is one timestep behind the code that reads the field. For instance, if one is at timestep i, the field data will correspond with timestep i but the blade data will correspond with the last timestep chosen. When a data set is opened initially both the blade and field items are set to the initial timestep. After one starts changing the timestep, the bug
Re: [Paraview] Interesting bug when using windbladereader and multiple time steps
Hi, There is actually another bug (which may or may not be related) which did not happen before with ParaView 3.8.0 but appears with the release branch of the ParaView that I am using. If I load a wind data set after connecting to 2-process localhost (mpirun -np 2 pvserver), the gradient of unstructured grid appears greyed out in the filters menu so I can't use that (which is necessary for me). If I run mpirun -np 1 pvserver, this is not an issue. It's only an issue with 2 processes or above. This doesn't happen if I load in a *foam file (i.e. openFOAM file); only with .wind and the most up-to-date release version of ParaView as opposed to 3.8.0. Sohail --- On Fri, 7/8/11, Andy Bauer andy.ba...@kitware.com wrote: From: Andy Bauer andy.ba...@kitware.com Subject: Re: [Paraview] Interesting bug when using windbladereader and multiple time steps To: Sohail Shafii sohailsha...@yahoo.com Cc: Patricia KFasel p...@lanl.gov, paraview@paraview.org Date: Friday, July 8, 2011, 10:22 AM 3 or 4 time steps should be fine. Thanks, Andy On Fri, Jul 8, 2011 at 1:18 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Thank you for looking into this. I only modified the windbladereader to handle a slight format change, so my knowledge here is a bit limited (in regards to its original design) but I'll help out when necessary. From what I see there are two outputs (one field, one blade). I guess the nice thing about having a separate output for the blade is that I can render it separately from the field. I can send you a data set with more timesteps (similar to what you are using now, actually), problem is my DSL connection has a poor upload. Would three or four timesteps do? Sohail --- On Fri, 7/8/11, Andy Bauer andy.ba...@kitware.com wrote: From: Andy Bauer andy.ba...@kitware.com Subject: Re: [Paraview] Interesting bug when using windbladereader and multiple time steps To: Sohail Shafii sohailsha...@yahoo.com, Fasel, Patricia K p...@lanl.gov Cc: paraview@paraview.org Date: Friday, July 8, 2011, 9:59 AM Ok, I'm able to replicate the bug. The test file I have (http://vtk.org/gitweb?p=VTKLargeData.git;a=tree;f=Data/WindBladeReader;h=3459e8ded59eb05796e922d90200bd4e83231df0;hb=HEAD) only has 2 time steps so I may need another set of files from you with more time steps. Playing around with the reader in ParaView, I think it may be worse than you realize. Only the time information for the field is getting updated properly and if you go back in time steps, the blade time step either still increases or stays the same if it has hit the last time step. I need to talk to our temporal pipeline expert to figure out the proper behavior. One question though -- is there a reason that the reader doesn't return a multiblock with the 3 grids in it instead of having 3 separate output ports? I'm thinking that if the output was changed to a multiblock then these time issues would go away. Andy On Thu, Jul 7, 2011 at 10:28 PM, Sohail Shafii sohailsha...@yahoo.com wrote: I believe that data set went you (related to the updates I made for the wind blade reader)...the one with ten timesteps, should do the trick. Do you still have it? Sohail --- On Thu, 7/7/11, Andy Bauer andy.ba...@kitware.com wrote: From: Andy Bauer andy.ba...@kitware.com Subject: Re: [Paraview] Interesting bug when using windbladereader and multiple time steps To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org Date: Thursday, July 7, 2011, 7:09 PM I'll take a look at it. I may need a time dependent data set if I can't figure it out though. Andy On Thu, Jul 7, 2011 at 9:13 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I've noticed a bug when using the windbladereader; the portion of the code that reads the blade file is one timestep behind the code that reads the field. For instance, if one is at timestep i, the field data will correspond with timestep i but the blade data will correspond with the last timestep chosen. When a data set is opened initially both the blade and field items are set to the initial timestep. After one starts changing the timestep, the bug makes an appearance. In the code, the blade and field specify their timesteps (individually) based on the bladeInfo and fieldInfo objects, respectively. These info objects are related to the outputVector of requestInformation -- I wonder if this is some kind of pipeline issue? Thanks, Sohail ___ 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
Re: [Paraview] Interesting bug when using windbladereader and multiple time steps
Thank you for looking into this. I only modified the windbladereader to handle a slight format change, so my knowledge here is a bit limited (in regards to its original design) but I'll help out when necessary. From what I see there are two outputs (one field, one blade). I guess the nice thing about having a separate output for the blade is that I can render it separately from the field. I can send you a data set with more timesteps (similar to what you are using now, actually), problem is my DSL connection has a poor upload. Would three or four timesteps do? Sohail --- On Fri, 7/8/11, Andy Bauer andy.ba...@kitware.com wrote: From: Andy Bauer andy.ba...@kitware.com Subject: Re: [Paraview] Interesting bug when using windbladereader and multiple time steps To: Sohail Shafii sohailsha...@yahoo.com, Fasel, Patricia K p...@lanl.gov Cc: paraview@paraview.org Date: Friday, July 8, 2011, 9:59 AM Ok, I'm able to replicate the bug. The test file I have (http://vtk.org/gitweb?p=VTKLargeData.git;a=tree;f=Data/WindBladeReader;h=3459e8ded59eb05796e922d90200bd4e83231df0;hb=HEAD) only has 2 time steps so I may need another set of files from you with more time steps. Playing around with the reader in ParaView, I think it may be worse than you realize. Only the time information for the field is getting updated properly and if you go back in time steps, the blade time step either still increases or stays the same if it has hit the last time step. I need to talk to our temporal pipeline expert to figure out the proper behavior. One question though -- is there a reason that the reader doesn't return a multiblock with the 3 grids in it instead of having 3 separate output ports? I'm thinking that if the output was changed to a multiblock then these time issues would go away. Andy On Thu, Jul 7, 2011 at 10:28 PM, Sohail Shafii sohailsha...@yahoo.com wrote: I believe that data set went you (related to the updates I made for the wind blade reader)...the one with ten timesteps, should do the trick. Do you still have it? Sohail --- On Thu, 7/7/11, Andy Bauer andy.ba...@kitware.com wrote: From: Andy Bauer andy.ba...@kitware.com Subject: Re: [Paraview] Interesting bug when using windbladereader and multiple time steps To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org Date: Thursday, July 7, 2011, 7:09 PM I'll take a look at it. I may need a time dependent data set if I can't figure it out though. Andy On Thu, Jul 7, 2011 at 9:13 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I've noticed a bug when using the windbladereader; the portion of the code that reads the blade file is one timestep behind the code that reads the field. For instance, if one is at timestep i, the field data will correspond with timestep i but the blade data will correspond with the last timestep chosen. When a data set is opened initially both the blade and field items are set to the initial timestep. After one starts changing the timestep, the bug makes an appearance. In the code, the blade and field specify their timesteps (individually) based on the bladeInfo and fieldInfo objects, respectively. These info objects are related to the outputVector of requestInformation -- I wonder if this is some kind of pipeline issue? Thanks, Sohail ___ 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
[Paraview] Interesting bug when using windbladereader and multiple time steps
Hi, I've noticed a bug when using the windbladereader; the portion of the code that reads the blade file is one timestep behind the code that reads the field. For instance, if one is at timestep i, the field data will correspond with timestep i but the blade data will correspond with the last timestep chosen. When a data set is opened initially both the blade and field items are set to the initial timestep. After one starts changing the timestep, the bug makes an appearance. In the code, the blade and field specify their timesteps (individually) based on the bladeInfo and fieldInfo objects, respectively. These info objects are related to the outputVector of requestInformation -- I wonder if this is some kind of pipeline issue? Thanks, Sohail ___ 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] Interesting bug when using windbladereader and multiple time steps
I believe that data set went you (related to the updates I made for the wind blade reader)...the one with ten timesteps, should do the trick. Do you still have it? Sohail --- On Thu, 7/7/11, Andy Bauer andy.ba...@kitware.com wrote: From: Andy Bauer andy.ba...@kitware.com Subject: Re: [Paraview] Interesting bug when using windbladereader and multiple time steps To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org Date: Thursday, July 7, 2011, 7:09 PM I'll take a look at it. I may need a time dependent data set if I can't figure it out though. Andy On Thu, Jul 7, 2011 at 9:13 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I've noticed a bug when using the windbladereader; the portion of the code that reads the blade file is one timestep behind the code that reads the field. For instance, if one is at timestep i, the field data will correspond with timestep i but the blade data will correspond with the last timestep chosen. When a data set is opened initially both the blade and field items are set to the initial timestep. After one starts changing the timestep, the bug makes an appearance. In the code, the blade and field specify their timesteps (individually) based on the bladeInfo and fieldInfo objects, respectively. These info objects are related to the outputVector of requestInformation -- I wonder if this is some kind of pipeline issue? Thanks, Sohail ___ 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] programmable sourece-programmable filter - streamtracer parallel
Have you tried to display each parts of the pipeline when in parallel? In other words, see if you can render something related to the programmable source, then the programmable filter, to see if the problem might occur earlier in the pipeline. Sohail --- On Wed, 4/6/11, Atanas Atanasov atana...@in.tum.de wrote: From: Atanas Atanasov atana...@in.tum.de Subject: [Paraview] programmable sourece-programmable filter - streamtracer parallel To: paraview@paraview.org Date: Wednesday, April 6, 2011, 2:26 PM Hi i have created the following pipeline on c++ based servermanager plugin(without qt) -programmable source(creates a structured grid)-programmable filter(adds velocities and pressure on the vertices) -ArbitrarySourceStreamTracer(Seeding thorugh points or thorugh line).This setup works fine with one process. In parallel the streamlines are not shown. Do you have experience with such kind of problem? Regards, A.Atanasov ___ 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] visualizing particles with different radius
It looks like your v, f and omega arrays (inputs to the glyph filter) cover the same set of points here. Each of them have 117 tuples. You want one glyph filter to apply to all points? Have you tried playing around with the Scale mode? I was able to get all the points to become spheres to become spheres that way. Sohail --- On Thu, 3/24/11, sarai anvari lamm...@gmail.com wrote: From: sarai anvari lamm...@gmail.com Subject: Re: [Paraview] visualizing particles with different radius To: Sohail Shafii sohailsha...@yahoo.com Date: Thursday, March 24, 2011, 8:23 AM thanks but as you sea in this vtk file there two types of atoms with different radius. i have tried applying two glyph filter but it only changes one type each time,,, On Wed, Mar 23, 2011 at 7:05 PM, Sohail Shafii sohailsha...@yahoo.com wrote: So in the glyph filter's Vectors drop-down you see the two types and you want it applied to both? If that's the case then I would add a second glyph filter into the pipeline, and have it applied to the other type. As far as I know each glyph filter can be applied to only one input vector from the pipeline. Sohail --- On Wed, 3/23/11, sarai anvari lamm...@gmail.com wrote: From: sarai anvari lamm...@gmail.com Subject: [Paraview] visualizing particles with different radius To: paraview@paraview.org Date: Wednesday, March 23, 2011, 7:00 AM I have produced *.vtk files which contain information about 6 coordinates in Cartesian space and radius of the particles. the particles have two different type. type 1 with radius 1 and type 2 with radius 5. each time i apply Glyph filter and choose sphere type only one type of the particles is changed. how is it possible to apply this filter to both types or even if i had three type with tree different radii how could i visualize these tree diffident radii? thanks. -Inline Attachment Follows- ___ 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] Getting a raytraced picture from Manta
Thank you...that version is actually way more stable, and I can see that the raytracer is working now. Phong shading does work (I was able to pick that successfully), and it works better than the rasterized Gouraud (as expected). I just don't see any shadows cast created on the polygons my filter created, I think that would be seen by default. Any ideas? Sohail --- On Tue, 3/22/11, David E DeMarle dave.dema...@kitware.com wrote: From: David E DeMarle dave.dema...@kitware.com Subject: Re: [Paraview] Getting a raytraced picture from Manta To: Sebastien Jourdain sebastien.jourd...@kitware.com Cc: Sohail Shafii sohailsha...@yahoo.com, paraview@paraview.org Date: Tuesday, March 22, 2011, 9:11 PM I highly recommend you use Manta r2439. As far as vtkManta's use of Manta goes, that version is the most stable and newer Manta versions are known to cause it to crash. David E DeMarle Kitware, Inc. RD Engineer 28 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-371-3971 x109 On Tue, Mar 22, 2011 at 3:17 PM, Sebastien Jourdain sebastien.jourd...@kitware.com wrote: Hi Sohail, if you are using 3.8.0, it is true that you will need to start a server (pvserver) and connect to it to allow manta to work properly with ParaView. You don't need to run in MPI, but it has to be in a different process than the client. (I mean no built-in mode) On the other hand, once you get connected to that pvserver (in localhost for example) you will have to load the plugin on both side (client+server). Once this is done, you should be able to close the 3D renderer and create a new one based on Manta. For the crash part, the best person will be Dave that I cc. But as he is traveling, he might not be able to answer to you right away. Seb On Tue, Mar 22, 2011 at 4:49 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Thanks. It turns out that I am working with an older version (3.8.0) since I am making preparations for a paper (due next week) and some code hasn't been ported to latest version yet. I basically have to add the plugins in the Remote plugins section o fthe Plugin manager. The plugin's a bit buggy and causes a crash, however, and I think that may be due to the fact that I am using an older version of Manta. Sohail --- On Tue, 3/22/11, Sebastien Jourdain sebastien.jourd...@kitware.com wrote: From: Sebastien Jourdain sebastien.jourd...@kitware.com Subject: Re: [Paraview] Getting a raytraced picture from Manta To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org Date: Tuesday, March 22, 2011, 1:32 PM Hi Sohail, you should just close the default 3D view and create the Manta one by clicking on its button. Then any object that you create will be shown in that raytracer view... And with the latest implementation, you don't need a remote server anymore. Seb On Tue, Mar 22, 2011 at 4:18 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I have been able to get Manta raytracer to compile with paraview and load in as a plugin. I even see the extra GUI options in the display tab of certain filters like the Tube filter. I'm just wondering how to get something raytraced as a result. I did a bit a research via the list and found a page that says that I must start the paraview server with the plugin loaded. The client would also have the plugin loaded and connect to the server. The instructions didn't say how to do that on the server side . Usually when I start the paraview server (pvserver), it does not launch the GUI but just waits for a client connection in the terminal. So I don't know how to load the plugin on the server...the command line options for pvserver didn't provide any hints. Please advise. Thanks, Sohail ___ 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] visualizing particles with different radius
So in the glyph filter's Vectors drop-down you see the two types and you want it applied to both? If that's the case then I would add a second glyph filter into the pipeline, and have it applied to the other type. As far as I know each glyph filter can be applied to only one input vector from the pipeline. Sohail --- On Wed, 3/23/11, sarai anvari lamm...@gmail.com wrote: From: sarai anvari lamm...@gmail.com Subject: [Paraview] visualizing particles with different radius To: paraview@paraview.org Date: Wednesday, March 23, 2011, 7:00 AM I have produced *.vtk files which contain information about 6 coordinates in Cartesian space and radius of the particles. the particles have two different type. type 1 with radius 1 and type 2 with radius 5. each time i apply Glyph filter and choose sphere type only one type of the particles is changed. how is it possible to apply this filter to both types or even if i had three type with tree different radii how could i visualize these tree diffident radii? thanks. -Inline Attachment Follows- ___ 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] Getting a raytraced picture from Manta
Thanks, that works. Sohail --- On Wed, 3/23/11, John Patchett patchett2...@gmail.com wrote: From: John Patchett patchett2...@gmail.com Subject: Re: [Paraview] Getting a raytraced picture from Manta To: Sohail Shafii sohailsha...@yahoo.com Cc: Sebastien Jourdain sebastien.jourd...@kitware.com, David E DeMarle dave.dema...@kitware.com, paraview@paraview.org Date: Wednesday, March 23, 2011, 2:23 PM Hi Sohail, Look in: Edit-Settings-Manta View There is check box for Shadows. --John On Wed, Mar 23, 2011 at 12:21 AM, Sohail Shafii sohailsha...@yahoo.com wrote: Thank you...that version is actually way more stable, and I can see that the raytracer is working now. Phong shading does work (I was able to pick that successfully), and it works better than the rasterized Gouraud (as expected). I just don't see any shadows cast created on the polygons my filter created, I think that would be seen by default. Any ideas? Sohail --- On Tue, 3/22/11, David E DeMarle dave.dema...@kitware.com wrote: From: David E DeMarle dave.dema...@kitware.com Subject: Re: [Paraview] Getting a raytraced picture from Manta To: Sebastien Jourdain sebastien.jourd...@kitware.com Cc: Sohail Shafii sohailsha...@yahoo.com, paraview@paraview.org Date: Tuesday, March 22, 2011, 9:11 PM I highly recommend you use Manta r2439. As far as vtkManta's use of Manta goes, that version is the most stable and newer Manta versions are known to cause it to crash. David E DeMarle Kitware, Inc. RD Engineer 28 Corporate Drive Clifton Park, NY 12065-8662 Phone: 518-371-3971 x109 On Tue, Mar 22, 2011 at 3:17 PM, Sebastien Jourdain sebastien.jourd...@kitware.com wrote: Hi Sohail, if you are using 3.8.0, it is true that you will need to start a server (pvserver) and connect to it to allow manta to work properly with ParaView. You don't need to run in MPI, but it has to be in a different process than the client. (I mean no built-in mode) On the other hand, once you get connected to that pvserver (in localhost for example) you will have to load the plugin on both side (client+server). Once this is done, you should be able to close the 3D renderer and create a new one based on Manta. For the crash part, the best person will be Dave that I cc. But as he is traveling, he might not be able to answer to you right away. Seb On Tue, Mar 22, 2011 at 4:49 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Thanks. It turns out that I am working with an older version (3.8.0) since I am making preparations for a paper (due next week) and some code hasn't been ported to latest version yet. I basically have to add the plugins in the Remote plugins section o fthe Plugin manager. The plugin's a bit buggy and causes a crash, however, and I think that may be due to the fact that I am using an older version of Manta. Sohail --- On Tue, 3/22/11, Sebastien Jourdain sebastien.jourd...@kitware.com wrote: From: Sebastien Jourdain sebastien.jourd...@kitware.com Subject: Re: [Paraview] Getting a raytraced picture from Manta To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org Date: Tuesday, March 22, 2011, 1:32 PM Hi Sohail, you should just close the default 3D view and create the Manta one by clicking on its button. Then any object that you create will be shown in that raytracer view... And with the latest implementation, you don't need a remote server anymore. Seb On Tue, Mar 22, 2011 at 4:18 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I have been able to get Manta raytracer to compile with paraview and load in as a plugin. I even see the extra GUI options in the display tab of certain filters like the Tube filter. I'm just wondering how to get something raytraced as a result. I did a bit a research via the list and found a page that says that I must start the paraview server with the plugin loaded. The client would also have the plugin loaded and connect to the server. The instructions didn't say how to do that on the server side . Usually when I start the paraview server (pvserver), it does not launch the GUI but just waits for a client connection in the terminal. So I don't know how to load the plugin on the server...the command line options for pvserver didn't provide any hints. Please advise. Thanks, Sohail ___ 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
[Paraview] Getting a raytraced picture from Manta
Hi, I have been able to get Manta raytracer to compile with paraview and load in as a plugin. I even see the extra GUI options in the display tab of certain filters like the Tube filter. I'm just wondering how to get something raytraced as a result. I did a bit a research via the list and found a page that says that I must start the paraview server with the plugin loaded. The client would also have the plugin loaded and connect to the server. The instructions didn't say how to do that on the server side . Usually when I start the paraview server (pvserver), it does not launch the GUI but just waits for a client connection in the terminal. So I don't know how to load the plugin on the server...the command line options for pvserver didn't provide any hints. Please advise. Thanks, Sohail ___ 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] Getting a raytraced picture from Manta
Thanks. It turns out that I am working with an older version (3.8.0) since I am making preparations for a paper (due next week) and some code hasn't been ported to latest version yet. I basically have to add the plugins in the Remote plugins section o fthe Plugin manager. The plugin's a bit buggy and causes a crash, however, and I think that may be due to the fact that I am using an older version of Manta. Sohail --- On Tue, 3/22/11, Sebastien Jourdain sebastien.jourd...@kitware.com wrote: From: Sebastien Jourdain sebastien.jourd...@kitware.com Subject: Re: [Paraview] Getting a raytraced picture from Manta To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org Date: Tuesday, March 22, 2011, 1:32 PM Hi Sohail, you should just close the default 3D view and create the Manta one by clicking on its button. Then any object that you create will be shown in that raytracer view... And with the latest implementation, you don't need a remote server anymore. Seb On Tue, Mar 22, 2011 at 4:18 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I have been able to get Manta raytracer to compile with paraview and load in as a plugin. I even see the extra GUI options in the display tab of certain filters like the Tube filter. I'm just wondering how to get something raytraced as a result. I did a bit a research via the list and found a page that says that I must start the paraview server with the plugin loaded. The client would also have the plugin loaded and connect to the server. The instructions didn't say how to do that on the server side . Usually when I start the paraview server (pvserver), it does not launch the GUI but just waits for a client connection in the terminal. So I don't know how to load the plugin on the server...the command line options for pvserver didn't provide any hints. Please advise. Thanks, Sohail ___ 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] Using streamtracer/gradient filters on multi-block data set
Hi, So the reason why my required filters are greyed out after loading this CAS data set is due to a few lines in the XML file. The streamtracer is an example: DataTypeDomain name=input_type composite_data_supported=1 DataType value=vtkDataSet/ /DataTypeDomain InputArrayDomain name=input_array attribute_type=point number_of_components=3/ Now I did introduce a vtkStructuredGrid into the pipeline, which inherits from vtkDtaSet, so that shouldn't be a problem with the DataType portion above. But what concerns me the InputArrayDomain requirement. By default the data set loads in a bunch of 1-component arrays. I thought that by creating the 3-component UVW array (see below), I wouldn't face any issues. But for some reason I still do -- maybe that requirement is unfulfilled? If anyone has any ideas, please let me know. Sohail --- On Wed, 3/2/11, Sohail Shafii sohailsha...@yahoo.com wrote: From: Sohail Shafii sohailsha...@yahoo.com Subject: Re: [Paraview] Using streamtracer/gradient filters on multi-block data set To: Berk Geveci berk.gev...@kitware.com Cc: paraview@paraview.org Date: Wednesday, March 2, 2011, 1:46 PM Not sure what you mean by that last statement -- do you mean transfer an array to pointdata? I tried to do this: input = self.GetInput() output = self.GetOuput() output = vtk.vtkStructuredGrid() output.ShallowCopy(input.GetBlock(0)) // create UVW array (3-component velocity array) output.GetPointData().AddArray(UVWArray) When I print output, I do see the added array in the point data section with the correct component and tuple count. However, the bounds of the grid are suspect (like -0.1 to 0.3 for x; -0.1 to 2.6e-18 for y) and the grid extent is 0, -1, 0, -1, 0, -1. I also tried: #yiv1702808980 p, #yiv1702808980 li {white-space:pre-wrap;}output.GetPointData().SetVectors(UVWArray) Didn't work either. Something is still not right here... Sohail --- On Wed, 3/2/11, Berk Geveci berk.gev...@kitware.com wrote: From: Berk Geveci berk.gev...@kitware.com Subject: Re: [Paraview] Using streamtracer/gradient filters on multi-block data set To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org Date: Wednesday, March 2, 2011, 3:40 AM I looks like you to apply cell data to point data. On Tue, Mar 1, 2011 at 9:06 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, Let me correct myself; the block contained inside is a vtkStructuredGrid. I did a printout what is in the pipeline (self.GetInput()): vtkMultiBlockDataSet (0x3a7c1d0) Debug: Off Modified Time: 164250 Reference Count: 2 Registered Events: (none) Source: (none) Information: 0x3a7c230 Data Released: False Global Release Data: Off UpdateTime: 164282 Release Data: Off UpdateExtent: Initialized Update Number Of Pieces: 1 Update Piece: 0 Update Ghost Level: 0 MaximumNumberOfPieces: -1 ExtentTranslator: (0x3a7df30) RequestExactExtent: Off Field Data: Debug: Off Modified Time: 155460 Reference Count: 1 Registered Events: (none) Number Of Arrays: 0 Number Of Components: 0 Number Of Tuples: 0 Number Of Children: 1 Child 0: vtkUnstructuredGrid Debug: Off Modified Time: 164208 Reference Count: 3 Registered Events: (none) Source: (none) Information: 0x4920790 Data Released: False Global Release Data: Off UpdateTime: 0 Release Data: Off UpdateExtent: Not Initialized MaximumNumberOfPieces: 1 RequestExactExtent: Off Field Data: Debug: Off Modified Time: 164167 Reference Count: 1 Registered Events: (none) Number Of Arrays: 0 Number Of Components: 0 Number Of Tuples: 0 Number Of Points: 34943 Number Of Cells: 31488 Cell Data: Debug: Off Modified Time: 164206 Reference Count: 1 Registered Events: (none) Number Of Arrays: 24 Array 0 name = PRESSURE Array 1 name = X_VELOCITY Array 2 name = X_VELOCITY_M1 Array 3 name = X_VELOCITY_M2 Array 4 name = Y_VELOCITY Array 5 name = Y_VELOCITY_M1 Array 6 name = Y_VELOCITY_M2 Array 7 name = Z_VELOCITY Array 8 name = Z_VELOCITY_M1 Array 9 name = Z_VELOCITY_M2 Array 10 name = TEMPERATURE Array 11 name = ENTHALPY Array 12 name = ENERGY_M1 Array 13 name = ENERGY_M2 Array 14 name = NUT Array 15 name = NUT_M1 Array 16 name = NUT_M2 Array 17 name = DENSITY Array 18 name = DENSITY_M1 Array 19 name = DENSITY_M2 Array 20 name = MU_LAM Array 21 name = MU_TURB Array 22 name = WALL_DIST Array 23 name = BODY_FORCES Number Of Components: 26 Number Of Tuples: 31488 Copy Tuple Flags: ( 1 1 1 1 1 0 1 ) Interpolate Flags: ( 1 1 1 1 1 0 0 ) Pass Through Flags: ( 1 1 1 1 1 1 1 ) Scalars: (none) Vectors: (none) Normals: (none
Re: [Paraview] Using streamtracer/gradient filters on multi-block data set
Not sure what you mean by that last statement -- do you mean transfer an array to pointdata? I tried to do this: input = self.GetInput() output = self.GetOuput() output = vtk.vtkStructuredGrid() output.ShallowCopy(input.GetBlock(0)) // create UVW array (3-component velocity array) output.GetPointData().AddArray(UVWArray) When I print output, I do see the added array in the point data section with the correct component and tuple count. However, the bounds of the grid are suspect (like -0.1 to 0.3 for x; -0.1 to 2.6e-18 for y) and the grid extent is 0, -1, 0, -1, 0, -1. I also tried: p, li { white-space: pre-wrap; }output.GetPointData().SetVectors(UVWArray) Didn't work either. Something is still not right here... Sohail --- On Wed, 3/2/11, Berk Geveci berk.gev...@kitware.com wrote: From: Berk Geveci berk.gev...@kitware.com Subject: Re: [Paraview] Using streamtracer/gradient filters on multi-block data set To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org Date: Wednesday, March 2, 2011, 3:40 AM I looks like you to apply cell data to point data. On Tue, Mar 1, 2011 at 9:06 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, Let me correct myself; the block contained inside is a vtkStructuredGrid. I did a printout what is in the pipeline (self.GetInput()): vtkMultiBlockDataSet (0x3a7c1d0) Debug: Off Modified Time: 164250 Reference Count: 2 Registered Events: (none) Source: (none) Information: 0x3a7c230 Data Released: False Global Release Data: Off UpdateTime: 164282 Release Data: Off UpdateExtent: Initialized Update Number Of Pieces: 1 Update Piece: 0 Update Ghost Level: 0 MaximumNumberOfPieces: -1 ExtentTranslator: (0x3a7df30) RequestExactExtent: Off Field Data: Debug: Off Modified Time: 155460 Reference Count: 1 Registered Events: (none) Number Of Arrays: 0 Number Of Components: 0 Number Of Tuples: 0 Number Of Children: 1 Child 0: vtkUnstructuredGrid Debug: Off Modified Time: 164208 Reference Count: 3 Registered Events: (none) Source: (none) Information: 0x4920790 Data Released: False Global Release Data: Off UpdateTime: 0 Release Data: Off UpdateExtent: Not Initialized MaximumNumberOfPieces: 1 RequestExactExtent: Off Field Data: Debug: Off Modified Time: 164167 Reference Count: 1 Registered Events: (none) Number Of Arrays: 0 Number Of Components: 0 Number Of Tuples: 0 Number Of Points: 34943 Number Of Cells: 31488 Cell Data: Debug: Off Modified Time: 164206 Reference Count: 1 Registered Events: (none) Number Of Arrays: 24 Array 0 name = PRESSURE Array 1 name = X_VELOCITY Array 2 name = X_VELOCITY_M1 Array 3 name = X_VELOCITY_M2 Array 4 name = Y_VELOCITY Array 5 name = Y_VELOCITY_M1 Array 6 name = Y_VELOCITY_M2 Array 7 name = Z_VELOCITY Array 8 name = Z_VELOCITY_M1 Array 9 name = Z_VELOCITY_M2 Array 10 name = TEMPERATURE Array 11 name = ENTHALPY Array 12 name = ENERGY_M1 Array 13 name = ENERGY_M2 Array 14 name = NUT Array 15 name = NUT_M1 Array 16 name = NUT_M2 Array 17 name = DENSITY Array 18 name = DENSITY_M1 Array 19 name = DENSITY_M2 Array 20 name = MU_LAM Array 21 name = MU_TURB Array 22 name = WALL_DIST Array 23 name = BODY_FORCES Number Of Components: 26 Number Of Tuples: 31488 Copy Tuple Flags: ( 1 1 1 1 1 0 1 ) Interpolate Flags: ( 1 1 1 1 1 0 0 ) Pass Through Flags: ( 1 1 1 1 1 1 1 ) Scalars: (none) Vectors: (none) Normals: (none) TCoords: (none) Tensors: (none) GlobalIds: (none) PedigreeIds: (none) Point Data: Debug: Off Modified Time: 164208 Reference Count: 1 Registered Events: (none) Number Of Arrays: 0 Number Of Components: 0 Number Of Tuples: 0 Copy Tuple Flags: ( 1 1 1 1 1 0 1 ) Interpolate Flags: ( 1 1 1 1 1 0 0 ) Pass Through Flags: ( 1 1 1 1 1 1 1 ) Scalars: (none) Vectors: (none) Normals: (none) TCoords: (none) Tensors: (none) GlobalIds: (none) PedigreeIds: (none) Bounds: Xmin,Xmax: (-0.1, 0.3) Ymin,Ymax: (-0.12, 2.60209e-18) Zmin,Zmax: (-0.12, 2.60209e-18) Compute Time: 379083 Number Of Points: 34943 Point Coordinates: 0x4aa82b0 Locator: 0 Number Of Pieces: 1 Piece: -1 Ghost Level: 0Maybe it's because the point data fields, like the scalars and vectors you mentioned, are empty? I could probably fill in the gaps. Sohail --- On Tue, 3/1/11, Berk Geveci berk.gev...@kitware.com wrote: From: Berk Geveci berk.gev...@kitware.com Subject: Re: [Paraview] Using streamtracer/gradient filters on multi-block data set To: Sohail Shafii sohailsha...@yahoo.com Cc
[Paraview] Using streamtracer/gradient filters on multi-block data set
Hi, I've read in a CAS file, which Paraview read in a Multi-block dataset, which contains within it a unstructured grid. Initially I had a problem where all of the data items were one-tuple arrays. So I created a UVW (3-tuple velocity) array and added it to the pipeline this way using a Python programmable filter: input = self.GetInput() output = self.GetOutput() output.ShallowCopy(input) create array here*** output.GetBlock(0).GetCellData().AddArray(UVWArray) Note that GetBlock(0) gives me access to the grid object. For some reason the gradients and streamtracer filters were greyed out and unavailable. So I decided to make the output variable a vtkUnstructuredGrid instead and transferred the UVW array to it this way: p, li { white-space: pre-wrap; }output.GetPointData().AddArray(UVWArray) Nope...didn't work. Any thoughts how I can modify the pipeline so that the streamtracer and gradients of unstructured data sets array can become usable? I'm guessing there is something in those classes which tells Paraview that they require certain input types. I've seen references to stuff like FillInputPortInformation. Thanks, Sohail p, li { white-space: pre-wrap; ___ 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] Using streamtracer/gradient filters on multi-block data set
Hi, Let me correct myself; the block contained inside is a vtkStructuredGrid. I did a printout what is in the pipeline (self.GetInput()): p, li { white-space: pre-wrap; } vtkMultiBlockDataSet (0x3a7c1d0) Debug: Off Modified Time: 164250 Reference Count: 2 Registered Events: (none) Source: (none) Information: 0x3a7c230 Data Released: False Global Release Data: Off UpdateTime: 164282 Release Data: Off UpdateExtent: Initialized Update Number Of Pieces: 1 Update Piece: 0 Update Ghost Level: 0 MaximumNumberOfPieces: -1 ExtentTranslator: (0x3a7df30) RequestExactExtent: Off Field Data: Debug: Off Modified Time: 155460 Reference Count: 1 Registered Events: (none) Number Of Arrays: 0 Number Of Components: 0 Number Of Tuples: 0 Number Of Children: 1 Child 0: vtkUnstructuredGrid Debug: Off Modified Time: 164208 Reference Count: 3 Registered Events: (none) Source: (none) Information: 0x4920790 Data Released: False Global Release Data: Off UpdateTime: 0 Release Data: Off UpdateExtent: Not Initialized MaximumNumberOfPieces: 1 RequestExactExtent: Off Field Data: Debug: Off Modified Time: 164167 Reference Count: 1 Registered Events: (none) Number Of Arrays: 0 Number Of Components: 0 Number Of Tuples: 0 Number Of Points: 34943 Number Of Cells: 31488 Cell Data: Debug: Off Modified Time: 164206 Reference Count: 1 Registered Events: (none) Number Of Arrays: 24 Array 0 name = PRESSURE Array 1 name = X_VELOCITY Array 2 name = X_VELOCITY_M1 Array 3 name = X_VELOCITY_M2 Array 4 name = Y_VELOCITY Array 5 name = Y_VELOCITY_M1 Array 6 name = Y_VELOCITY_M2 Array 7 name = Z_VELOCITY Array 8 name = Z_VELOCITY_M1 Array 9 name = Z_VELOCITY_M2 Array 10 name = TEMPERATURE Array 11 name = ENTHALPY Array 12 name = ENERGY_M1 Array 13 name = ENERGY_M2 Array 14 name = NUT Array 15 name = NUT_M1 Array 16 name = NUT_M2 Array 17 name = DENSITY Array 18 name = DENSITY_M1 Array 19 name = DENSITY_M2 Array 20 name = MU_LAM Array 21 name = MU_TURB Array 22 name = WALL_DIST Array 23 name = BODY_FORCES Number Of Components: 26 Number Of Tuples: 31488 Copy Tuple Flags: ( 1 1 1 1 1 0 1 ) Interpolate Flags: ( 1 1 1 1 1 0 0 ) Pass Through Flags: ( 1 1 1 1 1 1 1 ) Scalars: (none) Vectors: (none) Normals: (none) TCoords: (none) Tensors: (none) GlobalIds: (none) PedigreeIds: (none) Point Data: Debug: Off Modified Time: 164208 Reference Count: 1 Registered Events: (none) Number Of Arrays: 0 Number Of Components: 0 Number Of Tuples: 0 Copy Tuple Flags: ( 1 1 1 1 1 0 1 ) Interpolate Flags: ( 1 1 1 1 1 0 0 ) Pass Through Flags: ( 1 1 1 1 1 1 1 ) Scalars: (none) Vectors: (none) Normals: (none) TCoords: (none) Tensors: (none) GlobalIds: (none) PedigreeIds: (none) Bounds: Xmin,Xmax: (-0.1, 0.3) Ymin,Ymax: (-0.12, 2.60209e-18) Zmin,Zmax: (-0.12, 2.60209e-18) Compute Time: 379083 Number Of Points: 34943 Point Coordinates: 0x4aa82b0 Locator: 0 Number Of Pieces: 1 Piece: -1 Ghost Level: 0Maybe it's because the point data fields, like the scalars and vectors you mentioned, are empty? I could probably fill in the gaps. Sohail --- On Tue, 3/1/11, Berk Geveci berk.gev...@kitware.com wrote: From: Berk Geveci berk.gev...@kitware.com Subject: Re: [Paraview] Using streamtracer/gradient filters on multi-block data set To: Sohail Shafii sohailsha...@yahoo.com Cc: paraview@paraview.org Date: Tuesday, March 1, 2011, 5:07 PM By any chance, does this multi-block dataset contain a block that does not have arrays (scalars and/or vectors)? That's probably why those filters are disabled. You can always use Extract Blocks to extract the block you are interested in and then apply the filter. As for the Programmable Filter, if you set the output type to vtkUnstructuredGrid, you should be able to do: self.GetOutput().ShallowCopy(self.GetInput().GetBlock(0)) assuming that block 0 is a uns. grid. -berk On Tue, Mar 1, 2011 at 5:29 PM, Sohail Shafii sohailsha...@yahoo.com wrote: Hi, I've read in a CAS file, which Paraview read in a Multi-block dataset, which contains within it a unstructured grid. Initially I had a problem where all of the data items were one-tuple arrays. So I created a UVW (3-tuple velocity) array and added it to the pipeline this way using a Python programmable filter: input = self.GetInput() output = self.GetOutput() output.ShallowCopy(input) create array here*** output.GetBlock(0).GetCellData().AddArray(UVWArray) Note that GetBlock(0) gives me access