Actual source code: ex40.c


  2: static char help[] = "Tests taking part of existing array to create a new vector.\n\n";

  4: #include <petscvec.h>

  6: int main(int argc, char **argv)
  7: {
  8:   PetscMPIInt size;
  9:   PetscInt    n = 10, i;
 10:   PetscScalar array[10];
 11:   Vec         x;

 14:   PetscInitialize(&argc, &argv, (char *)0, help);
 15:   MPI_Comm_size(PETSC_COMM_WORLD, &size);

 18:   /* create vector */
 19:   for (i = 0; i < n; i++) array[i] = i;
 20:   n = n - 1;

 22:   VecCreateSeqWithArray(PETSC_COMM_SELF, 1, n, array + 1, &x);
 23:   VecView(x, PETSC_VIEWER_STDOUT_SELF);
 24:   VecDestroy(&x);
 25:   PetscFinalize();
 26:   return 0;
 27: }