Skip to content
Snippets Groups Projects
Commit 3338e62f authored by Oliver Sander's avatar Oliver Sander Committed by sander@PCPOOL.MI.FU-BERLIN.DE
Browse files

New local copy of LocalStiffness needs less initialization

[[Imported from SVN: r5561]]
parent 69f5d654
No related branches found
No related tags found
No related merge requests found
...@@ -5,9 +5,8 @@ ...@@ -5,9 +5,8 @@
#include <dune/common/fmatrix.hh> #include <dune/common/fmatrix.hh>
#include <dune/istl/matrixindexset.hh> #include <dune/istl/matrixindexset.hh>
#include <dune/istl/matrix.hh> #include <dune/istl/matrix.hh>
#include <dune/disc/operators/localstiffness.hh>
#include<dune/disc/operators/boundaryconditions.hh>
#include "localstiffness.hh"
#include "rigidbodymotion.hh" #include "rigidbodymotion.hh"
template<class GridView, class TargetSpace> template<class GridView, class TargetSpace>
...@@ -145,12 +144,6 @@ assemble(const Entity& element, ...@@ -145,12 +144,6 @@ assemble(const Entity& element,
this->A = 0; 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; double eps = 1e-4;
typedef typename Dune::Matrix<Dune::FieldMatrix<double,blocksize,blocksize> >::row_type::iterator ColumnIterator; typedef typename Dune::Matrix<Dune::FieldMatrix<double,blocksize,blocksize> >::row_type::iterator ColumnIterator;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment