<div dir="ltr"><br><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Jun 12, 2013 at 3:25 PM, Moroney, Catherine M (398D) <span dir="ltr"><<a href="mailto:Catherine.M.Moroney@jpl.nasa.gov" target="_blank">Catherine.M.Moroney@jpl.nasa.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello,<br>
<br>
I've got two arrays of the same shape that I read in from a file, and I'm trying to<br>
difference them.  Very simple stuff, but I'm getting weird answers.<br>
<br>
Here is the code:<br>
<br>
>>> counts1 = hfile1.read_grid_field("CFbA", "TerrainReferencedRCCMFraction_Num")<br>
>>> counts2 = hfile2.read_grid_field("CFbA", "TerrainReferencedRCCMFraction_Num")<br>
>>> counts1.max(), counts2.max()<br>
(13, 13)<br>
>>> counts1.min(), counts2.min()<br>
(0, 0)<br>
>>> numpy.all(counts1 == counts2)<br>
False<br>
>>> diff = counts1 - counts2<br>
>>> diff.max()<br>
4294967295      !! WHAT IS HAPPENING HERE ??<br>
>>> sum = counts1 + counts2<br>
>>> sum.max()<br>
26<br>
<br>
As you can see, the range of values in both arrays is 0 to 13, and the sum<br>
behaves normally, but the difference gives this weird number.<br>
<br>
When I create dummy arrays, the subtraction works fine.  So there must be some funny value<br>
lurking in either the counts1 or counts2 array, but the numpy.isnan() test returns False.<br>
<br>
Any ideas for how I debug this?<br>
<br>
Catherine<br>
<br></blockquote><div><br></div><div>Check the dtype of the arrays.  They are probably unsigned integers, and the subtraction leads to wrap-around in some cases.<br><br></div><div>For example:<br><span style="font-family:courier new,monospace"><br>
In [1]: x = np.array([0, 1, 2], dtype=np.uint32)<br><br>In [2]: y = np.array([1, 1, 1], dtype=np.uint32)<br><br>In [3]: x - y<br>Out[3]: array([4294967295,          0,          1], dtype=uint32)</span><br><br><br></div><div>
Warren<br><br><br></div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
_______________________________________________<br>
NumPy-Discussion mailing list<br>
<a href="mailto:NumPy-Discussion@scipy.org">NumPy-Discussion@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/numpy-discussion" target="_blank">http://mail.scipy.org/mailman/listinfo/numpy-discussion</a><br>
</blockquote></div><br></div></div>