Dear Nathan,
Many thanks for your reply.
Actually I took the simplest case in which this can be done, and usually, I
did it in this way.
I was thinking that the provided test script may lead to an alternative
solution. But I actually need to pass the value of the field to the
function for the current task.
I try to illustrate my task, actually, I have two distribution functions.
One is Maxwell-Juttner (MJ) distribution and power-law distribution in
momentum space. I have to find the intersection point (say p_int) between
these two distributions on each cell. Now the thing is that MJ distribution
needs the value of Temperature and density field of the cell and power-law
distribution needs the value of density and Mach number field of the same
cell to compute the intersection point of the said distribution functions
on that cell which I will select based on some criteria.
Hope this will illustrate you what exactly I am looking for.
Many thanks and regards,
-Prateek
On Fri, Nov 1, 2019 at 1:44 AM Nathan
Hi Prateek,
Instead of using your "check" function and trying to set values one-by-one, you can use numpy's fancy indexing to do the conditional setting of values in your _norm function. I've modified your script to do what I think you were trying to do here: https://gist.github.com/ngoldbaum/b3ed70139faf6da89f37ae56cb965faf
This script produces this image: https://i.imgur.com/6Fgpf53.png
Hope that helps!
If there's something else that you were trying to do that your script wasn't illustrating and my solution won't work for, please do write back with additional details.
-Nathan
On Thu, Oct 31, 2019 at 4:05 AM Prateek Gupta
wrote: Dear Yt-users,
I need to pass the value of the available/derived field corresponds to the index selected based on some criterion (using numpy.where()) to a scalar function. It should list out the indexes wheres the condition satisfies. It happens when this has been done outside the field but instead of getting 3D indexes it gives 1D YT_Array indexes I think. But I need to pass this value of field corresponds to the selected index to scalar function inside the field, which may give some problem (please find the attachment to see what exactly I am trying to do).
The attached script can run on the ENZO test data set, enzo_tiny_cosmology.
Many thanks and regards -Prateek
_______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org
_______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org