Applying Jacobian-Free Newton-Krylov methods in FeniCS

General
1

Dear. Researchers
I am solving a multiphysics problem using legacy fenics. Due to the ill-conditioning of Jacobain, I intend to utilize Jacobian-Free Newton-Krylov method for treating systems of nonlinear equations. However I do not know how to apply this methods inside fenics. I would be appreciated for any assistance to resolve this issue.
best regards.

2

Depending on what solver you would like for the linearized system, you can follow along the lines of: Matrix-free solvers in DOLFINx using PETSc — DOLFINx Python 0.11.0.dev0

3

I am so grateful. I will assess the page to find out the procedures.
best regards.