From 45aa63a662556bc9ed0f6018f9f3f580586f38a9 Mon Sep 17 00:00:00 2001 From: Michael Orlitzky Date: Sun, 13 Apr 2014 23:20:11 -0400 Subject: [PATCH] Remove some commented code from FEM.R1. --- src/FEM/R1.hs | 28 +++------------------------- 1 file changed, 3 insertions(+), 25 deletions(-) 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) -- 2.43.2