From 01744e1df9448a405696d91af709a7416db0cb78 Mon Sep 17 00:00:00 2001 From: Alberto Ruiz Date: Wed, 15 Oct 2008 12:41:28 +0000 Subject: expose-bug flag --- lib/Numeric/LinearAlgebra/Algorithms.hs | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) (limited to 'lib/Numeric/LinearAlgebra/Algorithms.hs') 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 @@ -142,9 +142,13 @@ 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" -- cgit v1.2.3