Use the 2-norm instead of the infty-norm for the fixed point error ratios.

[numerical-analysis.git] / src / Roots / Simple.hs

diff --git a/src/Roots/Simple.hs b/src/Roots/Simple.hs
index 3237d60217aeb9e94bafe6983cf8f6a5dd353727..f9b36fad6578733fe8cad0308e573a77aa37640b 100644 (file)
--- a/src/Roots/Simple.hs
+++ b/src/Roots/Simple.hs
@@ -260,6 +260,6 @@ fixed_point_error_ratios f x0 x_star p =
   zipWith (/) en_plus_one en_exp
   where
     xn = F.fixed_point_iterations f x0
-    en = map (\x -> norm (x_star - x)) xn
+    en = map (\x -> norm_2 (x_star - x)) xn
     en_plus_one = tail en
     en_exp = map (^p) en