[Numpy-discussion] Apparently non-deterministic behaviour of complex array multiplication
kneil
magnetotellurics at gmail.com
Fri Dec 2 22:10:51 EST 2011
OK - here is the status update on this problem:
1. To test for bad RAM, I saved the code plus data onto a usb drive and
manually transferred it to a colleague's desktop machine with same specs as
mine. The NaN values continued to appear at random, so it is unlikely to be
bad RAM - unless its contagious.
2. Here is how I am temporarily working around the problem: Right before
performing the multiplication of X*X.H where X=asmatrix(Y), I save X to file
using np.save('X',X). Then I reload it via X=np.load('X.npy'). Cast as
matrix: X=asmatrix(X); and carry on: S=X*X.H. I have not seen any NaN
although I have only run it a few times... but it seems to work.
I'd be very curious to hear any ideas on why this problem exists in the
first place, and why save/load gets me around it.
For the record I am running ubuntu 11.04, have 16GB RAM (not 12) and use
python2.7.
cheers,
Karl
Hi Pierre,
I confirmed with the guy who put together the machine that it is non-ECC
RAM. You know, now that i think about it, this machine seems to crash a
fair amount more often than its identical twin which sits on a desk near me.
I researched memtest a bit... downloaded and compiled it, but I do not quite
understand the finer points of using it... it seems that I want to remove my
RAM cards and test them one at a time. Do you know a good reference for
using it.
I think at this point the best thing to do will be to dump my data/code to
portable HDD and load it on the other computer with same specs as this one.
If it runs without generating any NaN then I will proceed to a full memtest.
Thanks for the advice.
-Karl
On Thu, Dec 1, 2011 at 2:47 PM, kneil <magnetotellurics at gmail.com> wrote:
Are you using non-ECC RAM, by chance? (Though if you have >4GB of ram, I
can't imagine that you wouldn't be using ECC...)
Alternately, have you run memtest lately? That sound suspiciously like bad
ram...
--
View this message in context: http://old.nabble.com/Apparently-non-deterministic-behaviour-of-complex-array-multiplication-tp32893004p32906553.html
Sent from the Numpy-discussion mailing list archive at Nabble.com.
More information about the NumPy-Discussion
mailing list