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

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