From a0b9bdde06e5cab1563bfbc989238b67a47c839d Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 5 Jan 2011 07:35:41 +0000
Subject: [PATCH] use BasisGridFunction instead of GenericGridFunction
[[Imported from SVN: r6668]]
---
dune/gfe/averageinterface.hh | 6 ++++--
1 file changed, 4 insertions(+), 2 deletions(-)
diff --git a/dune/gfe/averageinterface.hh b/dune/gfe/averageinterface.hh
index d897f571..afc67ca6 100644
--- a/dune/gfe/averageinterface.hh
+++ b/dune/gfe/averageinterface.hh
@@ -7,6 +7,8 @@
#include <dune/fufem/dgindexset.hh>
#include <dune/fufem/crossproduct.hh>
#include <dune/fufem/surfmassmatrix.hh>
+#include <dune/fufem/functions/basisgridfunction.hh>
+
#include <dune/solvers/numproc.hh>
#include "svd.hh"
@@ -425,7 +427,7 @@ void computeTotalForceAndTorque(const BoundaryPatchBase<GridView>& interface,
////////////////////////////////////////////////////////////////
typedef P1NodalBasis<GridView,double> P1Basis;
P1Basis p1Basis(interface.gridView());
- GenericGridFunction<P1Basis, dim> neumannFunction(p1Basis, boundaryStress);
+ BasisGridFunction<P1Basis, Dune::BlockVector<Dune::FieldVector<double, GridView::dimension> > > neumannFunction(p1Basis, boundaryStress);
// ///////////////////////////////////////////
// Loop and integrate over the interface
@@ -451,7 +453,7 @@ void computeTotalForceAndTorque(const BoundaryPatchBase<GridView>& interface,
const double integrationElement = segmentGeometry.integrationElement(quad[ip].position());
Dune::FieldVector<double,dim> value;
- neumannFunction.evalalllocal(*it->inside(), quadPos, value);
+ neumannFunction.evaluateLocal(*it->inside(), quadPos, value);
totalForce.axpy(quad[ip].weight() * integrationElement, value);
totalTorque.axpy(quad[ip].weight() * integrationElement, crossProduct(worldPos-center,value));
--
GitLab