Commit 29045c03 authored by Thomas Witkowski's avatar Thomas Witkowski

Work on more general and faster parallelization.

parent cf3fc63f
......@@ -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());
......
......@@ -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) {
......
......@@ -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);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment