expose-bug flag

1 files changed, 7 insertions, 1 deletions

diff --git a/lib/Numeric/LinearAlgebra/Algorithms.hs b/lib/Numeric/LinearAlgebra/Algorithms.hs
index 2c5f44c..00a0ab0 100644
--- a/lib/Numeric/LinearAlgebra/Algorithms.hs
+++ b/lib/Numeric/LinearAlgebra/Algorithms.hs
@@ -143,8 +143,12 @@ instance Field (Complex Double) where
 #if defined(WORKAROUND)
     multiply = mulCW
 #else
+#if defined(EXPOSEBUG)
+    multiply = multiplyC
+#else
     multiply = multiplyC3
 #endif
+#endif
 
 -- | Eigenvalues and Eigenvectors of a complex hermitian or real symmetric matrix using lapack's dsyev or zheev.
 --
@@ -664,7 +668,9 @@ cmultiply f st a b
       = unsafePerformIO $ do
         s <- createMatrix RowMajor (rows a) (cols b)
         app3 f mat a mat b mat s st
-        if toLists s== toLists s then return s else error $ "BRUTAL " ++ (show (toLists s))
+        if toLists s== toLists s
+            then return s
+            else error $ "NaN FOUND!! " ++ (show (toLists s))
         -- return s
 --    | otherwise = error $ st ++ " (matrix product) of nonconformant matrices"