From b1b445697db31b1603a31747ca31151f97ee7263 Mon Sep 17 00:00:00 2001
From: Alberto Ruiz <aruiz@um.es>
Date: Sat, 23 May 2015 12:31:32 +0200
Subject: join SContainer and Container using Fractional contexts

---
 packages/base/src/Data/Packed/Internal/Numeric.hs | 53 ++++++++---------------
 packages/base/src/Data/Packed/Numeric.hs          | 13 +++---
 2 files changed, 24 insertions(+), 42 deletions(-)

(limited to 'packages/base/src/Data/Packed')

diff --git a/packages/base/src/Data/Packed/Internal/Numeric.hs b/packages/base/src/Data/Packed/Internal/Numeric.hs
index 2e36de2..00ec70c 100644
--- a/packages/base/src/Data/Packed/Internal/Numeric.hs
+++ b/packages/base/src/Data/Packed/Internal/Numeric.hs
@@ -20,7 +20,7 @@ module Data.Packed.Internal.Numeric (
     -- * Basic functions
     ident, diag, ctrans,
     -- * Generic operations
-    SContainer(..), Container(..),
+    Container(..),
     scalar, conj, scale, arctan2, cmap,
     atIndex, minIndex, maxIndex, minElement, maxElement,
     sumElements, prodElements,
@@ -115,7 +115,7 @@ m ¿¿ ec = trans (trans m ?? ec)
 
 
 -- | Basic element-by-element functions for numeric containers
-class Element e => SContainer c e
+class Element e => Container c e
   where
     conj'        :: c e -> c e
     size'        :: c e -> IndexOf c
@@ -155,24 +155,20 @@ class Element e => SContainer c e
           -> [(IndexOf c, e)] -- ^ association list
           -> c e              -- ^ result
 
-
--- | Basic element-by-element functions for numeric containers
-class (Fractional e, SContainer c e) => Container c e
-  where
     -- | scale the element by element reciprocal of the object:
     --
     -- @scaleRecip 2 (fromList [5,i]) == 2 |> [0.4 :+ 0.0,0.0 :+ (-2.0)]@
-    scaleRecip  :: e -> c e -> c e
+    scaleRecip  :: Fractional e => e -> c e -> c e
     -- | element by element division
-    divide      :: c e -> c e -> c e
+    divide      :: Fractional e => c e -> c e -> c e
     --
     -- element by element inverse tangent
-    arctan2'     :: c e -> c e -> c e
+    arctan2'     :: Fractional e => c e -> c e -> c e
 
 
 --------------------------------------------------------------------------
 
-instance SContainer Vector CInt
+instance Container Vector CInt
   where
     conj' = id
     size' = dim
@@ -198,9 +194,11 @@ instance SContainer Vector CInt
     assoc' = assocV
     accum' = accumV
 --    cond' = condV condI
+    scaleRecip = undefined -- cannot match
+    divide = undefined
+    arctan2' = undefined
 
-
-instance SContainer Vector Float
+instance Container Vector Float
   where
     conj' = id
     size' = dim
@@ -226,16 +224,13 @@ instance SContainer Vector Float
     assoc' = assocV
     accum' = accumV
     cond' = condV condF
-
-instance Container Vector Float
-  where
     scaleRecip = vectorMapValF Recip
     divide = vectorZipF Div
     arctan2' = vectorZipF ATan2
 
 
 
-instance SContainer Vector Double
+instance Container Vector Double
   where
     conj' = id
     size' = dim
@@ -261,15 +256,12 @@ instance SContainer Vector Double
     assoc' = assocV
     accum' = accumV
     cond' = condV condD
-
-instance Container Vector Double
-  where
     scaleRecip = vectorMapValR Recip
     divide = vectorZipR Div
     arctan2' = vectorZipR ATan2
 
 
-instance SContainer Vector (Complex Double)
+instance Container Vector (Complex Double)
   where
     conj' = conjugateC
     size' = dim
@@ -295,15 +287,11 @@ instance SContainer Vector (Complex Double)
     assoc' = assocV
     accum' = accumV
     cond' = undefined -- cannot match
-
-
-instance Container Vector (Complex Double)
-  where
     scaleRecip = vectorMapValC Recip
     divide = vectorZipC Div
     arctan2' = vectorZipC ATan2
 
-instance SContainer Vector (Complex Float)
+instance Container Vector (Complex Float)
   where
     conj' = conjugateQ
     size' = dim
@@ -329,16 +317,13 @@ instance SContainer Vector (Complex Float)
     assoc' = assocV
     accum' = accumV
     cond' = undefined -- cannot match
-
-instance Container Vector (Complex Float)
-  where
     scaleRecip = vectorMapValQ Recip
     divide = vectorZipQ Div
     arctan2' = vectorZipQ ATan2
 
 ---------------------------------------------------------------
 
-instance (Num a, Element a, SContainer Vector a) => SContainer Matrix a
+instance (Num a, Element a, Container Vector a) => Container Matrix a
   where
     conj' = liftMatrix conj'
     size' = size
@@ -366,10 +351,6 @@ instance (Num a, Element a, SContainer Vector a) => SContainer Matrix a
     assoc' = assocM
     accum' = accumM
     cond' = condM
-
-
-instance (Fractional a, Container Vector a) => Container Matrix a
-  where
     scaleRecip x = liftMatrix (scaleRecip x)
     divide = liftMatrix2 divide
     arctan2' = liftMatrix2 arctan2'
@@ -404,7 +385,7 @@ conj = conj'
 scale :: Container c e => e -> c e -> c e
 scale = scale'
 
-arctan2 :: Container c e => c e -> c e -> c e
+arctan2 :: (Fractional e, Container c e) => c e -> c e -> c e
 arctan2 = arctan2'
 
 -- | like 'fmap' (cannot implement instance Functor because of Element class constraint)
@@ -754,7 +735,7 @@ buildV n f = fromList [f k | k <- ks]
 
 --------------------------------------------------------
 -- | conjugate transpose
-ctrans :: (SContainer Vector e, Element e) => Matrix e -> Matrix e
+ctrans :: (Container Vector e, Element e) => Matrix e -> Matrix e
 ctrans = liftMatrix conj' . trans
 
 -- | Creates a square matrix with a given diagonal.
@@ -810,7 +791,7 @@ class Transposable m mt | m -> mt, mt -> m
     -- | (conjugate) transpose
     tr :: m -> mt
 
-instance (SContainer Vector t) => Transposable (Matrix t) (Matrix t)
+instance (Container Vector t) => Transposable (Matrix t) (Matrix t)
   where
     tr = ctrans
 
diff --git a/packages/base/src/Data/Packed/Numeric.hs b/packages/base/src/Data/Packed/Numeric.hs
index 4d66f27..d11ecf9 100644
--- a/packages/base/src/Data/Packed/Numeric.hs
+++ b/packages/base/src/Data/Packed/Numeric.hs
@@ -31,7 +31,7 @@ module Data.Packed.Numeric (
     diag, ident,
     ctrans,
     -- * Generic operations
-    SContainer(..), Container(..), Numeric,
+    Container(..), Numeric,
     -- add, mul, sub, divide, equal, scaleRecip, addConstant,
     scalar, conj, scale, arctan2, cmap,
     atIndex, minIndex, maxIndex, minElement, maxElement,
@@ -88,7 +88,7 @@ Logarithmic spacing can be defined as follows:
 
 @logspace n (a,b) = 10 ** linspace n (a,b)@
 -}
-linspace :: (Container Vector e) => Int -> (e, e) -> Vector e
+linspace :: (Fractional e, Container Vector e) => Int -> (e, e) -> Vector e
 linspace 0 _     = fromList[]
 linspace 1 (a,b) = fromList[(a+b)/2]
 linspace n (a,b) = addConstant a $ scale s $ fromList $ map fromIntegral [0 .. n-1]
@@ -219,11 +219,11 @@ class Konst e d c | d -> c, c -> d
     --
     konst :: e -> d -> c e
 
-instance SContainer Vector e => Konst e Int Vector
+instance Container Vector e => Konst e Int Vector
   where
     konst = konst'
 
-instance (Num e, SContainer Vector e) => Konst e (Int,Int) Matrix
+instance (Num e, Container Vector e) => Konst e (Int,Int) Matrix
   where
     konst = konst'
 
@@ -246,11 +246,11 @@ class Build d f c e | d -> c, c -> d, f -> e, f -> d, f -> c, c e -> f, d e -> f
     --
     build :: d -> f -> c e
 
-instance SContainer Vector e => Build Int (e -> e) Vector e
+instance Container Vector e => Build Int (e -> e) Vector e
   where
     build = build'
 
-instance SContainer Matrix e => Build (Int,Int) (e -> e -> e) Matrix e
+instance Container Matrix e => Build (Int,Int) (e -> e -> e) Matrix e
   where
     build = build'
 
@@ -299,5 +299,6 @@ instance Numeric Double
 instance Numeric (Complex Double)
 instance Numeric Float
 instance Numeric (Complex Float)
+instance Numeric CInt
 
 
-- 
cgit v1.2.3