3 files changed, 36 insertions, 13 deletions

diff --git a/examples/tests.hs b/examples/experiments/bigtests.hs
index cd923cd..a3dfad2 100644
--- a/examples/tests.hs
+++ b/examples/experiments/bigtests.hs
@@ -11,15 +11,41 @@ pseudorandomR seed (n,m) = reshape m $ fromList $ take (n*m) $ randomRs (-100,10
 
 pseudorandomC seed (n,m) = toComplex (pseudorandomR seed (n,m), pseudorandomR (seed+1) (n,m))
 
-bigmat = m + trans m :: RM
-    where m = pseudorandomR 18 (1000,1000)
-bigmatc = mc + ctrans mc ::CM
-    where mc = pseudorandomC 19 (1000,1000)
+bigmat = m + trans m
+    where m = pseudorandomR 18 (1000,1000) :: Matrix Double
+bigmatc = mc + ctrans mc
+    where mc = pseudorandomC 19 (1000,1000) :: Matrix (Complex Double)
 
 utest str b = TestCase $ assertBool str b
 
 feye n = flipud (ident n) :: Matrix Double
 
+infixl 4 |~|
+a |~| b = dist a b < 10^^(-10)
+
+dist a b = r
+    where norm = pnorm Infinity
+          na = norm a
+          nb = norm b
+          nab = norm (a-b)
+          mx = max na nb
+          mn = min na nb
+          r = if mn < eps
+                then mx
+                else nab/mx
+
+square m = rows m == cols m
+
+unitary m = square m && m <> ctrans m |~| ident (rows m)
+
+eigProp m = complex m <> v |~| v <> diag s
+    where (s, v) = eig m
+
+eigSHProp m = m <> v |~| v <> real (diag s)
+              && unitary v
+              && m |~| v <> real (diag s) <> ctrans v
+    where (s, v) = eigSH m
+
 bigtests = do
     putStrLn "--------- big matrices -----"
     runTestTT $ TestList
@@ -31,8 +57,4 @@ bigtests = do
      ]
     return ()
 
-main = do
-    args <- getArgs
-    if "--big" `elem` args
-        then bigtests
-        else runTests 20
+main = bigtests
diff --git a/lib/Numeric/LinearAlgebra/Tests.hs b/lib/Numeric/LinearAlgebra/Tests.hs
index 96280fb..31ba2ff 100644
--- a/lib/Numeric/LinearAlgebra/Tests.hs
+++ b/lib/Numeric/LinearAlgebra/Tests.hs
@@ -13,8 +13,8 @@ Some tests.
 -}
 
 module Numeric.LinearAlgebra.Tests(
-  module Numeric.LinearAlgebra.Tests.Instances,
-  module Numeric.LinearAlgebra.Tests.Properties,
+--  module Numeric.LinearAlgebra.Tests.Instances,
+--  module Numeric.LinearAlgebra.Tests.Properties,
   qCheck, runTests
 --, runBigTests
 ) where
@@ -116,7 +116,8 @@ rotTest = fun (10^5) :~12~: rot 5E4
               where angles = toList $ linspace n (0,1)
 
 
--- | It runs all the tests.
+-- | All tests must pass with a maximum dimension of about 20
+--  (some tests may fail with bigger sizes due to precision loss).
 runTests :: Int  -- ^ maximum dimension
          -> IO ()
 runTests n = do
diff --git a/lib/Numeric/LinearAlgebra/Tests/Properties.hs b/lib/Numeric/LinearAlgebra/Tests/Properties.hs
index 0563e62..566d038 100644
--- a/lib/Numeric/LinearAlgebra/Tests/Properties.hs
+++ b/lib/Numeric/LinearAlgebra/Tests/Properties.hs
@@ -144,6 +144,6 @@ cholProp m = m |~| ctrans c <> c && upperTriang c
     where c = chol m
           pos = positiveDefinite m
 
-expmDiagProp m = expm (logm m) |~| complex m
+expmDiagProp m = expm (logm m) :~ 7 ~: complex m
     where logm m = matFunc log m