diff options
Diffstat (limited to 'lib/Numeric/LinearAlgebra/LAPACK.hs')
-rw-r--r-- | lib/Numeric/LinearAlgebra/LAPACK.hs | 8 |
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 | ||
262 | fixeig1 s = toComplex (subVector 0 r (asReal s), subVector r r (asReal s)) | 262 | fixeig1 s = toComplex' (subVector 0 r (asReal s), subVector r r (asReal s)) |
263 | where r = dim s | 263 | where r = dim s |
264 | 264 | ||
265 | fixeig [] _ = [] | 265 | fixeig [] _ = [] |
266 | fixeig [_] [v] = [comp v] | 266 | fixeig [_] [v] = [comp' v] |
267 | fixeig ((r1:+i1):(r2:+i2):r) (v1:v2:vs) | 267 | fixeig ((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 |
271 | fixeig _ _ = error "fixeig with impossible inputs" | 271 | fixeig _ _ = error "fixeig with impossible inputs" |
272 | 272 | ||