1 2 static char help[] = "Tests taking part of existing array to create a new vector.\n\n"; 3 4 #include <petscvec.h> 5 6 int main(int argc,char **argv) 7 { 8 PetscMPIInt size; 9 PetscInt n = 10,i; 10 PetscScalar array[10]; 11 Vec x; 12 13 PetscFunctionBeginUser; 14 PetscCall(PetscInitialize(&argc,&argv,(char*)0,help)); 15 PetscCallMPI(MPI_Comm_size(PETSC_COMM_WORLD,&size)); 16 PetscCheck(size == 1,PETSC_COMM_WORLD,PETSC_ERR_WRONG_MPI_SIZE,"This is a uniprocessor example only!"); 17 18 /* create vector */ 19 for (i=0; i<n; i++) array[i] = i; 20 n = n-1; 21 22 PetscCall(VecCreateSeqWithArray(PETSC_COMM_SELF,1,n,array+1,&x)); 23 PetscCall(VecView(x,PETSC_VIEWER_STDOUT_SELF)); 24 PetscCall(VecDestroy(&x)); 25 PetscCall(PetscFinalize()); 26 return 0; 27 } 28