On Wed, Nov 9, 2022 at 9:17 AM Edoardo alinovi <edoardo.alin...@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/>