From 150029d0ecd4bf05373fab10a06bcd4e77d5509f Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sat, 22 Mar 2014 18:42:08 +0000
Subject: [PATCH] Switch to second-order finite elements
[[Imported from SVN: r9667]]
---
cosserat-continuum.cc | 7 +++++++
1 file changed, 7 insertions(+)
diff --git a/cosserat-continuum.cc b/cosserat-continuum.cc
index 30be33a4..c611292c 100644
--- a/cosserat-continuum.cc
+++ b/cosserat-continuum.cc
@@ -1,5 +1,7 @@
#include <config.h>
+#define SECOND_ORDER
+
#include <fenv.h>
// Includes for the ADOL-C automatic differentiation library
@@ -25,6 +27,7 @@
#include <dune/fufem/functiontools/boundarydofs.hh>
#include <dune/fufem/functiontools/basisinterpolator.hh>
#include <dune/fufem/functionspacebases/p1nodalbasis.hh>
+#include <dune/fufem/functionspacebases/p2nodalbasis.hh>
#include <dune/solvers/solvers/iterativesolver.hh>
#include <dune/solvers/norms/energynorm.hh>
@@ -191,7 +194,11 @@ int main (int argc, char *argv[]) try
typedef GridType::LeafGridView GridView;
GridView gridView = grid->leafGridView();
+#ifdef SECOND_ORDER
+ typedef P2NodalBasis<GridView,double> FEBasis;
+#else
typedef P1NodalBasis<GridView,double> FEBasis;
+#endif
FEBasis feBasis(gridView);
// /////////////////////////////////////////
--
GitLab