Fixed bug in intialization procedure of parallel mesh distribution.

AMDiS/src/parallel/ParallelDomainBase.cc

@@ -72,11 +72,13 @@ namespace AMDiS {
     if (deserialized)
       return;
     
+
     // Test, if the mesh is the macro mesh only! Paritioning of the mesh is supported
     // only for macro meshes, so it will not work yet if the mesh is already refined
     // in some way.
     testForMacroMesh();
 
+
     // create an initial partitioning of the mesh
     partitioner->createPartitionData();
     // set the element weights, which are 1 at the very first begin
@@ -84,6 +86,7 @@ namespace AMDiS {
     // and now partition the mesh
     partitionMesh(adaptInfo);   
 
+
 #if (DEBUG != 0)
     debug::ElementIdxToDofs elMap;
     debug::createSortedDofs(mesh, elMap);
@@ -100,6 +103,7 @@ namespace AMDiS {
     ParallelDomainDbg::testAllElements(*this);
 #endif
 
+
     // === Create interior boundary information. ===
 
     createInteriorBoundaryInfo();
@@ -108,6 +112,7 @@ namespace AMDiS {
     ParallelDomainDbg::printBoundaryInfo(*this);
 #endif
 
+
     // === Create new global and local DOF numbering. ===
 
     createLocalGlobalNumbering();
@@ -172,7 +177,6 @@ namespace AMDiS {
       }
     }
 
-
     // === Remove periodic boundary conditions in sequential problem definition. ===
 
     // Remove periodic boundaries in boundary manager on matrices and vectors.
@@ -189,7 +193,7 @@ namespace AMDiS {
 	if (probStat[i]->getSolution()->getDOFVector(j)->getBoundaryManager())
 	  removeBoundaryCondition(const_cast<BoundaryIndexMap&>(probStat[i]->getSolution()->getDOFVector(j)->getBoundaryManager()->getBoundaryConditionMap()));
 	
-	if (probStat[i]->getRhs()->getDOFVector(i)->getBoundaryManager())
+	if (probStat[i]->getRhs()->getDOFVector(j)->getBoundaryManager())
 	  removeBoundaryCondition(const_cast<BoundaryIndexMap&>(probStat[i]->getRhs()->getDOFVector(j)->getBoundaryManager()->getBoundaryConditionMap()));
       }
     }