[petsc-users] how to reate an IS of size 0 on all the other processes and a seq Vec of the size zero on all other processes

[Barry Smith]

  just do it.

    MPI_Comm_rank(comm,&rank);

    if (rank) {
        ISCreateGeneral(comm,0….)
        VecCreateSeq(comm,0….
     } else {
        ISCreateGeneral(comm,N,…..)
        VecCreateSeq(comm,N…..)
     }

    Barry

Note for scalability reasons this is generally not a wise thing to do.

  
On Aug 22, 2013, at 5:39 PM, 丁老师 <ztdepyahoo at 163.com> wrote:

> 
> 
>