<div dir="ltr"><div><div>Can I get one of you to open an issue?<br><br><a href="https://github.com/yt-project/yt/issues/new">https://github.com/yt-project/yt/issues/new</a><br><br></div>If you can, please include a script to reproduce the problem, making use of one of the public datasets on <a href="http://yt-project.org/data">yt-project.org/data</a>.<br><br></div>-Nathan<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Nov 14, 2017 at 12:34 PM, Andrew James Emerick <span dir="ltr"><<a href="mailto:aje2123@columbia.edu" target="_blank">aje2123@columbia.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello Matias and Nathan,<div><br></div><div>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:</div><div><br></div><div><div style="font-size:12.8px">region = ds.region([0.5,0.5,0.5], [0.0, 0.0, 0.0], [1.0, 1.0, 1.0])</div><div style="font-size:12.8px">proj=ds.proj('density',axis='z<wbr>',weight_field='density',data_<wbr>source=region,method='integrat<wbr>e')</div></div><div><br></div><div>but this does not:</div><div><br></div><div><div style="font-size:12.8px">region = ds.region([0.5,0.5,0.5], [0.1,0.1,0.1], [0.9,0.9,0.9] )</div><div style="font-size:12.8px">proj =ds.proj('density',axis='z',<wbr>weight_field='density',data_<wbr>source=region,method='integrat<wbr>e')</div></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">In short, this indeed appears to be a bug, and possibly not specific to a single dataset.</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Andrew</div><div><br></div></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">On Tue, Nov 14, 2017 at 12:40 PM, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">This sounds like a bug. I don't think you should ever get back NaNs from a projection operation.<br><br>Is there any chance you can share the dataset that's triggering this behavior?<br><div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_8817100826249856841h5">On Tue, Nov 14, 2017 at 11:34 AM, Matias Enrique Suazo Campos <span dir="ltr"><<a href="mailto:matias.suazo@ug.uchile.cl" target="_blank">matias.suazo@ug.uchile.cl</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_8817100826249856841h5"><div dir="ltr">Dear YT users,<div><br></div><div>I'm trying to make a density projection, in the following lines I show you the lines that I have.</div><div><br></div><div>#Ramses output</div><div><div>ds = yt.load(path+"output_%05d"%sna<wbr>p+"/info_%05d.txt" %snap)</div><div>disk = ds.disk('max',[0,0,1],(20,'pc'<wbr>),(5,'pc'))</div><div>proj=ds.proj('density',axis='z<wbr>',weight_field='density',data_<wbr>source=disk,method='integrate'<wbr>)</div></div><div><br></div><div>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.</div><div><br></div><div>Thanks,</div><div>Matías Suazo</div><div>MSc student at Universidad de Chile</div></div>
<br></div></div>______________________________<wbr>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/lis<wbr>tinfo.cgi/yt-users-spacepope.o<wbr>rg</a><br>
<br></blockquote></div><br></div></div></div>
<br>______________________________<wbr>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/lis<wbr>tinfo.cgi/yt-users-spacepope.<wbr>org</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div></div></div><span class="HOEnZb"><font color="#888888">-- <br><div class="m_8817100826249856841gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div style="font-family:arial;font-size:small">NSF Graduate Fellow</div><div style="font-family:arial;font-size:small">Columbia University<br></div><div style="font-family:arial;font-size:small">Department of Astronomy</div></div></div></div></div>
</font></span></div>
<br>______________________________<wbr>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
<br></blockquote></div><br></div>