On Tue, Nov 8, 2022 at 12:05 PM Edoardo alinovi <[email protected]> wrote:

Hello Guys,

Thanks to your suggestions on the block matrices, my fully coupled solver is proceeding very well!

I am now about to take advantage of the block structure of the matrix using PCFIELDSPLIT. I have learned a bit from the user manual and followed with interest this discussion in the mailing list: https://lists.mcs.anl.gov/pipermail/petsc-users/2015-February/024154.html which is actually the exact same situation I am in, so I guess most of the command line options will be copy and paste.

I would like however to code them in fortran, as I usually provide some default implementation alongside the command line options.

While coding some of the options I got an error here in PCFieldSplitSetFields() which looks to be undefined. I am importing petscksp, do I need to import something else maybe?

Since it uses arrays, we will have to write the Fortran wrapper by hand. I will see if I can do it soon. Thanks, Matt

Thank you!

-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>