Fixed small parallel debug in deserialization code.

AMDiS/src/parallel/MeshDistributor.cc

@@ -171,15 +171,14 @@ namespace AMDiS {
       
       updateMacroElementInfo();
 
-      setRankDofs();
-
       removePeriodicBoundaryConditions();
 
       elObjDb.createMacroElementInfo(allMacroElements);
 
-      //      createBoundaryDofs();
       updateLocalGlobalNumbering();
 
+      setRankDofs();
+
 #if (DEBUG != 0)
       ParallelDebug::writeDebugFile(feSpaces[feSpaces.size() - 1], dofMap, 
 				    debugOutputDir + "mpi-dbg", "dat");
@@ -1688,11 +1687,6 @@ namespace AMDiS {
   
     lastMeshChangeIndex = mesh->getChangeIndex();
 
-    MSG("--------> REMOVE THIS LINE <------------\n");
-    debug::writeElementIndexMesh(mesh, debugOutputDir + "elementIndex-" + 
-				 lexical_cast<string>(mpiRank) + ".vtu");
-
-
 #if (DEBUG != 0)
     ParallelDebug::testDofContainerCommunication(*this);
 

AMDiS/src/parallel/PetscSolverFeti.cc

@@ -1103,7 +1103,6 @@ namespace AMDiS {
   {
     FUNCNAME("PetscSolverFeti::dbgMatrix()");
 
-#if (DEBUG != 0)
     int writeInteriorMatrix = -1;
     Parameters::get("parallel->debug->write interior matrix", 
 		    writeInteriorMatrix);
@@ -1184,6 +1183,11 @@ namespace AMDiS {
 
 	MatMult(mat_feti, unitVector, resultVector);
 
+	if (fetiPreconditioner != FETI_NONE) {
+	  PCApply(precon_feti, resultVector, unitVector);
+	  VecCopy(unitVector, resultVector);
+	}
+
 	PetscScalar *vals;
 	VecGetArray(resultVector, &vals);
 
@@ -1215,7 +1219,6 @@ namespace AMDiS {
       MSG("FETI-DP matrix written: %d x %d mat with %d nnz\n", 
 	  lagrangeMap.getOverallDofs(), lagrangeMap.getOverallDofs(), nnzCounter);
     }
-#endif
   }