Finite-differences approximation of the Jacobian-vector product

nekStab / nekStab

An open source toolbox based on Krylov methods for highly parallel global linear stability analysis of incompressible flow problems for Nek5000.

BSD 3-Clause "New" or "Revised" License

11 stars 4 forks source link

Finite-differences approximation of the Jacobian-vector product #37

Open loiseaujc opened 1 year ago

loiseaujc commented 1 year ago

Second and fourth-order finite-difference approximations of the Jacobian-vector product have been implemented in the dedicated branch. So far, it has been validated for the two-dimensional cylinder flow, both for the Newton-Krylov solver and the stability analysis. It doesn't work on the thermosyphon though. The Krylov vectors generated do not seem to make any sense. I am not sure where it comes from.

loiseaujc commented 1 year ago

@ricardofrantz : I seem to remember you made some changes to the finite difference branch on your local copy. Is that correct ? If so, could you commit ? I am running a few other tests and would like to merge to the main soon.