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;
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!");
18: /* create vector */
19: for (i = 0; i < n; i++) array[i] = i;
20: n = n - 1;
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: }