X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FRoots%2FSimple.hs;h=414b77787262e608ceef4a66981a5dc5ab82cc50;hb=2f54e89d36e835c58efcc281741632d457859b20;hp=0a1debff0a0b1db653efabdf2be00ce9a4b91ed4;hpb=59c49750fd2455574fe4e67ddd7e67a20321c8a8;p=numerical-analysis.git diff --git a/src/Roots/Simple.hs b/src/Roots/Simple.hs index 0a1debf..414b777 100644 --- a/src/Roots/Simple.hs +++ b/src/Roots/Simple.hs @@ -18,8 +18,9 @@ import Normed import qualified Roots.Fast as F import NumericPrelude hiding (abs) -import qualified Algebra.Absolute as Absolute import Algebra.Absolute (abs) +import qualified Algebra.Additive as Additive +import qualified Algebra.Algebraic as Algebraic import qualified Algebra.Field as Field import qualified Algebra.RealField as RealField import qualified Algebra.RealRing as RealRing @@ -87,7 +88,7 @@ bisect f a b epsilon = -- at x0. We delegate to the version that returns the number of -- iterations and simply discard the number of iterations. -- -fixed_point :: (Normed a, RealField.C b) +fixed_point :: (Normed a, Additive.C a, Algebraic.C b, RealField.C b) => (a -> a) -- ^ The function @f@ to iterate. -> b -- ^ The tolerance, @epsilon@. -> a -- ^ The initial value @x0@. @@ -100,7 +101,10 @@ fixed_point f epsilon x0 = -- the function @f@ with the search starting at x0 and tolerance -- @epsilon@. We delegate to the version that returns the number of -- iterations and simply discard the fixed point. -fixed_point_iteration_count :: (Normed a, RealField.C b) +fixed_point_iteration_count :: (Normed a, + Additive.C a, + RealField.C b, + Algebraic.C b) => (a -> a) -- ^ The function @f@ to iterate. -> b -- ^ The tolerance, @epsilon@. -> a -- ^ The initial value @x0@. @@ -118,7 +122,10 @@ fixed_point_iteration_count f epsilon x0 = -- -- This is used to determine the rate of convergence. -- -fixed_point_error_ratios :: (Normed a, RealField.C b) +fixed_point_error_ratios :: (Normed a, + Additive.C a, + RealField.C b, + Algebraic.C b) => (a -> a) -- ^ The function @f@ to iterate. -> a -- ^ The initial value @x0@. -> a -- ^ The true solution, @x_star@.