<div dir="ltr">Do you guys think this is related to the problem we encountered a few months ago with SPH datasets where regions of "zero" density were just filtered out of the LightRay to avoid NaNs and such? It seemed like our fix was too kludgy, but perhaps this is unrelated.</div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Feb 25, 2016 at 6:39 AM, Britton Smith <span dir="ltr"><<a href="mailto:brittonsmith@gmail.com" target="_blank">brittonsmith@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">Hi Bili,<div><br></div><div>I have looked into this and it seems there are some problems with generating rays with particle datasets. I cannot tell whether yours is a special case, but what seems to be happening is that certain fields are being returned with different sizes. The error happens in the LightRay code because it first tries to access the "t" field, which comes back with size 40. After that, it tries to access one of the coordinate fields ("x", "y", or "z") and that is being returned with size 81, which doesn't fit into the arrays made to hold it that already have size 40.</div><div><br></div><div>Anyway, you're not seeing this when you make the ray yourself if you first try to access a field other than the "t" field. This results in the field arrays being made with the larger size. If you do this and then try to access the "t" array, you will find that the first 40 values have real numbers but that the rest are just uninitialized values.</div><div><br></div><div>Unfortunately, I haven't gotten much past this as I'm having trouble seeing where the fcoords for this particular object are getting generated. I have created an issue, which you can follow here:</div><div><a href="https://bitbucket.org/yt_analysis/yt/issues/1175/fcoords-and-tcoords-not-the-same-shape-for" target="_blank">https://bitbucket.org/yt_analysis/yt/issues/1175/fcoords-and-tcoords-not-the-same-shape-for</a><br></div><div><br></div><div>That's about all I can do for now, but if anyone has any thoughts, help would be appreciated.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Britton</div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 24, 2016 at 7:21 PM, Bili Dong - Gmail <span dir="ltr"><<a href="mailto:qobilidop@gmail.com" target="_blank">qobilidop@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">Hi Britton,<div><br></div><div>I was using the yt sample dataset 'GadgetDiskGalaxy/snapshot_200.hdf5' to test the functionality.</div><div><br></div><div>Thanks,</div><div>Bili</div></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Feb 24, 2016 at 2:39 AM, Britton Smith <span dir="ltr"><<a href="mailto:brittonsmith@gmail.com" target="_blank">brittonsmith@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">Hi Bili,<div><br></div><div>My apologies for the long delay in replying to this thread. By now, I think you may be working with Trident. If things are working, that's great. However, Trident relies on much of the machinery that you were trying to use, so it's still possible that you see this error. The problem is likely related to how the LightRay generator is calculating the trajectory for the ray and not with the underlying ray machinery which it relies on. If it's possible for you to give me access to this dataset, I would be happy to take a look at this.</div><span><font color="#888888"><div><br></div><div>Britton</div></font></span></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 15, 2016 at 9:30 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">I think Cameron, Britton, and Devin have done this for their Trident project - you might want to contact them.</div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 15, 2016 at 3:23 PM, Bili Dong - Gmail <span dir="ltr"><<a href="mailto:qobilidop@gmail.com" target="_blank">qobilidop@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">Hi Matt,<div><br></div><div>I've tried to construct the ray and ortho_ray objects, neither of which can reproduce the error. They work just fine. I haven't succeeded in tracing down the error yet, but I'll keep working on that.</div><div><br></div><div>Besides, I'm wondering if there are examples of successfully making light rays in a particle dataset? That may help figuring out what might go wrong.</div><div><br></div><div>Thanks,</div><div>Bili</div></div><div><div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Feb 12, 2016 at 12:42 PM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com" target="_blank">matthewturk@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">Hi Bili,<div><br></div><div>Hm, I'm honestly not sure. This looks like an oddity in the way the octree is being constructed in spatial dimensions. Any chance you can reproduce this with a ray or ortho_ray object, which should simplify the number of moving parts?</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Tue, Feb 9, 2016 at 2:04 PM, Bili Dong - Gmail <span dir="ltr"><<a href="mailto:qobilidop@gmail.com" target="_blank">qobilidop@gmail.com</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><div dir="ltr">Hi yt users,<div><br></div><div>I'm trying to make a light ray for a gadget dataset using this script: <a href="http://paste.yt-project.org/show/6234/" target="_blank">http://paste.yt-project.org/show/6234/</a></div><div><br></div><div>But I get the following error:</div><div><br></div><div>```</div><div><pre style="overflow:auto;font-size:14px;padding:0px;margin-top:0px;margin-bottom:0px;line-height:inherit;word-wrap:break-word;color:rgb(0,0,0);background-color:transparent;border:0px;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline"><span style="color:rgb(70,130,180)"><span style="color:rgb(0,100,0)"><span style="color:rgb(70,130,180)"><span style="color:rgb(0,100,0)"><span style="color:rgb(70,130,180)"><span style="color:rgb(0,100,0)"><span style="color:rgb(70,130,180)"><span style="color:rgb(0,100,0)"><span style="color:rgb(70,130,180)"><span style="color:rgb(0,100,0)"><span style="color:rgb(70,130,180)"><span style="color:rgb(0,100,0)"><span style="color:rgb(70,130,180)"><span style="color:rgb(0,100,0)"><span style="color:rgb(70,130,180)"><span style="color:rgb(0,100,0)"><span style="color:rgb(70,130,180)"><span style="color:rgb(0,100,0)"><span style="color:rgb(70,130,180)"><span style="color:rgb(0,100,0)"><span style="color:rgb(139,0,0)"><span style="color:rgb(70,130,180)"><span style="color:rgb(0,100,0)"><span>/Users/qobilidop/repos/yt/yt/geometry/geometry_handler.py</span> in <span style="color:rgb(70,130,180)">fcoords<span style="color:rgb(0,0,139)">(self)</span>
<span style="color:rgb(0,100,0)"> 320</span> c <span style="color:rgb(0,0,139)">=</span> obj<span style="color:rgb(0,0,139)">.</span>select_fcoords<span style="color:rgb(0,0,139)">(</span>self<span style="color:rgb(0,0,139)">.</span>dobj<span style="color:rgb(0,0,139)">)</span><span style="color:rgb(0,0,139)"></span>
<span style="color:rgb(0,100,0)"> 321</span> <span style="color:rgb(0,100,0)">if</span> c<span style="color:rgb(0,0,139)">.</span>shape<span style="color:rgb(0,0,139)">[</span><span>0</span><span style="color:rgb(0,0,139)">]</span> <span style="color:rgb(0,0,139)">==</span> <span>0</span><span style="color:rgb(0,0,139)">:</span> <span style="color:rgb(0,100,0)">continue</span><span style="color:rgb(0,0,139)"></span>
<span style="color:rgb(0,100,0)">--> 322<span style="color:rgb(139,0,0)"> </span>ci<span style="color:rgb(0,0,139)">[</span>ind<span style="color:rgb(0,0,139)">:</span>ind<span style="color:rgb(0,0,139)">+</span>c<span style="color:rgb(0,0,139)">.</span>shape<span style="color:rgb(0,0,139)">[</span><span style="color:rgb(70,130,180)">0</span><span style="color:rgb(0,0,139)">]</span><span style="color:rgb(0,0,139)">,</span> <span style="color:rgb(0,0,139)">:</span><span style="color:rgb(0,0,139)">]</span> <span style="color:rgb(0,0,139)">=</span> c<span style="color:rgb(0,0,139)"></span>
<span> 323</span> ind <span style="color:rgb(0,0,139)">+=</span> c<span style="color:rgb(0,0,139)">.</span>shape<span style="color:rgb(0,0,139)">[</span><span style="color:rgb(70,130,180)">0</span><span style="color:rgb(0,0,139)">]</span><span style="color:rgb(0,0,139)"></span>
<span> 324</span> <span>return</span> ci<span style="color:rgb(0,0,139)"></span>
<span style="color:rgb(139,0,0)">ValueError</span>: could not broadcast input array from shape (81,3) into shape (40,3)</span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></span></pre><pre style="overflow:auto;font-size:14px;padding:0px;margin-top:0px;margin-bottom:0px;line-height:inherit;word-wrap:break-word;color:rgb(0,0,0);background-color:transparent;border:0px;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;white-space:normal">```</span><br></pre><pre style="overflow:auto;font-size:14px;padding:0px;margin-top:0px;margin-bottom:0px;line-height:inherit;word-wrap:break-word;color:rgb(0,0,0);background-color:transparent;border:0px;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline"><br></pre><pre style="overflow:auto;font-size:14px;padding:0px;margin-top:0px;margin-bottom:0px;line-height:inherit;word-wrap:break-word;color:rgb(0,0,0);background-color:transparent;border:0px;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;white-space:pre-wrap;vertical-align:baseline">M<font face="arial, sans-serif" style="color:rgb(34,34,34)"><span style="white-space:normal">y yt is on the </span><span style="white-space:normal">changeset </span></font><span style="color:rgb(34,34,34);font-family:arial,sans-serif">d47150dfbde6.</span></pre><pre style="overflow:auto;padding:0px;margin-top:0px;margin-bottom:0px;line-height:inherit;word-wrap:break-word;background-color:transparent;border:0px;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;vertical-align:baseline"><font face="arial, sans-serif"><span style="white-space:normal"><br></span></font></pre><pre style="overflow:auto;padding:0px;margin-top:0px;margin-bottom:0px;line-height:inherit;word-wrap:break-word;background-color:transparent;border:0px;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;vertical-align:baseline"><font face="arial, sans-serif"><span style="white-space:normal">I'm wondering if I'm doing it the correct way? And how could I solve the problem?</span></font></pre><pre style="overflow:auto;padding:0px;margin-top:0px;margin-bottom:0px;line-height:inherit;word-wrap:break-word;background-color:transparent;border:0px;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;vertical-align:baseline"><span style="font-family:arial,sans-serif"><br></span></pre><pre style="overflow:auto;padding:0px;margin-top:0px;margin-bottom:0px;line-height:inherit;word-wrap:break-word;background-color:transparent;border:0px;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;vertical-align:baseline"><span style="font-family:arial,sans-serif">Thanks,</span></pre><pre style="overflow:auto;padding:0px;margin-top:0px;margin-bottom:0px;line-height:inherit;word-wrap:break-word;background-color:transparent;border:0px;border-top-left-radius:0px;border-top-right-radius:0px;border-bottom-right-radius:0px;border-bottom-left-radius:0px;vertical-align:baseline"><span style="font-family:arial,sans-serif">Bili</span></pre></div></div>
<br></div></div>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>
<br>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>
</div></div><br>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>
</div></div><br>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>
</div></div><br>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>
</div></div><br>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>
</div></div><br>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div>Cameron Hummels<div><span style="font-size:12.8000001907349px">NSF Postdoctoral Fellow</span></div><div><span style="font-size:12.8000001907349px">Department of Astronomy</span></div><div><span style="font-size:12.8000001907349px">California Institute of Technology</span><br></div><div><a href="http://chummels.org" style="font-size:12.8000001907349px" target="_blank">http://chummels.org</a><br></div></div></div></div>
</div>