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

remove method getNeighborsPerVertex() which is already implemented in the base class

[[Imported from SVN: r5084]]
parent 945ba340
No related branches found
No related tags found
No related merge requests found
......@@ -10,41 +10,6 @@
#include "src/rodlocalstiffness.hh"
template <class GridType>
void RodAssembler<GridType>::
getNeighborsPerVertex(Dune::MatrixIndexSet& nb) const
{
const int gridDim = GridType::dimension;
const typename GridType::Traits::LevelIndexSet& indexSet = grid_->levelIndexSet(grid_->maxLevel());
int i, j;
int n = grid_->size(grid_->maxLevel(), gridDim);
nb.resize(n, n);
ElementIterator it = grid_->template lbegin<0>( grid_->maxLevel() );
ElementIterator endit = grid_->template lend<0> ( grid_->maxLevel() );
for (; it!=endit; ++it) {
for (i=0; i<it->template count<gridDim>(); i++) {
for (j=0; j<it->template count<gridDim>(); j++) {
int iIdx = indexSet.subIndex(*it,i,gridDim);
int jIdx = indexSet.subIndex(*it,j,gridDim);
nb.add(iIdx, jIdx);
}
}
}
}
template <class GridType>
void RodAssembler<GridType>::
assembleGradient(const std::vector<RigidBodyMotion<3> >& sol,
......
......@@ -89,8 +89,6 @@ class RodAssembler : public GeodesicFEAssembler<typename GridType::LeafGridView,
/** \brief Compute the energy of a deformation state */
double computeEnergy(const std::vector<RigidBodyMotion<3> >& sol) const;
void getNeighborsPerVertex(Dune::MatrixIndexSet& nb) const;
void getStrain(const std::vector<RigidBodyMotion<3> >& sol,
Dune::BlockVector<Dune::FieldVector<double, blocksize> >& strain) const;
......
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