<div dir="ltr">Hi Matthew,<div><br></div><div>First of all, thank you for the detailed report. That made it much easier to track down the issue. Just in case you run into issues in the future and want to share larger data files, you can use the yt curldrop for that:</div><div><br></div><div><a href="https://docs.hub.yt/services.html#curldrop">https://docs.hub.yt/services.html#curldrop</a><br></div><div><br></div><div>As Ben mentioned, your issue is eerily reminiscent of Issue 1084:</div><div><br></div><div><a href="https://bitbucket.org/yt_analysis/yt/issues/1084">https://bitbucket.org/yt_analysis/yt/issues/1084</a><br></div><div><br></div><div>This was fixed last year by Ben Keller in PR 1735:</div><div><br></div><div><a href="https://bitbucket.org/yt_analysis/yt/pull-requests/1735">https://bitbucket.org/yt_analysis/yt/pull-requests/1735</a><br></div><div><br></div><div>However, it turns out that the implementation in PR 1735 had a bug that your test script revealed. I think there's a one-liner fix:</div><div><br></div><div><a href="https://bitbucket.org/yt_analysis/yt/pull-requests/2430">https://bitbucket.org/yt_analysis/yt/pull-requests/2430</a><br></div><div><br></div><div>If you want to test out that fix, that would be very much appreciated. The easiest way to do so is to modify your current yt installation by hand (you can find it's location by looking at yt.__file__). That said, if you have time I'd recommend compiling yt from source:</div><div><br></div><div>$ # get rid of old yt installation via conda or pip</div><div>$ pip uninstall yt</div><div>$ conda uninstall yt</div><div><br></div><div>$ hg clone <a href="https://bitbucket.org/yt_analysis/yt">https://bitbucket.org/yt_analysis/yt</a> ./yt-hg</div><div>$ cd yt-hg</div><div><br></div><div>$ # pull bugfix from Nathan's fork of yt and update to bugfix changeset</div><div>$ hg pull -u -r 1a290fab65e4 <a href="https://bitbucket.org/ngoldbaum/yt">https://bitbucket.org/ngoldbaum/yt</a></div><div><br></div><div>$ # compile yt</div><div>$ python setup.py develop</div><div><br></div><div>This will build yt from source. It may take a few minutes because it will need to compile yt's C extensions. Once that's done, try rerunning your test script.</div><div><br></div><div>-Nathan</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 31, 2016 at 10:27 AM, B.W. Keller <span dir="ltr"><<a href="mailto:kellerbw@mcmaster.ca" target="_blank">kellerbw@mcmaster.ca</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>This sounds similar to a bug that I (thought)  I had fixed a while back.  I'll take a look at your issue, see if I can reproduce it and put in a fix.</div><div><br></div><div>Ben</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Mon, Oct 31, 2016 at 11:06 AM, Matthew Robert Brinson <span dir="ltr"><<a href="mailto:mbrinson@uwm.edu" target="_blank">mbrinson@uwm.edu</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="h5">Hi yt-users,<br>
<br>
Let’s try this again. The file I sent the first time was too large.<br>
<br>
I’m using yt version 3.3.1, and I’ve recently encountered an issue loading data correctly. I’m running SPH simulations of white dwarf mergers and want to know what the mass fractions of different elements are at different times. I have two files that contain the data: 1) a tipsy file, 2) and a text file that contains the mass fraction values of helium-4. The tipsy file loads just fine and I’m able to access the information from the helium file once the tipsy file is loaded. For some reason though, when I have more than 262144 particles in my simulation, yt does not load the helium mass fraction values correctly anymore. I’ve initialized all particles to have a 4% helium mass fraction, and all of the particles up to and including the particle indexed by 262143 have this value, but after that the particles all have 0% according to yt. I’ve checked the helium files directly and all values are at 4%. I’m confused as to why yt would stop loading values correctly after a certain point.<br>
<br>
I’ve included two tipsy files as well as their corresponding helium mass fraction files. I have also written a script called visualization_test_script.py that makes scatter plots of the white dwarfs showing particles that have correct and incorrect values of helium. I can send more tipsy files if necessary, but they are too large to send any more here.<br>
<br>
Thanks,<br>
Matt<br>
<br>
<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.<wbr>org</a><br>
<br></blockquote></div><br></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>