X-Git-Url: http://gitweb.michael.orlitzky.com/?p=numerical-analysis.git;a=blobdiff_plain;f=src%2FLinear%2FIteration.hs;h=6fb279c78bbb64fc5122070c82bce6f956e16a14;hp=06d30a802aa55717b6a9ec45b42c69a8b8f8b473;hb=f8fe5b7210c35c84193cfcf6c26fd37eac7a1734;hpb=e91c536e40c56cd12a2fbb06ec86c729253017a8

diff --git a/src/Linear/Iteration.hs b/src/Linear/Iteration.hs
index 06d30a8..6fb279c 100644
--- 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