<div dir="ltr">Hi,<div><br></div><div>I am using GRMES to solve the linear systems within Newton Raphson. </div><div>I am using finite element method to solve fracture problem with phase field method.(the formulation gives symmetric but indefinite hessian so can't use CG) </div><div><br></div><div>The total iteration number to solve one linear system is often very large( > 2000) even for a small sized problem(120K degree of freedom). I've tried different preconditioner, but it doesn't help a lot.  I have checked my hessian and right hand side, and sure they are correct. </div><div><br></div><div>I think it might stem from the great difference in element size(max/min=100) and damage stiffness around crack. </div><div><br></div><div>Anyone has idea to reduce the number of iteration?</div><div><br></div><div>Now, I get away with it by using LU, which is not ideal for large size problem (I will have to run it in the future)</div><div><br></div><div><br></div><div>Thanks,</div><div>Josh</div><div><br></div></div>