Re: [Paraview] Newbie question on Python Programmable Filters

[Andy Bauer]

It's even simpler than setting CCFLAGS to turn debug information on. Just
set CMAKE_BUILD_TYPE to "Debug". It may already be set to that.

If you're creating your VTK files from Python code, you may want to look at
PyEVTK (https://bitbucket.org/pauloh/pyevtk). I haven't tried the
unstructured grid/vtu writer from that but have tried the image data writer
from it.

Andy

On Thu, Sep 6, 2012 at 1:19 PM, Frank Horowitz
wrote:

> Thanks for the replies!
>
> At the moment, I need to finish a poster for a presentation on Saturday,
> but promise to get to this Sunday or Monday. (I'll need to re-build a
> debugging version of paraview. Any hints on how to accomplish this within
> ccmake/cmake will be gratefully received, since that system is new to me
> too. Is it as simple as setting a CCFLAGS environment variable before
> invoking make -- or something similar??)
>
> I also found that the 200**3 wavelet worked fine on my
> paraview-compiled-from-source Linux installation. Let's ignore the problem
> I found on the OSX box with the binary installer for the moment. (One
> problem at a time!)
>
> Because of that, I'm now suspicious of the upstream code that writes .vtu
> files.
>
> Also, responding here to your later comment from your next email:
>
> > I am going to disagree with Andy about this being a better solution. I
> > don't even think that it addresses what Frank is trying to do. But beyond
> > that Array Calculator is cumbersome and slow compared to the Python
> > Calculator or Python Programmable Filter.
> >
> > As far as I understand, Frank is not trying to combine components from
> > different files. He is trying to combine components written as different
> > arrays in the same file hence the need for sqrt(x*x + y*y + z*z). I think
> > that Andy is confused by inputs[0] - that's simple the way to access the
> > input (even if there is only one) of the Python Programmable Filter.
>
> Actually, there is one file per component per timeslice.
>
> I've had a quick look at the file format written by my upstream simulator.
> It's very simple.  Maybe I should just write the .vtu files from within my
> own script that reads the upstream meshes to avoid any other potential
> problems with the upstream codes.
>
> Does paraview.simple have an easy way to construct a field vector at the
> location of arbitrary points? If so, I think the Python script to do *that*
> will be faster than me trying to chase all of this down…
>
> Thanks again for the advice!
>
> Cheers,
> Frank Horowitz
>
>
>
> On 06/09/2012, at 1:00 PM, Berk Geveci wrote:
>
> > I just ran this with a 500^3 wavelet on my machine without any problems.
> > The whole thing (including a slice filter) takes 20-30 seconds but the
> > programmable filter is actually pretty much instantaneous (<0.1 seconds)
> > using a debug build from Git master. Is there any way you can attach to
> the
> > stuck ParaView with a debugger and look at the stack trace?
> >
> > On Wed, Sep 5, 2012 at 9:13 PM, Frank Horowitz
> > wrote:
> >
> >> Still more checking things out. A problem (THE problem???) can be
> >> demonstrated with your code snippet from a few messages back if you
> >> construct the Wavelet like so rather than using the default constructor:
> >>
> >> Wavelet(WholeExtent=(-100,100,-100,100,-100,100))
> >>
> >> N.B. this is on a different machine, with a stock Mac OSX  binary
> version
> >> 3.14.1 installed.
> >>
> >>
> >> On 05/09/2012, at 5:49 PM, Frank Horowitz wrote:
> >>
> >>> OK, more tracking the problem:
> >>>
> >>> The assumption underlying my "close up" description below is not
> >> correct.  The problem appears to be the actual _size_ of the dataset I'm
> >> feeding in. Even a single time step fed in to the filter (with some
> 866,000
> >> points) hangs at the x*x computation, even though that appeared to work
> >> perfectly fine with your 3D Wavelet as a source (only some 9,000
> points).
> >>>
> >>> I'm guessing a memory leak somewhere?  When was the last time paraview
> >> was run through valgrind?
> >>>
> >>> Cheers (and sorry about the red herring earlier),
> >>>  Frank
> >>>
> >>>
> >>>
> >>> On 05/09/2012, at 4:42 PM, Frank Horowitz wrote:
> >>>
>  OK, the 30,000 foot view:
> 
>  I have a time series of about 20 slices (currently) of values of each
> >> component (3 components total) of a vector displacement 3d spatial field
> >> from an elastic wave propagation simulation. At the moment, I'm simply
> >> trying to visualise the vector length to see the propagating phases of
> the
> >> wavefield. Because my upstream codes and scripts only write individual
> >> components of the vector field, I need to combine the components and
> then
> >> calculate the vector lengths *somewhere*.  The programmable filter
> seemed
> >> like a good place to try that, since once I am past the learning curve
> for
> >> it, I'll have the full flexibilty of numpy for doing almost anything
> else I
> >> want with these and similar data from my upstream

Re: [Paraview] Manipulating Big Data through Python Shell

Today I uploaded two videos to Youtube, showing the basics of what is 
happening. They can be found via the following links:

1.   http://www.youtube.com/watch?v=68yeT8t0TUc&feature=youtu.be

2.   http://www.youtube.com/watch?v=Linns_iSBro&feature=youtu.be


The init.py file used now uses "renView.InteractiveRender()" instead of just 
"Render()", as suggested by Utkarsh. (Thank you again for that suggestion by 
the way!)

The first video is the mouse-based manipulation of 20MB of earthquake data. As 
shown, LOD is not active and my computer has no problem panning, rotating, etc. 
in near real-time.

Now in the second video you see how ParaView acts when using the Wii Nunchuk. 
The first part, without LOD active, you can see the extreme delay after I let 
go. After turning on LOD, the rendering speed does improve quite a bit but it 
is still not nearly as instantaneous as the mouse is without LOD. Also the data 
doesn't fully re-render after releasing the joystick; a mouse click is still 
required. I think I can remedy that just by sending "Render()" when the 
joystick is released.

So my current goal is to get the Wii Nunchuk application running the ParaView 
rendering at least at 90% of the capacity of what the mouse can, only using LOD 
if it is required when using the mouse also. I'm clueless as what else I can do 
with sending commands to the Python Shell to speed things up.  Are there any 
ideas of how to reach this goal from where I am now? Whether it be improving 
the current process of sending Python commands or using a completely different 
process of manipulating the rendering from an external application.

Also here is some more insight as to what the Wii application is actually doing 
in terms of interacting with ParaView. Just in case it could be useful 
information:


1.   The application is written in C++. During initialization, it connects 
to the server socket and loads the init.py file through the following commands:

svr.sin_family = AF_INET;
svr.sin_port = htons(9000);
inet_aton("127.0.0.1", &svr.sin_addr);
s = socket(AF_INET, SOCK_STREAM, 0);
connect(s, (const sockaddr*)&svr, sizeof(sockaddr_in));
write(s, "execfile('init.py')", sizeof("execfile('init.py')"));


2.   From there, every time the joystick is moved outside of its "deadzone" 
the following lines are ran:

stringstream ss;
ss.str(std::string());
ss << "command(" << x_value << "," << y_value << ");";
const std::string temp = ss.str();
const char* cstr = temp.c_str();
write(s, cstr, strlen(cstr));


3.   I'm pretty sure that it is not the C++ side of my Wii application 
slowing things down. Using 'printf()'s, all the 'write()'s happen and finish 
without any noticeable delay.

I'm still open to any suggestions out there! All help is very appreciated :)

Thanks,

Travis J. Bueter
--
Missouri University of Science and Technology - Junior
B.S. Computer Engineering/Computer Science
tjb...@mail.mst.edu
(573)-238-5843

IT RSS
Treasurer - MS&T Robotics Competition Team
--


___
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] socket communicating error

[LUO Jian]

Hi All,

I tried to use a binary 3.14.1 version paraview client on windows 7 to
connect to a remote server which is a compiled 3.14.1 version in
linux.
After the "Client connected." message, the following error will occur
on the remote server:

ERROR: In 
/home/luoj/ParaView-3.14.1-Source/VTK/Parallel/vtkSocketCommunicator.cxx,
line 812
vtkSocketCommunicator (0x1347e8a0): Could not receive tag. 1

And the client on the window side will also stop working.

So how can I solve this problem?

Best,
Jian
___
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 on a Retina display

[Pierre-Olivier Dallaire]

Robert,

I have looked into this and the QT team are making adjustments to support HiDPI:
https://bugreports.qt-project.org/browse/QTBUG-23870

QT 5 beta is not "HiDPI" ready yet...

I don't know how difficult it would be after to update Paraview to support the 
same capabilities.
The temporary workaround affects only the fonts; icons and main window remain 
blurry.

Thanks

PO


On 2012-09-05, at 3:11 PM, Pierre-Olivier Dallaire 
 wrote:

> I'm compiling QT 5 beta to see what happens, I'll see if I can get ParaView 
> 3.14 compiled against it.
> 
> On 2012-09-05, at 1:27 PM, Robert Maynard  wrote:
> 
>> I do believe that is correct. I don't know what version of Qt will have 
>> retina support by default.
>> 
>> On Wed, Sep 5, 2012 at 10:26 AM, Pierre-Olivier Dallaire 
>>  wrote:
>> Thanks / yes, that takes care of rendering the fonts correctly, icons are 
>> still blurry but overall it's better.
>> My guess is that to get the full retina support, that means QT needs to be 
>> updated first … ?
>> 
>> 
>> On 2012-09-05, at 10:15 AM, Robert Maynard  
>> wrote:
>> 
>>> As a temporary workaround you can manually modify the paraview application 
>>> to enable retina support.
>>> I haven't verified that all the icons layout properly.
>>> 
>>> You need to edit the info.plist inside the Paraview.app/Contents/ folder 
>>> and add the following two lines as the
>>> first items in the :
>>> 
>>> NSPrincipalClass
>>> NSApplication
>>> 
>>> 
>>> On Wed, Sep 5, 2012 at 12:50 AM, Pierre-Olivier Dallaire 
>>>  wrote:
>>> Hi,
>>> 
>>> I'm testing ParaView on a new MacBook Pro Retina and was wondering if there 
>>> is
>>> any plan to support that hardware in a near future ? The application works 
>>> great
>>> but the GUI is not as "sharp" as it could be in on this display.
>>> 
>>> Regards,
>>> 
>>> PO
>>> ___
>>> 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
>>> 
>>> 
>>> 
>>> -- 
>>> Robert Maynard
>> 
>> 
>> 
>> 
>> -- 
>> Robert Maynard
> 
> ___
> 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] Newbie question on Python Programmable Filters

[Berk Geveci]

Just to make sure that I understand this: you have some code that reads
some files in a proprietary format and spits out 1 vtu file per component?
What is that format like? What is the script written in? It should be
pretty easy to change that script to combine components together and write
them out as a vector. You could also calculate derived quantities there and
write them to the vtu file. Or you could still calculate that using the
Programmable Filter.

On Thu, Sep 6, 2012 at 1:19 PM, Frank Horowitz
wrote:

> Thanks for the replies!
>
> At the moment, I need to finish a poster for a presentation on Saturday,
> but promise to get to this Sunday or Monday. (I'll need to re-build a
> debugging version of paraview. Any hints on how to accomplish this within
> ccmake/cmake will be gratefully received, since that system is new to me
> too. Is it as simple as setting a CCFLAGS environment variable before
> invoking make -- or something similar??)
>
> I also found that the 200**3 wavelet worked fine on my
> paraview-compiled-from-source Linux installation. Let's ignore the problem
> I found on the OSX box with the binary installer for the moment. (One
> problem at a time!)
>
> Because of that, I'm now suspicious of the upstream code that writes .vtu
> files.
>
> Also, responding here to your later comment from your next email:
>
> > I am going to disagree with Andy about this being a better solution. I
> > don't even think that it addresses what Frank is trying to do. But beyond
> > that Array Calculator is cumbersome and slow compared to the Python
> > Calculator or Python Programmable Filter.
> >
> > As far as I understand, Frank is not trying to combine components from
> > different files. He is trying to combine components written as different
> > arrays in the same file hence the need for sqrt(x*x + y*y + z*z). I think
> > that Andy is confused by inputs[0] - that's simple the way to access the
> > input (even if there is only one) of the Python Programmable Filter.
>
> Actually, there is one file per component per timeslice.
>
> I've had a quick look at the file format written by my upstream simulator.
> It's very simple.  Maybe I should just write the .vtu files from within my
> own script that reads the upstream meshes to avoid any other potential
> problems with the upstream codes.
>
> Does paraview.simple have an easy way to construct a field vector at the
> location of arbitrary points? If so, I think the Python script to do *that*
> will be faster than me trying to chase all of this down…
>
> Thanks again for the advice!
>
> Cheers,
> Frank Horowitz
>
>
>
> On 06/09/2012, at 1:00 PM, Berk Geveci wrote:
>
> > I just ran this with a 500^3 wavelet on my machine without any problems.
> > The whole thing (including a slice filter) takes 20-30 seconds but the
> > programmable filter is actually pretty much instantaneous (<0.1 seconds)
> > using a debug build from Git master. Is there any way you can attach to
> the
> > stuck ParaView with a debugger and look at the stack trace?
> >
> > On Wed, Sep 5, 2012 at 9:13 PM, Frank Horowitz
> > wrote:
> >
> >> Still more checking things out. A problem (THE problem???) can be
> >> demonstrated with your code snippet from a few messages back if you
> >> construct the Wavelet like so rather than using the default constructor:
> >>
> >> Wavelet(WholeExtent=(-100,100,-100,100,-100,100))
> >>
> >> N.B. this is on a different machine, with a stock Mac OSX  binary
> version
> >> 3.14.1 installed.
> >>
> >>
> >> On 05/09/2012, at 5:49 PM, Frank Horowitz wrote:
> >>
> >>> OK, more tracking the problem:
> >>>
> >>> The assumption underlying my "close up" description below is not
> >> correct.  The problem appears to be the actual _size_ of the dataset I'm
> >> feeding in. Even a single time step fed in to the filter (with some
> 866,000
> >> points) hangs at the x*x computation, even though that appeared to work
> >> perfectly fine with your 3D Wavelet as a source (only some 9,000
> points).
> >>>
> >>> I'm guessing a memory leak somewhere?  When was the last time paraview
> >> was run through valgrind?
> >>>
> >>> Cheers (and sorry about the red herring earlier),
> >>>  Frank
> >>>
> >>>
> >>>
> >>> On 05/09/2012, at 4:42 PM, Frank Horowitz wrote:
> >>>
>  OK, the 30,000 foot view:
> 
>  I have a time series of about 20 slices (currently) of values of each
> >> component (3 components total) of a vector displacement 3d spatial field
> >> from an elastic wave propagation simulation. At the moment, I'm simply
> >> trying to visualise the vector length to see the propagating phases of
> the
> >> wavefield. Because my upstream codes and scripts only write individual
> >> components of the vector field, I need to combine the components and
> then
> >> calculate the vector lengths *somewhere*.  The programmable filter
> seemed
> >> like a good place to try that, since once I am past the learning curve
> for
> >> it, I'll have the full flexibilty of numpy 

Re: [Paraview] Getting String arrays out of Field Data

[Joshua Murphy]

Thank you! I was using GetArray(), as I had no idea that
GetAbstractArray() even existed! (tab complete in the Python Shell doesn't
recognize it, and that is how I find most things as I learn how various
things work).

Again, Thank you very much! This solves a big problems for us!

-Josh
 

On 9/6/12 11:27 AM, "David Thompson"  wrote:

>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=H
>EAD
>
>
>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

[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] Newbie question on Python Programmable Filters

[Frank Horowitz]

Thanks for the replies!

At the moment, I need to finish a poster for a presentation on Saturday, but 
promise to get to this Sunday or Monday. (I'll need to re-build a debugging 
version of paraview. Any hints on how to accomplish this within ccmake/cmake 
will be gratefully received, since that system is new to me too. Is it as 
simple as setting a CCFLAGS environment variable before invoking make -- or 
something similar??)

I also found that the 200**3 wavelet worked fine on my 
paraview-compiled-from-source Linux installation. Let's ignore the problem I 
found on the OSX box with the binary installer for the moment. (One problem at 
a time!)

Because of that, I'm now suspicious of the upstream code that writes .vtu files.

Also, responding here to your later comment from your next email:

> I am going to disagree with Andy about this being a better solution. I
> don't even think that it addresses what Frank is trying to do. But beyond
> that Array Calculator is cumbersome and slow compared to the Python
> Calculator or Python Programmable Filter.
> 
> As far as I understand, Frank is not trying to combine components from
> different files. He is trying to combine components written as different
> arrays in the same file hence the need for sqrt(x*x + y*y + z*z). I think
> that Andy is confused by inputs[0] - that's simple the way to access the
> input (even if there is only one) of the Python Programmable Filter.

Actually, there is one file per component per timeslice.

I've had a quick look at the file format written by my upstream simulator. It's 
very simple.  Maybe I should just write the .vtu files from within my own 
script that reads the upstream meshes to avoid any other potential problems 
with the upstream codes. 

Does paraview.simple have an easy way to construct a field vector at the 
location of arbitrary points? If so, I think the Python script to do *that* 
will be faster than me trying to chase all of this down…

Thanks again for the advice!

Cheers,
Frank Horowitz


 
On 06/09/2012, at 1:00 PM, Berk Geveci wrote:

> I just ran this with a 500^3 wavelet on my machine without any problems.
> The whole thing (including a slice filter) takes 20-30 seconds but the
> programmable filter is actually pretty much instantaneous (<0.1 seconds)
> using a debug build from Git master. Is there any way you can attach to the
> stuck ParaView with a debugger and look at the stack trace?
> 
> On Wed, Sep 5, 2012 at 9:13 PM, Frank Horowitz
> wrote:
> 
>> Still more checking things out. A problem (THE problem???) can be
>> demonstrated with your code snippet from a few messages back if you
>> construct the Wavelet like so rather than using the default constructor:
>> 
>> Wavelet(WholeExtent=(-100,100,-100,100,-100,100))
>> 
>> N.B. this is on a different machine, with a stock Mac OSX  binary version
>> 3.14.1 installed.
>> 
>> 
>> On 05/09/2012, at 5:49 PM, Frank Horowitz wrote:
>> 
>>> OK, more tracking the problem:
>>> 
>>> The assumption underlying my "close up" description below is not
>> correct.  The problem appears to be the actual _size_ of the dataset I'm
>> feeding in. Even a single time step fed in to the filter (with some 866,000
>> points) hangs at the x*x computation, even though that appeared to work
>> perfectly fine with your 3D Wavelet as a source (only some 9,000 points).
>>> 
>>> I'm guessing a memory leak somewhere?  When was the last time paraview
>> was run through valgrind?
>>> 
>>> Cheers (and sorry about the red herring earlier),
>>>  Frank
>>> 
>>> 
>>> 
>>> On 05/09/2012, at 4:42 PM, Frank Horowitz wrote:
>>> 
 OK, the 30,000 foot view:
 
 I have a time series of about 20 slices (currently) of values of each
>> component (3 components total) of a vector displacement 3d spatial field
>> from an elastic wave propagation simulation. At the moment, I'm simply
>> trying to visualise the vector length to see the propagating phases of the
>> wavefield. Because my upstream codes and scripts only write individual
>> components of the vector field, I need to combine the components and then
>> calculate the vector lengths *somewhere*.  The programmable filter seemed
>> like a good place to try that, since once I am past the learning curve for
>> it, I'll have the full flexibilty of numpy for doing almost anything else I
>> want with these and similar data from my upstream simulation codes.
 
 Now the 10,000 foot view:
 I've been trying to simply calculate the sqrt(x*x + y*y + z*z) scalar
>> field (where x, y, and z are the respective components as a function of
>> position)  and let that flow on down the pipeline as a new spatial field to
>> be visualised. I assumed (correctly if I'm interpreting your answer
>> correctly) that if I built the filter to deal with the spatial aspects of
>> the calculation, the time series would be taken care of by paraview
>> iterating over all files pointed to by inputs[0] (a set of files containing
>> timesteps of the 

Re: [Paraview] Newbie question on Python Programmable Filters

[Berk Geveci]

I am going to disagree with Andy about this being a better solution. I
don't even think that it addresses what Frank is trying to do. But beyond
that Array Calculator is cumbersome and slow compared to the Python
Calculator or Python Programmable Filter.

As far as I understand, Frank is not trying to combine components from
different files. He is trying to combine components written as different
arrays in the same file hence the need for sqrt(x*x + y*y + z*z). I think
that Andy is confused by inputs[0] - that's simple the way to access the
input (even if there is only one) of the Python Programmable Filter.

On Thu, Sep 6, 2012 at 11:56 AM, Andy Bauer  wrote:

> Hmm, I tried with the WholeExtent of (-100, 100, -100, 100, -100, 100) on
> my machine and still didn't have problems. How much memory is on your
> system? I'd hate to think that ParaView is that poor with memory
> management. If I remember correctly, we used to have nightly dashboards
> running valgrind. I know that Bill Lorensen's project last fall was going
> through with valgrind to get rid of a bunch of memory issues.
>
> Well, the good news is that I think I have a better solution for you
> anyways. I get the feeling from your original script that you have a
> separate file for each time step and each component of your vector field.
> If all components of your vector field are in the same file you can ignore
> the first part of the steps below. Here's the solution:
>
> 1) load all of the files for the first time step (e.g. load the
> ..
> 2) use control and select all the file readers
> 3) use the append attributes filter to combine the fields from each input
> data set into a single output dataset (it assumes that the input data sets
> all have the same number of points and cells and the desired fields all
> have different names -- use the calculator filter to change the field names
> if needed)
>  file>
> 4) create a calculator filter
> 4a) compute the magnitude by inputting in "mag( name>*iHat + *jHat +  name>*kHat)" for the widget right above the calculator buttons
> 4b) if you want a specific name for the newly created field, set it in
> Result Array Name
>
> Let me know if this accomplishes your goal.
>
> Andy
>
>
> On Wed, Sep 5, 2012 at 9:13 PM, Frank Horowitz  > wrote:
>
>> Still more checking things out. A problem (THE problem???) can be
>> demonstrated with your code snippet from a few messages back if you
>> construct the Wavelet like so rather than using the default constructor:
>>
>> Wavelet(WholeExtent=(-100,100,-100,100,-100,100))
>>
>> N.B. this is on a different machine, with a stock Mac OSX  binary version
>> 3.14.1 installed.
>>
>>
>> On 05/09/2012, at 5:49 PM, Frank Horowitz wrote:
>>
>> > OK, more tracking the problem:
>> >
>> > The assumption underlying my "close up" description below is not
>> correct.  The problem appears to be the actual _size_ of the dataset I'm
>> feeding in. Even a single time step fed in to the filter (with some 866,000
>> points) hangs at the x*x computation, even though that appeared to work
>> perfectly fine with your 3D Wavelet as a source (only some 9,000 points).
>> >
>> > I'm guessing a memory leak somewhere?  When was the last time paraview
>> was run through valgrind?
>> >
>> > Cheers (and sorry about the red herring earlier),
>> >   Frank
>> >
>> >
>> >
>> > On 05/09/2012, at 4:42 PM, Frank Horowitz wrote:
>> >
>> >> OK, the 30,000 foot view:
>> >>
>> >> I have a time series of about 20 slices (currently) of values of each
>> component (3 components total) of a vector displacement 3d spatial field
>> from an elastic wave propagation simulation. At the moment, I'm simply
>> trying to visualise the vector length to see the propagating phases of the
>> wavefield. Because my upstream codes and scripts only write individual
>> components of the vector field, I need to combine the components and then
>> calculate the vector lengths *somewhere*.  The programmable filter seemed
>> like a good place to try that, since once I am past the learning curve for
>> it, I'll have the full flexibilty of numpy for doing almost anything else I
>> want with these and similar data from my upstream simulation codes.
>> >>
>> >> Now the 10,000 foot view:
>> >> I've been trying to simply calculate the sqrt(x*x + y*y + z*z) scalar
>> field (where x, y, and z are the respective components as a function of
>> position)  and let that flow on down the pipeline as a new spatial field to
>> be visualised. I assumed (correctly if I'm interpreting your answer
>> correctly) that if I built the filter to deal with the spatial aspects of
>> the calculation, the time series would be taken care of by paraview
>> iterating over all files pointed to by inputs[0] (a set of files containing
>> timesteps of the x components), inputs[1] (ditto for the y components), and
>> inputs[2] (ditto for the z components).  The set of *.vtu files for each
>> component have filenames specifying component and time step (e.g.: x001,
>> x00

Re: [Paraview] Getting String arrays out of Field Data

[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] Newbie question on Python Programmable Filters

[Berk Geveci]

I just ran this with a 500^3 wavelet on my machine without any problems.
The whole thing (including a slice filter) takes 20-30 seconds but the
programmable filter is actually pretty much instantaneous (<0.1 seconds)
using a debug build from Git master. Is there any way you can attach to the
stuck ParaView with a debugger and look at the stack trace?

On Wed, Sep 5, 2012 at 9:13 PM, Frank Horowitz
wrote:

> Still more checking things out. A problem (THE problem???) can be
> demonstrated with your code snippet from a few messages back if you
> construct the Wavelet like so rather than using the default constructor:
>
> Wavelet(WholeExtent=(-100,100,-100,100,-100,100))
>
> N.B. this is on a different machine, with a stock Mac OSX  binary version
> 3.14.1 installed.
>
>
> On 05/09/2012, at 5:49 PM, Frank Horowitz wrote:
>
> > OK, more tracking the problem:
> >
> > The assumption underlying my "close up" description below is not
> correct.  The problem appears to be the actual _size_ of the dataset I'm
> feeding in. Even a single time step fed in to the filter (with some 866,000
> points) hangs at the x*x computation, even though that appeared to work
> perfectly fine with your 3D Wavelet as a source (only some 9,000 points).
> >
> > I'm guessing a memory leak somewhere?  When was the last time paraview
> was run through valgrind?
> >
> > Cheers (and sorry about the red herring earlier),
> >   Frank
> >
> >
> >
> > On 05/09/2012, at 4:42 PM, Frank Horowitz wrote:
> >
> >> OK, the 30,000 foot view:
> >>
> >> I have a time series of about 20 slices (currently) of values of each
> component (3 components total) of a vector displacement 3d spatial field
> from an elastic wave propagation simulation. At the moment, I'm simply
> trying to visualise the vector length to see the propagating phases of the
> wavefield. Because my upstream codes and scripts only write individual
> components of the vector field, I need to combine the components and then
> calculate the vector lengths *somewhere*.  The programmable filter seemed
> like a good place to try that, since once I am past the learning curve for
> it, I'll have the full flexibilty of numpy for doing almost anything else I
> want with these and similar data from my upstream simulation codes.
> >>
> >> Now the 10,000 foot view:
> >> I've been trying to simply calculate the sqrt(x*x + y*y + z*z) scalar
> field (where x, y, and z are the respective components as a function of
> position)  and let that flow on down the pipeline as a new spatial field to
> be visualised. I assumed (correctly if I'm interpreting your answer
> correctly) that if I built the filter to deal with the spatial aspects of
> the calculation, the time series would be taken care of by paraview
> iterating over all files pointed to by inputs[0] (a set of files containing
> timesteps of the x components), inputs[1] (ditto for the y components), and
> inputs[2] (ditto for the z components).  The set of *.vtu files for each
> component have filenames specifying component and time step (e.g.: x001,
> x002, x003, etc., y001, y002, y003, etc., and z001, z002, z003 etc.).
> >>
> >> Now the close-up:
> >> As far as I can tell, after compiling from source from the git
> "release" branch (on an up-to-date Linux Mint distro -- which tracks Ubuntu
> Precise), my approach will work for a set of 3 single timestep arrays fed
> into the "ports" (or whatever paraview's terminology is for the inputs to a
> filter in the pipeline). In other words, a single timeslice
> (x001,y001,z001) appears on inputs[0:3] (numpy slice notation) and is
> processed correctly. (I'll go away and verify that shortly.)  I *think* the
> "hangs/race-conditions/whatever" start when there are multiple timesteps
> presented to each input port.  I was unsure whether the 'inputs[0],
> inputs[1], inputs[2]' syntax needed to be extended to deal with time series
> as above, and could not find documentation on the inputs[] array after
> extensive searching -- just a few snippets of example code that were all
> consistent with paraview dealing with the time series on its own.
> >>
> >> I'm inexperienced enough with paraview to not know if I have a
> conceptual error with the way I'm approaching this problem, or if I have
> simply stumbled over a bug. That's why I've posted to the mailing list!
> >>
> >> As always, your (and the rest of the list's) help is greatly
> appreciated!
> >>
> >> Cheers,
> >>  Frank
> >>
> >>
> >> On 05/09/2012, at 4:01 PM, Andy Bauer wrote:
> >>
> >>> I'm not sure I completely understand about your inputs[0]. As for
> iterating
> >>> over time with the python programmable filter, inside the filter it
> will
> >>> only have access to values from a single time step at each invocation
> of
> >>> RequestData(). You can think of the order of operations if you have 5
> time
> >>> steps as:
> >>> 1) iterate over each of the 5 time steps (i.e. load a file and update
> the
> >>> pipeline with that data set)
> >>> 2) update 

Re: [Paraview] Maximize ParaView across 3 monitors; NVIDIA

[Aashish Chaudhary]

On Thu, Sep 6, 2012 at 12:38 PM, David Zemon  wrote:

>  I'm looking into X extensions now.
>
> Thank you for the link to the 590. It's very promising. I dug into their
> website a bit more and it looks like 3D Vision Surround is only supported
> on Windows. It's good to know it's out there as an option, but hopefully X
> extensions will get us going without switching OSes and buying a new card.
>

Great. Also have a look at this link:
http://awesome.naquadah.org/wiki/Using_Multiple_Screens . I forgot all the
details but in our setup we have six node cluster each having a independent
card and we were able to get it working with x extensions but again it was
couple of years ago and I forgot all the details.

Thanks


>
> Thanks,
> David
>
>
> On 09/06/2012 11:15 AM, Aashish Chaudhary wrote:
>
> Hi David,
>
>  I have heard that GTK 590 supports up to three monitors.
>
>  http://www.geforce.com/hardware/desktop-gpus/geforce-gtx-590/features
>
>  There is an alternative via X extensions (with second GPU), it seemed to
> work for us couple of years ago but I am not sure the current status of it.
>
>  Thanks,
>
>
> On Thu, Sep 6, 2012 at 11:27 AM, David Zemon  wrote:
>
>>  Software: Ubuntu 12.04.1 running Xfce4
>> Hardware: Microway 1U rack, 1 x NVidia Quadro 5000, 1 x Tesla C2050
>> (access to more teslas), 3 x Dells S500 projectors, 1 x totally awesome
>> custom built immersive screen :)
>>
>> *Goal:* Maximize ParaView across all three monitors for seamless
>> integration.
>>
>> *Problem:* the Quadro 5000 supports only 2 monitors which means we need
>> a second GPU. Second GPU means second X screen and, as far as I know, a
>> window can not be maximized across two X screens.
>>
>> *Help:* What are my options for maximizing ParaView across all three
>> monitors?
>>
>> *Known solutions:* The only solution I've found so far is to use an AMD
>> card. Two problems with that: 1) the driver setup is supposed to be a real
>> PITA, 2) we don't have a FirePro sitting around to try out, so that would
>> mean buying one.
>>
>> I'm open to changing up the software stack if you know of something that
>> will work.
>>
>> Thanks,
>> David Zemon
>>
>> ___
>> 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
>>
>>
>
>
>  --
> | Aashish Chaudhary
> | R&D Engineer
> | Kitware Inc.
> | www.kitware.com
>
>
>


-- 
| Aashish Chaudhary
| R&D Engineer
| Kitware Inc.
| www.kitware.com
___
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] Maximize ParaView across 3 monitors; NVIDIA

[David Zemon]

I'm looking into X extensions now.

Thank you for the link to the 590. It's very promising. I dug into their 
website a bit more and it looks like 3D Vision Surround is only 
supported on Windows. It's good to know it's out there as an option, but 
hopefully X extensions will get us going without switching OSes and 
buying a new card.


Thanks,
David

On 09/06/2012 11:15 AM, Aashish Chaudhary wrote:

Hi David,

I have heard that GTK 590 supports up to three monitors.

http://www.geforce.com/hardware/desktop-gpus/geforce-gtx-590/features

There is an alternative via X extensions (with second GPU), it seemed 
to work for us couple of years ago but I am not sure the current 
status of it.


Thanks,


On Thu, Sep 6, 2012 at 11:27 AM, David Zemon > wrote:


Software: Ubuntu 12.04.1 running Xfce4
Hardware: Microway 1U rack, 1 x NVidia Quadro 5000, 1 x Tesla
C2050 (access to more teslas), 3 x Dells S500 projectors, 1 x
totally awesome custom built immersive screen :)

*Goal:* Maximize ParaView across all three monitors for seamless
integration.

*Problem:* the Quadro 5000 supports only 2 monitors which means we
need a second GPU. Second GPU means second X screen and, as far as
I know, a window can not be maximized across two X screens.

*Help:* What are my options for maximizing ParaView across all
three monitors?

*Known solutions:* The only solution I've found so far is to use
an AMD card. Two problems with that: 1) the driver setup is
supposed to be a real PITA, 2) we don't have a FirePro sitting
around to try out, so that would mean buying one.

I'm open to changing up the software stack if you know of
something that will work.

Thanks,
David Zemon

___
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




--
| Aashish Chaudhary
| R&D Engineer
| Kitware Inc.
| www.kitware.com 


___
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] Plugin License question

[Felipe Bordeu]

Hello,

In our lab we use a home made plugin to read the solutions files of ours 
computations.
The plugin is written in c++ (using Qt for the UI, xdmf for the storage 
and blas for the calculations).

So nothing else that the libraries downloaded during a Superbuild and Qt

The question is, can I put only the binary of the plugin in a webpage 
for free ???


We don't want to put the source of the plugin because is in constant 
development.

I don't know a lot about the Qt licenses (commercial LGPL).

Any advice ?


--

Felipe Bordeu Weldt
Ingénieur de Recherche
-
Tél. : 33 (0)2 40 37 16 57
Fax. : 33 (0)2 40 74 74 06
felipe.bor...@ec-nantes.fr
Institut GeM - UMR CNRS 6183
École Centrale Nantes
1 Rue de La Noë, 44321 Nantes, FRANCE
-

___
Powered by www.kitware.com

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

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

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

Re: [Paraview] Newbie question on Python Programmable Filters

[Berk Geveci]

Given the trivial amount of data involved, I can't imagine this being a
memory leak. Maybe if it hung after 10 time steps... Do you have a way of
verifying that your numpy installation works properly? So for example, what
happens if you do this in the programmable filter:

import numpy
a = numpy.arange(100)
m = a * a

-berk

On Wed, Sep 5, 2012 at 5:49 PM, Frank Horowitz
wrote:

> OK, more tracking the problem:
>
> The assumption underlying my "close up" description below is not correct.
>  The problem appears to be the actual _size_ of the dataset I'm feeding in.
> Even a single time step fed in to the filter (with some 866,000 points)
> hangs at the x*x computation, even though that appeared to work perfectly
> fine with your 3D Wavelet as a source (only some 9,000 points).
>
> I'm guessing a memory leak somewhere?  When was the last time paraview was
> run through valgrind?
>
> Cheers (and sorry about the red herring earlier),
> Frank
>
>
>
> On 05/09/2012, at 4:42 PM, Frank Horowitz wrote:
>
> > OK, the 30,000 foot view:
> >
> > I have a time series of about 20 slices (currently) of values of each
> component (3 components total) of a vector displacement 3d spatial field
> from an elastic wave propagation simulation. At the moment, I'm simply
> trying to visualise the vector length to see the propagating phases of the
> wavefield. Because my upstream codes and scripts only write individual
> components of the vector field, I need to combine the components and then
> calculate the vector lengths *somewhere*.  The programmable filter seemed
> like a good place to try that, since once I am past the learning curve for
> it, I'll have the full flexibilty of numpy for doing almost anything else I
> want with these and similar data from my upstream simulation codes.
> >
> > Now the 10,000 foot view:
> > I've been trying to simply calculate the sqrt(x*x + y*y + z*z) scalar
> field (where x, y, and z are the respective components as a function of
> position)  and let that flow on down the pipeline as a new spatial field to
> be visualised. I assumed (correctly if I'm interpreting your answer
> correctly) that if I built the filter to deal with the spatial aspects of
> the calculation, the time series would be taken care of by paraview
> iterating over all files pointed to by inputs[0] (a set of files containing
> timesteps of the x components), inputs[1] (ditto for the y components), and
> inputs[2] (ditto for the z components).  The set of *.vtu files for each
> component have filenames specifying component and time step (e.g.: x001,
> x002, x003, etc., y001, y002, y003, etc., and z001, z002, z003 etc.).
> >
> > Now the close-up:
> > As far as I can tell, after compiling from source from the git "release"
> branch (on an up-to-date Linux Mint distro -- which tracks Ubuntu Precise),
> my approach will work for a set of 3 single timestep arrays fed into the
> "ports" (or whatever paraview's terminology is for the inputs to a filter
> in the pipeline). In other words, a single timeslice (x001,y001,z001)
> appears on inputs[0:3] (numpy slice notation) and is processed correctly.
> (I'll go away and verify that shortly.)  I *think* the
> "hangs/race-conditions/whatever" start when there are multiple timesteps
> presented to each input port.  I was unsure whether the 'inputs[0],
> inputs[1], inputs[2]' syntax needed to be extended to deal with time series
> as above, and could not find documentation on the inputs[] array after
> extensive searching -- just a few snippets of example code that were all
> consistent with paraview dealing with the time series on its own.
> >
> > I'm inexperienced enough with paraview to not know if I have a
> conceptual error with the way I'm approaching this problem, or if I have
> simply stumbled over a bug. That's why I've posted to the mailing list!
> >
> > As always, your (and the rest of the list's) help is greatly appreciated!
> >
> > Cheers,
> >   Frank
> >
> >
> > On 05/09/2012, at 4:01 PM, Andy Bauer wrote:
> >
> >> I'm not sure I completely understand about your inputs[0]. As for
> iterating
> >> over time with the python programmable filter, inside the filter it will
> >> only have access to values from a single time step at each invocation of
> >> RequestData(). You can think of the order of operations if you have 5
> time
> >> steps as:
> >> 1) iterate over each of the 5 time steps (i.e. load a file and update
> the
> >> pipeline with that data set)
> >> 2) update the python programmable filter output for the current data set
> >> resulting from the loaded file
> >>
> >> If you want time statistics for your variables I think you probably
> want to
> >> use the calculator filter to square the value at each point and then use
> >> the temporal statistics filter to get an average of those values. After
> >> that you can just multiply the results by the number of time steps.
> >>
> >> If I'm misunderstanding what you're trying to do, I'd suggest giving a
> >> d

Re: [Paraview] Newbie question on Python Programmable Filters

[Berk Geveci]

When you have a time series, the programmable filter will work on the
current time step only. So if you animate or change time through the GUI,
the filter should update properly. But each execution of the filter would
work only on 1 timestep.

On Wed, Sep 5, 2012 at 2:56 PM, Frank Horowitz
wrote:

>  OK Andy,
>
> Thanks for that!
>
> I've compiled and installed Paraview from git's release branch, and (a
> slight modification) of your code runs correctly. I think that means that
> my paraview/python integration is semi-sane.
>
> However, my problem still occurs. Recall that I have a time series of .vtu
> files (indexed by integers built-in to their filenames) being fed into
> inputs[0] (at least as I understand it).
>
> In my original code snippet, I assumed that the "m = x*x" statement would
> iterate not only over the point values of the spatial array, but also over
> all timeslices.  Is that a correct assumption? Is that the source of my
> problems???
>
> Thanks again for your help,
> Frank Horowitz
>
>  --
> *From:* Andy Bauer [andy.ba...@kitware.com]
> *Sent:* Tuesday, September 04, 2012 7:57 PM
> *To:* Frank Horowitz
> *Cc:* paraview@paraview.org
> *Subject:* Re: [Paraview] Newbie question on Python Programmable Filters
>
>  Hi Frank,
>
> I didn't have a problem with the following on my windows machine using the
> ParaView 3.14.1 installer.
> ===
> try: paraview.simple
>
> except: from paraview.simple import *
>
> paraview.simple._DisableFirstRenderCameraReset()
>
> Wavelet1 = Wavelet()
>
> RenderView1 = GetRenderView()
>
> ProgrammableFilter1 = ProgrammableFilter()
>
>  ProgrammableFilter1.PythonPath = ''
>
> ProgrammableFilter1.RequestInformationScript = ''
>
> ProgrammableFilter1.Script = 'x = inputs[0].PointData[\'RTData\']\nprint
> "x shape =",x.shape\n\nm = x*x # I _thought_ this should run at numpy
> speeds?\nprint "m shape=",m.shape\n\n# more debugging code to finish
> computation of the norm omitted\n\noutput.PointData.append(m,"Displacement
> Norm")'
>
> Show()
>
> Render()
>
> =
>
> Does this work for you?  If it doesn't, I wonder if it's an issue with
> having multiple pythons and numpy or system environments causing problems.
> What version of paraview are you using and how did you get it
> built/installed on your machine?
>
> Andy
>
> On Tue, Sep 4, 2012 at 6:04 PM, Frank Horowitz  > wrote:
>
>> Hal Canary wrote on Tue Sep 4 15:22:01 EDT 2012:
>> > On 09/04/2012 02:20 PM, Frank Horowitz wrote:
>> > >
>> >  x = inputs[0].PointData['Scalars_']
>> >
>> >
>> >
>> > I thought one needs to convert a vtkarray to a numpy array with
>> >
>> >   x = numpy.array(inputs[0].PointData['Scalars_'])
>> >
>> > before doing anything with it.
>>
>>
>>  Computations after that style of construction hang too.  To my eye, the
>> results of such an expression should be 100% standard numpy, obeying 100%
>> standard numpy semantics.
>>
>> It appears that there is a bug in numpy integration to the Python
>> Programmable Filter.
>>
>> I'll re-compile from source and report back on the results.  Is the bug
>> tracking system evident from the homepage?
>>
>> Cheers,
>> Frank Horowitz
>>
>>
>>
>>
>> ___
>> 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] Newbie question on Python Programmable Filters

[Berk Geveci]

I am just starting to read this thread so I wanted to answer various
questions here separately.

Your assumption is correct. m = x*x will be 100% evaluated by NumPy. In
theory, this should be pretty much instantaneous for ~1M elements, even for
20 timeslices.

On Tue, Sep 4, 2012 at 2:20 PM, Frank Horowitz
wrote:

> Hi All,
>
> I'm trying to visualize some vector component field time series written by
> a seismic wave propagation code, and converted into .vtu files by a set of
> upstream scripts.
>
> I'm trying to compute the norm of the vector at each timeslice and then
> visualize the resulting scalar field. The data are all PointData. There are
> approximately 866,000 values for each compoment per timeslice, and about 20
> time slices.
>
> I'm trying to do this computation via the Python Programmable Filter
> (mostly so I can learn its tricks for future reference). However, I'm
> running into an apparent hang. The following code (complete with many more
> debugging print statements than I would use in production) is the simplest
> code I can make fail reliably:
>
> --- begin code snippet ---
>
> x = inputs[0].PointData['Scalars_']
> print "x shape =",x.shape
>
> y = inputs[1].PointData['Scalars_']
> print "y shape =",y.shape
>
> z = inputs[2].PointData['Scalars_']
> print "z shape =",z.shape
>
> m = x*x  # I _thought_ this should run at numpy speeds…
> print "m shape=",m.shape
>
> # more debugging code to finish computation of the norm omitted
>
> output.PointData.append(m,"Displacement Norm")
>
> --- end code snippet ---
>
> I can run everything -- cutting and pasting initial subsets of the code
> into the program text box -- with the print output I expected all the way
> to the "m = x*x" statement.  If I include the "m=x*x" statement the process
> hangs. It does not even print out the preceding print statements.
>
> I was under the impression that such syntax as my "x*x" used above would
> be evaluated by numpy, and the (implicit) loop should be executed at C
> speeds.  Am I correct in that understanding?  If not, how should I go about
> coding an equivalent computation in the Programmable Filter that would
> execute at C (numpy) speeds?
>
> Thanks in advance for any help you might be able to provide,
> Frank Horowitz
> Cornell, Earth and Atmospheric Sciences
>
>
>
>
>
> ___
> 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] Maximize ParaView across 3 monitors; NVIDIA

[Aashish Chaudhary]

Hi David,

I have heard that GTK 590 supports up to three monitors.

http://www.geforce.com/hardware/desktop-gpus/geforce-gtx-590/features

There is an alternative via X extensions (with second GPU), it seemed to
work for us couple of years ago but I am not sure the current status of it.

Thanks,


On Thu, Sep 6, 2012 at 11:27 AM, David Zemon  wrote:

>  Software: Ubuntu 12.04.1 running Xfce4
> Hardware: Microway 1U rack, 1 x NVidia Quadro 5000, 1 x Tesla C2050
> (access to more teslas), 3 x Dells S500 projectors, 1 x totally awesome
> custom built immersive screen :)
>
> *Goal:* Maximize ParaView across all three monitors for seamless
> integration.
>
> *Problem:* the Quadro 5000 supports only 2 monitors which means we need a
> second GPU. Second GPU means second X screen and, as far as I know, a
> window can not be maximized across two X screens.
>
> *Help:* What are my options for maximizing ParaView across all three
> monitors?
>
> *Known solutions:* The only solution I've found so far is to use an AMD
> card. Two problems with that: 1) the driver setup is supposed to be a real
> PITA, 2) we don't have a FirePro sitting around to try out, so that would
> mean buying one.
>
> I'm open to changing up the software stack if you know of something that
> will work.
>
> Thanks,
> David Zemon
>
> ___
> 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
>
>


-- 
| Aashish Chaudhary
| R&D Engineer
| Kitware Inc.
| www.kitware.com
___
Powered by www.kitware.com

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

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

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

Re: [Paraview] Newbie question on Python Programmable Filters

[Andy Bauer]

Hmm, I tried with the WholeExtent of (-100, 100, -100, 100, -100, 100) on
my machine and still didn't have problems. How much memory is on your
system? I'd hate to think that ParaView is that poor with memory
management. If I remember correctly, we used to have nightly dashboards
running valgrind. I know that Bill Lorensen's project last fall was going
through with valgrind to get rid of a bunch of memory issues.

Well, the good news is that I think I have a better solution for you
anyways. I get the feeling from your original script that you have a
separate file for each time step and each component of your vector field.
If all components of your vector field are in the same file you can ignore
the first part of the steps below. Here's the solution:

1) load all of the files for the first time step (e.g. load the
..
2) use control and select all the file readers
3) use the append attributes filter to combine the fields from each input
data set into a single output dataset (it assumes that the input data sets
all have the same number of points and cells and the desired fields all
have different names -- use the calculator filter to change the field names
if needed)

4) create a calculator filter
4a) compute the magnitude by inputting in "mag(*iHat + *jHat + *kHat)" for the widget right above the calculator buttons
4b) if you want a specific name for the newly created field, set it in
Result Array Name

Let me know if this accomplishes your goal.

Andy

On Wed, Sep 5, 2012 at 9:13 PM, Frank Horowitz
wrote:

> Still more checking things out. A problem (THE problem???) can be
> demonstrated with your code snippet from a few messages back if you
> construct the Wavelet like so rather than using the default constructor:
>
> Wavelet(WholeExtent=(-100,100,-100,100,-100,100))
>
> N.B. this is on a different machine, with a stock Mac OSX  binary version
> 3.14.1 installed.
>
>
> On 05/09/2012, at 5:49 PM, Frank Horowitz wrote:
>
> > OK, more tracking the problem:
> >
> > The assumption underlying my "close up" description below is not
> correct.  The problem appears to be the actual _size_ of the dataset I'm
> feeding in. Even a single time step fed in to the filter (with some 866,000
> points) hangs at the x*x computation, even though that appeared to work
> perfectly fine with your 3D Wavelet as a source (only some 9,000 points).
> >
> > I'm guessing a memory leak somewhere?  When was the last time paraview
> was run through valgrind?
> >
> > Cheers (and sorry about the red herring earlier),
> >   Frank
> >
> >
> >
> > On 05/09/2012, at 4:42 PM, Frank Horowitz wrote:
> >
> >> OK, the 30,000 foot view:
> >>
> >> I have a time series of about 20 slices (currently) of values of each
> component (3 components total) of a vector displacement 3d spatial field
> from an elastic wave propagation simulation. At the moment, I'm simply
> trying to visualise the vector length to see the propagating phases of the
> wavefield. Because my upstream codes and scripts only write individual
> components of the vector field, I need to combine the components and then
> calculate the vector lengths *somewhere*.  The programmable filter seemed
> like a good place to try that, since once I am past the learning curve for
> it, I'll have the full flexibilty of numpy for doing almost anything else I
> want with these and similar data from my upstream simulation codes.
> >>
> >> Now the 10,000 foot view:
> >> I've been trying to simply calculate the sqrt(x*x + y*y + z*z) scalar
> field (where x, y, and z are the respective components as a function of
> position)  and let that flow on down the pipeline as a new spatial field to
> be visualised. I assumed (correctly if I'm interpreting your answer
> correctly) that if I built the filter to deal with the spatial aspects of
> the calculation, the time series would be taken care of by paraview
> iterating over all files pointed to by inputs[0] (a set of files containing
> timesteps of the x components), inputs[1] (ditto for the y components), and
> inputs[2] (ditto for the z components).  The set of *.vtu files for each
> component have filenames specifying component and time step (e.g.: x001,
> x002, x003, etc., y001, y002, y003, etc., and z001, z002, z003 etc.).
> >>
> >> Now the close-up:
> >> As far as I can tell, after compiling from source from the git
> "release" branch (on an up-to-date Linux Mint distro -- which tracks Ubuntu
> Precise), my approach will work for a set of 3 single timestep arrays fed
> into the "ports" (or whatever paraview's terminology is for the inputs to a
> filter in the pipeline). In other words, a single timeslice
> (x001,y001,z001) appears on inputs[0:3] (numpy slice notation) and is
> processed correctly. (I'll go away and verify that shortly.)  I *think* the
> "hangs/race-conditions/whatever" start when there are multiple timesteps
> presented to each input port.  I was unsure whether the 'inputs[0],
> inputs[1], inputs[2]' syntax needed to be extended to 

[Paraview] Maximize ParaView across 3 monitors; NVIDIA

[David Zemon]

Software: Ubuntu 12.04.1 running Xfce4
Hardware: Microway 1U rack, 1 x NVidia Quadro 5000, 1 x Tesla C2050 
(access to more teslas), 3 x Dells S500 projectors, 1 x totally awesome 
custom built immersive screen :)


*Goal:* Maximize ParaView across all three monitors for seamless 
integration.


*Problem:* the Quadro 5000 supports only 2 monitors which means we need 
a second GPU. Second GPU means second X screen and, as far as I know, a 
window can not be maximized across two X screens.


*Help:* What are my options for maximizing ParaView across all three 
monitors?


*Known solutions:* The only solution I've found so far is to use an AMD 
card. Two problems with that: 1) the driver setup is supposed to be a 
real PITA, 2) we don't have a FirePro sitting around to try out, so that 
would mean buying one.


I'm open to changing up the software stack if you know of something that 
will work.


Thanks,
David Zemon
___
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] Debugging reverse connections

[Utkarsh Ayachit]

Glad to help.

On Thu, Sep 6, 2012 at 11:05 AM, Renato Elias  wrote:
> Hi Utkarsh,
>
> I guess we've found the problem. The connection was really getting lost in
> our local network. After tweaking some servers in our lab, the things are
> now working.
>
> Thanks for your help!
>
> Renato.
>
>
> On Thu, Sep 6, 2012 at 11:58 AM, Utkarsh Ayachit
>  wrote:
>>
>> I think that's you first clue. The fact that even telnet cannot
>> connect back to your client machine/port indicates that something is
>> either blocking the connection or it's just getting lost in the
>> network. Does ping "my.ip.external.address" work from the server host?
>>
>> Utkarsh
>>
>> On Thu, Sep 6, 2012 at 10:03 AM, Renato Elias  wrote:
>> > Hi Utkarsh,
>> >
>> > neither an error message I get from pvserver. It just keep waiting for
>> > connection until exit due to time out... the client stays there with the
>> > "waiting for server connection" screen open and nothing happens.
>> >
>> > I'm not familiar with telnet commands but tried the following:
>> >
>> > telnet> open my.ip.external.address 1
>> > Trying my.ip.external.address...
>> >
>> > and it just hangs trying to connect.
>> >
>> > The client in this test is a Macbook running Lion. I've tried with the
>> > firewall open and closed as well.
>> >
>> > Any more hints?
>> >
>> > Renato.
>> >
>> > On Thu, Sep 6, 2012 at 10:46 AM, Utkarsh Ayachit
>> >  wrote:
>> >>
>> >> Renato,
>> >>
>> >> Here's what I get the server spewing when I run "pvserver --rc"
>> >> without any client already running:
>> >>
>> >> --
>> >> Connecting to client (reverse connection requested)...
>> >> Connection URL: csrc://miranda:1
>> >> ERROR: In
>> >>
>> >> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkSocket.cxx,
>> >> line 481
>> >> vtkClientSocket (0xfbc520): Socket error in call to connect. Connection
>> >> refused.
>> >>
>> >> ERROR: In
>> >>
>> >> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkClientSocket.cxx,
>> >> line 53
>> >> vtkClientSocket (0xfbc520): Failed to connect to server miranda:1
>> >>
>> >> Warning: In
>> >>
>> >> /home/utkarsh/Kitware/ParaView3/ParaView/ParaViewCore/ClientServerCore/Core/vtkTCPNetworkAccessManager.cxx,
>> >> line 302
>> >> vtkTCPNetworkAccessManager (0x8cdc50): Connect failed.  Retrying for
>> >> 59.9997 more seconds.
>> >>
>> >> ERROR: In
>> >>
>> >> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkSocket.cxx,
>> >> line 481
>> >> vtkClientSocket (0xfbc520): Socket error in call to connect. Connection
>> >> refused.
>> >>
>> >> ERROR: In
>> >>
>> >> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkClientSocket.cxx,
>> >> line 53
>> >> vtkClientSocket (0xfbc520): Failed to connect to server miranda:1
>> >> --
>> >>
>> >> Thus if the socket connection failed, it should start reporting error
>> >> and attempting again. Since the server totally just freezes, even when
>> >> no client is waiting for it at the other end, I am getting somehow the
>> >> BSD socket "connect" call is never returning success or failure.
>> >>
>> >> Try this: from the server node, try to telnet to you client host/port.
>> >> Does that work?
>> >>
>> >> Utkarsh
>> >>
>> >>
>> >> On Thu, Sep 6, 2012 at 9:34 AM, Renato Elias  wrote:
>> >> > Hi everybody,
>> >> >
>> >> > I'd like to understand what is happening with a client-server
>> >> > connection
>> >> > I'm
>> >> > trying to create in PV-3.14.1.
>> >> >
>> >> > pvserver is running on a SGI Altix ICE compiled with OpenMPI 1.6 and
>> >> > Mesa
>> >> > library for offscreen rendering. Everything run smoothly from the
>> >> > compilation and installation but when I try to connect with a client
>> >> > (same
>> >> > version), they can't get connected. The machine runs Torque batch
>> >> > system
>> >> > but
>> >> > I certified that pvserver was running after the client. Tests I've
>> >> > done
>> >> > to
>> >> > remove some common problems.
>> >> >
>> >> > 1). Client and server in the same network
>> >> >
>> >> > case a). Server running in the computing nodes under Torque batch
>> >> > system.
>> >> > Nothing happens. pvserver just sits with the following message:
>> >> >
>> >> > Connecting to client (reverse connection requested)...
>> >> > Connection URL: csrc://my.client.local.ip:1
>> >> >
>> >> > case b). Server running in the fron-end interactively (launched from
>> >> > command
>> >> > line). It would remove problems due to batch system and NAT not
>> >> > running
>> >> > in
>> >> > the computing nodes but the same message posted before appears and
>> >> > nothing
>> >> > happens.
>> >> >
>> >> > In both cases I also did tests with client's firewall turned off with
>> >> > similar result (nothing happens)
>> >> >
>> >> > 2). Client running remotely.
>> >> >
>> >> > Same behavior mentioned before.
>> >> >
>> >> > The question is: Is there a "verbose mode reverse connection" that
>> >> > could
>> >> > help me underst

Re: [Paraview] Debugging reverse connections

[Renato Elias]

Hi Utkarsh,

I guess we've found the problem. The connection was really getting lost in
our local network. After tweaking some servers in our lab, the things are
now working.

Thanks for your help!

Renato.

On Thu, Sep 6, 2012 at 11:58 AM, Utkarsh Ayachit <
utkarsh.ayac...@kitware.com> wrote:

> I think that's you first clue. The fact that even telnet cannot
> connect back to your client machine/port indicates that something is
> either blocking the connection or it's just getting lost in the
> network. Does ping "my.ip.external.address" work from the server host?
>
> Utkarsh
>
> On Thu, Sep 6, 2012 at 10:03 AM, Renato Elias  wrote:
> > Hi Utkarsh,
> >
> > neither an error message I get from pvserver. It just keep waiting for
> > connection until exit due to time out... the client stays there with the
> > "waiting for server connection" screen open and nothing happens.
> >
> > I'm not familiar with telnet commands but tried the following:
> >
> > telnet> open my.ip.external.address 1
> > Trying my.ip.external.address...
> >
> > and it just hangs trying to connect.
> >
> > The client in this test is a Macbook running Lion. I've tried with the
> > firewall open and closed as well.
> >
> > Any more hints?
> >
> > Renato.
> >
> > On Thu, Sep 6, 2012 at 10:46 AM, Utkarsh Ayachit
> >  wrote:
> >>
> >> Renato,
> >>
> >> Here's what I get the server spewing when I run "pvserver --rc"
> >> without any client already running:
> >>
> >> --
> >> Connecting to client (reverse connection requested)...
> >> Connection URL: csrc://miranda:1
> >> ERROR: In
> >>
> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkSocket.cxx,
> >> line 481
> >> vtkClientSocket (0xfbc520): Socket error in call to connect. Connection
> >> refused.
> >>
> >> ERROR: In
> >>
> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkClientSocket.cxx,
> >> line 53
> >> vtkClientSocket (0xfbc520): Failed to connect to server miranda:1
> >>
> >> Warning: In
> >>
> /home/utkarsh/Kitware/ParaView3/ParaView/ParaViewCore/ClientServerCore/Core/vtkTCPNetworkAccessManager.cxx,
> >> line 302
> >> vtkTCPNetworkAccessManager (0x8cdc50): Connect failed.  Retrying for
> >> 59.9997 more seconds.
> >>
> >> ERROR: In
> >>
> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkSocket.cxx,
> >> line 481
> >> vtkClientSocket (0xfbc520): Socket error in call to connect. Connection
> >> refused.
> >>
> >> ERROR: In
> >>
> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkClientSocket.cxx,
> >> line 53
> >> vtkClientSocket (0xfbc520): Failed to connect to server miranda:1
> >> --
> >>
> >> Thus if the socket connection failed, it should start reporting error
> >> and attempting again. Since the server totally just freezes, even when
> >> no client is waiting for it at the other end, I am getting somehow the
> >> BSD socket "connect" call is never returning success or failure.
> >>
> >> Try this: from the server node, try to telnet to you client host/port.
> >> Does that work?
> >>
> >> Utkarsh
> >>
> >>
> >> On Thu, Sep 6, 2012 at 9:34 AM, Renato Elias  wrote:
> >> > Hi everybody,
> >> >
> >> > I'd like to understand what is happening with a client-server
> connection
> >> > I'm
> >> > trying to create in PV-3.14.1.
> >> >
> >> > pvserver is running on a SGI Altix ICE compiled with OpenMPI 1.6 and
> >> > Mesa
> >> > library for offscreen rendering. Everything run smoothly from the
> >> > compilation and installation but when I try to connect with a client
> >> > (same
> >> > version), they can't get connected. The machine runs Torque batch
> system
> >> > but
> >> > I certified that pvserver was running after the client. Tests I've
> done
> >> > to
> >> > remove some common problems.
> >> >
> >> > 1). Client and server in the same network
> >> >
> >> > case a). Server running in the computing nodes under Torque batch
> >> > system.
> >> > Nothing happens. pvserver just sits with the following message:
> >> >
> >> > Connecting to client (reverse connection requested)...
> >> > Connection URL: csrc://my.client.local.ip:1
> >> >
> >> > case b). Server running in the fron-end interactively (launched from
> >> > command
> >> > line). It would remove problems due to batch system and NAT not
> running
> >> > in
> >> > the computing nodes but the same message posted before appears and
> >> > nothing
> >> > happens.
> >> >
> >> > In both cases I also did tests with client's firewall turned off with
> >> > similar result (nothing happens)
> >> >
> >> > 2). Client running remotely.
> >> >
> >> > Same behavior mentioned before.
> >> >
> >> > The question is: Is there a "verbose mode reverse connection" that
> could
> >> > help me understand what is happening?
> >> >
> >> > Any help would be appreciated
> >> >
> >> > Thanks
> >> >
> >> > --
> >> > Renato N. Elias
> >> > =
> >> > Professor, at Department of Civil Engineering
> 

Re: [Paraview] Debugging reverse connections

[Utkarsh Ayachit]

I think that's you first clue. The fact that even telnet cannot
connect back to your client machine/port indicates that something is
either blocking the connection or it's just getting lost in the
network. Does ping "my.ip.external.address" work from the server host?

Utkarsh

On Thu, Sep 6, 2012 at 10:03 AM, Renato Elias  wrote:
> Hi Utkarsh,
>
> neither an error message I get from pvserver. It just keep waiting for
> connection until exit due to time out... the client stays there with the
> "waiting for server connection" screen open and nothing happens.
>
> I'm not familiar with telnet commands but tried the following:
>
> telnet> open my.ip.external.address 1
> Trying my.ip.external.address...
>
> and it just hangs trying to connect.
>
> The client in this test is a Macbook running Lion. I've tried with the
> firewall open and closed as well.
>
> Any more hints?
>
> Renato.
>
> On Thu, Sep 6, 2012 at 10:46 AM, Utkarsh Ayachit
>  wrote:
>>
>> Renato,
>>
>> Here's what I get the server spewing when I run "pvserver --rc"
>> without any client already running:
>>
>> --
>> Connecting to client (reverse connection requested)...
>> Connection URL: csrc://miranda:1
>> ERROR: In
>> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkSocket.cxx,
>> line 481
>> vtkClientSocket (0xfbc520): Socket error in call to connect. Connection
>> refused.
>>
>> ERROR: In
>> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkClientSocket.cxx,
>> line 53
>> vtkClientSocket (0xfbc520): Failed to connect to server miranda:1
>>
>> Warning: In
>> /home/utkarsh/Kitware/ParaView3/ParaView/ParaViewCore/ClientServerCore/Core/vtkTCPNetworkAccessManager.cxx,
>> line 302
>> vtkTCPNetworkAccessManager (0x8cdc50): Connect failed.  Retrying for
>> 59.9997 more seconds.
>>
>> ERROR: In
>> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkSocket.cxx,
>> line 481
>> vtkClientSocket (0xfbc520): Socket error in call to connect. Connection
>> refused.
>>
>> ERROR: In
>> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkClientSocket.cxx,
>> line 53
>> vtkClientSocket (0xfbc520): Failed to connect to server miranda:1
>> --
>>
>> Thus if the socket connection failed, it should start reporting error
>> and attempting again. Since the server totally just freezes, even when
>> no client is waiting for it at the other end, I am getting somehow the
>> BSD socket "connect" call is never returning success or failure.
>>
>> Try this: from the server node, try to telnet to you client host/port.
>> Does that work?
>>
>> Utkarsh
>>
>>
>> On Thu, Sep 6, 2012 at 9:34 AM, Renato Elias  wrote:
>> > Hi everybody,
>> >
>> > I'd like to understand what is happening with a client-server connection
>> > I'm
>> > trying to create in PV-3.14.1.
>> >
>> > pvserver is running on a SGI Altix ICE compiled with OpenMPI 1.6 and
>> > Mesa
>> > library for offscreen rendering. Everything run smoothly from the
>> > compilation and installation but when I try to connect with a client
>> > (same
>> > version), they can't get connected. The machine runs Torque batch system
>> > but
>> > I certified that pvserver was running after the client. Tests I've done
>> > to
>> > remove some common problems.
>> >
>> > 1). Client and server in the same network
>> >
>> > case a). Server running in the computing nodes under Torque batch
>> > system.
>> > Nothing happens. pvserver just sits with the following message:
>> >
>> > Connecting to client (reverse connection requested)...
>> > Connection URL: csrc://my.client.local.ip:1
>> >
>> > case b). Server running in the fron-end interactively (launched from
>> > command
>> > line). It would remove problems due to batch system and NAT not running
>> > in
>> > the computing nodes but the same message posted before appears and
>> > nothing
>> > happens.
>> >
>> > In both cases I also did tests with client's firewall turned off with
>> > similar result (nothing happens)
>> >
>> > 2). Client running remotely.
>> >
>> > Same behavior mentioned before.
>> >
>> > The question is: Is there a "verbose mode reverse connection" that could
>> > help me understand what is happening?
>> >
>> > Any help would be appreciated
>> >
>> > Thanks
>> >
>> > --
>> > Renato N. Elias
>> > =
>> > Professor, at Department of Civil Engineering
>> > COPPE/Federal University of Rio de Janeiro
>> > PO Box 68506, Rio de Janeiro, RJ 21945-970, Brazil
>> >
>> >
>> > ___
>> > 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
>> >
>
>
>
>
> --
> Renato N. Elias
> ===

[Paraview] Getting String arrays out of Field Data

[Joshua Murphy]

Hello,

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?

Thanks,

Joshua Murphy
Laboratory for Atmospheric and Space Physics
University of Colorado, Boulder


___
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] Debugging reverse connections

[Renato Elias]

Hi Utkarsh,

neither an error message I get from pvserver. It just keep waiting for
connection until exit due to time out... the client stays there with the
"waiting for server connection" screen open and nothing happens.

I'm not familiar with telnet commands but tried the following:

telnet> open my.ip.external.address 1
Trying my.ip.external.address...

and it just hangs trying to connect.

The client in this test is a Macbook running Lion. I've tried with the
firewall open and closed as well.

Any more hints?

Renato.

On Thu, Sep 6, 2012 at 10:46 AM, Utkarsh Ayachit <
utkarsh.ayac...@kitware.com> wrote:

> Renato,
>
> Here's what I get the server spewing when I run "pvserver --rc"
> without any client already running:
>
> --
> Connecting to client (reverse connection requested)...
> Connection URL: csrc://miranda:1
> ERROR: In
> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkSocket.cxx,
> line 481
> vtkClientSocket (0xfbc520): Socket error in call to connect. Connection
> refused.
>
> ERROR: In
> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkClientSocket.cxx,
> line 53
> vtkClientSocket (0xfbc520): Failed to connect to server miranda:1
>
> Warning: In
> /home/utkarsh/Kitware/ParaView3/ParaView/ParaViewCore/ClientServerCore/Core/vtkTCPNetworkAccessManager.cxx,
> line 302
> vtkTCPNetworkAccessManager (0x8cdc50): Connect failed.  Retrying for
> 59.9997 more seconds.
>
> ERROR: In
> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkSocket.cxx,
> line 481
> vtkClientSocket (0xfbc520): Socket error in call to connect. Connection
> refused.
>
> ERROR: In
> /home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkClientSocket.cxx,
> line 53
> vtkClientSocket (0xfbc520): Failed to connect to server miranda:1
> --
>
> Thus if the socket connection failed, it should start reporting error
> and attempting again. Since the server totally just freezes, even when
> no client is waiting for it at the other end, I am getting somehow the
> BSD socket "connect" call is never returning success or failure.
>
> Try this: from the server node, try to telnet to you client host/port.
> Does that work?
>
> Utkarsh
>
>
> On Thu, Sep 6, 2012 at 9:34 AM, Renato Elias  wrote:
> > Hi everybody,
> >
> > I'd like to understand what is happening with a client-server connection
> I'm
> > trying to create in PV-3.14.1.
> >
> > pvserver is running on a SGI Altix ICE compiled with OpenMPI 1.6 and Mesa
> > library for offscreen rendering. Everything run smoothly from the
> > compilation and installation but when I try to connect with a client
> (same
> > version), they can't get connected. The machine runs Torque batch system
> but
> > I certified that pvserver was running after the client. Tests I've done
> to
> > remove some common problems.
> >
> > 1). Client and server in the same network
> >
> > case a). Server running in the computing nodes under Torque batch system.
> > Nothing happens. pvserver just sits with the following message:
> >
> > Connecting to client (reverse connection requested)...
> > Connection URL: csrc://my.client.local.ip:1
> >
> > case b). Server running in the fron-end interactively (launched from
> command
> > line). It would remove problems due to batch system and NAT not running
> in
> > the computing nodes but the same message posted before appears and
> nothing
> > happens.
> >
> > In both cases I also did tests with client's firewall turned off with
> > similar result (nothing happens)
> >
> > 2). Client running remotely.
> >
> > Same behavior mentioned before.
> >
> > The question is: Is there a "verbose mode reverse connection" that could
> > help me understand what is happening?
> >
> > Any help would be appreciated
> >
> > Thanks
> >
> > --
> > Renato N. Elias
> > =
> > Professor, at Department of Civil Engineering
> > COPPE/Federal University of Rio de Janeiro
> > PO Box 68506, Rio de Janeiro, RJ 21945-970, Brazil
> >
> >
> > ___
> > 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
> >
>



-- 
Renato N. Elias
=
Professor, at Department of Civil Engineering
COPPE/Federal University of Rio de Janeiro
PO Box 68506, Rio de Janeiro, RJ 21945-970, Brazil
___
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:
htt

Re: [Paraview] Debugging reverse connections

[Utkarsh Ayachit]

Renato,

Here's what I get the server spewing when I run "pvserver --rc"
without any client already running:

--
Connecting to client (reverse connection requested)...
Connection URL: csrc://miranda:1
ERROR: In 
/home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkSocket.cxx,
line 481
vtkClientSocket (0xfbc520): Socket error in call to connect. Connection refused.

ERROR: In 
/home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkClientSocket.cxx,
line 53
vtkClientSocket (0xfbc520): Failed to connect to server miranda:1

Warning: In 
/home/utkarsh/Kitware/ParaView3/ParaView/ParaViewCore/ClientServerCore/Core/vtkTCPNetworkAccessManager.cxx,
line 302
vtkTCPNetworkAccessManager (0x8cdc50): Connect failed.  Retrying for
59.9997 more seconds.

ERROR: In 
/home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkSocket.cxx,
line 481
vtkClientSocket (0xfbc520): Socket error in call to connect. Connection refused.

ERROR: In 
/home/utkarsh/Kitware/ParaView3/ParaView/VTK/Common/System/vtkClientSocket.cxx,
line 53
vtkClientSocket (0xfbc520): Failed to connect to server miranda:1
--

Thus if the socket connection failed, it should start reporting error
and attempting again. Since the server totally just freezes, even when
no client is waiting for it at the other end, I am getting somehow the
BSD socket "connect" call is never returning success or failure.

Try this: from the server node, try to telnet to you client host/port.
Does that work?

Utkarsh


On Thu, Sep 6, 2012 at 9:34 AM, Renato Elias  wrote:
> Hi everybody,
>
> I'd like to understand what is happening with a client-server connection I'm
> trying to create in PV-3.14.1.
>
> pvserver is running on a SGI Altix ICE compiled with OpenMPI 1.6 and Mesa
> library for offscreen rendering. Everything run smoothly from the
> compilation and installation but when I try to connect with a client (same
> version), they can't get connected. The machine runs Torque batch system but
> I certified that pvserver was running after the client. Tests I've done to
> remove some common problems.
>
> 1). Client and server in the same network
>
> case a). Server running in the computing nodes under Torque batch system.
> Nothing happens. pvserver just sits with the following message:
>
> Connecting to client (reverse connection requested)...
> Connection URL: csrc://my.client.local.ip:1
>
> case b). Server running in the fron-end interactively (launched from command
> line). It would remove problems due to batch system and NAT not running in
> the computing nodes but the same message posted before appears and nothing
> happens.
>
> In both cases I also did tests with client's firewall turned off with
> similar result (nothing happens)
>
> 2). Client running remotely.
>
> Same behavior mentioned before.
>
> The question is: Is there a "verbose mode reverse connection" that could
> help me understand what is happening?
>
> Any help would be appreciated
>
> Thanks
>
> --
> Renato N. Elias
> =
> Professor, at Department of Civil Engineering
> COPPE/Federal University of Rio de Janeiro
> PO Box 68506, Rio de Janeiro, RJ 21945-970, Brazil
>
>
> ___
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Please keep messages on-topic and check the ParaView Wiki at:
> http://paraview.org/Wiki/ParaView
>
> Follow this link to subscribe/unsubscribe:
> http://www.paraview.org/mailman/listinfo/paraview
>
___
Powered by www.kitware.com

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

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

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

[Paraview] Debugging reverse connections

[Renato Elias]

Hi everybody,

I'd like to understand what is happening with a client-server connection
I'm trying to create in PV-3.14.1.

pvserver is running on a SGI Altix ICE compiled with OpenMPI 1.6 and Mesa
library for offscreen rendering. Everything run smoothly from the
compilation and installation but when I try to connect with a client (same
version), they can't get connected. The machine runs Torque batch system
but I certified that pvserver was running after the client. Tests I've done
to remove some common problems.

1). Client and server in the same network

case a). Server running in the computing nodes under Torque batch system.
Nothing happens. pvserver just sits with the following message:

Connecting to client (reverse connection requested)...
Connection URL: csrc://my.client.local.ip:1

case b). Server running in the fron-end interactively (launched from
command line). It would remove problems due to batch system and NAT not
running in the computing nodes but the same message posted before appears
and nothing happens.

In both cases I also did tests with client's firewall turned off with
similar result (nothing happens)

2). Client running remotely.

Same behavior mentioned before.

The question is: Is there a "verbose mode reverse connection" that could
help me understand what is happening?

Any help would be appreciated

Thanks

-- 
Renato N. Elias
=
Professor, at Department of Civil Engineering
COPPE/Federal University of Rio de Janeiro
PO Box 68506, Rio de Janeiro, RJ 21945-970, Brazil
___
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