From 1e564ccf064a3486532aeb32d689e5ca4fd2e260 Mon Sep 17 00:00:00 2001 From: Alberto Ruiz Date: Fri, 20 Jun 2014 08:58:19 +0200 Subject: don't compute eigenvectors in eigOnlyH! --- packages/base/src/Numeric/LinearAlgebra/LAPACK.hs | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'packages/base/src/Numeric/LinearAlgebra/LAPACK.hs') diff --git a/packages/base/src/Numeric/LinearAlgebra/LAPACK.hs b/packages/base/src/Numeric/LinearAlgebra/LAPACK.hs index b32b67f..e088fdc 100644 --- a/packages/base/src/Numeric/LinearAlgebra/LAPACK.hs +++ b/packages/base/src/Numeric/LinearAlgebra/LAPACK.hs @@ -320,7 +320,7 @@ eigOnlyS = vrev . fst. eigSHAux (dsyev 0) "eigS'" . fmat -- | Eigenvalues of a hermitian complex matrix, using LAPACK's /zheev/ with jobz == \'N\'. -- The eigenvalues are sorted in descending order. eigOnlyH :: Matrix (Complex Double) -> Vector Double -eigOnlyH = vrev . fst. eigSHAux (zheev 1) "eigH'" . fmat +eigOnlyH = vrev . fst. eigSHAux (zheev 0) "eigH'" . fmat vrev = flatten . flipud . reshape 1 -- cgit v1.2.3