Add Eq constraints throughout Linear.Iteration.

diff --git a/src/Linear/Iteration.hs b/src/Linear/Iteration.hs
index 06d30a802aa55717b6a9ec45b42c69a8b8f8b473..6fb279c78bbb64fc5122070c82bce6f956e16a14 100644 (file)
--- a/src/Linear/Iteration.hs
+++ b/src/Linear/Iteration.hs
@@ -41,7 +41,7 @@ import Normed ( Normed(..) )
 -- | A generalized implementation for Jacobi, Gauss-Seidel, etc. All
 --   that we really need to know is how to construct the matrix M, so we
 --   take a function that does it as an argument.
-classical_iteration :: (Field.C a, Arity m)
+classical_iteration :: (Eq a, Field.C a, Arity m)
                  => (Mat m m  a -> Mat m m  a)
                  -> Mat m m  a
                  -> Mat m N1 a
@@ -60,7 +60,7 @@ classical_iteration m_function matrix b x_current =
 -- | Perform one iteration of successive over-relaxation.
 --
 sor_iteration :: forall m a.
-                 (Field.C a, Arity m)
+                 (Eq a, Field.C a, Arity m)
               => a -- ^ Omega
               -> Mat m m  a -- ^ Matrix A
               -> Mat m N1 a -- ^ Vector b
@@ -79,7 +79,7 @@ sor_iteration omega =
 
 -- | Compute an infinite list of SOR iterations starting with the
 --   vector x0.
-sor_iterations :: (Field.C a, Arity m)
+sor_iterations :: (Eq a, Field.C a, Arity m)
                => a
                -> Mat m m  a
                -> Mat m N1 a
@@ -90,7 +90,7 @@ sor_iterations omega matrix b =
 
 
 -- | Perform one iteration of Gauss-Seidel.
-gauss_seidel_iteration :: (Field.C a, Arity m)
+gauss_seidel_iteration :: (Eq a, Field.C a, Arity m)
                        => Mat m m  a
                        -> Mat m N1 a
                        -> Mat m N1 a
@@ -100,7 +100,7 @@ gauss_seidel_iteration = sor_iteration one
 
 -- | Compute an infinite list of Gauss-Seidel iterations starting with
 --   the vector x0.
-gauss_seidel_iterations :: (Field.C a, Arity m)
+gauss_seidel_iterations :: (Eq a, Field.C a, Arity m)
                         => Mat m m  a
                         -> Mat m N1 a
                         -> Mat m N1 a
@@ -126,7 +126,7 @@ gauss_seidel_iterations matrix b =
 --   >>> jacobi_iteration m b x1
 --   ((0.0),(0.25))
 --
-jacobi_iteration :: (Field.C a, Arity m)
+jacobi_iteration :: (Eq a, Field.C a, Arity m)
                  => Mat m m  a
                  -> Mat m N1 a
                  -> Mat m N1 a
@@ -137,7 +137,7 @@ jacobi_iteration =
 
 -- | Compute an infinite list of Jacobi iterations starting with the
 --   vector x0.
-jacobi_iterations :: (Field.C a, Arity m)
+jacobi_iterations :: (Eq a, Field.C a, Arity m)
                   => Mat m m  a
                   -> Mat m N1 a
                   -> Mat m N1 a