summaryrefslogtreecommitdiff
path: root/packages/base/src/Internal/LAPACK.hs
diff options
context:
space:
mode:
authorAlberto Ruiz <aruiz@um.es>2015-06-29 11:43:48 +0200
committerAlberto Ruiz <aruiz@um.es>2015-06-29 11:43:48 +0200
commit18dad63040ded686187204f549d38100f62fe388 (patch)
treeec2bee42d7f5566370f289df2b5a248bbefe6e74 /packages/base/src/Internal/LAPACK.hs
parent982b24442018c14510ce9bcf4d0e402613fcbea2 (diff)
use lda in lus
Diffstat (limited to 'packages/base/src/Internal/LAPACK.hs')
-rw-r--r--packages/base/src/Internal/LAPACK.hs8
1 files changed, 4 insertions, 4 deletions
diff --git a/packages/base/src/Internal/LAPACK.hs b/packages/base/src/Internal/LAPACK.hs
index 049c11e..124e353 100644
--- a/packages/base/src/Internal/LAPACK.hs
+++ b/packages/base/src/Internal/LAPACK.hs
@@ -544,11 +544,11 @@ foreign import ccall unsafe "lu_l_C" zgetrf :: R :> C ::> Ok
544 544
545-- | LU factorization of a general real matrix, using LAPACK's /dgetrf/. 545-- | LU factorization of a general real matrix, using LAPACK's /dgetrf/.
546luR :: Matrix Double -> (Matrix Double, [Int]) 546luR :: Matrix Double -> (Matrix Double, [Int])
547luR = luAux dgetrf "luR" . fmat 547luR = luAux dgetrf "luR"
548 548
549-- | LU factorization of a general complex matrix, using LAPACK's /zgetrf/. 549-- | LU factorization of a general complex matrix, using LAPACK's /zgetrf/.
550luC :: Matrix (Complex Double) -> (Matrix (Complex Double), [Int]) 550luC :: Matrix (Complex Double) -> (Matrix (Complex Double), [Int])
551luC = luAux zgetrf "luC" . fmat 551luC = luAux zgetrf "luC"
552 552
553luAux f st a = unsafePerformIO $ do 553luAux f st a = unsafePerformIO $ do
554 lu <- copy ColumnMajor a 554 lu <- copy ColumnMajor a
@@ -566,11 +566,11 @@ foreign import ccall unsafe "luS_l_C" zgetrs :: C ::> R :> C ::> Ok
566 566
567-- | Solve a real linear system from a precomputed LU decomposition ('luR'), using LAPACK's /dgetrs/. 567-- | Solve a real linear system from a precomputed LU decomposition ('luR'), using LAPACK's /dgetrs/.
568lusR :: Matrix Double -> [Int] -> Matrix Double -> Matrix Double 568lusR :: Matrix Double -> [Int] -> Matrix Double -> Matrix Double
569lusR a piv b = lusAux dgetrs "lusR" (fmat a) piv (fmat b) 569lusR a piv b = lusAux dgetrs "lusR" (fmat a) piv b
570 570
571-- | Solve a real linear system from a precomputed LU decomposition ('luC'), using LAPACK's /zgetrs/. 571-- | Solve a real linear system from a precomputed LU decomposition ('luC'), using LAPACK's /zgetrs/.
572lusC :: Matrix (Complex Double) -> [Int] -> Matrix (Complex Double) -> Matrix (Complex Double) 572lusC :: Matrix (Complex Double) -> [Int] -> Matrix (Complex Double) -> Matrix (Complex Double)
573lusC a piv b = lusAux zgetrs "lusC" (fmat a) piv (fmat b) 573lusC a piv b = lusAux zgetrs "lusC" (fmat a) piv b
574 574
575lusAux f st a piv b 575lusAux f st a piv b
576 | n1==n2 && n2==n =unsafePerformIO $ do 576 | n1==n2 && n2==n =unsafePerformIO $ do