#include "petscao.h" PetscErrorCode AOCreateMappingIS(IS isapp, IS ispetsc, AO *aoout)
comm | - MPI communicator that is to share AO | |
isapp | - index set that defines an ordering | |
ispetsc | - index set that defines another ordering, maybe NULL for identity IS |
Notes: the index sets isapp and ispetsc need NOT contain the all the integers 0 to N-1, that is there CAN be "holes" in the indices. Use AOCreateBasic() or AOCreateBasicIS() if they do not have holes for better performance.