I'm pretty certain now that this is a numeric precision problem. The discrepancy only shows up for very small numbers. For now I'll ignore it. Rob. Rob wrote:

I found a slight discrepancy in my Python EM code output when I use the invert() function. It may be that the C routine that I am comparing it to actually has the bug. For now I don't know. It would be cool to have a large two dimensional array that has a known inverse. I think NIST keeps an array collection, but I don't know if they have the inverses. Rob.

-- The Numeric Python EM Project

www.members.home.net/europax

_______________________________________________ Numpy-discussion mailing list Numpy-discussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpy-discussion

-- The Numeric Python EM Project www.members.home.net/europax