move conjugate* to Container

The benefit of this patch is that Data.Packed.Internal.Matrix now has no numeric dependency.

There is no 'optimisation' class for conjugateD, instead, different instances are directly specified in each Container instance.

Unfortunately, an unsafeIO call now exists in Numeric.Container. 

1 files changed, 6 insertions, 19 deletions

diff --git a/lib/Data/Packed/Internal/Matrix.hs b/lib/Data/Packed/Internal/Matrix.hs
index f3bab76..7b823de 100644
--- a/lib/Data/Packed/Internal/Matrix.hs
+++ b/lib/Data/Packed/Internal/Matrix.hs
@@ -257,28 +257,30 @@ class (Storable a) => Element a where
     transdata = transdataP -- transdata'
     constantD  :: a -> Int -> Vector a
     constantD = constantP -- constant'
+{-
     conjugateD :: Vector a -> Vector a
     conjugateD = id
+-}
 
 instance Element Float where
     transdata  = transdataAux ctransF
     constantD  = constantAux cconstantF
-    conjugateD = id
+--    conjugateD = id
 
 instance Element Double where
     transdata  = transdataAux ctransR
     constantD  = constantAux cconstantR
-    conjugateD = id
+--    conjugateD = id
 
 instance Element (Complex Float) where
     transdata  = transdataAux ctransQ
     constantD  = constantAux cconstantQ
-    conjugateD = conjugateQ
+--    conjugateD = conjugateQ
 
 instance Element (Complex Double) where
     transdata  = transdataAux ctransC
     constantD  = constantAux cconstantC
-    conjugateD = conjugateC
+--    conjugateD = conjugateC
 
 -------------------------------------------------------------------
 
@@ -390,21 +392,6 @@ constantP a n = unsafePerformIO $ do
     return v
 foreign import ccall "constantP" cconstantP :: Ptr () -> CInt -> Ptr () -> CInt -> IO CInt
 
----------------------------------------
-
-conjugateAux fun x = unsafePerformIO $ do
-    v <- createVector (dim x)
-    app2 fun vec x vec v "conjugateAux"
-    return v
-
-conjugateQ :: Vector (Complex Float) -> Vector (Complex Float)
-conjugateQ = conjugateAux c_conjugateQ
-foreign import ccall "conjugateQ" c_conjugateQ :: TQVQV
-
-conjugateC :: Vector (Complex Double) -> Vector (Complex Double)
-conjugateC = conjugateAux c_conjugateC
-foreign import ccall "conjugateC" c_conjugateC :: TCVCV
-
 ----------------------------------------------------------------------
 
 -- | Extracts a submatrix from a matrix.