diff options
Diffstat (limited to 'packages/base')
-rw-r--r-- | packages/base/src/Internal/LAPACK.hs | 4 |
1 files changed, 3 insertions, 1 deletions
diff --git a/packages/base/src/Internal/LAPACK.hs b/packages/base/src/Internal/LAPACK.hs index ff55688..27d1f95 100644 --- a/packages/base/src/Internal/LAPACK.hs +++ b/packages/base/src/Internal/LAPACK.hs | |||
@@ -506,8 +506,10 @@ linearSolveGTAux2 g f st dl d du b | |||
506 | | ndl == nd - 1 && | 506 | | ndl == nd - 1 && |
507 | ndu == nd - 1 && | 507 | ndu == nd - 1 && |
508 | nd == r = unsafePerformIO . g $ do | 508 | nd == r = unsafePerformIO . g $ do |
509 | dl' <- head . toRows <$> copy ColumnMajor (fromRows [dl]) | ||
510 | du' <- head . toRows <$> copy ColumnMajor (fromRows [du]) | ||
509 | s <- copy ColumnMajor b | 511 | s <- copy ColumnMajor b |
510 | (dl # d # du #! s) f #| st | 512 | (dl' # d # du' #! s) f #| st |
511 | return s | 513 | return s |
512 | | otherwise = error $ st ++ " of nonsquare matrix" | 514 | | otherwise = error $ st ++ " of nonsquare matrix" |
513 | where | 515 | where |