From 1f4aa4e88d8e4131ce681e8d0aa308f38334299e Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sun, 20 Nov 2011 11:55:57 +0000
Subject: [PATCH] choose the correct quadrature order for the matrices that
create the norms to measure the discretization error
[[Imported from SVN: r8243]]
---
harmonicmaps-eoc.cc | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/harmonicmaps-eoc.cc b/harmonicmaps-eoc.cc
index 5b8470ce..b432f74f 100644
--- a/harmonicmaps-eoc.cc
+++ b/harmonicmaps-eoc.cc
@@ -4,6 +4,7 @@
//#define HARMONIC_ENERGY_FD_INNER_GRADIENT
#define THIRD_ORDER
//#define SECOND_ORDER
+const int order = 3;
#include <dune/common/bitsetvector.hh>
#include <dune/common/parametertree.hh>
@@ -226,8 +227,8 @@ int main (int argc, char *argv[]) try
FEBasis referenceBasis(referenceGrid->leafView());
OperatorAssembler<FEBasis,FEBasis> operatorAssembler(referenceBasis, referenceBasis);
- LaplaceAssembler<GridType, FEBasis::LocalFiniteElement, FEBasis::LocalFiniteElement> laplaceLocalAssembler;
- MassAssembler<GridType, FEBasis::LocalFiniteElement, FEBasis::LocalFiniteElement> massMatrixLocalAssembler;
+ LaplaceAssembler<GridType, FEBasis::LocalFiniteElement, FEBasis::LocalFiniteElement> laplaceLocalAssembler(2*(order-1));
+ MassAssembler<GridType, FEBasis::LocalFiniteElement, FEBasis::LocalFiniteElement> massMatrixLocalAssembler(2*order);
typedef Dune::BCRSMatrix<Dune::FieldMatrix<double,1,1> > ScalarMatrixType;
ScalarMatrixType laplace, massMatrix;
--
GitLab