summaryrefslogtreecommitdiff
path: root/lib/Numeric/LinearAlgebra/Algorithms.hs
diff options
context:
space:
mode:
Diffstat (limited to 'lib/Numeric/LinearAlgebra/Algorithms.hs')
-rw-r--r--lib/Numeric/LinearAlgebra/Algorithms.hs6
1 files changed, 3 insertions, 3 deletions
diff --git a/lib/Numeric/LinearAlgebra/Algorithms.hs b/lib/Numeric/LinearAlgebra/Algorithms.hs
index a3f541b..7223cd9 100644
--- a/lib/Numeric/LinearAlgebra/Algorithms.hs
+++ b/lib/Numeric/LinearAlgebra/Algorithms.hs
@@ -484,7 +484,7 @@ zh k v = fromList $ replicate (k-1) 0 ++ (1:drop k xs)
484 where xs = toList v 484 where xs = toList v
485 485
486zt 0 v = v 486zt 0 v = v
487zt k v = vjoin [subVector 0 (dim v - k) v, konst 0 k] 487zt k v = vjoin [subVector 0 (dim v - k) v, konst' 0 k]
488 488
489 489
490unpackQR :: (Field t) => (Matrix t, Vector t) -> (Matrix t, Matrix t) 490unpackQR :: (Field t) => (Matrix t, Vector t) -> (Matrix t, Matrix t)
@@ -640,10 +640,10 @@ luFact (l_u,perm) | r <= c = (l ,u ,p, s)
640 c = cols l_u 640 c = cols l_u
641 tu = triang r c 0 1 641 tu = triang r c 0 1
642 tl = triang r c 0 0 642 tl = triang r c 0 0
643 l = takeColumns r (l_u |*| tl) |+| diagRect 0 (konst 1 r) r r 643 l = takeColumns r (l_u |*| tl) |+| diagRect 0 (konst' 1 r) r r
644 u = l_u |*| tu 644 u = l_u |*| tu
645 (p,s) = fixPerm r perm 645 (p,s) = fixPerm r perm
646 l' = (l_u |*| tl) |+| diagRect 0 (konst 1 c) r c 646 l' = (l_u |*| tl) |+| diagRect 0 (konst' 1 c) r c
647 u' = takeRows c (l_u |*| tu) 647 u' = takeRows c (l_u |*| tu)
648 (|+|) = add 648 (|+|) = add
649 (|*|) = mul 649 (|*|) = mul