Skip to content
Snippets Groups Projects
Commit 42e5f9d7 authored by Oliver Sander's avatar Oliver Sander
Browse files

Remove trailing whitespace

parent e4fb27a8
Branches
Tags
No related merge requests found
......@@ -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)
{
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment