From 42e5f9d77112956b89de1323b55f06a78a577143 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 12 Feb 2015 22:41:05 +0100
Subject: [PATCH] Remove trailing whitespace
---
dune/fufem/discretizationerror.hh | 26 +++++++++++++-------------
1 file changed, 13 insertions(+), 13 deletions(-)
diff --git a/dune/fufem/discretizationerror.hh b/dune/fufem/discretizationerror.hh
index b2bd1e5f..d21b4d09 100644
--- a/dune/fufem/discretizationerror.hh
+++ b/dune/fufem/discretizationerror.hh
@@ -24,7 +24,7 @@ public:
QuadratureRuleKey quadKey);
static double computeH1HalfNormDifferenceSquared(const GridView& gridView,
- const VirtualDifferentiableFunction<Dune::FieldVector<ctype,dim>, Dune::FieldVector<double,1> >* u,
+ const VirtualDifferentiableFunction<Dune::FieldVector<ctype,dim>, Dune::FieldVector<double,1> >* u,
const VirtualDifferentiableFunction<Dune::FieldVector<ctype,dim>, Dune::FieldVector<double,1> >* v,
QuadratureRuleKey quadKey);
@@ -37,7 +37,7 @@ computeL2Error(const VirtualGridViewFunction<GridView,Dune::FieldVector<double,1
const Dune::VirtualFunction<Dune::FieldVector<ctype,dim>, Dune::FieldVector<double,1> >* b,
QuadratureRuleKey quadKey)
{
-
+
// The error to be computed
double error = 0;
@@ -47,42 +47,42 @@ computeL2Error(const VirtualGridViewFunction<GridView,Dune::FieldVector<double,1
typename GridView::template Codim<0>::Iterator eEndIt = gridView.template end<0>();
for (; eIt!=eEndIt; ++eIt) {
-
+
// Get quadrature formula
quadKey.setGeometryType(eIt->type());
const Dune::QuadratureRule<ctype,dim>& quad = QuadratureRuleCache<double, dim>::rule(quadKey);
-
+
for (size_t i=0; i<quad.size(); i++) {
-
+
// Evaluate function a
Dune::FieldVector<double,1> aValue;
a->evaluateLocal(*eIt, quad[i].position(),aValue);
-
+
// Evaluate function b. If it is a grid function use that to speed up the evaluation
Dune::FieldVector<double,1> bValue;
if (dynamic_cast<const VirtualGridViewFunction<GridView,Dune::FieldVector<double,1> >*>(b))
- dynamic_cast<const VirtualGridViewFunction<GridView,Dune::FieldVector<double,1> >*>(b)->evaluateLocal(*eIt,
+ dynamic_cast<const VirtualGridViewFunction<GridView,Dune::FieldVector<double,1> >*>(b)->evaluateLocal(*eIt,
quad[i].position(),
bValue
);
else
b->evaluate(eIt->geometry().global(quad[i].position()), bValue);
-
+
// integrate error
error += (aValue - bValue) * (aValue - bValue) * quad[i].weight() * eIt->geometry().integrationElement(quad[i].position());
-
+
}
-
+
}
-
- return std::sqrt(error);
+
+ return std::sqrt(error);
}
template <class GridView>
double DiscretizationError<GridView>::
computeH1HalfNormDifferenceSquared(const GridView& gridView,
- const VirtualDifferentiableFunction<Dune::FieldVector<ctype,dim>, Dune::FieldVector<double,1> >* u,
+ const VirtualDifferentiableFunction<Dune::FieldVector<ctype,dim>, Dune::FieldVector<double,1> >* u,
const VirtualDifferentiableFunction<Dune::FieldVector<ctype,dim>, Dune::FieldVector<double,1> >* v,
QuadratureRuleKey quadKey)
{
--
GitLab