[petsc-users] AOCreateBasic & AOCreateBasicIS

[Hui Zhang]

Is the following correct? 

From my understanding before, 

AOCreateBasic(MPI_Comm comm,PetscInt napp,const PetscInt myapp[],const PetscInt mypetsc[],AO *aoout)

should be called by all the processors in 'comm' and each involved processor provides only part of all the indices.
For example, 'comm' contains two processors proc0, 

on proc0:  napp= 3, myapp[]= {0,1,2}, mypetsc[]= {4,3,2}

on proc1:  napp= 2, myapp[]= {3,4}, mypetsc[]= {0,1}

so that union of myapp[] of proc0 and proc1 gives 0,..,4, and union of myapp[] of proc0 and proc1 also gives 0,..,4.

This seemed to work even if I fed AOApplicationToPetsc with indices on proc0 larger than 2.
But now I suspect I was wrong after reading 

	http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/AO/AOCreateBasic.html

The question is whether AOCreateBasic will do union of the input indices automatically, or the AO is only for the input indices?