X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=src%2FFEM%2FR1.hs;h=dd5c11b023b9f267ddf2ac70479e2b4c6a004f85;hb=3274b18a2fb48f26e6582d9dcee11b4067230911;hp=28317a6faecabaf0a079215a192f961a7d5e5be3;hpb=45aa63a662556bc9ed0f6018f9f3f580586f38a9;p=numerical-analysis.git

diff --git a/src/FEM/R1.hs b/src/FEM/R1.hs
index 28317a6..dd5c11b 100644
--- a/src/FEM/R1.hs
+++ b/src/FEM/R1.hs
@@ -35,7 +35,7 @@ import qualified Algebra.RealField as RealField ( C )
 import qualified Algebra.ToRational as ToRational ( C )
 import Data.Vector.Fixed ( Arity, S )
 import NumericPrelude hiding ( abs )
-import qualified Prelude as P
+import Prelude ()
 
 import Integration.Gaussian ( gaussian )
 import Linear.Matrix (
@@ -615,8 +615,7 @@ relative_error_pointwise :: forall m n l a.
                          -> a -- ^ The point @x@ at which to compute the error.
                          -> a
 relative_error_pointwise pde params u x =
-  cent * ( u_exact - u_fem ) / u_exact
+  cent * ( abs $ (u x) - u_fem ) / ( abs $ u x )
   where
-    u_exact = abs $ u x
     u_fem = evaluate' (solution pde params) x
     cent = fromInteger 100