Actual source code: ex21.c
1: /*$Id: ex21.c,v 1.13 2001/09/11 16:32:10 bsmith Exp $*/
3: static char help[] = "Tests VecMax() with index.\n\
4: -n <length> : vector length\n\n";
6: #include petscvec.h
7: #include petscsys.h
11: int main(int argc,char **argv)
12: {
13: int n = 5,ierr,idx;
14: PetscScalar value;
15: Vec x;
16: PetscRandom rand;
18: PetscInitialize(&argc,&argv,(char*)0,help);
19: PetscOptionsGetInt(PETSC_NULL,"-n",&n,PETSC_NULL);
21: /* create vector */
22: VecCreate(PETSC_COMM_WORLD,&x);
23: VecSetSizes(x,PETSC_DECIDE,n);
24: VecSetFromOptions(x);
26: PetscRandomCreate(PETSC_COMM_WORLD,RANDOM_DEFAULT_REAL,&rand);
27: VecSetRandom(rand,x);
28: PetscRandomDestroy(rand);
30: VecView(x,PETSC_VIEWER_STDOUT_WORLD);
32: VecMax(x,&idx,&value);
33: PetscPrintf(PETSC_COMM_WORLD,"Maximum value %g index %d\n",value,idx);
34: VecMin(x,&idx,&value);
35: PetscPrintf(PETSC_COMM_WORLD,"Minimum value %g index %d\n",value,idx);
37: VecDestroy(x);
39: PetscFinalize();
40: return 0;
41: }
42: