[yt-users] Projection

[Nathan Goldbaum]

Can I get one of you to open an issue?

https://github.com/yt-project/yt/issues/new

If you can, please include a script to reproduce the problem, making use of
one of the public datasets on yt-project.org/data.

-Nathan

On Tue, Nov 14, 2017 at 12:34 PM, Andrew James Emerick <aje2123 at columbia.edu
> wrote:

> Hello Matias and Nathan,
>
> I was able to reproduce this error on my own Enzo data set(s) using
> similar code. I did a little more playing around and it looks like it only
> breaks when having a weight_field. weight_field = None seems to be O.K. .
> It also occurs for the "region" objects, in addition to the disk.
> Interestingly, as Matias pointed out, doing something like this works:
>
> region = ds.region([0.5,0.5,0.5], [0.0, 0.0, 0.0], [1.0, 1.0, 1.0])
> proj=ds.proj('density',axis='z',weight_field='density',data_
> source=region,method='integrate')
>
> but this does not:
>
> region = ds.region([0.5,0.5,0.5], [0.1,0.1,0.1], [0.9,0.9,0.9] )
> proj =ds.proj('density',axis='z',weight_field='density',data_
> source=region,method='integrate')
>
> In short, this indeed appears to be a bug, and possibly not specific to a
> single dataset.
>
>
>
>
> Andrew
>
>
> On Tue, Nov 14, 2017 at 12:40 PM, Nathan Goldbaum <nathan12343 at gmail.com>
> wrote:
>
>> This sounds like a bug. I don't think you should ever get back NaNs from
>> a projection operation.
>>
>> Is there any chance you can share the dataset that's triggering this
>> behavior?
>>
>> On Tue, Nov 14, 2017 at 11:34 AM, Matias Enrique Suazo Campos <
>> matias.suazo at ug.uchile.cl> wrote:
>>
>>> Dear YT users,
>>>
>>> I'm trying to make a density projection, in the following lines I show
>>> you the lines that I have.
>>>
>>> #Ramses output
>>> ds = yt.load(path+"output_%05d"%snap+"/info_%05d.txt" %snap)
>>> disk = ds.disk('max',[0,0,1],(20,'pc'),(5,'pc'))
>>> proj=ds.proj('density',axis='z',weight_field='density',data_
>>> source=disk,method='integrate')
>>>
>>> However when I print proj["density"] I get a list of nan, however when I
>>> use the whole domain as my data_source I get reasonable numbers.
>>>
>>> Thanks,
>>> Matías Suazo
>>> MSc student at Universidad de Chile
>>>
>>> _______________________________________________
>>> yt-users mailing list
>>> yt-users at lists.spacepope.org
>>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>>
>>>
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
>
>
> --
> NSF Graduate Fellow
> Columbia University
> Department of Astronomy
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20171114/0fa7907f/attachment-0001.html>