The issue is resolved.<div><br></div><div>It turns out that PlotCollection uses the min/max values of the x,y axes as bin edges so the edge values may sometimes not get included correctly.</div><div><br></div><div>By using BinnedProfile2D and add_field with weight=None, I can get consistent answers as compared to using TotalQuantity by either:</div>
<div><br></div><div>1) expanding the min/max by a tiny fraction<br>2) switch end_collect=True (which places values outside the bins to the edges)</div><div><br></div><div>Thanks Matt, mystery solved.</div><div><br></div><div>
From</div><div>G.S.</div><div><br><div class="gmail_quote">On Wed, Mar 6, 2013 at 4:41 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">Hi Geoffrey<br>
<div class="im"><br>
On Wed, Mar 6, 2013 at 4:38 PM, Geoffrey So <<a href="mailto:gsiisg@gmail.com">gsiisg@gmail.com</a>> wrote:<br>
> I'm trying to use add_phase_object to get a 2D distribution and summing over<br>
> the y axis to see the total quantity in each x bins.  I am using the old<br>
> plot collection to do the phase distribution because it is parallelized, the<br>
> binned2D is serial I believe.<br>
<br>
</div>If you specify lazy_reader = True it will be run in parallel.<br>
<div class="im"><br>
><br>
> I'm checking the numbers by using TotalQuantity which sums all the cells in<br>
> the simulation, and I should get the same number when I sum over the color<br>
> axis of the phase plot (I'm assuming this is true, can someone also verify<br>
> the logic?), but I am getting slightly different number.  More curious, is<br>
> when the y axis is the same as the color axis, I get different numbers as<br>
> when the axis names are different.  I first saw the discrepancy using other<br>
> fields, but in this script I am demonstrating with fields that everyone<br>
> should have access to in an Enzo simulation.<br>
><br>
> <a href="http://paste.yt-project.org/show/3222/" target="_blank">http://paste.yt-project.org/show/3222/</a><br>
><br>
> Can someone run the script very quick and verify that they get different<br>
> numbers printed on screen?<br>
<br>
</div>I can't right now, but what if you specify a dtype in the call to<br>
sum() ?  For instance:<br>
<br>
.sum(dtype="float64")<br>
<br>
-Matt<br>
<div class="im"><br>
><br>
> Is there an easier way to do what I'm trying to do in parallel?<br>
><br>
> From<br>
> G.S.<br>
><br>
> PS: This issue is in both dev-yt I have on my laptop and on Gordon<br>
><br>
> my laptop uses:<br>
> changeset:   6566:6794b2798b5a<br>
> branch:      yt<br>
> tag:         tip<br>
> user:        Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>><br>
> date:        Wed Oct 24 15:58:34 2012 -0400<br>
> summary:     Swap out get_smallest_dx for 8192.  Should fix the buildbot.<br>
> Gordon I have:<br>
> changeset:   7124:4151cc9ff9b2<br>
> branch:      yt<br>
> tag:         tip<br>
> user:        Kacper Kowalik <<a href="mailto:xarthisius.kk@gmail.com">xarthisius.kk@gmail.com</a>><br>
> date:        Sat Feb 02 11:03:07 2013 +0100<br>
> files:       yt/utilities/logger.py<br>
> description:<br>
> Add `NullHandler` to `ytLogger` when stream logging is suppresed. Fixes #492<br>
><br>
</div>> _______________________________________________<br>
> yt-dev mailing list<br>
> <a href="mailto:yt-dev@lists.spacepope.org">yt-dev@lists.spacepope.org</a><br>
> <a href="http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org</a><br>
><br>
_______________________________________________<br>
yt-dev mailing list<br>
<a href="mailto:yt-dev@lists.spacepope.org">yt-dev@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org</a><br>
</blockquote></div><br></div>