ldl factorization

2 files changed, 7 insertions, 1 deletions

diff --git a/packages/tests/src/Numeric/LinearAlgebra/Tests.hs b/packages/tests/src/Numeric/LinearAlgebra/Tests.hs
index d9bc9a0..2ff1580 100644
--- a/packages/tests/src/Numeric/LinearAlgebra/Tests.hs
+++ b/packages/tests/src/Numeric/LinearAlgebra/Tests.hs
@@ -587,6 +587,9 @@ runTests n = do
     putStrLn "------ luSolve"
     test (linearSolveProp (luSolve.luPacked) . rSqWC)
     test (linearSolveProp (luSolve.luPacked) . cSqWC)
+    putStrLn "------ ldlSolve"
+    test (linearSolveProp (ldlSolve.ldlPacked) . rSymWC)
+    test (linearSolveProp (ldlSolve.ldlPacked) . cSymWC)
     putStrLn "------ cholSolve"
     test (linearSolveProp (cholSolve.chol) . rPosDef)
     test (linearSolveProp (cholSolve.chol) . cPosDef)
diff --git a/packages/tests/src/Numeric/LinearAlgebra/Tests/Instances.hs b/packages/tests/src/Numeric/LinearAlgebra/Tests/Instances.hs
index 904ae05..7c54535 100644
--- a/packages/tests/src/Numeric/LinearAlgebra/Tests/Instances.hs
+++ b/packages/tests/src/Numeric/LinearAlgebra/Tests/Instances.hs
@@ -16,7 +16,7 @@ module Numeric.LinearAlgebra.Tests.Instances(
     Rot(..),    rRot,cRot,
     Her(..),    rHer,cHer,
     WC(..),     rWC,cWC,
-    SqWC(..),   rSqWC, cSqWC,
+    SqWC(..),   rSqWC, cSqWC, rSymWC, cSymWC,
     PosDef(..), rPosDef, cPosDef,
     Consistent(..), rConsist, cConsist,
     RM,CM, rM,cM,
@@ -176,6 +176,9 @@ cWC (WC m) = m :: CM
 rSqWC (SqWC m) = m :: RM
 cSqWC (SqWC m) = m :: CM
 
+rSymWC (SqWC m) = m + tr m :: RM
+cSymWC (SqWC m) = m + tr m :: CM
+
 rPosDef (PosDef m) = m :: RM
 cPosDef (PosDef m) = m :: CM