From: Michael Orlitzky Date: Mon, 14 Apr 2014 03:20:11 +0000 (-0400) Subject: Remove some commented code from FEM.R1. X-Git-Url: http://gitweb.michael.orlitzky.com/?p=numerical-analysis.git;a=commitdiff_plain;h=45aa63a662556bc9ed0f6018f9f3f580586f38a9 Remove some commented code from FEM.R1. --- diff --git a/src/FEM/R1.hs b/src/FEM/R1.hs index f99b3dc..28317a6 100644 --- a/src/FEM/R1.hs +++ b/src/FEM/R1.hs @@ -583,33 +583,11 @@ solution pde params = combine ci ni = ci*(ni xi) --- energy_true :: (Arity m, Arity n, Arity l, --- Algebraic.C a, Eq a, RealField.C a, ToRational.C a) --- => PDE a --- -> Params m n (S l) a --- -> (a -> a) -- ^ True solution @u@ --- -> (a -> a) -- ^ Derivative of true solution @u'@ --- -> a --- energy_true pde params u u' = --- case (bdy pde) of --- Left (Dirichlet (x1,x2)) -> --- sqrt $ bilinear_form u u' u u' --- where --- two = fromInteger 2 --- q = affine_inv (x1,x2) --- bilinear_form w w' v v' = (x2 - x1)*(gaussian integrand)/two --- where --- integrand x = ((big_A pde) (q x))*(w' (q x))*(v' (q x)) --- + ((c pde) (q x))*(w (q x))*(v (q x)) - --- _ -> error "Neumann BCs not implemented." - - energy_fem :: (Arity m, Arity n, Arity l, Algebraic.C a, Eq a, RealField.C a, ToRational.C a) - => PDE a - -> Params m n (S l) a - -> a + => PDE a + -> Params m n (S l) a + -> a energy_fem pde params = (coefficients pde params) `dot` (big_F pde params)