3 files changed, 19 insertions, 1 deletions

diff --git a/THANKS b/THANKS
index d8a9e0e..df47b12 100644
--- a/THANKS
+++ b/THANKS
@@ -76,3 +76,4 @@ and all the people in the Haskell mailing lists for their help.
   case that /tmp is not writable. He also contributed Binary instances for
   Vector and Matrix.
 
+- Vivian McPhail added support for Float and Complex Float elements.
diff --git a/lib/Data/Packed/Vector.hs b/lib/Data/Packed/Vector.hs
index bb89506..c5b8ba7 100644
--- a/lib/Data/Packed/Vector.hs
+++ b/lib/Data/Packed/Vector.hs
@@ -21,7 +21,8 @@ module Data.Packed.Vector (
     vecdisp,
 --    moved to Numeric.LinearAlgebra.Interface typeclass
 --    vectorFMax, vectorFMin, vectorFMaxIndex, vectorFMinIndex,
---    vectorMax, vectorMin, vectorMaxIndex, vectorMinIndex,
+--    vectorMax, vectorMin,
+    vectorMaxIndex, vectorMinIndex,
     mapVector, zipVector,
     fscanfVector, fprintfVector, freadVector, fwriteVector,
     foldLoop, foldVector, foldVectorG
diff --git a/lib/Numeric/LinearAlgebra/Tests.hs b/lib/Numeric/LinearAlgebra/Tests.hs
index 65b0e4c..24b5bb5 100644
--- a/lib/Numeric/LinearAlgebra/Tests.hs
+++ b/lib/Numeric/LinearAlgebra/Tests.hs
@@ -356,6 +356,7 @@ makeUnitary v | realPart n > 1    = v / scalar n
 -- | Performance measurements.
 runBenchmarks :: IO ()
 runBenchmarks = do
+  --cholBench
     solveBench
     subBench
     multBench
@@ -455,3 +456,18 @@ solveBench = do
     solveBenchN 500
     solveBenchN 1000
     -- solveBenchN 1500
+
+--------------------------------
+
+cholBenchN n = do
+    let x = uniformSample 777 (2*n) (replicate n (-1,1))
+        a = trans x <> x
+    a `seq` putStrLn ""
+    time ("chol " ++ show n) (chol a)
+
+cholBench = do
+    cholBenchN 1200
+    cholBenchN 600
+    cholBenchN 300
+--    cholBenchN 150
+--    cholBenchN 50