Skip to content
GitLab

Fixing a bug in Process::setInitialConditions under PETSc

Merged wenqing requested to merge wenqing/ogs:fix_petsc into master

In Process::setInitialConditions, MathLib::LinAlg::setLocalAccessibleVector(x) is missed call before x is used for the computation the initial secondary variables. The bug leads to incorrect initial secondary variable assignment for some cases when PETSc is used.

This bug fixing fixes the issue addressed in #3217 (closed). As a verification, the example in #3217 (closed) is added as a benchmark. The new PETSc related benchmark use the existing reference file. The following figures show the comparison of the results solved by Eigen and PETSc linear solver, respectively: TH1

This MR closes #3217 (closed).

  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