Skip to content
Snippets Groups Projects
Commit 96f78d22 authored by Sebastien Blaise's avatar Sebastien Blaise
Browse files

Repaires system petsc bug I introduced

parent 104bf0dd
No related branches found
No related tags found
No related merge requests found
......@@ -206,10 +206,8 @@ int linearSystemPETScBlockDouble::systemSolve()
KSPSetOperators(_ksp, _a, _a, SAME_NONZERO_PATTERN);
else
KSPSetOperators(_ksp, _a, _a, DIFFERENT_NONZERO_PATTERN);
if (_matrixModified && _parameters["matrix_reuse"]!="same_matrix"){
MatAssemblyBegin(_a, MAT_FINAL_ASSEMBLY);
MatAssemblyEnd(_a, MAT_FINAL_ASSEMBLY);
}
MatAssemblyBegin(_a, MAT_FINAL_ASSEMBLY);
MatAssemblyEnd(_a, MAT_FINAL_ASSEMBLY);
_matrixModified=false;
VecAssemblyBegin(_b);
VecAssemblyEnd(_b);
......
......@@ -280,10 +280,8 @@ int linearSystemPETSc<scalar>::systemSolve()
_try(KSPSetOperators(_ksp, _a, _a, SAME_NONZERO_PATTERN));
else
_try(KSPSetOperators(_ksp, _a, _a, DIFFERENT_NONZERO_PATTERN));
if (_matrixModified && linearSystem<scalar>::_parameters["matrix_reuse"]!="same_matrix"){
_try(MatAssemblyBegin(_a, MAT_FINAL_ASSEMBLY));
_try(MatAssemblyEnd(_a, MAT_FINAL_ASSEMBLY));
}
_try(MatAssemblyBegin(_a, MAT_FINAL_ASSEMBLY));
_try(MatAssemblyEnd(_a, MAT_FINAL_ASSEMBLY));
_matrixModified=false;
/*MatInfo info;
MatGetInfo(_a, MAT_LOCAL, &info);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment