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:
This MR closes #3217 (closed).
-
Feature description was added to the changelog -
Tests covering your feature were added? -
Any new feature or behavior change was documented?