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 lengthnn";
6: #include petscvec.h
7: #include petscsys.h
9: int main(int argc,char **argv)
10: {
11: int n = 5,ierr,idx;
12: PetscScalar value;
13: Vec x;
14: PetscRandom rand;
16: PetscInitialize(&argc,&argv,(char*)0,help);
17: PetscOptionsGetInt(PETSC_NULL,"-n",&n,PETSC_NULL);
19: /* create vector */
20: VecCreate(PETSC_COMM_WORLD,&x);
21: VecSetSizes(x,PETSC_DECIDE,n);
22: VecSetFromOptions(x);
24: PetscRandomCreate(PETSC_COMM_WORLD,RANDOM_DEFAULT_REAL,&rand);
25: VecSetRandom(rand,x);
26: PetscRandomDestroy(rand);
28: VecView(x,PETSC_VIEWER_STDOUT_WORLD);
30: VecMax(x,&idx,&value);
31: PetscPrintf(PETSC_COMM_WORLD,"Maximum value %g index %dn",value,idx);
32: VecMin(x,&idx,&value);
33: PetscPrintf(PETSC_COMM_WORLD,"Minimum value %g index %dn",value,idx);
35: VecDestroy(x);
37: PetscFinalize();
38: return 0;
39: }
40: