Thanks Burlen. I think I am finally done with this. I really appreciate all 
yours and others help in getting me here. thanks...

________________________________
From: Burlen Loring [[email protected]]
Sent: Wednesday, May 21, 2014 12:05 PM
To: Su, Simon M CTR USARMY ARL (US); [email protected]
Subject: Re: [Paraview] ParaView programmable filter -> what datatype to use to 
store complex number and complex number visualization

So you have 3 numpy arrays each with 3 complex values?

In any case you'll use vtk double  array for this, and you'll want to explictly 
handle real and imaginary parts. PV can iso surface the real and imaginary 
parts separately. There are a couple of ways you could do that. 1) each array 
of complex values could go into a 2 component vtkDoubleArray, first component 
is the real part and second component is the imaginary part. 2) each array of 
complex values could go into 2 1 component double arrays, the first double 
array for the real part and the second double array for the imaginary part. 
Since you mentioned iso surfacing I'm guessing option 2 will be the easiest 
path for you.

Burlen

On 05/21/2014 07:20 AM, Su, Simon M CTR USARMY ARL (US) wrote:
Hello,

I am calculating the roots using numpy in the python programmable filter. 
numpy.roots(coeff) is returning an array of type ndarray with complex128 array 
elements.


<type 'numpy.ndarray'>

[-0.77184451+1.11514251j -0.77184451-1.11514251j 0.54368901+0.j ]

array size = 3

(-0.771844506346+1.11514250804j)

<type 'numpy.complex128'>

(-0.771844506346-1.11514250804j)

<type 'numpy.complex128'>

(0.543689012692+0j)

<type 'numpy.complex128'>


In my calculation I will have 3 roots so I will be creating three variables 
root1, root2, root3, and  add them to the output data via


output.GetPointData().AddArray(root1)
output.GetPointData().AddArray(root2)
output.GetPointData().AddArray(root3)



what would the vtk type that I need to use to create root1, root2, and root3?

Can ParaView draw an isosurface using root1, or root2, or root3?

Any help on this is much appreciated.

thanks
-simon







_______________________________________________
Powered by www.kitware.com<http://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

Reply via email to