<div dir="ltr">Hi all!<div><br></div><div>Recently I faced with necessity of iterative solver use for magnetodynamics. I started with this case:</div><div><br></div><div><a href="https://gitlab.onelab.info/doc/models/wikis/Magnetodynamics-with-cohomology-conditions">https://gitlab.onelab.info/doc/models/wikis/Magnetodynamics-with-cohomology-conditions</a><br></div><div><br></div><div>I used getdp with petsc, and I could not get to convergence.</div><div>I have tried:</div><div><br></div><div>getdp <a href="http://indheat.pro">indheat.pro</a> -msh indheat.msh -solve MagDynAVComplex -solver solver.par<br></div><div><br></div><div>solver.par:</div><div><br></div><div><div>-pc_type icc</div><div>-ksp_type gmres</div><div>-ksp_max_it 1000<br></div></div><div><div>-ksp_rtol 1e-12</div></div><div><br></div><div>I have tried also other combinations of preconditioner and solver (icc + cg, ilu + gmres, bjacobi + gmres, bjacobi + cg ...).</div><div>So far the best convergence was with icc+gmres, which was only to 1e-9 after 1000 iterations, but the results were unphysical. Direct solver gives proper results.</div><div><br></div><div>Do you have any clue why is there such behavior and is there a way to get iterative solver work?</div><div><br></div><div>Best Regards, Vadims</div></div>