creating parallel vectors

Barry Smith bsmith at mcs.anl.gov
Thu Apr 2 22:04:19 CDT 2009


On Apr 2, 2009, at 9:19 PM, Khan, Irfan wrote:

> Thank you for the prompt help. It have benefited me a lot.
>
> I have another question about creating parallel vectors. Is it  
> possible to create parallel vectors by specifying the global indices  
> of the entries?
> For instance in
> VecCreateGhost(MPI_Comm comm,PetscInt n,PetscInt N,PetscInt  
> nghost,const PetscInt global_ordering[],Vec *vv);
>
> is it possible to give the global ordering of all the local+ghost  
> entries of the vector through "global_ordering". Thus the vector  
> will not create a new petsc ordering.

   No, indexing for PETSc parallel vectors is ALWAYS with the  
numbering where the first process has 0 to n1-1 second process n1 to  
n2_1 etc.

    Barry

>
>
>
> Thank you
> Irfan



More information about the petsc-users mailing list