Skip to content

Draft: [MathLib/PETSc] Call MatZeroRowsColumns instead of MatZeroRows for one compute thread

wenqing requested to merge wenqing/ogs:update_PETScMatrix into master

OGS with PETSc enabled can be run in the serial manner, i.e. with one compute thread by MPI launcher as mpiexec -n 1 ogs .... To prescribe the Dirichlet boundary condition, MatZeroRows is called to zero the off-diagonal entries of the rows corresponding to the known entries of the solution vector. MatZeroRows leads to an unsymmetrical matrix but with less communication load. For the computation with one compute thread, this MR replaces MatZeroRows with MatZeroRowsColumns, which keeps matrix symmetry for symmetry matrix.

  1. Feature description was added to the changelog
  2. Tests covering your feature were added?
  3. Any new feature or behavior change was documented?
Edited by wenqing

Merge request reports