From 99bb590e8cc49df4bf31c72fc5fa067380f44688 Mon Sep 17 00:00:00 2001
From: Thomas Witkowski <thomas.witkowski@gmx.de>
Date: Tue, 23 Oct 2012 12:26:14 +0000
Subject: [PATCH] bo eh, part iii
---
AMDiS/src/parallel/PetscSolverFeti.cc | 5 +----
1 file changed, 1 insertion(+), 4 deletions(-)
diff --git a/AMDiS/src/parallel/PetscSolverFeti.cc b/AMDiS/src/parallel/PetscSolverFeti.cc
index 48a9ed03..72402ac3 100644
--- a/AMDiS/src/parallel/PetscSolverFeti.cc
+++ b/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();
--
GitLab