From e102da3cbe01987e648aa1280f417cee71757c6f Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 21 Oct 2009 18:58:09 +0000
Subject: [PATCH] remove method getNeighborsPerVertex() which is already
 implemented in the base class

[[Imported from SVN: r5084]]
---
 src/rodassembler.cc | 35 -----------------------------------
 src/rodassembler.hh |  2 --
 2 files changed, 37 deletions(-)

diff --git a/src/rodassembler.cc b/src/rodassembler.cc
index e64da341..561a1af3 100644
--- a/src/rodassembler.cc
+++ b/src/rodassembler.cc
@@ -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,
diff --git a/src/rodassembler.hh b/src/rodassembler.hh
index 5d9c7140..0a76b32b 100644
--- a/src/rodassembler.hh
+++ b/src/rodassembler.hh
@@ -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;
 
-- 
GitLab