summaryrefslogtreecommitdiff
path: root/lib/Numeric/LinearAlgebra/LAPACK.hs
diff options
context:
space:
mode:
Diffstat (limited to 'lib/Numeric/LinearAlgebra/LAPACK.hs')
-rw-r--r--lib/Numeric/LinearAlgebra/LAPACK.hs8
1 files changed, 4 insertions, 4 deletions
diff --git a/lib/Numeric/LinearAlgebra/LAPACK.hs b/lib/Numeric/LinearAlgebra/LAPACK.hs
index cb48571..fbc5460 100644
--- a/lib/Numeric/LinearAlgebra/LAPACK.hs
+++ b/lib/Numeric/LinearAlgebra/LAPACK.hs
@@ -259,14 +259,14 @@ eigRaux m = unsafePerformIO $ do
259 where r = rows m 259 where r = rows m
260 g ra ca pa = dgeev ra ca pa 0 0 nullPtr 260 g ra ca pa = dgeev ra ca pa 0 0 nullPtr
261 261
262fixeig1 s = toComplex (subVector 0 r (asReal s), subVector r r (asReal s)) 262fixeig1 s = toComplex' (subVector 0 r (asReal s), subVector r r (asReal s))
263 where r = dim s 263 where r = dim s
264 264
265fixeig [] _ = [] 265fixeig [] _ = []
266fixeig [_] [v] = [comp v] 266fixeig [_] [v] = [comp' v]
267fixeig ((r1:+i1):(r2:+i2):r) (v1:v2:vs) 267fixeig ((r1:+i1):(r2:+i2):r) (v1:v2:vs)
268 | r1 == r2 && i1 == (-i2) = toComplex (v1,v2) : toComplex (v1,scale (-1) v2) : fixeig r vs 268 | r1 == r2 && i1 == (-i2) = toComplex' (v1,v2) : toComplex' (v1,scale (-1) v2) : fixeig r vs
269 | otherwise = comp v1 : fixeig ((r2:+i2):r) (v2:vs) 269 | otherwise = comp' v1 : fixeig ((r2:+i2):r) (v2:vs)
270 where scale = vectorMapValR Scale 270 where scale = vectorMapValR Scale
271fixeig _ _ = error "fixeig with impossible inputs" 271fixeig _ _ = error "fixeig with impossible inputs"
272 272