bo eh, part iii

AMDiS/src/parallel/PetscSolverFeti.cc

@@ -1281,12 +1281,12 @@ namespace AMDiS {
 	  // === Create mass matrix solver ===
 	  
 	  if (!massMatrixSolver) {
-	    MSG("START CREATE MASS MATRIX!\n");
 	    ParallelDofMapping *massMapping = new ParallelDofMapping(COMPONENT_WISE);
 	    massMapping->init(meshDistributor->getMeshLevelData(), 
 			      pressureFeSpace, pressureFeSpace);
 	    massMapping->setDofComm(meshDistributor->getDofComm());
 	    massMapping->setMpiComm(meshDistributor->getMeshLevelData().getMpiComm(0), 0);
+	    massMapping->setComputeMatIndex(true);
 	    (*massMapping)[0] = interfaceDofMap[pressureComponent];
 	    massMapping->update();	    
 
@@ -1302,9 +1302,7 @@ namespace AMDiS {
 						 mpiCommGlobal,
 						 mpiCommLocal);
 	    massMatrixSolver->setDofMapping(massMapping);
-	    MSG("START ASM ===\n");
 	    massMatrixSolver->fillPetscMatrix(&massMatrix);
-	    MSG("END ASM ===\n");
 
 	    int r, c;
 	    MatGetSize(massMatrixSolver->getMatInterior(), &r, &c);
@@ -1313,7 +1311,6 @@ namespace AMDiS {
 	    MSG("MASS MAT INFO:  size = %d x %d   nnz = %d\n",
 		r, c,
 		static_cast<int>(info.nz_used));
-	    MSG("END CREATE MASS MATRIX!\n");
 
 	    fetiInterfaceLumpedPreconData.ksp_mass = 
 	      massMatrixSolver->getSolver();