<div dir="ltr">Hello Petsc Users,<div><br clear="all"><div>I am solving a fortran framework to simulate phase-field fracture problem. Upto the failure point the code is working fine. But as the crack-phase starts to propagate, there is abrupt changes in the Stiffness Matrix, so determinant becomes nearly equal to zero.  Hence, the code crashes.</div><div> </div><div>I tried with mumps package with pc lu. Also tried, superlu, gmres and fgmres with jacobi, bjacobi, lu. But the problem still persist.</div><div><br></div><div>Can anyone suggest, Which solver and preconditioner should i use?<br></div><div><br></div><div>Thank you</div><div> </div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Avtar Singh</div><div>Research Scholar </div><div>Multiscale Mechanics and Multiphysics Lab</div><div>IIT Roorkee, Roorkee</div></div></div></div></div>
</div></div>