VecGetValues(Vec x,PetscInt ni,const PetscInt ix[],PetscScalar y[]) The documentation says ix is global index, isn't it the local index? Billy