SimVascular / svFSIplus

This repository contains a line-by-line conversion of the Fortran svFSI solver to C++.
Other
6 stars 19 forks source link

Possible bug in BICG linear solver #183

Open aabrown100-git opened 3 months ago

aabrown100-git commented 3 months ago

Description

When running the LV_Guccione_passive test case (https://github.com/SimVascular/svFSIplus/tree/main/tests/cases/struct/LV_Guccione_passive) for 1000 timesteps, the simulation encounters NaNs at timestep 62. When I change the linear solver from "BICG" to "GMRES", the simulation runs to completion.

Reproduction

Expected behavior

Both BICG and GMRES should produce similar results.

Additional context

No response

Code of Conduct