From 3338e62f8833c6e9e30613a47e86eb3ee7ebcc51 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sat, 13 Feb 2010 10:14:39 +0000
Subject: [PATCH] New local copy of LocalStiffness needs less initialization
[[Imported from SVN: r5561]]
---
src/localgeodesicfestiffness.hh | 9 +--------
1 file changed, 1 insertion(+), 8 deletions(-)
diff --git a/src/localgeodesicfestiffness.hh b/src/localgeodesicfestiffness.hh
index 42aa5744..9b084c4a 100644
--- a/src/localgeodesicfestiffness.hh
+++ b/src/localgeodesicfestiffness.hh
@@ -5,9 +5,8 @@
#include <dune/common/fmatrix.hh>
#include <dune/istl/matrixindexset.hh>
#include <dune/istl/matrix.hh>
-#include <dune/disc/operators/localstiffness.hh>
-#include<dune/disc/operators/boundaryconditions.hh>
+#include "localstiffness.hh"
#include "rigidbodymotion.hh"
template<class GridView, class TargetSpace>
@@ -145,12 +144,6 @@ assemble(const Entity& element,
this->A = 0;
- for (int i=0; i<nDofs; i++) {
- this->b[i] = 0;
- for (int j=0; j<this->bctype[i].size(); j++)
- this->bctype[i][j] = Dune::BoundaryConditions::neumann;
- }
-
double eps = 1e-4;
typedef typename Dune::Matrix<Dune::FieldMatrix<double,blocksize,blocksize> >::row_type::iterator ColumnIterator;
--
GitLab