Re: [petsc-users] Matrix-Free J*v in PETSc

Thank you very much for your help and the examples provided. I was able to implement the Jacobian-times-vector (Jv) approach as recommended, and it works great!. Art El mar, 8 jul 2025 a las 13:36, Zhang, Hong () escribió: > For a python example, please take a look at > src/binding/pe

[petsc-users] Matrix-Free J*v in PETSc

matrix-free Newton-Krylov method to approximate the Jacobian and have adjusted the tolerances to achieve convergence, as recommended by Barry. In that case, I obtain similar integration times with scikits.odes CVODE without using the Jacobian times vector. Best regards, Art

Re: [petsc-users] High Memory Usage with PETSc BDF Solver Compared to scikits.odes CVODE

pport thanks Art El lun, 16 jun 2025 a las 20:51, Barry Smith () escribió: > >Can you please post the entire code so we can run it ourselves to > reproduce the problem. > > I suspect the code maybe using a dense matrix to represent the > Jacobian. This ill require

[petsc-users] High Memory Usage with PETSc BDF Solver Compared to scikits.odes CVODE

lar results to the rk solver in PETSC. I suspect the issue might be related to the way the Jacobian is built in PETSC, but scikits.odes.sundials.cvode works perfectly without requiring the Jacobian. I would greatly appreciate any suggestions or examples on how to properly set up the BDF solver with PETSc. Thanks Art