used an ugly const_cast to recalculate Galerkin product.

[[Imported from SVN: r2771]]

istl/paamg/hierarchy.hh

@@ -554,10 +554,10 @@ namespace Dune
       GalerkinProduct productBuilder;
 
       for(Iterator level = matrices_.finest(), coarsest=matrices_.coarsest(); level!=coarsest; ++amap) {
-        Matrix& fine = level->matrix();
+        const Matrix& fine = level->matrix();
         ++level;
-        Matrix& coarse = level->matrix();
-        productBuilder.calculate(fine, *(*amap), coarse);
+        //Matrix& coarse(level->matrix());
+        productBuilder.calculate(fine, *(*amap), const_cast<Matrix&>(level->matrix()));
 
       }
     }

istl/paamg/pmatrix.hh

@@ -40,16 +40,15 @@ namespace Dune
        * @brief Get the locally stored matrix.
        * @return The locally stored matrix.
        */
-      const Matrix& matrix()
+      const Matrix& matrix() const
       {
         return *matrix_;
       }
-
       /**
        * @brief Get the index set that maps global indices to matrix rows.
        *  @return The index set.
        */
-      const IndexSet& indexSet()
+      const IndexSet& indexSet() const
       {
         return *indices_;
       }
@@ -58,7 +57,7 @@ namespace Dune
        * @brief Get the information about remote indices also present locally.
        * @return The remote index information.
        */
-      const RemoteIndices& remoteIndices()
+      const RemoteIndices& remoteIndices() const
       {
         return *rIndices_;
       }

istl/paamg/test/hierarchytest.cc

@@ -41,6 +41,9 @@ int main(int argc, char** argv)
 
   hierarchy.build(criterion);
 
+  hierarchy.recalculateGalerkin();
+
+
   MPI_Finalize();
 
 }