KVSlab / turtleFSI

Monolithic Fluid-Structure Interaction (FSI) solver
https://turtlefsi2.readthedocs.io/en/latest/
GNU General Public License v3.0
63 stars 23 forks source link

added check when reading checkpoint #28

Closed keiyamamo closed 1 year ago

keiyamamo commented 1 year ago

This PR adds checking system when restarting from the previous simulation. I noticed that you get

Error:   Unable to solve linear system using PETSc Krylov solver.
Reason:  Solution failed to converge in 0 iterations (PETSc reason DIVERGED_PC_FAILED, residual norm ||r|| = 0.000000e+00).
 Where:   This error was encountered inside PETScKrylovSolver.cpp.

even when the simulation actually diverged. With the current implementation, it is not possible to see if this error comes from the actual divergence or mumps. By adding this check, we can at least tell if our checkpoint contains reasonable value or not.