[petsc-users] On PCFIELDSPLIT and its implementation

[Edoardo alinovi]

Hello,

I have tried a different way to create the splitting:

    ui(1) = 0
    ui(2) = 1
    pi(1) = 2


    call ISCreateGeneral(PETSC_COMM_WORLD, 2, ui, PETSC_COPY_VALUES, isu,
ierr)
    call ISCreateGeneral(PETSC_COMM_WORLD, 1, pi, PETSC_COPY_VALUES, isp,
ierr)


    call PCFieldSplitSetIS(mypc, "0", isu, ierr)
    call PCFieldSplitSetIS(mypc, "1", isp, ierr)

However I get even worst erros!

*[0]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------*
*[0]PETSC ERROR: Nonconforming object sizes*
*[0]PETSC ERROR: Local column sizes 0 do not add up to total number of
columns 1*
*[0]PETSC ERROR: See https://petsc.org/release/faq/
<https://petsc.org/release/faq/> for trouble shooting.*


ISes looks like:

IS Object: 1 MPI process *---> isu*
  type: general
Number of indices in set 2
0 0
1 1
IS Object: 1 MPI process  * ---> isp*
  type: general
Number of indices in set 1
0 2

I really have no idea on how to deal with this thing... 😫 I guess I am not
setting ok the index of the splitting, but cannot really figure out what I
am doing wrong.

Any help is much appreciated :)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221110/7e462df8/attachment-0001.html>