[petsc-users] On PCFIELDSPLIT and its implementation

[Matthew Knepley]

On Wed, Nov 9, 2022 at 9:17 AM Edoardo alinovi <edoardo.alinovi at gmail.com>
wrote:

> Thanks,
>
> the stuff I am doing is within my code, so I am not sure you can reproduce
> it.
>

How about just making a small code that fills those nonzeros with 1s. We
just want to figure out why your
sparsity pattern is not working. We have lots of working PCFIELDSPLIT
examples, but I don't think it would
be useful to show you those.

  Thanks,

     Matt


> I am just doing this:
>
>             call PCSetType(mypc, PCFIELDSPLIT, ierr)
>
>             call PCFieldSplitSetBlockSize(mypc, 4-bdim, ierr)
>
>                 ufields(1) = 0
>                 ufields(2) = 1
>                 pfields(1) = 2
>                 call PCFieldSplitSetFields(mypc, "u", 2, ufields, ufields,
> ierr)
>                 call PCFieldSplitSetFields(mypc, "p", 1, pfields, pfields,
> ierr)
>
> On an MPIBAIJ matrix with bs = 3.
>
>
>
>
>

-- 
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/20221109/f34f4752/attachment.html>