[petsc-users] Using matrix-free with KSP

[Zhang, Hong]

You can refer to the example src/ts/tutorials/advection-diffusion-reaction/ex5adj_mf.c
which uses a matrix-free approach for a 2D multi-component reaction-diffusion PDE.
In particular, you only need to look at MyIMatMult() and ignore the transposed version
which is intended for adjoint calculations.

Hong (Mr.)

On Sep 6, 2022, at 10:00 PM, Tu, Jiannan <Jiannan_Tu at uml.edu<mailto:Jiannan_Tu at uml.edu>> wrote:

I am using TS IMEX to solve a large DAE system. The DAE is obtained by applying finite FV method to 3-D multi-specie ion/neutral fluids equations with magnetic induction equation. The Jacobian for stiff part is formed by using MatSetValuesStencil(). The Jacobian matrix is very sparse, no more than 10 non-zeros on each row. MatSetValuesStencil requires local to global mapping by calling ISLocaltoGlobalMapping(). Could you please instruct me how to use local to global mapping?

I also tried using DMCreateMatrix() to create the Jacobian. While local to global mapping is not necessary, the matrix takes too much memory and requires 64-bit indices. I would prefer to take the advantage of sparsity of the Jacobian, pre-allocate the matrix to use as less as possible memory so that the code can be run on a multi-core desktop.

Thank you very much for your advice.

Jiannan

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220907/2e860064/attachment.html>