[Numpy-discussion] NumPy/SciPy LAPACK version
robert.kern at gmail.com
Mon Jul 16 12:41:00 EDT 2007
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.
"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
More information about the NumPy-Discussion