import qualified Algebra.ToRational as ToRational ( C )
import Data.Vector.Fixed ( Arity, S )
import NumericPrelude hiding ( abs )
import qualified Algebra.ToRational as ToRational ( C )
import Data.Vector.Fixed ( Arity, S )
import NumericPrelude hiding ( abs )
(x1,x2) = (mesh params) !!! (i,0)
q = affine_inv (x1,x2)
integrand x = ((f pde) (q x)) * (this_N x)
(x1,x2) = (mesh params) !!! (i,0)
q = affine_inv (x1,x2)
integrand x = ((f pde) (q x)) * (this_N x)
(x1,x2) = (mesh params) !!! (k,0)
q = affine_inv (x1,x2)
integrand x = ((big_A pde) (q x)) * (these_N's x)
(x1,x2) = (mesh params) !!! (k,0)
q = affine_inv (x1,x2)
integrand x = ((big_A pde) (q x)) * (these_N's x)
(x1,x2) = (mesh params) !!! (k,0)
q = affine_inv (x1,x2)
integrand x = ((c pde) (q x)) * (these_Ns x)
(x1,x2) = (mesh params) !!! (k,0)
q = affine_inv (x1,x2)
integrand x = ((c pde) (q x)) * (these_Ns x)
relative_error pde params energy_true =
cent * sqrt(energy_true - (energy_fem pde params)/energy_true)
where
relative_error pde params energy_true =
cent * sqrt(energy_true - (energy_fem pde params)/energy_true)
where