Commit d32efd78 authored by Thomas Witkowski's avatar Thomas Witkowski
Browse files

Small bugfix for periodic boundary conditions.

parent 8f7b82ff
...@@ -760,14 +760,14 @@ namespace AMDiS { ...@@ -760,14 +760,14 @@ namespace AMDiS {
ERROR_EXIT("Not yet implemented!\n"); ERROR_EXIT("Not yet implemented!\n");
} }
// TODO: ExitMatrix should be called after finishInsertion!
if (assembleMatrix && matrix->getBoundaryManager())
matrix->getBoundaryManager()->exitMatrix(matrix);
assembledMatrix[i][j] = true; assembledMatrix[i][j] = true;
if (assembleMatrix) if (assembleMatrix)
matrix->finishInsertion(); matrix->finishInsertion();
// TODO: ExitMatrix should be called after finishInsertion!
if (assembleMatrix && matrix->getBoundaryManager())
matrix->getBoundaryManager()->exitMatrix(matrix);
if (matrix) if (matrix)
nnz += matrix->getBaseMatrix().nnz(); nnz += matrix->getBaseMatrix().nnz();
...@@ -966,7 +966,7 @@ namespace AMDiS { ...@@ -966,7 +966,7 @@ namespace AMDiS {
if (systemMatrix && (*systemMatrix)[row][col]) if (systemMatrix && (*systemMatrix)[row][col])
(*systemMatrix)[row][col]->getBoundaryManager()->addBoundaryCondition(periodic); (*systemMatrix)[row][col]->getBoundaryManager()->addBoundaryCondition(periodic);
if (rhs) if (rhs && row == col)
rhs->getDOFVector(row)->getBoundaryManager()->addBoundaryCondition(periodic); rhs->getDOFVector(row)->getBoundaryManager()->addBoundaryCondition(periodic);
} }
......
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