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 (
-> 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