Re: [Paraview] exodusII reader & writer in parallel

2017-01-13 Thread David Thompson
Hi Andy,

>> ... a basic assumption of the Exodus file format: each rank is assigned a 
>> subset of the cells and nodes in the simulation domain, since otherwise 
>> there is no state to store. ...
>> 
> My worry here is that with the processes that have data possibly changing 
> between time steps I'm not sure how well it would be read back in since there 
> probably isn't a regular way to deal with that. ...

Exodus assumes the node IDs and cell-connectivity are time-constant within a 
single file. Each timestep the number of nodes and/or cells changes, a new 
Exodus file should be created. There's a convention for how filenames should 
change at adaptation steps as opposed to process rank.

>>> > At this point I'm thinking it's better to modify the reader than the 
>>> > writer...
>> 
>> I would be very wary of this.
> 
> Yeah, this scares me as well with all of the blocks that ExodusII wants to 
> have in there but my thought was to basically skip reading in files that had 
> the num_nodes dimension not specified. Maybe there just isn't a good solution 
> here based on the format requirements. I'll test out modifying the reader 
> though and see what happens...

Be aware that the vtkPExodusIIReader behaves differently on rank 0; during 
UpdateInformation, I believe it is the only rank to do some tasks in order to 
avoid slamming parallel filesystems with tons of identical requests. I know one 
of the tasks it does it discover the list of files to split among the ranks. 
I'm less sure about transmitting lists of array names to all processes, but it 
may do that from its own list of arrays, too. There was a problem long ago 
(that got fixed) where element blocks with no cells were missing arrays -- I 
don't think we've ever had processes with no points/nodes, though.

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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] exodusII reader & writer in parallel

2017-01-13 Thread David Thompson
Hi all,

> It's looking like there could be a decent amount of hacking in the ExodusII 
> code to get the writer to work properly for this (as well as the reader to 
> read it back in without issues). The main issue appears to be that NetCDF 
> doesn't allow creating a dimension of 0 length/size. Thus the stuff in 
> VTK/ThirdParty/exodusII/vtkexodusII would need to be changed. 
> 
> Both the reader and writer make the assumption that if the num_nodes 
> dimension exists that there are points to be read in. 

Yes, that is a basic assumption of the Exodus file format: each rank is 
assigned a subset of the cells and nodes in the simulation domain, since 
otherwise there is no state to store. If rank 0 of a simulation has no nodes 
(and thus cannot have cells either), then it should probably not participate in 
writing the Exodus dataset. You could create an MPI subcontroller without the 
rank 0 process and have the writer use that.

> At this point I'm thinking it's better to modify the reader than the writer...

I would be very wary of this.

David

> 
> On Fri, Jan 13, 2017 at 2:59 PM, Andy Bauer  wrote:
> Hi Ken,
> 
> Thanks for the input. 
> 
> There is an explicit check in the writer to see if there are any points 
> before writing out point data. I took that check out and am hitting a NetCDF 
> error that the "num_nodes" dimension isn't specified. It will probably take a 
> bit of investigating to fix this but at least you've helped me go down the 
> correct path.
> 
> Thanks,
> Andy
> 
> On Fri, Jan 13, 2017 at 2:44 PM, Moreland, Kenneth  wrote:
> Andy,
> 
>  
> 
> That sounds like a bug in our Exodus writer to me. I’m not positive, but I’m 
> pretty sure that you can specify in Exodus point and cell arrays if no grid 
> points or cells exist. The writer is probably making a shortcut and skipping 
> that if there is no actual data.
> 
>  
> 
> -Ken
> 
>  
> 
> From: ParaView [mailto:paraview-boun...@paraview.org] On Behalf Of Andy Bauer
> Sent: Friday, January 13, 2017 12:00 PM
> To: paraview@paraview.org
> Subject: [EXTERNAL] [Paraview] exodusII reader & writer in parallel
> 
>  
> 
> Hi,
> 
> I'm trying out the ExodusII writer in parallel and had some questions about 
> how it should work in general when there isn't any data on process 0. 
> Currently what happens in ParaView is that a file for each process is written 
> out but the file corresponding to process 0 doesn't have any grid information 
> or field data (an example is attached). When I read this back into ParaView 
> using the built-in server (i.e. in serial) the ExodusII reader gets the 
> proper points and cells but no field data is read in.
> 
> So my question is this, does the ExodusII file format allow specifying point 
> and cell arrays in a file if no grid points or cells exist?
> 
> I'm tempted to fix this issue by modifying the reader to properly pass the 
> field data information to process 0 from another one that has data but wanted 
> to get the community's thoughts on this before I go through the 
> implementation. The alternative would be to modify the writer to include 
> point and cell data information and that's probably a better solution, 
> assuming that the ExodusII format allows for that.
> 
> Thanks,
> 
> Andy
> 
> ps. Attached is a sample set of ExodusII files that doesn't have any data on 
> parallel.ex2.4.0 in case someone wants to play around with it.
> 
> 
> 
> ___
> 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
> 
> Search the list archives at: http://markmail.org/search/?q=ParaView
> 
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] ParaView 5.2 Release Candidate 2 available for download

2016-10-21 Thread David Thompson
Hi Fabian,

Your attempt to set CMAKE_OSX_SDK=10.12 was probably on the superbuild project, 
and it is not passed to the Qt subproject as QT_MAC_SDK. Try passing the proper 
value to CMAKE_OSX_SDK (i.e., the full path which is usually 
/Applications/Xcode.app/Contents/Developer/Platforms/MacOSX.platform/Developer/SDKs/MacOSX10.12.sdk)
 and see if the superbuild properly determines the value of QMAKE_MAC_SDK from 
that.

David

> But your guess was right:
> 
> xcodebuild -showsdks
> iOS SDKs:
>   iOS 10.0-sdk iphoneos10.0
> 
> iOS Simulator SDKs:
>   Simulator - iOS 10.0-sdk iphonesimulator10.0
> 
> macOS SDKs:
>   macOS 10.12 -sdk macosx10.12
> 
> tvOS SDKs:
>   tvOS 10.0   -sdk appletvos10.0
> 
> tvOS Simulator SDKs:
>   Simulator - tvOS 10.0   -sdk appletvsimulator10.0
> 
> watchOS SDKs:
>   watchOS 3.0 -sdk watchos3.0
> 
> watchOS Simulator SDKs:
>   Simulator - watchOS 3.0 -sdk watchsimulator3.0
> 
> 
> 
>> Am 21.10.2016 um 16:11 schrieb Fabian Wein :
>> 
>> Did not work:
>> 
>> -DCMAKE_OSX_SDK=macosx10.12
>> 
>> in cmake ../paraview-superbuild -DENABLE_cxx11=ON -DENABLE_hdf5=ON 
>> -DENABLE_qt5=ON -Dsuperbuild_download_location=~/tmp/pvsb_downloads 
>> -DCMAKE_OSX_ARCHITECTURES=x86_64 -DCMAKE_OSX_DEPLOYMENT_TARGET=10.12 
>> -DCMAKE_OSX_SYSROOT=/Applications/Xcode.app/Contents/Developer/Platforms/MacOSX.platform/Developer/SDKs/MacOSX10.12.sdk
>>  -DCMAKE_OSX_SDK=macosx10.12
>> 
>> results in 
>> 
>> erforming build step for 'qt5'
>> Info: creating cache file 
>> /Users/fwein/tmp/pvsb/superbuild/qt5/build/qtbase/.qmake.cache
>> /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/ranlib:
>>  file: ../../../lib/libQt5Bootstrap.a(qvector.o) has no symbols
>> /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/bin/ranlib:
>>  file: ../../../lib/libQt5Bootstrap.a(qvector.o) has no symbols
>> clang: warning: libstdc++ is deprecated; move to libc++ with a minimum 
>> deployment target of OS X 10.9
>> clang: warning: libstdc++ is deprecated; move to libc++ with a minimum 
>> deployment target of OS X 10.9
>> Makefile:13164: warning: overriding commands for target 
>> `.moc/qeventdispatcher_cf.moc'
>> Makefile:12994: warning: ignoring old commands for target 
>> `.moc/qeventdispatcher_cf.moc'
>> clang: warning: libstdc++ is deprecated; move to libc++ with a minimum 
>> deployment target of OS X 10.9
>> Undefined symbols for architecture x86_64:
>> "std::__1::__vector_base_common::__throw_length_error() const", 
>> referenced from:
>> void std::__1::vector> std::__1::allocator 
>> >::__push_back_slow_path(QMetaMethodBuilderPrivate&&)
>>  in qmetaobjectbuilder.o
>> void std::__1::vector> std::__1::allocator 
>> >::__push_back_slow_path(QMetaPropertyBuilderPrivate&&)
>>  in qmetaobjectbuilder.o
>> void std::__1::vector> std::__1::allocator 
>> >::__push_back_slow_path(QMetaEnumBuilderPrivate&&) 
>> in qmetaobjectbuilder.o
>> "void std::__1::__sort&, int*>(int*, int*, 
>> std::__1::__less&)", referenced from:
>> QString::multiArg(int, QString const**) const in qstring.o
>> 
>> QSortFilterProxyModelPrivate::_q_sourceHeaderDataChanged(Qt::Orientation, 
>> int, int) in qsortfilterproxymodel.o
>> QSortFilterProxyModel::removeRows(int, int, QModelIndex const&) in 
>> qsortfilterproxymodel.o
>> "std::bad_alloc::bad_alloc()", referenced from:
>> qBadAlloc() in qglobal.o
>> QVector >::reallocData(int, int, 
>> QFlags) in qglobal.o
>> ld: symbol(s) not found for architecture x86_64
>> clang: error: linker command failed with exit code 1 (use -v to see 
>> invocation)
>> make[6]: *** [../../lib/QtCore.framework/QtCore] Error 1
>> make[5]: *** [sub-corelib-make_first] Error 2
>> make[4]: *** [sub-src-make_first] Error 2
>> make[3]: *** [module-qtbase-make_first] Error 2
>> CMake Error at /Users/fwein/tmp/pvsb/superbuild/sb-qt5-build.cmake:43 
>> (message):
>> Failed with exit code 2
>> 
>> I have the impressions that my other issues and comments on the gitlab issue 
>> tracker are removed, so I stop posting reports of attempts there.
>> 
>> 
>> 
>>> CMAKE_OSX_SDK needs to be something like `macosx10.9`.  I forget the 
>>> command to make XCode list the valid sdk versions available.  For your case 
>>> I would guess `macosx10.12` is the correct value.  Qt looks up the path to 
>>> the SDK from that string.
>>> 
>>> HTH,
>>> Shawn
>>> 
>>> On Fri, Oct 21, 2016 at 7:07 AM, Fabian Wein  wrote:
>>> I added a lot of the stuff to the issue tracker.
>>> 
>>> I’m currently lost with building qt5 on macOS 10.12. It’s a qmake issue, I 
>>> tracked it down to
>>> 
>>> Creating qmake... Running configuration tests... Failed to process makespec 
>>> for platform ‚macx-clang' Info: creating stash file 
>>> /Users/fwein/tmp/pvsb/superbuild/qt5/build/qtbase/.qmake.stash Project 
>>> ERROR: QMAKE_MAC_SD

Re: [Paraview] wall shear stress

2016-09-18 Thread David Thompson
Hi Fande,

You should be able to use the gradient filter to compute the velocity gradient 
and examine the component of the gradient normal to the wall. That will be 
proportional to the shear stress at the wall. Is the wall aligned with a 
coordinate axis? If not, you will need to extract the wall surface and generate 
normal vectors, then dot the gradient with the normal vectors.

David

> On Sep 18, 2016, at 16:01, Fande Kong  wrote:
> 
> Hi Developers,
> 
> I am solving an incompressible Navier-Stokes equations, and I have solutions 
> for the velocity (u)  and the pressure (p). Is it possible to show the wall 
> shear stress in paraview using the velocity?  
> 
> 
> Fande Kong,
> ___
> 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
> 
> Search the list archives at: http://markmail.org/search/?q=ParaView
> 
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Finding cell coordinates in Python

2016-06-18 Thread David Thompson
> We have been successful at implementing python scripts and routines to 
> analyze much of our data, but we would like for the ability to obtain the 
> [X,Y,Z] coordinates for each cell. Is there currently a way to do this? This 
> is CTH AMR data.

Perhaps they have been writing client-side scripts to get at the metadata but 
now want access to the full dataset?

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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] [EXTERNAL] python calculator confusion

2016-06-16 Thread David Thompson
Hi all,

> ...
> 1. Make sure you set "Array Association" on "Python Calculator" to "Cell 
> Data".
> 2. The volume() code internally uses vtkCellQuality filter. Now, when
> I debugged into it, I see that is doesn;t support voxel cell types --
> the cells produced by vtkImageData.

Yes, that is an artifact of using the verdict library to compute volume. The 
cell quality filter is being used to compute a "quality" metric that is rarely 
used as a primary indicator of mesh quality, so this is not really the intended 
use of the filter. The cell quality of axis aligned boxes is not usually of 
interest to anyone, so when writing the quality filter I did not consider that 
cell type... only those types directly supported by verdict. While volume can 
be obtained for other types by tetrahedralizing and summing, that does not work 
for other metrics like Frobenius norm. Many element shapes do not have 
well-characterized quality metrics.

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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Using Paraview for modeling orbital motion of many objects/database connectivity?

2016-02-02 Thread David Thompson
Hi Cynthia,

It's not something ParaView supports without any modifications, but it is not 
impossible. VTK, the toolkit ParaView is built upon, provides ways to populate 
tables from SQL database queries and also ways to turn those tables into 
polygonal datasets that can be rendered. It would take a small amount of 
programming to connect those two together into a ParaView pipeline object.

Depending on what's in the database, it would probably also require some 
filters to turn the satellite ephemera into a time-series of point coordinates 
for animations. I'm less certain about how much work that would be; I'm not 
familiar enough with orbital mechanics to know how much work that is.

Hope that helps,
David

> On Feb 2, 2016, at 12:26 PM, Hart, Cynthia Joyce - (cyndyb) 
>  wrote:
> 
> Hi folks -
> 
> Forgive me if this question is entirely out of scope of Paraview’s 
> capabilities, but I’m new to the software, and still learning what it can and 
> cannot do.  I have a rather large scale project that I have been asked to 
> look into the feasibility of.  It involves a simulation of the display of a 
> multitude of satellite ‘and other’ objects orbiting the Earth, with the 
> correct orbital motions, and in an animated display, no less.  Due to the 
> number of objects, I have the idea in mind that a database of some sort is 
> the best way to store the information of the objects, and have whatever is 
> driving the simulation call on that database.  I have no idea if Paraview 
> could be - or has been - used with a database.  Is this at all do-able with 
> Paraview?  And if not, what do you think I ought to be considering, for 
> something of this scope?  
> 
> Thanks much,
> 
> Cynthia
> 
> 
> ___
> 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
> 
> Search the list archives at: http://markmail.org/search/?q=ParaView
> 
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] 3D widget in branded app panel?

2015-11-09 Thread David Thompson
Hi Utkarsh,

> I believe you're missing a pq3DWIdget::select() call after the showWigdte().

That did it.

Thanks!
David

> 
> On Mon, Nov 9, 2015 at 10:15 AM, David Thompson
>  wrote:
>> Hi all,
>> 
>> I need to add some 3D widgets to a custom ParaView application. I've 
>> succeeded in creating a pqImplicitPlaneWidget and pqSphereWidget and adding 
>> their Qt widgets to the UI. However, the 3D widget never appears in the 
>> render view; the plane/sphere show up, but not the widgets. It's unclear 
>> from the documentation/code what proxy I should pass to the widget 
>> constructors for the "controlled" proxy. This is what I'm currently doing:
>> 
>>  pqServer* server = paraViewApp->getActiveServer();
>>  pqObjectBuilder* builder = paraViewApp->getObjectBuilder();
>>  QPointer source =
>>builder->createSource("sources", "PlaneSource", server);
>>  vtkSMProxy* ctrlProxy = source->getProxy();
>>  vtkSMProxy* refProxy = ...; // proxy of the pqPipelineSource whose bounds 
>> should be used by the widget
>>  pq3DWidget* pvwidget = new pqImplicitPlaneWidget(refProxy, ctrlProxy);
>>  pqActiveObjects& actives(pqActiveObjects::instance());
>>  pvwidget->setView(actives.activeView());
>>  pqDataRepresentation* rep =
>>builder->createDataRepresentation(
>>  source->getOutputPort(0), actives.activeView())
>>  panelLayout->addWidget(pvwidget);
>>  pvwidget->showWidget();
>> 
>> I've even tried
>> 
>>  double bds[6] = { -85, 85, -202, -34, -216, -47 };
>>  pvwidget->resetBounds(bds);
>> 
>> to fix the bounds to something meaningful, but this has no effect.
>> 
>> Any hints would be useful.
>> 
>>Thanks,
>>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
>> 
>> Search the list archives at: http://markmail.org/search/?q=ParaView
>> 
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


[Paraview] 3D widget in branded app panel?

2015-11-09 Thread David Thompson
Hi all,

I need to add some 3D widgets to a custom ParaView application. I've succeeded 
in creating a pqImplicitPlaneWidget and pqSphereWidget and adding their Qt 
widgets to the UI. However, the 3D widget never appears in the render view; the 
plane/sphere show up, but not the widgets. It's unclear from the 
documentation/code what proxy I should pass to the widget constructors for the 
"controlled" proxy. This is what I'm currently doing:

  pqServer* server = paraViewApp->getActiveServer();
  pqObjectBuilder* builder = paraViewApp->getObjectBuilder();
  QPointer source =
builder->createSource("sources", "PlaneSource", server);
  vtkSMProxy* ctrlProxy = source->getProxy();
  vtkSMProxy* refProxy = ...; // proxy of the pqPipelineSource whose bounds 
should be used by the widget
  pq3DWidget* pvwidget = new pqImplicitPlaneWidget(refProxy, ctrlProxy);
  pqActiveObjects& actives(pqActiveObjects::instance());
  pvwidget->setView(actives.activeView());
  pqDataRepresentation* rep =
builder->createDataRepresentation(
  source->getOutputPort(0), actives.activeView())
  panelLayout->addWidget(pvwidget);
  pvwidget->showWidget();

I've even tried

  double bds[6] = { -85, 85, -202, -34, -216, -47 };
  pvwidget->resetBounds(bds);

to fix the bounds to something meaningful, but this has no effect.

Any hints would be useful.

Thanks,
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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] [EXT] Trusses to Tubes

2015-10-16 Thread David Thompson
Hi Dennis,

> ... running extract surface and then tube is working very well for me.  I 
> guess I never thought about a "line" element having a surface but apparently 
> that's a line!

Yes, the "Extract Surface" filter's name is a little misleading. What it really 
does is convert a dataset from an unstructured grid (which may contain 
volumetric elements) to a polydata (which may only contain 0-, 1-, or 2-d 
primitives). In the process, any volume elements are represented by the 
portions of their bounding surface not shared with other volume elements.

The Tube filter only runs on polydata. So it is not exactly extracting the 
surface of a line so much as converting representations.

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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Trusses to Tubes

2015-10-15 Thread David Thompson
Hi Dennis,

> I have a Exodus multi-block dataset with some blocks that are truss elements 
> (unstructured grid).Each truss represents a cord and I would like to run 
> the Tube filter on them to display their actual diameter so that I can check 
> how close they come to touching each other.  Of course, the Tube only runs on 
> lines, not trusses. ...

I can imagine 2 issues at play:

1. One issue may be that the Tube filter expects polydata, not unstructured 
grids (which is what the Exodus reader generates by default). Running "Extract 
Surface" on the reader output should let you run the Tube filter on the result.

2. Another issue is that the Exodus reader does not currently read element 
block attributes, which is normally where truss cross-sectional area is stored. 
Thus you would have to know the proper radius to pass to the Tube filter.

It would not be a huge amount of work to add Exodus attributes as field data 
nor to have the reader run a tube filter on truss blocks, but it is more than I 
have time to do right now. If you want to write a patch, I can point you to 
where changes would need to be made.

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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] [EXT] Re: Exodus Block Names

2015-08-31 Thread David Thompson
Hi all,

The Exodus reader creates 8 top-level blocks to identify different types of 
data present in each file (corresponding to the names you see for 0-7). Each of 
these blocks contains child blocks for the actual data loaded from the file 
(which may be a subset of the data present in the file). It is these sub-blocks 
whose metadata will hold the names from the Exodus file.

David

> On Aug 31, 2015, at 10:40 AM, Dennis Conklin  
> wrote:
> 
> Shawn,
>  
> Ok, I’m dense.
>  
> When I run your code inside the Programmable Filter, I see some confusing 
> things:
> mbi.GetNumberOfBlocks returns 8, which is NOT the number of blocks in my 
> model but IS the number of MetaData blocks.
>  
> The mbi.GetMetaData[i].Get(vtk.vtkCompositeDataSet.NAME()) then is
> Value of i Name
> 0Element Blocks
> 1Face Blocks
> 2Edge Blocks
> 3Element Sets
> 4Side Sets
> 5Face Sets
> 6Edge Sets
> 7Node Sets
> 8
>  
> I am poking around mbi.GetMetaData[0] (Element Blocks) but I still haven’t 
> found any Block Names there.I feel like I am completely missing something 
> here, but I have no idea what it is.
>  
> Dennis
>  
>  
> From: Shawn Waldon [mailto:shawn.wal...@kitware.com] 
> Sent: Monday, August 31, 2015 9:46 AM
> To: Dennis Conklin
> Cc: Paraview (paraview@paraview.org)
> Subject: Re: [EXT] Re: [Paraview] Exodus Block Names
>  
> Hi Dennis,
> 
> The metadata is on the reader's output, which is a vtkMultiBlockDataSet.  
> reader.GetOutput() should get you the dataset in your python script.  Inside 
> the programmable filter you will need to get the input dataset 
> (self.GetInput() should get you the input dataset and self.GetOutput() should 
> get you the output dataset).  So something like the following for your 
> programmable filter:
> 
> mbi = self.GetInput()
> mbo = self.GetOutput()
> 
> mbo.ShallowCopy(mbi)
> 
> for i in range(mbo.GetNumberOfBlocks()):
>   metadata = mbo.GetMetaData(i)
>   name = metadata.Get(vtk.vtkCompositeDataSet.NAME())
>   # do something with the name
>  
>  
> HTH,
> Shawn
>  
> On Mon, Aug 31, 2015 at 8:31 AM, Dennis Conklin  
> wrote:
> Shawn,
>  
> Thanks for that tip but I can’t seem to access this structure.   Really, I 
> need it within the Programmable Filter, but even when I run a Python script 
> and try to find it directly in an Exodus reader, I can’t seem to locate this 
> metadata.
>  
> If I print dir(ExodusReader), there doesn’t seem to be anything about 
> metadata.
> Dennis
>  
> From: Shawn Waldon [mailto:shawn.wal...@kitware.com] 
> Sent: Friday, August 28, 2015 4:44 PM
> To: Dennis Conklin
> Cc: Paraview (paraview@paraview.org)
> Subject: [EXT] Re: [Paraview] Exodus Block Names
>  
> Hi Dennis,
> 
> The block name is in the block metadata, which is not where I looked the 
> first time either.  Here is a code snippet that shows how to access it.
> 
> mb = vtk.vtkMultiBlockDataSet()
> ...
> for i in range(mb.GetNumberOfBlocks):
>   metadata = mb.GetMetaData(i)
>   name = metadata.Get(vtk.vtkCompositeDataSet.NAME())
> 
> HTH,
> Shawn
>  
> On Fri, Aug 28, 2015 at 4:13 PM, Dennis Conklin  
> wrote:
> All,
>  
> Once again, I need help.  We are starting to assign names to our Exodus 
> blocks (outside of Paraview).   This is very useful for ease and clarity of 
> post processing and allows convenient reference to actual components.   The 
> good news is that Paraview seems perfectly happy to read in the Block Names 
> and use them in the Properties Panel and the Multi-block Inspector and the 
> Find Data screen if they exist (no action or programming required for this).  
> So, for instance, instead of Paraview generating the non-useful name of 
> “Unnamed block ID: 13 Type: hex” when the blockname is empty, it will 
> automagically use the more useful Blockname of “Tread”  if that is in the 
> Exodus file.  
>  
> However, there has been an unintended consequence of this change.   If I 
> select some elements and examine them in Spreadsheet View, I see Block Number 
> 14( the +1 offset from Block_ID and Block Number is NOT the problem!).   And 
> now I can no longer associate the element (Block Number 14) with it’s block 
> in the Properties Panel or the Multi-block Inspector or the Find Data panel – 
> where the same block is called “Tread”.
>  
> So, to remedy this situation I would like to add an element variable of 
> BlockName, which would contain “Tread” for all the elements in BlockID 13, 
> etc., etc. for all the blocks.
>  
> My only problem is I can’t seem to find the list of Block Names.   I’m pretty 
> weak at vtk, so it must be available but I surrender – I can’t seem to find 
> it.   If someone can help me find the list of Block Names from inside the 
> Programmable Filter, I’ll be happy t

Re: [Paraview] coords not available

2015-05-01 Thread David Thompson
> Remember that image data and rectilinear grids don't explicitly store point 
> locations. Can the cut filter be used instead? He can do the crinkle cut if 
> he wants the entire cell.

Now that we have data array subclasses that provide only iterator-based access, 
it might be a good idea to revisit this decision. It would be easy to create 
new arrays that present just one component of explicit coordinates or to 
present implicit coordinates without allocating memory. It certainly would make 
a lot of things easier.

David

> 
> On Fri, May 1, 2015 at 12:51 PM, Scott, W Alan  wrote:
> Why aren’t coords available to filters in ParaView?  I have a user that wants 
> to threshold on one of the components of coord, and every time he has to use 
> the calculator to get at this variable.
>  
> Thanks,
>  
> Alan
>  
>  
> 
> ___
> 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
> 
> Search the list archives at: http://markmail.org/search/?q=ParaView
> 
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/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
> 
> Search the list archives at: http://markmail.org/search/?q=ParaView
> 
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] how not to scale the thickness of the Arrow Glyph?

2015-03-25 Thread David Thompson
Hi Kazuyoshi,

I do not think the glyph filter will do what you want it to. You might use a 
Python programmable filter[1] to turn each point into a line segment (i.e., a 
line from (x,y,z) to (x+dx, y+dy, z+dz)). Then run "Cell Centers" on the 
resulting line segment data and add oriented, unscaled cone glyphs at each cell 
center to indicate the direction along the line segments. That would draw lines 
along the path plus cones pointing along the direction of the path.

David

[1]: http://www.paraview.org/Wiki/Python_Programmable_Filter


> On Mar 25, 2015, at 8:47 PM, Kazuyoshi Furutaka 
>  wrote:
> 
> Dear ParaView experts,
> 
> I'm visualizing paths of neutrons in some media.
> Now I prepared the data as point data with displacement
> vector to the next point attached to each point, and
> visualize them using the Glyph filter with the following
> setting:
>  Glyph Type: Arrow
>  Orientation: Orient
>  Scale Mode "vector"
>  Scale Factor = 1.0
>  Glyph Mode: All Points
> 
> The problem is that some arrows are too thick...
> 
> So here are my questions.
> 
> (1)How can I make the thickness of the arrows not
>   to scale (only the lengths do)?
> 
> (2)It may be better to use "2D Glyph" instead of
>   "Arrow", but the 2D Glyph is centered at a
>   data point, instead of starting from the point...
>   How can I make the 2D Glyphs start from the data
>   points?
> 
> (3)Do you have any recommendation on how to visualize
>   such a particle paths?
> 
> Thanks in advance,
> 
> Kazuyoshi
> --
> Kazuyoshi Furutaka
> furutaka _dot_ kazuyoshi _at_ jaea _dot_ go _dot_ jp
> ___
> 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
> 
> Search the list archives at: http://markmail.org/search/?q=ParaView
> 
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] colormap nan opacity

2015-03-20 Thread David Thompson
Hi Ryan,

> I use paraview to render ocean model data. My datasets are usually 
> rectilinear grids ...
> ... The colormap editor allows you to choose your NaN color, but this dialog 
> doesn't include an alpha / opacity option. ...
> Instead, what I have to do is use a clip filter to mask the NaNs, which ends 
> up turning my rectilinear grid data into unstructure grid data, ...

I haven't worked with rectilinear grids in a while... but at least one 
structured-data PV representation accepts a blanking array. Would it be 
possible to use that to cull cells instead of the NaN color?

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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Writer plugin?

2015-01-29 Thread David Thompson
Hi David

>   
> 
>   
>   
> 
> 
> ?

Nope.

Thanks,
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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Writer plugin?

2015-01-29 Thread David Thompson
Hi John,

> I used to make my writers the same as filters ...
> Adding the writer as a filter is an easy option!

Thanks! Changing the ProxyGroup to "filters" and the WriterProxy element to a 
SourceProxy seems to have worked. I still want to understand why I can't expose 
it as a writer, but at least I can move on to the actual work.

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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Writer plugin?

2015-01-28 Thread David Thompson
Hi Burlen,

> The WriterProxy may have a special "file_name_method" attribute that needs to 
> be set. I'm not sure if that's required, optional, or long since obsolete 
> since I haven't used our plugin in quite a while.

I did see the "file_name_method" (and "supports_parallel" and "parallel_only") 
attribute, but not all of the writers have it. When I add it, nothing happens.

Thanks,
David

> 
> 1314 class="vtkSQBOVWriter"
> 15 file_name_method="SetFileName"
> 16 supports_parallel="1">
> 
> On 1/28/2015 8:54 AM, David Thompson wrote:
>> Hi Andy,
>> 
>>> I think you may be missing the FileName property. I thought that PV looked 
>>> for that attribute to know that it should be classified as a writer. Also, 
>>> you'll want to hide that property.
>> Adding the FileName property and hints did not work.
>> 
>>  Thanks,
>>  David
>> 
>>> So:
>>>   >> command="SetFileName"
>>> number_of_elements="1"
>>>     label="File Name">
>>> 
>>>   The file name.
>>> 
>>>   
>>> 
>>> 
>>> and:
>>>   
>>> 
>>> 
>>> 
>>>   
>>> 
>>> 
>>> 
>>> 
>>> On Wed, Jan 28, 2015 at 11:34 AM, David Thompson 
>>>  wrote:
>>> Hi all,
>>> 
>>> I have a ParaView plugin that contains a new writer which I would like to 
>>> expose. However, while the plugin does properly expose a new view and 
>>> representation, the writer does not appear to get added (when I try to save 
>>> a vtkTable such as a statistics filter output, it doesn't appear as an 
>>> option).
>>> 
>>> 1. Is there any documentation for adding writers in plugins? This page has 
>>> a comment or two but doesn't appear complete:
>>> 
>>> http://www.paraview.org/Wiki/ParaView/Plugin_HowTo
>>> 
>>> 2. If you're willing to look over the XML, I've attached it below. See any 
>>> mistakes?
>>> 
>>> Thanks,
>>> David
>>> 
>>> 
>>>   
>>> 
>>>   
>>> >>   class="vtkPJSONTableWriter"
>>>   name="PJSONTableWriter">
>>>   
>>> Write tabular data to a JSON file.
>>>   
>>>   >> command="SetInputConnection"
>>> name="Input">
>>> >>   composite_data_supported="0"
>>>   name="input_type">
>>>   
>>> 
>>> The input filter/source whose output dataset is to
>>> written to the file.
>>>   
>>>   
>>> 
>>>   
>>> 
>>> 
>>>   
>>> 
>>> 
>>> ___
>>> 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
>>> 
>>> Search the list archives at: http://markmail.org/search/?q=ParaView
>>> 
>>> Follow this link to subscribe/unsubscribe:
>>> http://public.kitware.com/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
>> 
>> Search the list archives at: http://markmail.org/search/?q=ParaView
>> 
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Writer plugin?

2015-01-28 Thread David Thompson
Hi Andy,

> I think you may be missing the FileName property. I thought that PV looked 
> for that attribute to know that it should be classified as a writer. Also, 
> you'll want to hide that property.

Adding the FileName property and hints did not work.

Thanks,
David

> So:
>command="SetFileName"
> number_of_elements="1"
> label="File Name">
> 
>   The file name.
> 
>   
> 
> 
> and:
>   
>     
>     
>  />
>   
> 
> 
> 
> 
> On Wed, Jan 28, 2015 at 11:34 AM, David Thompson  
> wrote:
> Hi all,
> 
> I have a ParaView plugin that contains a new writer which I would like to 
> expose. However, while the plugin does properly expose a new view and 
> representation, the writer does not appear to get added (when I try to save a 
> vtkTable such as a statistics filter output, it doesn't appear as an option).
> 
> 1. Is there any documentation for adding writers in plugins? This page has a 
> comment or two but doesn't appear complete:
> 
> http://www.paraview.org/Wiki/ParaView/Plugin_HowTo
> 
> 2. If you're willing to look over the XML, I've attached it below. See any 
> mistakes?
> 
> Thanks,
> David
> 
> 
>   
> 
>   
>class="vtkPJSONTableWriter"
>   name="PJSONTableWriter">
>   
> Write tabular data to a JSON file.
>   
>command="SetInputConnection"
> name="Input">
>composite_data_supported="0"
>   name="input_type">
>   
> 
> The input filter/source whose output dataset is to
> written to the file.
>   
>   
>  />
>   
> 
> 
>   
> 
> 
> ___
> 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
> 
> Search the list archives at: http://markmail.org/search/?q=ParaView
> 
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


[Paraview] Writer plugin?

2015-01-28 Thread David Thompson
Hi all,

I have a ParaView plugin that contains a new writer which I would like to 
expose. However, while the plugin does properly expose a new view and 
representation, the writer does not appear to get added (when I try to save a 
vtkTable such as a statistics filter output, it doesn't appear as an option).

1. Is there any documentation for adding writers in plugins? This page has a 
comment or two but doesn't appear complete:

http://www.paraview.org/Wiki/ParaView/Plugin_HowTo

2. If you're willing to look over the XML, I've attached it below. See any 
mistakes?

Thanks,
David


  

  

  
Write tabular data to a JSON file.
  
  

  

The input filter/source whose output dataset is to
written to the file.
  
  

  


  


___
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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] [EXTERNAL] Re: temporal threshold

2015-01-26 Thread David Thompson
> You may be right.  I am just dredging up old memories from a few years ago.  
> David Thompson should remember – he helped write the filter.  David?

I am sure that at one point it worked in parallel, with the caveat that samples 
of points or point-data will over-represent values on process boundaries 
(because ParaView does not have knowledge of how many processes hold any given 
process-boundary point).

Err, and as I see Ken just mentioned, I only worked on the statistics filters 
that Temporal Statistics uses under the hood.

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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] [EXTERNAL] Opacity crashing ParaView and Value Range field for creating multiple slices is missing

2015-01-14 Thread David Thompson
Hi Mary,

Try clicking the small gear icon in the property panel. That should display the 
advanced settings like value range.

David



> On Jan 14, 2015, at 23:30, Mary Draney Blomme  wrote:
> 
> Alan,
> 
> Thanks for your reply. I'm out of the office until Friday, so I'll have to 
> wait to try upgrading my graphics drivers on the desktop until then. The 
> laptop problem, however, continues. I'm attaching a screen snapshot of what I 
> see. Am I missing something?
> 
> Thanks again,
> 
> - Mary
> 
>> On Wed, Jan 14, 2015 at 8:57 PM, Scott, W Alan  wrote:
>> Mary,
>> 
>> Since one computer works and the other doesn’t, sounds like a graphics card 
>> issue to me.  Try updating your graphics drivers?  Also, try PV 4.3.0.  I 
>> know there was a problem with some volume rendering for older versions of 
>> OpenGL.  This should be fixed for PV 4.3.0.  Note that the next version of 
>> ParaView will require a somewhat newer version of OpenGL than PV 4.3.0, so 
>> if this is the problem, you are on borrowed  time for volume rendering and 
>> opacity rendering.
>> 
>>  
>> 
>> I just tried, and got multiple slices to work just fine.  The trick is to 
>> click on the blue plus sign, slightly right of the Value Range window.  You 
>> can then edit the value that appears.  Does that work?
>> 
>>  
>> 
>> Alan
>> 
>>  
>> 
>>  
>> 
>>  
>> 
>>  
>> 
>> From: ParaView [mailto:paraview-boun...@paraview.org] On Behalf Of Mary 
>> Draney Blomme
>> Sent: Wednesday, January 14, 2015 1:16 PM
>> To: paraview@paraview.org
>> Subject: [EXTERNAL] [Paraview] Opacity crashing ParaView and Value Range 
>> field for creating multiple slices is missing
>> 
>>  
>> 
>> Hi,
>> 
>>  
>> 
>> I'm new to ParaView and going through the first webinar, using the demo data 
>> set disk_out.ref. I'm using ParaView 4.2.0 on two different Windows 64-bit 
>> computers.
>> 
>>  
>> 
>> I'm having two problems, one on each computer.
>> 
>>  
>> 
>> 1. On the first computer, a desktop, when I go to change opacity of the main 
>> disk_out_ref.ex2 object, ParaView crashes. Is this a problem with my 
>> graphics card?
>> 
>>  
>> 
>> 2. On the second computer, a laptop, the opacity changes just fine, however, 
>> for creating multiple slices, the Value Range field is missing (under the 
>> Properties tab for Slice 1, immediately under the "Triangulate the slice" 
>> option). How do I get it to show up? I have tried quitting and restarting, 
>> and also tried ParaView 4.3.0 and had no success.
>> 
>>  
>> 
>> Thanks,
>> 
>>  
>> 
>> - Mary
>> 
> 
> 
> ___
> 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
> 
> Search the list archives at: http://markmail.org/search/?q=ParaView
> 
> Follow this link to subscribe/unsubscribe:
> http://public.kitware.com/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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] ANN: ParaView 4.3.0-RC1 now available for download

2015-01-12 Thread David Thompson
Hi Utkarsh,

> To get the year started, we have the ParaView 4.3.0, Release Candidate 1 
> binaries available for download [1]. ...


Will it be possible to build branded ParaView applications from an *installed* 
version of ParaView 4.3 instead of requiring a build directory? I haven't been 
able to find a bug in the tracker for this.

Thanks,
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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] FW: sideset ids in paraview

2014-11-12 Thread David Thompson
Hi Alan,

> We are working with a large mesh that has side sets defined by string names 
> and 
> not integer values.  Does Paraview assign an integer if value to the side set 
> and if so,
> how does it assign the id and is there a way to extract that from paraview?
> Thanks!

The Exodus reader generates a multiblock dataset with one toplevel block 
containing all the side sets. The index of a side set inside that toplevel 
block is an integer and can be obtained in a Python programmable filter. The 
advantage of this is that it's cheap, storage-wise. The disadvantage is that 
the index will vary depending on what side sets you choose to load.

The Exodus reader can also provide a cell array identifying the side set ID (an 
integer stored in the Exodus file). See

  
http://vtk.org/gitweb?p=VTK.git;a=blob;f=IO/Exodus/vtkExodusIIReader.h;h=10de0366f262732d333ba45187917228d428c21e;hb=HEAD#l112

for more information. I believe this is exposed in ParaView, but forget what 
the button is called... it used to be "Generate ObjectIds" but someone  decided 
that was unclear and was going to change it to something that *I* think is 
unclear. :-)

The advantage of the generated IDs is that they do not vary depending on what 
you choose to load. The disadvantage is that it is not stored on the block, but 
on each cell in the block, so it takes up a lot of memory.

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://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] FW: (ParaView) Lofting Shells (INC000000978412)

2014-10-28 Thread David Thompson
> My user wants to know if ParaView can "loft" shells?  Any ideas?

Not as far as I know. One issue is that different FEM codes parameterize shell 
elements differently:

+ Some assume that the surface is a true mid-surface. Others allow different 
"top" and "bottom" heights. Others "loft" only in one direction.

+ Some treat neighboring shells as constraining the offset nodes for problems 
where contact may be an issue. Others do not do contact.

The latter in particular can lead to tricky problems when angles between 
neighboring shell elements are acute. Doing something akin to the "tube" 
filter, but for quads+tris would be easy as a Python programmable filter. Doing 
something where neighboring elements were taken into account would be 
significantly more work.

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://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Labeling specific points in dataset

2014-10-16 Thread David Thompson
Hi Kit,

If you can find a way to create a selection of the points, the selected points 
can be labeled. Is that feasible?

David

On Oct 16, 2014, at 7:29, Kit Chambers  wrote:

> Hi Everyone,
> 
> I have a dataset of about 1000 or so points and I want to label say every 
> 100'th point with the point id. I can use the DataLabelRepresentation if I 
> want to label every point, but as far as I can tell this cannot label only 
> certain points it is everything or nothing. I could also create a new dataset 
> with every 100th point in it and use that with the DataLabelRepresentation 
> but then the point ids would be wrong.  
> 
> Any ideas or suggestions are welcome.
> 
> Cheers
> 
> Kit
> 
> 
> 
> 
> ___
> 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://public.kitware.com/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://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] [ParaView] Unofficial ParaView 4.1.0 SDK for Windows 32 and 64-bit (UNCLASSIFIED)

2014-08-12 Thread David Thompson
> Thanks. And per Robert, ParaView binaries are done with vs2008 due to Python
> build system. I tried to build python 2.7.8 (the very last dependency
> library) with vs2013 express but failed miserably (even with what Gerson
> Kurz how to list to build python 3.3.4 with Visual Studio 2013). Thankfully
> I was able to find the link to download vs2008 express iso.

David Sansome and JC Fillion-Robin have created a CMake-ified version of Python 
that will build with other versions of VS:

  https://github.com/davidsansome/python-cmake-buildsystem/

but currently (as of a month ago) it does not support debug builds -- which is 
something we most probably require before moving ParaView to use it.

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://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] problem on the max field value over time

2014-07-28 Thread David Thompson
Hi,

> When I run the descriptive statistics filter, I do not find the option to 
> subsample the data, so I think it is applied to the whole dataset.

The option is called "Training Fraction" and should be 1.0 if you want the 
description to include all of the data.

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://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] problem on the max field value over time

2014-07-25 Thread David Thompson
> Uh, oh. That's not right.

When you run the descriptive statistics filter, are you running on all of the 
data or a subsample (I believe a subsample is the default)?

David

> 
> I've submitted a bug report about this problem with descriptive statistics 
> (http://www.paraview.org/Bug/view.php?id=14886). Thanks for reporting the 
> problem.
> 
> -Ken
> 
> From: linzhenhua 
> Date: Wednesday, July 23, 2014 3:05 AM
> To: "Paraview@paraview.org" 
> Subject: [EXTERNAL] [Paraview] problem on the max field value over time
> 
> Hi,
> 
> I want to plot the temporal variation of the maximum field value over time. I 
> use two methods and the results are different in paraview.
> 
> The first method is to use the Descriptive statistics filter (Green line). 
> The second method is to use the Find data tool(orange line). It seems to me 
> that the latter one is more reasonable as the max field value should decrease 
> with time. The result from the find data tool also agree with the scale data 
> range.
> 
> Does this mean that there is a bug in Descriptive Statistics filter?
> 
> 
> 
> 
> ___
> 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://public.kitware.com/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://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] [EXTERNAL] Re: Rotational velocities

2014-07-08 Thread David Thompson
> I've never done this, and maybe there is an existing established way to do 
> it, but a simple Google search did not reveal it.

The proper Google search phrase to find mechanical engineering sites is 
"instant center of rotation". Most of them deal with 2-D though, as that's 
what's covered in undergraduate mechanics. If you find an advanced mechanics 
book (or an aerospace mechanics book), you should find the 3-D case.

> ... I'm guessing that he would probably be interested in rotations around the 
> center of mass.

I think it's more likely for someone to be interested in the rotation about the 
instant center than the center of mass, but the angular velocity about the 
center of mass is also a useful thing to know. (It is also nice to know the 
moment of inertia about whatever axis you are dealing with because then you can 
estimate rotational kinetic energy as I\omega^2 and the linear kinetic energy 
as m*v^2 and compare the two.)

> An average of all the point coordinates might be a close enough approximation 
> to the center of mass.

I would be nervous about using that approximation in a general-purpose tool 
because the angular velocities that result can vary significantly for small 
changes in the position of the center of mass when velocity is large.

You could use the mesh quality filter and cell-center filters to avoid bias 
from the mesh. Or do a least squares fit to the instant center of rotation and 
the angular velocity all at once.

> Next, find the translation of the rigid body by looking at the displacement 
> at the center of mass. That displacement is the translation. Now subtract 
> that translation vector from all the displacements to get what I will call 
> the rotation displacement (a term I just made up).

If the body is deforming (or the element block contains independently-moving 
disjoint bodies), both the estimated displacement and your "rotation 
displacement" include variations due to non-rigid deflections.

> Now you can use any of these rotation displacements to find the rotation. 
> Technically any one displacement will give you the rotation (because it is a 
> rigid body), but you should pick one far away from the axis of rotation.

The problem I have with this is that -- assuming non-rigid deflections are 
significant -- the ones far from the axis of rotation are likely the ones most 
affected by non-rigid-body deflections.

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://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Rotational velocities

2014-07-08 Thread David Thompson
>> 2)  I want to calculate the rotational velocities of a component (block 
>> of elements) in paraview using the nodal translational velocities. How do I 
>> do that?

> It would be good to have some more clear mathematical definition of (2). I am 
> assuming it involves some sort of integration but it is not clear.

Assuming that the elements are (mostly) undergoing rigid-body motion, you could 
estimate the instant center of rotation and decompose the velocity into 
rotational and linear components. It should only take 3 points and their 
velocities (in 3-D, or 2 pts in 2-D) to estimate the instant center. For the 
2-D case, see

  http://emweb.unl.edu/NEGAHBAN/EM373/note16/note.htm

In 3-D, you get an axis, not a point. Having more than 3 (or 2 in 2-D) points 
generates a distribution of instant centers. For a deformable body, the 
distribution might not be anything like a point/line in space. You can imagine 
a pancake being flipped and see that the instantaneous angular velocity at the 
edges (far from the "axis of flip") would oscillate above and below that of the 
portion nearest the "axis of flip" -- because pancakes can curl (taco?) a good 
deal while flipping. How should ParaView report that? As a single angular 
velocity for the pancake, averaged over the entire cake? What if there are 
multiple pancakes? One average for each connected component? One for the full 
stack? :-)

For each set of nodes N_i you want to use in an estimate, you could deflect 
nodal coordinates by their velocity N_i' = N_i + delta_t * V_i and then compute 
the optimal rotation matrix using the Kabsch algorithm:

  http://en.wikipedia.org/wiki/Kabsch_algorithm

>From the optimal rotation matrix, you can choose the eigenvector with the 
>largest eigenvalue and use it as an estimate of the instant center of 
>rotation. With the instant center, you can then decompose each nodal velocity 
>into an angular velocity magnitude and linear velocity magnitude -- and 
>average over all the nodes in a connected component/block/full stack.

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://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Color legend and log scaling

2014-07-01 Thread David Thompson
Hi Alan,

> I would like to propose a few changes in our log scaling algorithm for 
> painting a dataset by a variable.   I discussed this with Utkarsh, and he 
> asked that I bounce it off the e-mail list.  So, here goes.
>  
> Currently, when a user log scales a variable, if all data is positive, 
> ParaView just uses the normal min and max.  There are times when this is not 
> proper – for instance when looking at the temperature or density of material 
> in a supernova, or velocity of outbound gas.  Another example is large data, 
> with noise around zero.  I would like to propose that we have a user 
> selectable option to set the minimum at maximum*10^-q, where q is user 
> defined but defaults to 4.  In other words, the minimum would be set to 
> 1*10^-4 of what the maximum is.

Avoiding a "window" around 0 in the initial view sounds good to me. However, I 
can imagine some cases where the data spans more than 4 orders of magnitude. 
One thing I've seen (debatably bad, but something ParaView must deal with) are 
simulations/datasets that use large numbers (i.e., 1e+38) to mark invalid 
values. (The LDAV climate data does this.) Showing a plot with the initial view 
set to [1e34,1e38] would not be useful, since it would only show invalid 
values. Another is chemical reaction simulations where concentrations span much 
more than 4 orders of magnitude (I've seen some span 11 or 12 orders of 
magnitude, but 5 or 6 can be common).

What choosing q to ensure that a significant fraction (say 90%?) of the data is 
actually on-screen? It not terribly hard (even in parallel) to extract a 
fixed-size sample that approximates a histogram to within a few percent. We 
could use that to determine where the bulk of the data resided and ensure that 
the q-value does not leave more than 10% offscreen.

> If all of the user’s data is negative, ParaView grumbles, and then seg faults 
> using a current master git pull – not optimal behavior.  In PV 4.1, it just 
> sets min and max to 0.  I would like to propose that ParaView calculate the 
> log of the data, as follows:  Index= -(log(abs(Var))).  Then, just draw the 
> color legend as normal – for instance, red at top, white in the middle and 
> blue at bottom.  Tick marks will be the reverse of positive log scaling – 
> with the dense numbers, more negative numbers at the bottom and less dense, 
> less negative numbers at the top.

I'm not sure I understand this, especially your use of "dense". It sounds like 
you have a particular dataset in mind where the probability density is low near 
zero. Are you saying you want the color scale to be different in the case of 
data that is all negative numbers? Or that log plots in general should have 
colors reversed?

> The problem arises with data that spans positive and negative numbers.  Since 
> the log of 0 is infinity, we have to deal with very small numbers in a 
> special way.  I propose that we find maxVal = max(maximum, abs(minimum)).  
> Then, we set the color bar to run from maxVal to –maxVal.  We log scale the 
> top half of the color legend, running from maxVal to maxVal*10^-4, and we 
> reverse log scale the bottom half of the color legend, running from 
> –maxVal*10^-4 to –maxVal.  We calculate this negative range the same as the 
> all negative data section above.  All data between maxVal*10^-4 and 
> -maxVal*10^-4 would remain white by default, or user selectable black.

I like the idea of having a custom color for the range [-10^q,10^q] and I think 
I understand and agree with the simplified bounds/ticks for the case where the 
data crosses zero -- at least when it firmly crosses the origin. However, whose 
chemistry simulations I mentioned above would occasionally have very slightly 
negative concentrations. Physically they can never be negative but 
floating-point precision meant that some were. In this case, we might have a 
distribution of concentrations that went from -1e-16 to 1 with the bulk in the 
range [1e-10,1]. Using the algorithm above for the axes would assign half of 
the color palette to the range [-1,-1e-4], half to [1e-4,1] and draw a 
significant fraction of the data as "black". If instead the histogram was used 
to choose q, then we might be able to decide that all of the "naughty" data 
were outliers and select a log color-scale of [1e-10,1].

David

> User selectable functionality would be as follows:
>   • to allow/ not allow negative numbers (default allow)
>   • to be able to change the q exponent (i.e., 4 above) (default 4)
>   • to be able to change the painting color that is too small (default 
> white)
>   • to clamp minimum to some number (such as q == 4 above).  (default on).
>  
> Thoughts?
>  
> Alan
>  
>  
> ___
> Powered by www.kitware.com
> 
> Visit other Kitware open-source projects 
> athttp://www.kitware.com/opensource/opensource.html
> 
> Please keep messages on-topic and check the ParaView Wiki 
> at:

Re: [Paraview] 3d printing

2014-06-19 Thread David Thompson
Hi Alan,

> Is there any functionality to allow a cool ParaView visualization to be 
> exported so that a 3D printer can print it?

If you triangulate the surface and save as STL, then you can 3-D print the 
result. (STL only accepts triangles and it is up to you to ensure that the 
surface is closed.)

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://public.kitware.com/mailman/listinfo/paraview


Re: [Paraview] Seeking advice for high-order and high-dimensional discontinuous Galerkin data

2014-04-11 Thread David Thompson
Hi Noah,

> ...
>   • The solution on any given element is smooth, represented by arbitrary 
> order polynomials and sampled at surface and interior node locations.  I tend 
> to use high order elements compared to typical practice, i.e. 7th, 8th, 
> 9th-order. I’d like the visualization to capture some of the element’s 
> internal solution structure rather than just a linear interpolation between 
> element vertices.

If you are dealing with meshes where regular sampling of higher order elements 
is feasible and the number of samples per element per axis is high (>8?), you 
might consider using composite datasets (AMR or multiblock datasets), with each 
block holding image data for a single element.

>   • The solution is multi-valued at faces where two elements are adjacent 
> and I want the visualization to well-capture the surface discontinuity.

That can be troublesome with precisely coincident grid points since rendering 
just those boundaries will exhibit z-fighting. If you isocontour, you will see 
cracks at element boundaries which reflect discontinuities. A better 
understanding of what you intend to visualize might help decide how to organize 
the data.

> ... Obviously ParaView is focused on 2D and 3D data. The underlying VTK looks 
> to only handle up to 3D vertex position vectors, for example.  My takeaway 
> from this, is it is not possible to describe a 4D or higher mesh to ParaView 
> directly.

No.

> ... My survey of the documentation makes me think that I’ll be able to 
> provide parameters to the reader plugin that control the dataset.  For 
> example, specify what dimensions to slice and where to slice.  I suppose I 
> need to similarly parameterize the load of the mesh.  Does this sound like a 
> workable approach?

Yes. You may also want to take a look at 
http://www.vtk.org/Wiki/VTK/InSituDataStructures and vtkTypedDataArray as a way 
to provide 3-d point coordinates from a higher-dimensional array (and the 
remaining dimensions as other independent arrays). But it's unclear to me -- 
assuming you regularly sample phase space the same way you sample positional 
space -- how you would handle places where the projection to a subset of those 
dimensions resulted in coincident elements in the projected space.

> For the data discontinuities at element faces, I’ve seen discussions in this 
> forum that suggest unstructured grids with duplicated/co-located vertices 
> allows for assigning two different data values at the same physical location. 
>  Since I am working with structured grids and data, I’m wondering if a 
> different technique might work.  How about inserting zero-volume collapsed 
> elements between each real element?

I think this would be problematic for many VTK filters.

> Lastly, for the high order visualization, ... the most straight forward 
> approach is to interpolate my structured high order elements... 
> vtkBridgeDataSet has been mentioned as a possible approach but requiring much 
> development effort.  That was several years ago.  Has there been any progress 
> by others with similar high-order or spectral data sets?

Sadly, no. Although we have had a lot of interest from academia, there haven't 
been paying customers for this.

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] How to indicate cell normals.

2013-09-23 Thread David Thompson
Hi Yusuke,

> I am able to get to point data of Normals.
> But if I select Normals, it is not displayed in paraview.( it means that 
> behavior of the normal is not observed in paraview)
> I analyze the axial symmetry model.
> I use Extract Surface filter and generate surface normal filter.
> Should  I use other filter??

If you want normal glyphs (arrows drawn at points showing the direction of the 
normals), then you should use the Glyph filter after Generate Surface Normals.

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] How to indicate cell normals.

2013-09-20 Thread David Thompson
Hi Yusuke,

Your ls-dyna dataset probably has volumetric cells instead of a surface mesh. 
Run the Extract Surface filter on your data first and then try generating 
normals.

David

On Sep 21, 2013, at 0:12, uverworld_1...@yahoo.co.jp wrote:

> Hi! All
> 
> I use ls-dyna and I open d3plot file in paraview.
> I want to indicate cell normals.
> So, I try to apply the generate surface normals.
> However, I can't apply it.
> 
> Please tell me how to indicate cell normals vectors and contour in paraview.
> 
> Thank you
> Yusuke
> ___
> 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] d3plot transform to vtk file

2013-09-16 Thread David Thompson
Hi Yusuke,

I'm not sure what you mean by "central axis". Do you mean the central axis of 
the simulation or the middle layer of a thick shell element? If the latter, 
note that the LS-DYNA file format only stores fluid pressure information at 
nodes or on SPH particles. If the shell elements do not have nodes on their 
middle layer, then no pressure information is available there and you will have 
to interpolate it.

Also, are you sure the pressure variable is being saved in the d3plot file? 
LS-DYNA may be configured to save pressure or not. 

Finally, there are other files (the time history database: d3thdt; the 
interface force database: intfor) that VTK's LS-DYNA does not support but that 
LS-PrePost is able to read. If the pressure information you want is stored in 
these files, VTK will not be able to read it.

David

On Sep 15, 2013, at 12:43 PM, uverworld_1...@yahoo.co.jp wrote:

> Dear all 
> 
> I want to observe pressure fluctuation on central axis in paraview.
> I use LS-DYNA.
> In ls-prepost, I can't observe pressure fluctuation on central axis. So, I 
> decided to use paraview.
> I can indicate pressure to use ls-prepost but I can't it in paraview.
> So, I  thought if d3plot transformed to vtk format, I can indicate pressure 
> in paraview. 
> I want to convert d3polt file to vtk format to indicate pressure in paraview.
> However I don't understand that methods.
> So, please tell me how to convert d3plot file to vtk.
> Should I use lsdynareader??
> 
> Thank you 
> Best regards
> Yusuke
> ___
> 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] d3plot to vtk

2013-09-13 Thread David Thompson
Hi,

> I use LS-DYNA.
> I want to convert d3polt file to vtk format.
> However I don't understand that methods.
> So, please tell me how to convert d3plot file to vtk.

I'm not sure I understand; are you having problems loading the data into 
ParaView or saving the results as a VTK dataset?

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] Tube Filter for ExodusII-BEAM

2013-09-12 Thread David Thompson
Hi Ken,

>> Hmm, I'm not sure you have access to the element attribute or block
>> properties within ParaView. ...
>> Dave, do you know if there is a way to get at this information from an
>> Exodus file?
> 
> Here is a small patch that will get the VTK reader loading data. A vtkpython 
> script is also attached that demonstrates it on the exEdgeFaceElem.exii test 
> file in vtkData.

I found some of my copious spare time lying around and submitted a patch to 
Gerrit (topic 3235). If you can review, I'll merge.

Thanks,
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] Tube Filter for ExodusII-BEAM

2013-08-26 Thread David Thompson
Hi Ken and Sven,

> Hmm, I'm not sure you have access to the element attribute or block
> properties within ParaView. ParaView is really optimized for reading
> results data.
> 
> Dave, do you know if there is a way to get at this information from an
> Exodus file?


Here is a small patch that will get the VTK reader loading data. A vtkpython 
script is also attached that demonstrates it on the exEdgeFaceElem.exii test 
file in vtkData.

However, I leave it as an exercise to the reader to expose things in 
ParaView... it was hard coming up with all the little icons for different types 
of data arrays and now you need 3 more (edge, face, and element block 
attributes)! Plus you may need to add some convenience methods to the 
vtkExodusIIReader class to expose the different types of attributes in a way 
that lets ParaView map them to properties.

David



0001-Read-requested-Exodus-block-attributes.patch
Description: Binary data
#!/usr/bin/env vtkpython
from vtk import *
rdr = vtkExodusIIReader()
rdr.SetFileName('/path/to/VTKData/Data/edgeFaceElem.exii')
rdr.UpdateInformation()
rdr.SetObjectAttributeStatus(rdr.ELEM_BLOCK,0,'SPAGHETTI',1)
rdr.SetObjectAttributeStatus(rdr.ELEM_BLOCK,0,'WESTERN',1)
rdr.Update()
na = rdr.GetOutputDataObject(0).GetBlock(0).GetBlock(0).GetCellData().GetNumberOfArrays()
for i in range(na):
  print rdr.GetOutputDataObject(0).GetBlock(0).GetBlock(0).GetCellData().GetArray(i).GetName()
___
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] ANN: ParaView 4.0.0-RC1 now available for download

2013-05-30 Thread David Thompson
Hi all,

> May be legend's font size? It's very small (default size is 12)

I have found the cause of this and have a fix, but it may not make it into the 
4.0 release.

David


> 
> On Thu May 23 2013 21:59:04 Utkarsh Ayachit wrote:
>>> minor, the font size for contours is extremely small, would it be
>>> possible to increase the default value closer to the source text font
>>> size.
>> 
>> I am not sure I follow, which font size are you referring to? On the
>> contours panel? Can you attach a screenshot to make things clearer,
>> please? Thanks.
>> 
>> Utkarsh
>> ___
>> 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


Re: [Paraview] segfault opening the color scale editor / vtkPVProminentValuesInformation?

2013-05-06 Thread David Thompson
Hi Burlen,

I'm working on it. It shouldn't take too long to fix. Thanks for the report.

David

On May 6, 2013, at 1:13 PM, Burlen Loring  wrote:

> Steps to reproduce:
> 
> Start pvserver running parallel 2 ranks
> Open the attached dataset, load b array
> Slice the dataset such that the slice contains no data from at least one rank.
> Color slice by b array magnitude (default)
> Open the color scale editor -> crash
> 
> I added the test dataset to the bug report: 
> http://paraview.org/Bug/view.php?id=14051
> 
> On 05/03/2013 04:12 PM, Burlen Loring wrote:
>> with git master from this afternoon, I'm experiencing segfaults when I try 
>> to open the color scale editor. seems to be related to 
>> vtkPVProminentValuesInformation. As soon as I get 4 or more ranks I have the 
>> segfault. Is this a known issue? 
>> 
>> Program received signal SIGSEGV, Segmentation fault. 
>> 0x7fde3f6dbee1 in vtkInformationVector::GetInformationObject (this=0x0, 
>> index=0) at 
>> /work/ext/ParaView/sqtk-pv/ParaView/VTK/Common/Core/vtkInformationVector.cxx:157
>> 157   if(index >= 0 && index < this->NumberOfInformationObjects) 
>> (gdb) where 
>> #0  0x7fde3f6dbee1 in vtkInformationVector::GetInformationObject 
>> (this=0x0, index=0) at 
>> /work/ext/ParaView/sqtk-pv/ParaView/VTK/Common/Core/vtkInformationVector.cxx:157
>>  
>> #1  0x7fde3f5c5950 in vtkAbstractArray::UpdateDiscreteValueSet 
>> (this=0x3ce0130, uncertainty=0, minimumProminence=0) at   
>> /work/ext/ParaView/sqtk-pv/ParaView/VTK/Common/Core/vtkAbstractArray.c 
>> #2  0x7fde3f5c4d81 in vtkAbstractArray::GetProminentComponentValues 
>> (this=0x3ce0130, comp=-1, values=0x3cd2f30, uncertainty=0, 
>> minimumProminence=0) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/VTK/Common/Core/vtkAbstractArray.cxx:548
>> #3  0x7fde383df72d in 
>> vtkPVProminentValuesInformation::CopyDistinctValuesFromObject 
>> (this=0x3ccfeb0, array=0x3ce0130) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/ParaViewCore/ClientServerCore/Rendering/vtkPVProminentValuesInformation.cxx:391
>> #4  0x7fde383df5b7 in 
>> vtkPVProminentValuesInformation::CopyFromLeafDataObject (this=0x3ccfeb0, 
>> dobj=0x4084620) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/ParaViewCore/ClientServerCore/Rendering/vtkPVProminentValuesInformation.cxx:360
>> #5  0x7fde383df2ad in 
>> vtkPVProminentValuesInformation::CopyFromCompositeDataSet (this=0x3ccbd80, 
>> cds=0x40bb800) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/ParaViewCore/ClientServerCore/Rendering/vtkPVProminentValuesInformation.cxx:309
>> #6  0x7fde383def5f in vtkPVProminentValuesInformation::CopyFromObject 
>> (this=0x3ccbd80, obj=0x405dd90) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/ParaViewCore/ClientServerCore/Rendering/vtkPVProminentValuesInformation.cxx:263
>> #7  0x7fde401554f0 in vtkPVSessionCore::GatherInformationInternal 
>> (this=0x23a3b90, information=0x3ccbd80, globalid=4483) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/ParaViewCore/ServerImplementation/Core/vtkPVSessionCore.cxx:838
>> #8  0x7fde401555cb in vtkPVSessionCore::GatherInformation 
>> (this=0x23a3b90, location=17, information=0x3ccbd80, globalid=4483) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/ParaViewCore/ServerImplementation/Core/vtkPVSessionCore.cxx:858
>> #9  0x7fde401512c0 in vtkPVSessionBase::GatherInformation 
>> (this=0x23a9ad0, location=17, information=0x3ccbd80, globalid=4483) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/ParaViewCore/ServerImplementation/Core/vtkPVSessionBase.cxx:258
>> #10 0x7fde40160b7e in vtkPVSessionServer::GatherInformationInternal 
>> (this=0x23a9ad0, location=17, classname=0x3cbb9f8 
>> "vtkPVProminentValuesInformation", globalid=4483, stream=...) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/ParaViewCore/ServerImplementation/Core/vtkPVSessionServer.cxx:567
>> #11 0x7fde401608f0 in vtkPVSessionServer::OnClientServerMessageRMI 
>> (this=0x23a9ad0, message=0x408dec0, message_length=135) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/ParaViewCore/ServerImplementation/Core/vtkPVSessionServer.cxx:527
>> #12 0x7fde4015ee0c in (anonymous namespace)::RMICallback 
>> (localArg=0x23a9ad0, remoteArg=0x408dec0, remoteArgLength=135)   
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/ParaViewCore/ServerImplementation/Core/vtkPVSessionServer.cxx:56
>> #13 0x7fde3fb0c4bf in vtkMultiProcessController::ProcessRMI 
>> (this=0x2cdfe60, remoteProcessId=1, arg=0x408dec0, argLength=135, 
>> rmiTag=55625) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/VTK/Parallel/Core/vtkMultiProcessController.cxx:650
>> #14 0x7fde3fb0c06e in vtkMultiProcessController::ProcessRMIs 
>> (this=0x2cdfe60, reportErrors=0, dont_loop=1) 
>> at 
>> /work/ext/ParaView/sqtk-pv/ParaView/VTK/Parallel/Core/vtkMultiProcessController.cxx:597
>> #15 0x7fde3fe64291 in vtkTCPNetworkAccessManager::ProcessEventsInternal 
>> (this=0x2394e90, timeout_msecs=

Re: [Paraview] writing volume to stl as topology

2013-04-15 Thread David Thompson
Hi Dmitry,

> I have a tetrahedral mesh in vtk unstructured grid. How can I write this 
> volume to stl (not just the surface) like a set of separated triangles? 
> Paraview display my data as a surface. Here is some screenshots of my model 
> and clipped model. My goal is just to handle this topology and then print it 
> by 3D printer.

The STL format only accepts triangles, not tetrahedra or other volumetric 
elements. Your 3D printer software will infer what regions are inside and 
outside the surface as long as the surface is manifold (i.e., a closed, 
orientable surface that is watertight). If you have a model that is not 
watertight, you can try using ParaView's "Clean" filter to fix it.

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] Using a programmable filter to apply custom rgb colors to a ExodusII data set

2013-03-04 Thread David Thompson
Hi Mark,

> I'm trying to use a programmable filter to apply colors to an ExodusII data 
> set.  ... From that what I assume I have to do is create a 
> vtkUnsignedCharArray with length equal to the number of cells in the model, 
> assign my colors, and then display.

Yes, you are correct.

> 1) I understand the ExodusIIReader is, well, a reader -- it will generate the 
> data blocks, it's not actually representing the data.

I'm not sure what you mean by "not actually representing the data." It *does* 
represent the data, segregated into a series of blocks (which is how it is 
stored on disk). The blocks are arranged into a 2-level hierarchy. The first 
level determines what the blocks represent (elements, faces, edges, or subsets 
thereof). The second level contains the actual data.

>  But I don't get how to access at the generated data.  Say I have a single 
> block of elements called 'hexes' with cell data 'stresses', how do I access 
> that cell data in the programmable filter?

I've put a script below that colors the first element block in the can.ex2 
dataset (available from the ParaViewData git repo: 
http://www.paraview.org/gitweb?p=ParaViewData.git;a=blob;f=Data/can.ex2;hb=HEAD 
). It shows how you can add a color array to one block (in the example below, 
ids.GetBlock(0).GetBlock(0) is the first "leaf" node of the tree containing the 
Exodus data). In reality, you probably want to loop over all leaf nodes in the 
dataset. You can use GetNumberOfBlocks() to set up a loop that calls GetBlock 
with each valid index.

David

ids = self.GetInput()
ods = self.GetOutput()

ods.ShallowCopy(ids)

cellData = ids.GetBlock(0).GetBlock(0).GetCellData()
eqps = cellData.GetArray('EQPS')
nt = eqps.GetNumberOfTuples()

color = vtk.vtkUnsignedCharArray()
color.SetName('Color')
color.SetNumberOfComponents(3)
color.SetNumberOfTuples(nt)
for i in range(nt):
   ival = eqps.GetTuple1(i)
   red = 255.0 * ival
   green = 255.0 * (1 - ival)
   blue= 255.0
   color.SetTuple3(i, red, green, blue)

cellData.SetScalars(color) 
___
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] Where's Display tab?

2013-01-12 Thread David Thompson

> How can I change display and Style settings such as line width, etc, in the 
> recent ParaView (v 3.98.0)?

The Display tab and filter Properties tab have been merged in v3.98. Also, 
infrequently-used settings only show up when the "advanced" button at the top 
right of the combined filter-properties/display-settings tab is clicked.

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] Categorical data in scalar bar

2012-12-31 Thread David Thompson
Hi Tim,

> We have a dataset that ranges from 0 to 7 and represents categorical data. Is 
> it possible through a python script to explicitly set the tick labels on a 
> scalar bar so I can manually map them to strings? For instance:
> 
> {0.0 -> "None", 1.0 -> "I", 2.0 -> "J",...} 

In version 3.98 of Paraview, you can do the following at the Python shell in 
ParaView:

# Unfortunately, there is not an easy way to access the Brewer color presets 
from Python;
# Here is one of the categorical schemes:
nanColor=[0.7019607843137254, 0.7019607843137254, 0.7019607843137254]
rgbPoints=[0, 0.4,0.7607843137254902, 0.6470588235294118, \
   1, 0.9882352941176471, 0.5529411764705883, 0.3843137254901961, \
   2, 0.5529411764705883, 0.6274509803921569, 0.796078431372549, \
   3, 0.9058823529411765, 0.5411764705882353, 0.7647058823529411, \
   4, 0.6509803921568628, 0.8470588235294118, 0.32941176470588235, \
   5, 1.0,0.8509803921568627, 0.1843137254901961, \
   6, 0.8980392156862745, 0.7686274509803922, 0.5803921568627451, \
   7, 0.7019607843137254, 0.7019607843137254, 0.7019607843137254]
# Create a vector of strings holding scalar values and their annotations:
annotations = ['0', 'None', '1', 'I', '2', 'J', '3', 'K', '4', 'L', '5', 'M', 
'6', 'N', '7', 'P']
# Create a categorical lookup table
lkup = CreateLookupTable(\
  RGBPoints=rgbPoints, NanColor=nanColor, IndexedLookup=1, 
Annotations=annotations, ColorSpace="HSV")
# Color by the new lookup table:
SetDisplayProperties(LookupTable = lkup)
SetDisplayProperties(ColorAttributeType = "CELL_DATA") # or "POINT_DATA"
SetDisplayProperties(ColorArrayName = "YourCategoricalArray")
Render()


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] How to Rename Blocks in Exodus File

2012-12-11 Thread David Thompson
Hi Dennis,

> I would like to rename the Exodus blocks to useful names, in a macro script, 
> reading the names from a file. ...

Before doing this in ParaView, does it make sense to change the Exodus file 
instead? That way the names would stay changed across ParaView sessions and any 
other applications that use the file will have access to them. You should be 
able to use netcdf4-python ( http://code.google.com/p/netcdf4-python/ ) to 
get/set the element block names (since Exodus files are just netCDF files that 
follow a particular naming convention). If the "eb_names" variable exists in 
the netCDF file, you can alter it. Otherwise, you can add it and save the file. 
See the VTKData repository ( http://vtk.org/VTKData.git ) for an example file 
(edgeFaceElem.exii) that has names in it. The ncdump and ncgen programs that 
come with NetCDF are also useful for editing small files; I would advise using 
them to examine and modify edgeFaceElem.exii before writing a Python script to 
change things.

> so using src=ExodusIIReader I can use SetData on src.ElementBlocks and change 
> that - which accomplishes nothing except de-activating all my blocks

I'm not sure how "src=ExodusIIReader; src.ElementBlocks" works at all. Are you 
using Tools->Python Shell to run Python commands? Or a Programmable Filter?

> However,  I need to change src.ElementBlocks.Available and I can't figure out 
> how to do that as there is no SetData and nothing else I tried works


Given my confusion above, I cannot be sure, but it looks like you are trying to 
modify the reader and not the data that it loads. But the reader in ParaView is 
part of a pipeline that creates the dataset and modifies it as the user 
requests changes (e.g., load variable "A", turn off block "B", and so on), so 
even if you were to modify the reader's output data, it would not stay changed.

If you want the names to be different in filters that connect to the reader, 
you'll have to create an intermediate Programmable Filter that runs a Python 
script to change block names. See 
http://www.paraview.org/Wiki/Python_Programmable_Filter for more information, 
but be aware that you'll need some knowledge of how VTK multiblock datasets are 
named to make this work.

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] Statistics/Mesh quality

2012-12-11 Thread David Thompson
Hi Matthias,

> I would be interested in learning more about mesh quality evaluation, and 
> would appreciate a pointer to some literature or other information sources.

This wiki page http://www.vtk.org/Wiki/VTK/mesh_quality has a link to user 
documentation for the quality library VTK uses (named Verdict). The 
documentation has references to papers on the different quality measures.

> Which formats can be used to load the mesh into ParaView?

The File->Open dialog has a pull-down list at the bottom describing the 
different file formats ParaView supports.

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] Statistics/Mesh quality

2012-12-10 Thread David Thompson
> I am using Paraview to evaluate mesh quality (scaled Jacobian measurement) 
> and want to know, how many elements have a larger quality measurement than, 
> e.g., 0.8, 0.6, 0.4 and 0.2. Is there an easy way to extract this 
> information? Thank you very much!

You could run the histogram filter on the dataset but that won't leave you with 
complete control over the bin cutoffs for quality values. It's probably what 
you want, though.

You can also create a Threshold filter for each quality (0.8, 0.6, 0.4, 0.2) of 
interest and look at the information tabs for a report of the number of cells. 
That is awkward but you get to pick exactly the quality values of interest, 
even if they are non-uniform.

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] Setting the colour/scale palete used for glyphs

2012-11-19 Thread David Thompson
Hi Kit,

> ... I would like to ensure that field is coloured (by default) using a 
> specific colour palette. More specifically I want to it to use 2 discrete 
> colours, rather than the rainbow scale. Is this possible?

You can set the number of colors in the lookup table by
+ clicking on the color scale editor button in the toolbar,
+ ensuring that the "Use Discrete Colors" checkbox is ON, and
+ changing the "Resolution" slider to 2.

As long as the field values are uniformly spaced, you can handle more than 2 
colors by changing the resolution appropriately (but note that you may wish to 
use the "Rescale Range" button to keep color band edges from falling on special 
values).

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] Mesh Quality

2012-11-09 Thread David Thompson
Hi Fastl,

>  I am using the Filter called 'Mesh Quality' with the 'Scaled Jacobian' 
> option for hexahedral elements. I have some elements which have only six 
> different vertices (instead of eight, which leads to a wedge) and the result 
> is a value like 1e+30. I want to analyze this result, because in general the 
> range for a 'Scaled Jacobian' is between -1 and +1. Can you provide me with 
> the equation how this 'Scaled Jacobian' is evaluated in Paraview? 

ParaView uses the Verdict mesh quality library to evaluate the scaled Jacobian 
quality function. The code for this is located here:

  
http://vtk.org/gitweb?p=VTK.git;a=blob;f=ThirdParty/verdict/vtkverdict/V_HexMetric.cpp;hb=HEAD#l1526

> ... I would expect to get an error, but instead I think I get extremely high 
> values.

Verdict substitutes a very large number (VERDICT_DBL_MAX) in the place of 
errors.

Hope this helps,
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] Converting an XML colour scale to a list for RGBPoints in batch

2012-11-07 Thread David Thompson
Hi Joe,

> In the GUI, this is easy, you just click "Edit color map" -> "Choose preset" 
> -> "Import", select the xml file and you're away.  I can't seem to be able to 
> do this in batch though.  Is there any way I can?

I don't know of a way offhand, but you might try using Tools->Start Trace and 
Tools->Stop Trace to load a dataset and import a colormap to color it. That's 
how I discover Python equivalents for many things in ParaView.

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] Can you use bands in colour maps?

2012-11-06 Thread David Thompson
Hi Joe,

> By default, all the colour maps in ParaView are continuous between points.  
> I'd like to have only one colour in between points I set in the XML map, is 
> this possible?

In the current release of ParaView, you can click on the "Edit color map" 
button in the toolbar and turn on "Use Discrete Colors" in the color editor. By 
changing the number of colors with the neighboring slider, uniform ranges of 
scalar values will be assigned a single color.

In an upcoming release, discrete color maps will be available. Discrete 
colormaps do not interpolate colors; instead they are looked up in a table, so 
a single large difference in values will not force any small variations to be 
lumped into a single color band. However, this requires your scalar field to 
take on a small number of distinct values (32 or fewer); colors are assigned to 
a single value, not a range.

To map non-uniform ranges of values to a single color each, you really should 
create a derived, integer-valued variable and use it for coloring -- regardless 
of whether you use continuous colormaps in the current release or plan to use 
categorical color maps later. The calculator filter's floor, ceil, and round 
functions might be worth a try in creating the derived field; if there's not a 
simple expression to convert each range of interest to an integer, the Python 
calculator will let you apply a tabular function.

Finally, it *may* also be possible to achieve what you desire with a custom XML 
color map by placing control points with different colors very close to each 
other to create a discontinuous color map. This would be difficult to adjust in 
the color editor should the control points need to be repositioned, but might 
be the easiest thing to do if the bands are fixed values or fixed fractions of 
the total range.

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] Getting String arrays out of Field Data

2012-09-06 Thread David Thompson
Hi Joshua,

Sorry, in my previous e-mail I forgot this was the ParaView and not VTK mailing 
list, so I posted a vtkpython script. What I said still applies if you are 
using the Python Programmable Filter, but of course you can use ParaView to 
load the dataset.

pdi = self.GetPolyDataInput()
pdo =  self.GetPolyDataOutput()
s=pdi.GetCellData().GetAbstractArray(0)
print s.GetVariantValue(1).ToString()
# Just copy the input to make ParaView happy:
pdo.DeepCopy( pdi )

This works for me with the political.vtp dataset

  http://vtk.org/gitweb?p=VTKData.git;a=blob_plain;f=Data/political.vtp;hb=HEAD


Does that address your problem? Or are you writing client-side python scripts 
using the simple interface to retrieve metadata from the server?

David

>> I have written a custom filter for Heliospheric data, and I am placing the 
>> model metadata into the Field Data section of ParaView.
>> 
>> Many of the metadata items are strings. I can view these strings just fine 
>> in Spread Sheet view, but anytime I try to get them out of Field Data via 
>> python, I get a null value return.  Doubles and Ints return just fine, but I 
>> cannot seem to figure out how to get the strings.
>> 
>> Can someone point me in the right direction for retrieving the strings from 
>> Field Data via Python?
> 
> What method are you using to get the array from the vtkDataSetAttributes 
> instance that holds the field data? You need to use GetAbstractArray instead 
> of GetArray. Below is a small example that gets string values from cell data 
> on the political.vtp file included in VTKData. If you have an example file 
> with strings in field data, I can take a look at that, too. It should be as 
> simple as replacing the GetCellData() method with GetFieldData().
> 
>   David
> 
> from vtk import *
> rdr = vtkXMLPolyDataReader()
> rdr.SetFileName( '/path/to/VTKData/Data/political.vtp' )
> rdr.Update();
> mesh = rdr.GetOutput()
> cd = mesh.GetCellData()
> print 'Should have 3 cell data arrays', cd.GetNumberOfArrays()
> cd.GetArray(0) # Returns None
> sa = cd.GetAbstractArray(0)
> print 'Array has %d tuples (should be 980)' % sa.GetNumberOfTuples()
> print 'Should be Albania:', sa.GetVariantValue(1).ToString()
> 

___
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 String arrays out of Field Data

2012-09-06 Thread David Thompson
Hi Joshua,

> I have written a custom filter for Heliospheric data, and I am placing the 
> model metadata into the Field Data section of ParaView.
> 
> Many of the metadata items are strings. I can view these strings just fine in 
> Spread Sheet view, but anytime I try to get them out of Field Data via 
> python, I get a null value return.  Doubles and Ints return just fine, but I 
> cannot seem to figure out how to get the strings.
> 
> Can someone point me in the right direction for retrieving the strings from 
> Field Data via Python?

What method are you using to get the array from the vtkDataSetAttributes 
instance that holds the field data? You need to use GetAbstractArray instead of 
GetArray. Below is a small example that gets string values from cell data on 
the political.vtp file included in VTKData. If you have an example file with 
strings in field data, I can take a look at that, too. It should be as simple 
as replacing the GetCellData() method with GetFieldData().

David

from vtk import *
rdr = vtkXMLPolyDataReader()
rdr.SetFileName( '/path/to/VTKData/Data/political.vtp' )
rdr.Update();
mesh = rdr.GetOutput()
cd = mesh.GetCellData()
print 'Should have 3 cell data arrays', cd.GetNumberOfArrays()
cd.GetArray(0) # Returns None
sa = cd.GetAbstractArray(0)
print 'Array has %d tuples (should be 980)' % sa.GetNumberOfTuples()
print 'Should be Albania:', sa.GetVariantValue(1).ToString()

___
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

2012-08-17 Thread David Thompson
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] Numpy masking (via programm filter) not quite working in parallel

2012-08-17 Thread David Thompson
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  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] Frustum selection display bug in RC3

2011-11-11 Thread David Thompson

Was anyone able to reproduce this? If so, I'll file a bug.

I can reproduce it, yep.

Thanks. I've filed bug 12721:
  http://paraview.org/Bug/view.php?id=12721

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] Frustum selection display bug in RC3

2011-11-10 Thread David Thompson

Was anyone able to reproduce this? If so, I'll file a bug.

David

On Nov 4, 2011, at 9:39 , Thompson, David C wrote:


Hi all,

I've been fiddling with ParaView 3.12 RC3 and found a small bug...  
when the frustum selection is being rendered as an outline and you  
delete the target of the selection, the frustum sticks around and  
cannot be "unshown". Attached is a test file illustrating the  
problem (created on Mac OS X 10.7).


   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] Concept: A world without the Apply button

2011-06-01 Thread David Thompson

Utkarsh,

I agree that composite readers (and probably other composite filters)  
can provide the dataset structure in the RequestDataObject pass, but I  
don't know that many currently do. Certainly the ExodusIIReader  
doesn't. It would be easy enough to provide an empty hierarchy of *all  
the available* blocks the file(s) contain, but doing that with all of  
the block/set selections plus populating empty arrays on each leaf  
(which I think would be useful) would be a lot of work.


David

On Jun 1, 2011, at 13:33 , Utkarsh Ayachit wrote:


Dave,

Output types become a problem only for composite-datasets, for others
data-types are known after RequestDataObject pass which happens before
RequestData(). I am still not sure how (or if) to deal with composite
datasets in this case. We can expect the reader to give meta-data for
composite datasets and leaves in  the dataset in more details as well,
but I am not sure we need to yet.

As far as the scalars/vectors/tensors go, no filter in paraview
requires active-scalars/tensors etc, it just checks for existence of
array with a give number of components. The meta-data should indeed
contain information about arrays and number of components.

Utkarsh

On Wed, Jun 1, 2011 at 4:09 PM, David Thompson   
wrote:
I like the idea. One possible problem that comes to mind is  
filters that
produce different output types based on the input. Will enough of  
the
metadata be passed through the pipeline before Execute() is called  
for those
filters to hint at their output type? Otherwise, you'll be forced  
to hit

Execute before adding any more filters.


Or, for that matter, what about filters which are disabled in the  
menu when

there is no scalar/vector/tensor data on their input?

   David


On Jun 1, 2011, at 12:30 , Utkarsh Ayachit wrote:


Folks,

I have been investigating the ability to create pipelines in  
ParaView

without having to hit Apply at every single stage in the pipeline.
This would make it much easier to deal with really large data in a
usage scenario where the scientist is fairly familiar with the
visualization pipeline he's interested in setting up.

Following the recent changes to the ParaView ServerManager, the
ServerManager level requirement that filters need to be updated  
before
connecting has disappeared. However, we still need to "Apply" to  
get

updated information about arrays etc.

After a couple on discussions with Berk and Dave, I've consolidated
the thoughts on a Wiki page. If any one has any thoughts to  
share, it

would be great.


http://www.paraview.org/ParaView3/index.php/No-Apply_Mechanism_for_Creating_Pipelines

Note that this is merely a concept. We have no plans of  
implementing

it in near future (3.12/4.0).

Utkarsh
___
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] Concept: A world without the Apply button

2011-06-01 Thread David Thompson
I like the idea. One possible problem that comes to mind is filters  
that produce different output types based on the input. Will enough  
of the metadata be passed through the pipeline before Execute() is  
called for those filters to hint at their output type? Otherwise,  
you'll be forced to hit Execute before adding any more filters.


Or, for that matter, what about filters which are disabled in the menu  
when there is no scalar/vector/tensor data on their input?


David


On Jun 1, 2011, at 12:30 , Utkarsh Ayachit wrote:


Folks,

I have been investigating the ability to create pipelines in ParaView
without having to hit Apply at every single stage in the pipeline.
This would make it much easier to deal with really large data in a
usage scenario where the scientist is fairly familiar with the
visualization pipeline he's interested in setting up.

Following the recent changes to the ParaView ServerManager, the
ServerManager level requirement that filters need to be updated  
before

connecting has disappeared. However, we still need to "Apply" to get
updated information about arrays etc.

After a couple on discussions with Berk and Dave, I've consolidated
the thoughts on a Wiki page. If any one has any thoughts to share, it
would be great.

http://www.paraview.org/ParaView3/index.php/No-Apply_Mechanism_for_Creating_Pipelines

Note that this is merely a concept. We have no plans of implementing
it in near future (3.12/4.0).

Utkarsh
___
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] Concept: A world without the Apply button

2011-06-01 Thread David Thompson

Hi Utkarsh,

I like the idea. One possible problem that comes to mind is filters  
that produce different output types based on the input. Will enough of  
the metadata be passed through the pipeline before Execute() is called  
for those filters to hint at their output type? Otherwise, you'll be  
forced to hit Execute before adding any more filters.


David

On Jun 1, 2011, at 12:30 , Utkarsh Ayachit wrote:


Folks,

I have been investigating the ability to create pipelines in ParaView
without having to hit Apply at every single stage in the pipeline.
This would make it much easier to deal with really large data in a
usage scenario where the scientist is fairly familiar with the
visualization pipeline he's interested in setting up.

Following the recent changes to the ParaView ServerManager, the
ServerManager level requirement that filters need to be updated before
connecting has disappeared. However, we still need to "Apply" to get
updated information about arrays etc.

After a couple on discussions with Berk and Dave, I've consolidated
the thoughts on a Wiki page. If any one has any thoughts to share, it
would be great.

http://www.paraview.org/ParaView3/index.php/No-Apply_Mechanism_for_Creating_Pipelines

Note that this is merely a concept. We have no plans of implementing
it in near future (3.12/4.0).

Utkarsh
___
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] Maybe a silly question?

2011-03-23 Thread David Thompson
Thanks again, and I took a look on the programmable filter code. In  
my case the block consists of more than 1 million tetrahedrons and  
it will be time consuming to run iterations. Is there a way that I  
can export the cells and their corresponding point information to a  
numpy array? The ultimate goal is to calculate the volume of each  
tetrahedron.


Your help is much appreciated!


There is a VTK-numpy bridge which is documented on the wiki. I don't  
recall exactly where but I'm sure a search engine will find it for you.


David

On Wed, Mar 23, 2011 at 4:45 PM, David Thompson   
wrote:


I have an Exodus II type file loaded into Paraview and trying to get  
the point ID for a group of cells. I know for vtk type object there  
are functions like GetCell(), GetPoint()


I assume you're talking about using Python from ParaView?

Yup, that's exactly what I am doing here.

There are 2 ways to use Python in ParaView: writing a programmable  
filter and scripting the user interface. The former is evaluated on  
the server where you can get access to (at least a portion of) the  
actual dataset. The latter is evaluated on the client, which does  
not have access to the datasets -- instead it has proxy objects that  
send messages to the servers, telling them what to do. So, you can  
access the points and cells from within a programmable filter, but  
not from within a client-side script.




and just wonder if there are such kind of routines available for  
Exodus II files.


Yes, the ExodusII reader outputs VTK objects in a multiblock  
dataset. You can use GetCell(), GetPoint(), and such on each  
individual block... you just have to decide on a particular block.  
The ExodusII reader separates things into blocks because each one  
may have different point and/or cell variables. There are actually 2  
levels of blocks; the top level segregates datasets into groups that  
correspond to exodus blocks vs. sets. The second level breaks  
individual blocks or sets into separate datasets. Does that help?



I am new in dealing with ex2 file, could you please be detailed?  
There is just one block in the dataset, so how to get the vtk object  
out from it so I can access the points and cells?



See the section of this page
 http://www.vtk.org/Wiki/Python_Programmable_Filter
named "Dealing with Composite Datasets". It describes how to iterate  
over all of the blocks in the dataset and run a function (in this  
case, one named "flatten") on each block.


   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] Maybe a silly question?

2011-03-23 Thread David Thompson


I have an Exodus II type file loaded into Paraview and trying to  
get the point ID for a group of cells. I know for vtk type object  
there are functions like GetCell(), GetPoint()



I assume you're talking about using Python from ParaView?


Yup, that's exactly what I am doing here.


There are 2 ways to use Python in ParaView: writing a programmable  
filter and scripting the user interface. The former is evaluated on  
the server where you can get access to (at least a portion of) the  
actual dataset. The latter is evaluated on the client, which does not  
have access to the datasets -- instead it has proxy objects that send  
messages to the servers, telling them what to do. So, you can access  
the points and cells from within a programmable filter, but not from  
within a client-side script.




and just wonder if there are such kind of routines available for  
Exodus II files.


Yes, the ExodusII reader outputs VTK objects in a multiblock  
dataset. You can use GetCell(), GetPoint(), and such on each  
individual block... you just have to decide on a particular block.  
The ExodusII reader separates things into blocks because each one  
may have different point and/or cell variables. There are actually  
2 levels of blocks; the top level segregates datasets into groups  
that correspond to exodus blocks vs. sets. The second level breaks  
individual blocks or sets into separate datasets. Does that help?




I am new in dealing with ex2 file, could you please be detailed?  
There is just one block in the dataset, so how to get the vtk object  
out from it so I can access the points and cells?



See the section of this page
  http://www.vtk.org/Wiki/Python_Programmable_Filter
named "Dealing with Composite Datasets". It describes how to iterate  
over all of the blocks in the dataset and run a function (in this  
case, one named "flatten") on each block.


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] Maybe a silly question?

2011-03-23 Thread David Thompson

Hi Ken,

I have an Exodus II type file loaded into Paraview and trying to get  
the point ID for a group of cells. I know for vtk type object there  
are functions like GetCell(), GetPoint()


I assume you're talking about using Python from ParaView?

and just wonder if there are such kind of routines available for  
Exodus II files.



Yes, the ExodusII reader outputs VTK objects in a multiblock dataset.  
You can use GetCell(), GetPoint(), and such on each individual  
block... you just have to decide on a particular block. The ExodusII  
reader separates things into blocks because each one may have  
different point and/or cell variables. There are actually 2 levels of  
blocks; the top level segregates datasets into groups that correspond  
to exodus blocks vs. sets. The second level breaks individual blocks  
or sets into separate datasets. Does that help?


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] Unstructured Grid Descriptive Statistics

2010-12-07 Thread David Thompson

Hi Kristjan,

You are correct that the statistics filters do not weight values by  
cell area. It's torturous, but it might be possible to use the "Mesh  
Quality" filter to compute the cell areas/volumes, a calculator filter  
to compute the product of the area and a cell variable, and then run  
descriptive statistics on the resulting field. But being able to  
perform area/volume-weighted statistics (at least for cell data)  
sounds like a good use case to support. Would you be willing to file a  
feature request at http://paraview.org/Bug/ and attach your example to  
it?


Thanks,
David

On Dec 7, 2010, at 7:31 , Kristjan Onu wrote:


Hi,

I'm having difficulty using the Descriptive Statistics filter with  
data

specified on unstructured grids.

The two files attached to this email demonstrate the problem. The mean
of y obtained with 2_cell_square.vtu is 0.5 whereas with
5_cell_square.vtu the mean is 0.417.

I think this happens because the mean is calculated over all the  
points
without taking into account the cell area. Can anyone explain what  
steps

I should follow to obtain a mean of 0.5 with 5_cell_square.vtu?

Thanks,

Kristjan
<2_cell_square.vtu><5_cell_square.vtu>



___
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] Descriptive statistics filter

2010-10-13 Thread David Thompson

My second question concerns the keyword 'Cardinality' in descriptive
statistics When I make a slice of my mesh and apply the
statistics to it, for example, the cardinality values are far bigger
than the number of nodes of the sliced plane or curve. Is cardinality
always the total number of nodes of the grid?


The reason that you see a large number of points for the cardinality
of a slice is that slicing creates a new mesh whose nodes are not
always nodes of the original mesh... they are intersections of mesh
edges with the slice plane.


Slicing MAY also convert hexahedral cells into triangles before  
intersecting edges with the slice plane (I just don't recall for sure).


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] Descriptive statistics filter

2010-10-13 Thread David Thompson

Hi Martin,


 I have an ouput file that comes from a cfd simulation. After I load
the data in paraview and apply the descriptive statistics filter, I
see a table presenting mean value, min., max. ... of my quantities.
I'd like to ask if it's possible to export this table as a text file
so that I can plot the values of certain columns with gnuplot, for
example.


Yes, clicking on the table so it becomes the active view and then  
choosing the File->Save Data menu will let you export it as a CSV file.



 My second question concerns the keyword 'Cardinality' in descriptive
statistics. I'd expect this to relate somehow to the number of degrees
of freedom of my mesh, but the numbers in 'Cardinality' column do not
always correspond. When I make a slice of my mesh and apply the
statistics to it, for example, the cardinality values are far bigger
than the number of nodes of the sliced plane or curve. Is cardinality
always the total number of nodes of the grid?


Cardinality is the number of observations encountered. For point data,  
it will be the number of nodes in the mesh (PLUS an extra entry for  
each time a node is shared on more than one process if you have  
ParaView running in parallel). For cell data, it will be the number of  
cells in mesh.


The reason that you see a large number of points for the cardinality  
of a slice is that slicing creates a new mesh whose nodes are not  
always nodes of the original mesh... they are intersections of mesh  
edges with the slice plane. If you really want properties of a subset  
of mesh nodes, you can use the Elevation filter plus a Threshold to  
select a subset of cells or points and pass that to Descriptive  
Statistics. You might also be able to use the selection mechanism to  
select nodes of the original mesh with your view's camera set to  
parallel projection (instead of perspective projection). Once you have  
the nodes selected, use the Extract Selection filter to create a new  
dataset from the selected nodes.


Hope this helps,
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] Using Threshold

2010-07-06 Thread David Thompson
>>>  I definitely want to break it up into a point cloud, probably using 
>>> the glyph filter to create a vertex cell at each point. How do I go about 
>>> this?
>> 
>> Load your data. Go to the Filters menu and choose "Glyph". When the filter 
>> options appear in the Object Inspector panel, set the "Glyph Type" to "2D 
>> Glyph". Then, a "2D Glyph" box will appear. Select "Vertex" in its drop-down 
>> menu and click the filter's Apply button.
>> 
>> Note that if you have more than 5000 points, you'll need to turn off the 
>> "Mask Points" option or it will randomly choose 5000 points and only output 
>> those (this was done so that rendering of large datasets would be possible, 
>> but since you plan to subset the output of the Glyph filter with a Threshold 
>> filter there is no need to mask points).
> 
> ... Also, when I do regular 2D glyphs, I find that they are flat and only 
> properly visible when the camera is looking down on them from the 
> z-direction.  How do I ensure that the glyphs are always properly visible to 
> the camera?

Michael,

As far as I know, there is no way to do this in ParaView. You can use a vector 
field to orient the 2-D glyphs at each point, but their orientation is fixed 
relative to the dataset not the camera.

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] Using Threshold

2010-07-06 Thread David Thompson
>>> I have been using the Threshold filter for some time and now I am 
>>> discovering it doesn't do what I want it to do. I have a surface made up of 
>>> a set of points, each point has a location in 2-space, and a function value 
>>> attached to it, -2, -1, 0 or 1. I would like to filter out all points that 
>>> don't have a particular value, say zero. I plot the full shape, and place a 
>>> Threshold filter on top of it, excluding everything above (say) 0.2,a nd 
>>> everything below -0.2. If I tick the box that says "All scalars", none of 
>>> my points show up, even though I can see by the colour coding that there 
>>> are some points with the zero value.  If I untick the box then those points 
>>> do appear, but so do neighbouring points with other values.  What's going 
>>> on here?
>> 
>> The threshold filter copies cells, not points, from the input dataset to the 
>> output dataset. So, since you have point data it must use a rule to decide 
>> whether or not a given input cell should be included in the output based on 
>> the points that make up the cells. When you click on "All scalars," it uses 
>> the rule that *all* the points of a given input cell must be inside the 
>> threshold range (i.e., between -0.2 and 0.2) in order for the cell to be 
>> copied. It sounds like none of your cells have all 0 values at their points, 
>> so the output is empty. When "All scalars" is unchecked, the rule is that 
>> *any* point of a given input cell matching the threshold criterion will 
>> result in the input cell being copied to the output. In this case, cells 
>> with *any single* point value between -0.2 and 0.2 will cause the cell (and 
>> *all* of its points, even those which don't meet the threshold) being copied 
>> to the output. So, what is it that you really want to do? You might consider 
>> clipping by 
 the scalar field (that will get the shape right but break cells up into 
pieces) or converting the input dataset from cells that connect points to just 
a point cloud by using the glyph filter to create a VERTEX cell at each input 
point. You could then threshold the point cloud. In that case you'll get just 
the points that match as output, but they won't be connected to other points as 
they were in the original mesh.
> 
> I understand now.  I definitely want to break it up into a point cloud, 
> probably using the glyph filter to create a vertex cell at each point. How do 
> I go about this?

Hi Michael,

Load your data. Go to the Filters menu and choose "Glyph". When the filter 
options appear in the Object Inspector panel, set the "Glyph Type" to "2D 
Glyph". Then, a "2D Glyph" box will appear. Select "Vertex" in its drop-down 
menu and click the filter's Apply button.

Note that if you have more than 5000 points, you'll need to turn off the "Mask 
Points" option or it will randomly choose 5000 points and only output those 
(this was done so that rendering of large datasets would be possible, but since 
you plan to subset the output of the Glyph filter with a Threshold filter there 
is no need to mask points).

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] Using Threshold

2010-07-05 Thread David Thompson
> I have been using the Threshold filter for some time and now I am discovering 
> it doesn't do what I want it to do.
> I have a surface made up of a set of points, each point has a location in 
> 2-space, and a function value attached to it, -2, -1, 0 or 1.
> I would like to filter out all points that don't have a particular value, say 
> zero.
> I plot the full shape, and place a Threshold filter on top of it, excluding 
> everything above (say) 0.2,a nd everything below -0.2.
> If I tick the box that says "All scalars", none of my points show up, even 
> though I can see by the colour coding that there are some points with the 
> zero value.  If I untick the box then those points do appear, but so do 
> neighbouring points with other values.  What's going on here?

Hi Michael,

The threshold filter copies cells, not points, from the input dataset to the 
output dataset. So, since you have point data it must use a rule to decide 
whether or not a given input cell should be included in the output based on the 
points that make up the cells.

When you click on "All scalars," it uses the rule that *all* the points of a 
given input cell must be inside the threshold range (i.e., between -0.2 and 
0.2) in order for the cell to be copied. It sounds like none of your cells have 
all 0 values at their points, so the output is empty.

When "All scalars" is unchecked, the rule is that *any* point of a given input 
cell matching the threshold criterion will result in the input cell being 
copied to the output. In this case, cells with *any single* point value between 
-0.2 and 0.2 will cause the cell (and *all* of its points, even those which 
don't meet the threshold) being copied to the output.

So, what is it that you really want to do? You might consider clipping by the 
scalar field (that will get the shape right but break cells up into pieces) or 
converting the input dataset from cells that connect points to just a point 
cloud by using the glyph filter to create a VERTEX cell at each input point. 
You could then threshold the point cloud. In that case you'll get just the 
points that match as output, but they won't be connected to other points as 
they were in the original mesh.

Hope this helps,
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] Indication of delay between interactive and still render?

2010-06-25 Thread David Thompson
> All are good suggestions but won't it confuse people into thinking
> that paraview is busy doing something which may result in people
> waiting to interact until the overlay disappears which till have
> totally reverse effect :).

Well, then we'll have to come up with some subtle psychological twist to let 
people know. Perhaps the overlay could be of a prone person who, after 2 
seconds, stands up and performs calisthenics until the still render is 
complete. ;-)

Or more seriously, the overlay could be present during all interactive 
rendering and just start "counting down" (I still prefer a clock-like wedge to 
a countdown, especially given that it's a short 2-second interval) when the 
mouse is released. That way the overlay's presence indicates interactive 
rendering, a modified form would indicate the still rendering is being 
performed, and its absence indicates that the image is a completed still 
rendering.

David


> On Thu, Jun 24, 2010 at 10:44 PM, Sebastien Jourdain
>  wrote:
>> Me too but as a count down...
>> 
>> On Thu, Jun 24, 2010 at 8:38 PM, Marcus D. Hanwell
>>  wrote:
>>> On Thu, Jun 24, 2010 at 8:31 PM, David Thompson  wrote:
>>>>> 
>>>>> That is really a good point. The more I use it, the more I realize
>>>>> this is absolutely essential. But how ? :). What could serve as a good
>>>>> indication?
>>>> 
>>>> 1. Different colors on the border surrounding the pane (but this competes
>>>> with the red/blue indicating which pane is "Active")
>>>> 2. An LED-like glyph in the pane's title bar?
>>>> 3. A small overlay in one corner of the pane? (Ken's favorite: a pie chart
>>>> showing either time remaining before a still render or time spent waiting
>>>> for a still render to finish. :-)
>>>> 
>>> You beat me to it - I was just writing we could add a small overlay in a
>>> corner of the display to give the user a subtle visual cue. Option 3 gets my
>>> vote too ;-)
>>> Marcus
>>>>> 
>>>>> On Tue, Jun 22, 2010 at 12:00 PM, Moreland, Kenneth 
>>>>> wrote:
>>>>>> 
>>>>>> I just got a chance to try out the new feature that pauses for 2 seconds
>>>>>> between the last interactive render and the still render to give you a
>>>>>> chance to make further movements
>>>>>> (http://paraview.org/Bug/view.php?id=10232).  Totally cool.
>>>>>> 
>>>>>> One thing I noticed though is that there is no indication of an
>>>>>> intentional
>>>>>> pause.  Is there a subtle way to let the user know when the interface is
>>>>>> pausing and when it is busy rendering?  Or is it fine the way it is?
>>>>>> 
>>>>>> -Ken
>>>>>> 
>>>>>>     Kenneth Moreland
>>>>>>***  Sandia National Laboratories
>>>>>> ***
>>>>>> *** *** ***  email: kmo...@sandia.gov
>>>>>> **  ***  **  phone: (505) 844-8919
>>>>>>***  web:   http://www.cs.unm.edu/~kmorel
>>>>>> 
>>>>>> 
>>>>>> ___
>>>>>> Powered by www.kitware.com
>>>>>> 
>>>>>> Visit other Kitware open-source projects at
>>>>>> http://www.kitware.com/opensource/opensource.html
>>>>>> 
>>>>>> Please keep messages on-topic and check the ParaView Wiki at:
>>>>>> http://paraview.org/Wiki/ParaView
>>>>>> 
>>>>>> Follow this link to subscribe/unsubscribe:
>>>>>> http://www.paraview.org/mailman/listinfo/paraview
>>>>>> 
>>>>>> 
>>>>> ___
>>>>> 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] Indication of delay between interactive and still render?

2010-06-24 Thread David Thompson

That is really a good point. The more I use it, the more I realize
this is absolutely essential. But how ? :). What could serve as a good
indication?


1. Different colors on the border surrounding the pane (but this  
competes with the red/blue indicating which pane is "Active")

2. An LED-like glyph in the pane's title bar?
3. A small overlay in one corner of the pane? (Ken's favorite: a pie  
chart showing either time remaining before a still render or time  
spent waiting for a still render to finish. :-)


David

On Tue, Jun 22, 2010 at 12:00 PM, Moreland, Kenneth  
 wrote:
I just got a chance to try out the new feature that pauses for 2  
seconds
between the last interactive render and the still render to give  
you a

chance to make further movements
(http://paraview.org/Bug/view.php?id=10232).  Totally cool.

One thing I noticed though is that there is no indication of an  
intentional
pause.  Is there a subtle way to let the user know when the  
interface is

pausing and when it is busy rendering?  Or is it fine the way it is?

-Ken

     Kenneth Moreland
***  Sandia National Laboratories
***
*** *** ***  email: kmo...@sandia.gov
**  ***  **  phone: (505) 844-8919
***  web:   http://www.cs.unm.edu/~kmorel


___
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at:
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview



___
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] Fwd: plotting additional variables in paraview

2010-03-18 Thread David Thompson
Hi all,

I have committed a fix for this issue and verified that ParaView can
read VStress from brick_with_hole2-out.e . I've cleared the continuous
dashboard and will watch the nightlies. It should be in ParaView 3.8.

Berk, please let me know if/when the contract w/ EDF goes through and I
can have the reader tag array components with names.

David

> While I'm not deep enough into Paraview to contribute specifics to this
> discussion, I appreciate your continuing to pursue a resolution.
> 
> However, as a software developer of analysis codes, I always appreciated
> having names associated with various columns of output data to avoid any
> confusion by the end user (or application programmer, like me) about the
> order of data, especially tensor components.  Or, heaven forbid, having
> the user specify a column number...
> 
> Thanks again,
> John
> 
> 
> On 2/12/2010 9:02 AM, Thompson, David C wrote:
> > Berk,
> >
> > I'll see if perhaps I can bend the new glommer to behave
> > like the previous one until component names appear in
> > ParaView... a few months out sounds like post-3.8 and I
> > wouldn't want to have an official release that misbehaved.
> > But I'm really excited about being able to name components,
> > even though it seems like a small detail -- I guess it's a peeve
> > of mine. :-)
> >
> >  David
> > 
> > From: Berk Geveci [berk.gev...@kitware.com]
> > Sent: Friday, February 12, 2010 06:54
> > To: Thompson, David C; Utkarsh Ayachit
> > Cc: John Schiermeier
> > Subject: Re: Fwd: [Paraview] plotting additional variables in paraview
> >
> > Now for the question. One thing I have noticed is that with the box.ex2
> > file in ParaViewData, only part of the stress tensor (sigxx, sigyy,
> > sigzz, sigxy, sigyz, sigzx) is getting glommed because the last 3
> > entries are not in the "canonical" order for a symmetric tensor. I
> > haven't had time to see how ParaView currently handles this but I'd like
> > it to permute the variables into canonical order and glom them into one
> > "sig" variable. The problem is that analysts expecting the order to
> > match the simulation output order will be mislead. What would be nice is
> > if I could *name* each component of a vector/tensor vtkDataArray and
> > have that show up in the toolbar instead of the {Magnitude, 0, 1, 2, 3,
> > 4, 5} combo-box entries that currently appear. It would be easy to do
> > store the names in vtkInformation objects but I don't know how easy it
> > would be to make the ParaView GUI use them (or have VTK copy them down
> > the pipeline). Any ideas? If this is too hard, I will just match the old
> > reader behavior. Otherwise I'll improve the reader's handling of
> > glomming.
> >
> > Funny that you should mention this. We have a deliverable to EDF that is 
> > supposed to do exactly this. We haven't thought about how to implement it 
> > yet - in the GUI vs. in the pipeline. Do you think it is OK to have a 
> > temporary solution until we add this feature for EDF? I can't remember the 
> > deliverable date but it is a few months out.
> >
> > -berk
> >
> 



___
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] vtkStdString

2010-01-13 Thread David Thompson
> Could someone provide me an example 'to lowcase' a vtkStdString ?
> I would like to apply a tolower() method but don't know how.
> 
>   vtkStdString units;
>   if (units.find(" since ") != vtkStdString::npos)
>  {
>  ...
> 
Take a look at VTK/Utilities/kwsys/SystemTools.hxx.in and
RegularExpression.hxx. SystemTools::LowerCase is a routine to downcase
strings and RegularExpression lets you define searches. To use them, you
need to
#include 
#include 
Note that the includes use vtksys, not kwsys.

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] Serial remote views of distributed data

2009-10-30 Thread David Thompson
Can’t use just use the existing 3D view by internally shutting off  
distributed rendering for that view, which will cause the geometry  
to be delivered to the client?


I want neither distributed rendering nor geometry rendered at the  
client... I'd like geometry rendered on the rank 0 pvserver and image  
delivery to the client. Is that possible? (Well, OK, I know it's  
possible but how hard will it be?)


David



On 10/29/09 1:35 PM, "David Thompson"  wrote:

Hi all,

I'm looking to convert a (currently serial) VTK application into a
ParaView-style client (one that talks to a distributed set of pvserver
processes) and have a few questions about pqView, distributed data,
and remote display. My situation is that I have
1. large-distributed datasets as sources
2. pipelines which extract very small subsets of the data from #1
3. several vtkView subclasses which do 3-d rendering/interaction of
the subsets from #2 in a non-parallel way, and
4. clients which will be run on machines without decent 3-D rendering
hardware.

I would like to create new pqView classes that collect data from the
distributed pvservers to the rank 0 process, do the rendering on the
rank 0 pvserver process (using the serial vtkViews), and use
ParaView's remote image delivery mechanism to present the images to
the client.

Are there existing views that do this? Is the view the place where all
of this should be performed or are there other places where the data
aggregation should take place? I see that some of the OverView plugins
are client-side views but really what I'm looking for is rendering on
the rank 0 pvserver process, not the client.

Thanks,
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




     Kenneth Moreland
***  Sandia National Laboratories
***
*** *** ***  email: kmo...@sandia.gov
**  ***  **  phone: (505) 844-8919
***  web:   http://www.cs.unm.edu/~kmorel




___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


[Paraview] Serial remote views of distributed data

2009-10-29 Thread David Thompson

Hi all,

I'm looking to convert a (currently serial) VTK application into a  
ParaView-style client (one that talks to a distributed set of pvserver  
processes) and have a few questions about pqView, distributed data,  
and remote display. My situation is that I have

1. large-distributed datasets as sources
2. pipelines which extract very small subsets of the data from #1
3. several vtkView subclasses which do 3-d rendering/interaction of  
the subsets from #2 in a non-parallel way, and
4. clients which will be run on machines without decent 3-D rendering  
hardware.


I would like to create new pqView classes that collect data from the  
distributed pvservers to the rank 0 process, do the rendering on the  
rank 0 pvserver process (using the serial vtkViews), and use  
ParaView's remote image delivery mechanism to present the images to  
the client.


Are there existing views that do this? Is the view the place where all  
of this should be performed or are there other places where the data  
aggregation should take place? I see that some of the OverView plugins  
are client-side views but really what I'm looking for is rendering on  
the rank 0 pvserver process, not the client.


Thanks,
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 3.7 from cvs on Mac OS X, x86_64, Qt 4.5.3 w/Cocoa

2009-10-28 Thread David Thompson

I've noticed that the filters menu is greyed out after loading a valid
data set on the version of paraview I built yesterday from cvs (also
updated yesterday).  I have the same version of paraview built on an
i386 mac using the Qt SDK instead of the Cocoa 32/64  framework, and
it does not have this problem. Also, the filters appear to work (e.g.
if I hit the isosurface button the filter runs correctly).  Any idea
what is happening here? I believe somebody mentioned that the QtDBus
framework is messed up in the Qt Cocoa package, maybe I just need to
build Qt from source.


I've run into this problem as well. I don't have a fix, but do have a
workaround... pressing the Alt+Space key will let you type in the name
of a filter to create it.

Hope this helps,
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


[Paraview] ANNOUNCE: New statistics features in ParaView

2009-09-15 Thread David Thompson

Hi all,

This is a friendly neighborhood feature announcement: there's now a  
new ParaView plugin (in the Plugins/Statistics directory) that  
presents the statistics algorithms recently added to VTK as part of  
the ParaView UI. You can
- compute descriptive statistics (mean, variance, min, max, skewness,  
kurtosis),

- compute contingency tables,
- perform k-means analysis,
- examine correlations between arrays, and
- perform principal component analysis on arrays.

Some more details on the plugin are at
  http://www.paraview.org/Wiki/Statistical_analysis

Please let me know if you have any problems.

David

--
David Thompson, MS 9159(925)294-2537
Sandia National Laboratories   dcth...@sandia.gov
Visualization and Scientific Computing Bldg/Rm 915/S136
PO Box 969, Livermore, CA 94551-0969


___
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] greek letters

2009-09-11 Thread David Thompson

On Sep 9, 2009, at 3:04  PDT, Ricardo Reis wrote:

On Tue, 8 Sep 2009, David Thompson wrote:

1. Rendering of character glyphs other than 7-bit ASCII
2. Lay out and rendering of mathematical equations (and possibly  
other formatted text like tables)

3. Easy input of equations (and other text for fancy layout).


Can't an approch like the one on this plugin, used by Pidgin (a  
instant messenger program) to render LaTeX, work?

http://sourceforge.net/projects/pidgin-latex


It's one alternative. I don't think the 3 tasks I outlined exclude it.  
However it does require you to have LaTeX installed on your system.


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


[Paraview] ParaView on exotic platforms

2009-09-08 Thread David Thompson

Hi all,

I've been working to get ParaView (with Python turned ON) compiled on  
an HPC platform that doesn't support shared libraries. It almost  
works, but there is at least one little glitch that I don't see an  
easy fix for: during CMake configuration, ParaView creates a file named

  Utilities/VTKPythonWrapping/Executable/pvpythonmodules.h
which defines a routine that initializes all of the python modules  
ParaView requires. Unfortunately, this file is written by a macro  
defined in FindPythonLibs.cmake which uses a list of the modules and  
assumes that the function to initialize each one will be named

'init' + PREFIX + moduleName + SUFFIX
This doesn't work for the mpi4py module because the build system has  
overridden the name of the dynamically-loaded module name. Instead of  
libMPIPython.so, it is simply MPIPython.so.


So far, I've manually edited the pvpythonmodules file but it would be  
nice to have a real solution committed. Should I
1. Un-override the non-standard library name when building static  
Python modules
2. Change PYTHON_ADD_MODULE and PYTHON_WRITE_MODULES_HEADER to keep a  
list of target names (and (a) how can I do this when the target name  
might be changed later in CMake parsing? (b) will VTK's  
FindPythonLibs.cmake always be used or might the version distributed  
with CMake also need to change?)

3. Do something else?

Thanks,
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] greek letters

2009-09-08 Thread David Thompson

I would change that a bit. The ability to use Latex is not the main
issue here. Rather, it is the ability to use symbols and other unicode
characters as well as (possibly) the ability to typeset equations.
Whether this is entered as Latex or MathML or something else is more
of an implementation detail, isn't it?


Actually, I believe there are 3 separate issues, but they are all  
important: MathML is horribly tedious to type in and GUI equation  
editors are frustratingly slow if you know LaTeX so please don't  
discount human interaction. The issues are:


1. Rendering of character glyphs other than 7-bit ASCII
2. Lay out and rendering of mathematical equations (and possibly other  
formatted text like tables)

3. Easy input of equations (and other text for fancy layout).

Also, Berk didn't say it, but at least in VTK you can now use a Qt  
renderer that will draw any Unicode string so #1 above is done but not  
taken advantage of in ParaView. However, this does have the  
disadvantage of requiring Qt on the server side since I believe labels  
are drawn by all the processes running pvserver.


Just my 2 cents,
David


On Mon, Sep 7, 2009 at 9:37 AM, Ricardo Reis  
wrote:

On Mon, 7 Sep 2009, Berk Geveci wrote:


I'd suggest putting this in paraview.uservoice.com. I would be
interested to see how many people are interested in this feature.


Added:

http://paraview.uservoice.com/pages/11350-general/suggestions/309511-add-ability-to-use-latex-in-text-and-axis-text?ref=title

please vote ;)

 Ricardo Reis

 'Non Serviam'

 PhD candidate @ Lasef
 Computational Fluid Dynamics, High Performance Computing, Turbulence
 http://www.lasef.ist.utl.pt

 Cultural Instigator @ Rádio Zero
 http://www.radiozero.pt

 Keep them Flying! Ajude/help a Aero Fénix!

 http://www.aeronauta.com/aero.fenix

 http://www.flickr.com/photos/rreis/

___
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] LSDYNA support in 3.6.1

2009-08-13 Thread David Thompson

I have another program that has supposedly generated an LSDYNA file
named *.k. ParaView 3.6.1 does recognize the file type as something it
thinks it can open but when I select the file and click "Open" from
the file dialog I get another window that asks me to specify the
reader to use and LSDYNA isn't in the list. Thoughts?
The *.k files are input decks. The LSDyna reader in VTK really only  
reads output decks (named d3plot, d3plot1, ...). If those are present,  
it will read the input deck to obtain part names but does not read the  
input deck for geometry or boundary condition information.


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] pvserver maxed out

2009-06-16 Thread David Thompson
> When I try that (on my MacBook Pro w/ 2 procs):
> mpirun --mca mpi_yield_when_idle 1 -np 2 pvserver
> 
> I still have one CPU pegged at 98%.  Are you saying that shouldn't be  
> the case?

It was my understanding that was what mpi_yield_when_idle was there to
fix -- since oversubscribed slots will otherwise have multiple processes
in busy-wait loops -- but I haven't fiddled with it in a while.

David

> On Jun 16, 2009, at 3:47 PM, David Thompson wrote:
> 
> >
> >> I confirm that with mpich2 version 1.0.7, using the default device
> >> (ch3:sock), the load of an idle pvserver is near 0. With ch3:nemesis,
> >> all but one of the processes use a full cpu. With OpenMPI, I think  
> >> all
> >> processes use a full cpu.
> >
> > This is configurable on the mpirun command line:
> >
> > http://www.open-mpi.org/faq/?category=running#force-aggressive- 
> > degraded
> >
> > If you set the MCA parameter mpi_yield_when_idle to a non-zero value,
> > OpenMPI will not aggressively poll inside MPI routines. For example:
> >
> >  mpirun --mca mpi_yield_when_idle 1 -np 4 pvserver
> >
> > David
> >
> >> In conclusion, if you are running pvserver on a Mac laptop, use  
> >> mpich2
> >> with default build options.
> >>
> >> -berk
> >>
> >> On Mon, Jun 15, 2009 at 9:11 PM, Berk  
> >> Geveci wrote:
> >>> For this reason, I use MPICH2 on my MacBook Pro and OpenMPI on my  
> >>> Mac
> >>> Pro. I can't remember which device I ended up settling on - I think
> >>> the default. If I remember correctly, ch3:nemesis was causing the  
> >>> same
> >>> issue as OpenMPI.
> >>>
> >>> -berk
> >>>
> >>> On Mon, Jun 15, 2009 at 1:44 PM, Randy  
> >>> Heiland wrote:
> >>>> Thanks Utkarsh.  Just to add to this thread, I asked some of the  
> >>>> OMPI
> >>>> developers and apparently it's not possible to avoid this behavior.
> >>>> -Randy
> >>>> On Jun 15, 2009, at 11:35 AM, Utkarsh Ayachit wrote:
> >>>>
> >>>> Yup, that's expected. That's because of your MPI implementation.  
> >>>> OpenMPI
> >>>> does a busy wait when waiting for mpi messages. I think there's an
> >>>> environment variable or something that you can set to disable  
> >>>> that, I am no
> >>>> sure. Maybe someone on the mailing list knows, or try checking  
> >>>> the OpenMPI
> >>>> documentation.
> >>>>
> >>>> Utkarsh
> >>>>
> >>>> On Mon, Jun 15, 2009 at 11:28 AM, Randy Heiland  >>>> > wrote:
> >>>>>
> >>>>> When I use 'mpirun -np 2 pvserver'  on my dual core MacBook  
> >>>>> (either
> >>>>> connecting to it via paraview client or standalone), it keeps  
> >>>>> one of my CPUs
> >>>>> pegged at nearly 100%.  Is this to be expected?  (Just running  
> >>>>> the serial
> >>>>> './pvserver' of course just generates 'Listening...Waiting for  
> >>>>> client...'
> >>>>> with no CPU load).
> >>>>>
> >>>>> % ./pvserver --version
> >>>>> ParaView3.7
> >>>>>
> >>>>> % ompi_info
> >>>>>   Open MPI: 1.3.1
> >>>>>
> >>>>> thanks, Randy
> >>>>> ___
> >>>>> 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] pvserver maxed out

2009-06-16 Thread David Thompson

> I confirm that with mpich2 version 1.0.7, using the default device
> (ch3:sock), the load of an idle pvserver is near 0. With ch3:nemesis,
> all but one of the processes use a full cpu. With OpenMPI, I think all
> processes use a full cpu.

This is configurable on the mpirun command line:

http://www.open-mpi.org/faq/?category=running#force-aggressive-degraded

If you set the MCA parameter mpi_yield_when_idle to a non-zero value,
OpenMPI will not aggressively poll inside MPI routines. For example:

  mpirun --mca mpi_yield_when_idle 1 -np 4 pvserver

David

> In conclusion, if you are running pvserver on a Mac laptop, use mpich2
> with default build options.
> 
> -berk
> 
> On Mon, Jun 15, 2009 at 9:11 PM, Berk Geveci wrote:
> > For this reason, I use MPICH2 on my MacBook Pro and OpenMPI on my Mac
> > Pro. I can't remember which device I ended up settling on - I think
> > the default. If I remember correctly, ch3:nemesis was causing the same
> > issue as OpenMPI.
> >
> > -berk
> >
> > On Mon, Jun 15, 2009 at 1:44 PM, Randy Heiland wrote:
> >> Thanks Utkarsh.  Just to add to this thread, I asked some of the OMPI
> >> developers and apparently it's not possible to avoid this behavior.
> >> -Randy
> >> On Jun 15, 2009, at 11:35 AM, Utkarsh Ayachit wrote:
> >>
> >> Yup, that's expected. That's because of your MPI implementation. OpenMPI
> >> does a busy wait when waiting for mpi messages. I think there's an
> >> environment variable or something that you can set to disable that, I am no
> >> sure. Maybe someone on the mailing list knows, or try checking the OpenMPI
> >> documentation.
> >>
> >> Utkarsh
> >>
> >> On Mon, Jun 15, 2009 at 11:28 AM, Randy Heiland  
> >> wrote:
> >>>
> >>> When I use 'mpirun -np 2 pvserver'  on my dual core MacBook (either
> >>> connecting to it via paraview client or standalone), it keeps one of my 
> >>> CPUs
> >>> pegged at nearly 100%.  Is this to be expected?  (Just running the serial
> >>> './pvserver' of course just generates 'Listening...Waiting for client...'
> >>> with no CPU load).
> >>>
> >>> % ./pvserver --version
> >>> ParaView3.7
> >>>
> >>> % ompi_info
> >>>Open MPI: 1.3.1
> >>>
> >>> thanks, Randy
> >>> ___
> >>> Powered by www.kitware.com
> >>>
> >>> Visit other Kitware open-source projects at
> >>> http://www.kitware.com/opensource/opensource.html
> >>>
> >>> Please keep messages on-topic and check the ParaView Wiki at:
> >>> http://paraview.org/Wiki/ParaView
> >>>
> >>> Follow this link to subscribe/unsubscribe:
> >>> http://www.paraview.org/mailman/listinfo/paraview
> >>
> >>
> >>
> >> ___
> >> Powered by www.kitware.com
> >>
> >> Visit other Kitware open-source projects at
> >> http://www.kitware.com/opensource/opensource.html
> >>
> >> Please keep messages on-topic and check the ParaView Wiki at:
> >> http://paraview.org/Wiki/ParaView
> >>
> >> Follow this link to subscribe/unsubscribe:
> >> http://www.paraview.org/mailman/listinfo/paraview
> >>
> >>
> >
> ___
> Powered by www.kitware.com
> 
> Visit other Kitware open-source projects at 
> http://www.kitware.com/opensource/opensource.html
> 
> Please keep messages on-topic and check the ParaView Wiki at: 
> http://paraview.org/Wiki/ParaView
> 
> Follow this link to subscribe/unsubscribe:
> http://www.paraview.org/mailman/listinfo/paraview
> 


___
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] error in X86_64: libvtkalglib.so

2008-12-16 Thread David Thompson
Ricardo,

As Berk mentions, install rules have been added. Can you say exactly
when your had a checkout that was working and when you updated and it
stopped working? I added vtkalglib to VTK 2008-12-04 15:57:36 GMT and
added the install rule 2008-12-15 23:12:04 GMT. Depending on when things
stopped working, I may have more to fix.

Thanks,
David


On Tue, 2008-12-16 at 07:11 -0700, Berk Geveci wrote:
> Ricardo,
> 
> Update your cvs and try again please. This is a new library that was
> missing install rules. It should be fixed.
> 
> -berk
> 
> 
> 
> On Tue, Dec 16, 2008 at 8:48 AM, Ricardo Reis  wrote:
> > On Tue, 16 Dec 2008, Olesen, Mark wrote:
> >
> >> Since bin/paraview finds paraview-real, it would seem that the
> >> lib/paraview-X.X directory is actually there.
> >>
> >> Did you attempt to install to the same directory as your build
> >> directory?  This is where I've encounter the same problem before.
> >
> > No, I'm using a diferent dir and when I do "make install" he puts stuff
> > there.
> >
> >> If you just edit the CMAKE_INSTALL_PREFIX:PATH in CMakeCache.txt (in the
> >> build dir) and try 'make install' again.
> >
> > I have this in my CMAKE_INSTALL_PREFIX
> >
> > CMAKE_INSTALL_PREFIX:PATH=/usr/local/paraviewcvs.shared
> >
> > so this is not the problem, it seems...
> >
> > Should I atach my CMakeCache.txt ?
> >
> >  Ricardo Reis
> >
> >  'Non Serviam'
> >
> >  PhD student @ Lasef
> >  Computational Fluid Dynamics, High Performance Computing, Turbulence
> >  http://www.lasef.ist.utl.pt
> >
> >  &
> >
> >  Cultural Instigator @ Rádio Zero
> >  http://www.radiozero.pt
> >
> >  http://www.flickr.com/photos/rreis/
> > ___
> > ParaView mailing list
> > ParaView@paraview.org
> > http://www.paraview.org/mailman/listinfo/paraview
> >
> >
> ___
> ParaView mailing list
> ParaView@paraview.org
> http://www.paraview.org/mailman/listinfo/paraview
> 


___
ParaView mailing list
ParaView@paraview.org
http://www.paraview.org/mailman/listinfo/paraview


Re: [Paraview] Paraview/Kitware site down

2008-12-15 Thread David Thompson
> I've been thinking about this for quite some time, but the website(s)
> being down prompted me to finally post it.
> 
> I know that Paraview and VTK have been using cvs for a very long time,
> but perhaps it might be time to consider migrating to a more modern
> source code management (SCM) system - for example, git.
This came up a while back with respect to CMake (another Kitware
project). You can see the thread here, along with some of Kitware's
concerns.
http://www.mail-archive.com/cm...@cmake.org/msg10880.html

David


___
ParaView mailing list
ParaView@paraview.org
http://www.paraview.org/mailman/listinfo/paraview