Hmm, I tried extracting out the components using the Calculator filter and then using the PythonCalculator filter to do that math.atan2 function. Underneath the covers I believe there is a conversion from a vtkDataArray to a numpy array and then back again. I don't think this is getting done for the math.atan2 function (it's not listed in the web page of python calculator functions) so you may be stuck here unless someone else has an idea.<br>
<br>Andy<br><br><div class="gmail_quote">On Wed, Feb 16, 2011 at 10:25 AM, Nico Schlömer <span dir="ltr"><<a href="mailto:nico.schloemer@gmail.com">nico.schloemer@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
That sounds *exactly* like what I need. Great!<br>
I'm having problems selecting the data sets though: Where<br>
"atan(psi__Y/psi__X)" works with the calculator filter, "math.atan2(<br>
psi__Y, psi__X)" fails in the PythonCalculator with<br>
<br>
============== *snip* ==============<br>
Traceback (most recent call last):<br>
File "<string>", line 19, in <module><br>
File "<string>", line 27, in RequestData<br>
NameError: name 'psi__Y' is not defined<br>
============== *snap* ==============<br>
<br>
I looked at <a href="http://www.paraview.org/Wiki/Python_Calculator" target="_blank">http://www.paraview.org/Wiki/Python_Calculator</a> but<br>
couldn't figure out what's wrong.<br>
<font color="#888888"><br>
--Nico<br>
</font><div><div></div><div class="h5"><br>
<br>
<br>
On Wed, Feb 16, 2011 at 2:43 PM, Tim Gallagher <<a href="mailto:tim.gallagher@gatech.edu">tim.gallagher@gatech.edu</a>> wrote:<br>
> You can use the Python Calculator, which will let you use any valid Python expression -- such as math.atan2().<br>
><br>
> Tim<br>
><br>
> ----- Original Message -----<br>
> From: "Nico Schlömer" <<a href="mailto:nico.schloemer@gmail.com">nico.schloemer@gmail.com</a>><br>
> To: "Andy Bauer" <<a href="mailto:andy.bauer@kitware.com">andy.bauer@kitware.com</a>><br>
> Cc: <a href="mailto:paraview@paraview.org">paraview@paraview.org</a><br>
> Sent: Wednesday, February 16, 2011 8:29:02 AM<br>
> Subject: Re: [Paraview] unstructured grid, programmable filters: extract, process array data<br>
><br>
> Almost there!<br>
><br>
>> atan(coordsX/coordsY)<br>
><br>
> If instead of coordsX I use psi__{X,Y}, things work. However, I need<br>
> something slightly different than atan: I'd like to get the complex<br>
> argument of a complex number X+iY (i.e., the angle with the positive<br>
> real axis). It's quite ~similar~ to what atan does, see for example<br>
> <<a href="http://gcc.gnu.org/onlinedocs/gcc-4.0.4/gfortran/ATAN2.html" target="_blank">http://gcc.gnu.org/onlinedocs/gcc-4.0.4/gfortran/ATAN2.html</a>>. In<br>
> Python, the function is called math.atan2() as well.<br>
> ParaView's calculator filter does however not support atan2(), and<br>
> that was my original motivation to look at programmable filters<br>
> actually. Well, it'd be much nicer if the same thing could be done<br>
> with the calculator of course!<br>
> Workarounds?<br>
><br>
> --Nico<br>
><br>
><br>
><br>
><br>
><br>
> On Tue, Feb 15, 2011 at 6:56 PM, Andy Bauer <<a href="mailto:andy.bauer@kitware.com">andy.bauer@kitware.com</a>> wrote:<br>
>> This sounds like a job for the calculator filter. "coordsX^2+coordsY^2" for<br>
>> your first example (or "psi__X^2+psi__Y^2" for the psi variables) and<br>
>> "atan(coordsX/coordsY)" for the second example.<br>
>><br>
>> It can operate on point or cell data (set in the Attribute Model drop-down<br>
>> menu), replace invalid numbers with a specified value, and you can name the<br>
>> resulting array as well.<br>
>><br>
>> Andy<br>
>><br>
>> 2011/2/15 Nico Schlömer <<a href="mailto:nico.schloemer@gmail.com">nico.schloemer@gmail.com</a>><br>
>>><br>
>>> Hi all,<br>
>>><br>
>>> I've been looking at<br>
>>> <<a href="http://www.paraview.org/Wiki/Python_Programmable_Filter" target="_blank">http://www.paraview.org/Wiki/Python_Programmable_Filter</a>> to create a<br>
>>> programmable filter, but I'm a bit lost as to whether or not I need<br>
>>> GetPolyData, CopyAttributes,...<br>
>>><br>
>>> Essentially, I have an unstructured grid in an ExodusII file (like the<br>
>>> attached file actually) with point data "psi" in it that has two<br>
>>> components (X,Y). I can plot psi:X, psi:Y, and psi:magnitude alright<br>
>>> with the default GUI options, but what I would really like to have<br>
>>> plotted is X^2+Y^2 (no sqrt()) and something along the lines of<br>
>>> arctan(Y/X).<br>
>>><br>
>>> Is there simple code somewhere some that takes array data off of an<br>
>>> unstructured grid, does something with it and returns the extra array<br>
>>> from the filter?<br>
>>><br>
>>> Cheers,<br>
>>> Nico<br>
>>><br>
>>> _______________________________________________<br>
>>> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
>>><br>
>>> Visit other Kitware open-source projects at<br>
>>> <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
>>><br>
>>> Please keep messages on-topic and check the ParaView Wiki at:<br>
>>> <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
>>><br>
>>> Follow this link to subscribe/unsubscribe:<br>
>>> <a href="http://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><br>
>>><br>
>><br>
>><br>
> _______________________________________________<br>
> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
><br>
> Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
><br>
> Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="http://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><br>
> _______________________________________________<br>
> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
><br>
> Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
><br>
> Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="http://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><br>
><br>
_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><br>
</div></div></blockquote></div><br>