[petsc-users] On PCFIELDSPLIT and its implementation

[Matthew Knepley]

On Tue, Nov 8, 2022 at 12:05 PM Edoardo alinovi <edoardo.alinovi at gmail.com>
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/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221108/4338c51a/attachment.html>