Small bugfix

AMDiS/src/parallel/PetscSolver.cc

@@ -43,7 +43,11 @@ namespace AMDiS {
   {
     interiorMap = interiorDofs;
     coarseSpaceMap = coarseDofs;
+  }
+
 
+  void PetscSolver::updateSubdomainData()
+  {
     if (mpiCommLocal.Get_size() == 1) {
       rStartInterior = 0;
       nGlobalOverallInterior = interiorMap->getOverallDofs();
@@ -63,7 +67,7 @@ namespace AMDiS {
     }
   }
 
-  
+ 
   void PetscSolver::solve(Vec &rhs, Vec &sol)
   {
     FUNCNAME("PetscSolver::solve()");

AMDiS/src/parallel/PetscSolver.h

@@ -202,6 +202,8 @@ namespace AMDiS {
     /// Returns a vector with the FE spaces of all components of a system vector.
     vector<const FiniteElemSpace*> getFeSpaces(SystemVector *vec);
 
+    void updateSubdomainData();
+
   protected:
     MeshDistributor *meshDistributor;
 

AMDiS/src/parallel/PetscSolverGlobalMatrix.cc

@@ -22,6 +22,7 @@ namespace AMDiS {
     FUNCNAME("PetscSolverGlobalMatrix::fillPetscMatrix()");
 
     if (coarseSpaceMap != NULL) {
+      updateSubdomainData();
       fillPetscMatrixWithCoarseSpace(mat);
       return;
     }
@@ -53,6 +54,7 @@ namespace AMDiS {
 	o_nnz = NULL;
       }
 
+      updateSubdomainData();
       createPetscNnzStructure(mat);
       lastMeshNnz = meshDistributor->getLastMeshChangeIndex();
     }
@@ -128,8 +130,6 @@ namespace AMDiS {
     int nRowsOverallInterior = interiorMap->getOverallDofs();
 
     if (subdomainLevel == 0) {
-      nGlobalOverallInterior = nRowsOverallInterior;
-
       MatCreateSeqAIJ(mpiCommLocal, nRowsRankInterior, nRowsRankInterior,
 		      60, PETSC_NULL, &matIntInt);
     } else {