[Numpy-discussion] NumPy/SciPy LAPACK version

[Robert Kern]

Kevin Jacobs <jacobs at bioinformed.com> wrote:
> Mea culpa on the msqrt example, however I still think it is wrong to get
> a complex square-root back when a real valued result is expected and exists.

No, in floating point you accumulate error. Those 1e-22j's are part of the
actual result. Some systems like MATLAB implicitly silent such small imaginary
components; we don't.

-- 
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless enigma
 that is made terrible by our own mad attempt to interpret it as though it had
 an underlying truth."
  -- Umberto Eco