Remove some commented code from FEM.R1.

diff --git a/src/FEM/R1.hs b/src/FEM/R1.hs
index f99b3dcace7be32ee7df3ba07a488dc012285ae2..28317a6faecabaf0a079215a192f961a7d5e5be3 100644 (file)
--- 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)