diff --git a/AMDiS/src/RefinementManager3d.cc b/AMDiS/src/RefinementManager3d.cc index 7a34ab91c4472e11123982870eeb9b0a6c45d87e..b14180c4c8fe3b59d491f237ad54232aef83042d 100644 --- a/AMDiS/src/RefinementManager3d.cc +++ b/AMDiS/src/RefinementManager3d.cc @@ -675,8 +675,6 @@ namespace AMDiS { bool foundEdge = false; while (elInfo2) { - MSG("TRY TO FIND ON EL %d\n", elInfo2->getElement()->getIndex()); - for (int i = 0; i < 6; i++) { DofEdge edge2 = elInfo2->getElement()->getEdge(i); if (edge2.first == *(edge[0]) && @@ -690,8 +688,6 @@ namespace AMDiS { // must be refined at least once to get a refinement edge. if (i == 0) { - MSG("AND FOUND!\n"); - // Edge is refinement edge, so add it to refine list. refineList.setElType(n_neigh, elInfo2->getType()); diff --git a/AMDiS/src/parallel/MeshManipulation.cc b/AMDiS/src/parallel/MeshManipulation.cc index b982e741766ebff95ac4ec8a9c7628716d63c780..256cfa2fd5bda66c027f30a51317e64c63491ffd 100644 --- a/AMDiS/src/parallel/MeshManipulation.cc +++ b/AMDiS/src/parallel/MeshManipulation.cc @@ -396,7 +396,6 @@ namespace AMDiS { // In this scenario, the element must be refined due to the structure // code, but the refinement does not introduce new DOFs on the face, // that should be adapted. Thus, we can ommit the refinement. - if (subObj == FACE) { if (s0 != -1 && s1 == -1 || s0 == -1 && s1 != -1) { if (ithObj <= 1 && code.lookAhead() == 0) { diff --git a/AMDiS/src/parallel/PetscSolverGlobalMatrix.cc b/AMDiS/src/parallel/PetscSolverGlobalMatrix.cc index 13608e10ff4ebc4e4ba097afc7ca3f450db77595..602ef981e4649fb01a6a6241ad36f1314e8e356c 100644 --- a/AMDiS/src/parallel/PetscSolverGlobalMatrix.cc +++ b/AMDiS/src/parallel/PetscSolverGlobalMatrix.cc @@ -72,8 +72,8 @@ namespace AMDiS { // === Create PETSc matrix with the computed nnz data structure. === MatCreateMPIAIJ(PETSC_COMM_WORLD, nRankRows, nRankRows, nOverallRows, nOverallRows, - 0, d_nnz, 0, o_nnz, &petscMatrix); - + 0, d_nnz, 0, o_nnz, &petscMatrix); + #if (DEBUG != 0) MSG("Fill petsc matrix 1 needed %.5f seconds\n", MPI::Wtime() - wtime); #endif @@ -150,7 +150,6 @@ namespace AMDiS { if (!zeroStartVector) KSPSetInitialGuessNonzero(solver, PETSC_TRUE); - // PETSc. KSPSolve(solver, petscRhsVec, petscSolVec);