[Numpy-discussion] why does eigvalsh return a complex array?
Charles R Harris
Fri Apr 6 21:04:02 CDT 2012
On Fri, Apr 6, 2012 at 1:56 AM, Christoph Groth <email@example.com> wrote:
> I noticed that numpy.linalg.eigvalsh returns a complex array, even
> though mathematically the resulting eigenvalues are guaranteed to be
> Looking at the source code, the underlying zheevd routine of LAPACK
> indeed returns an array of real numbers which is than converted to
> complex in the numpy wrapper.
> Does numpy policy require the type of the result to be the same as the
> type of input? Copying an array twice to arrive at the original result
> seems pointless to me.
I think this should be fixed, the problem is the wrapper from a complex
array. Not sure what the easiest fix is. I expect eigh is similar in
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the NumPy-Discussion