[Numpy-discussion] NumPy/SciPy LAPACK version
Mon Jul 16 11:41:00 CDT 2007
Kevin Jacobs <firstname.lastname@example.org> 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