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