xref: /petsc/src/sys/tests/ex58.c (revision 3c859ba3a04a72e8efcb87bc7ffc046a6cbab413)
1 
2 static char help[] = "Tests PetscGlobalMinMax\n\n";
3 
4 #include <petscsys.h>
5 
6 int main(int argc,char **argv)
7 {
8   PetscErrorCode ierr;
9   PetscMPIInt    size,rank;
10   PetscInt       li[2],gi[2] = {-1, -1};
11   PetscReal      lr[2],gr[2] = {-1., -1.};
12 
13   ierr = PetscInitialize(&argc,&argv,(char*)0,help);if (ierr) return ierr;
14   ierr = MPI_Comm_size(PETSC_COMM_WORLD,&size);CHKERRMPI(ierr);
15   ierr = MPI_Comm_rank(PETSC_COMM_WORLD,&rank);CHKERRMPI(ierr);
16 
17   li[0] = 4 + rank;
18   li[1] = -3 + size - rank;
19   ierr = PetscGlobalMinMaxInt(PETSC_COMM_WORLD,li,gi);CHKERRQ(ierr);
20   if (gi[0] != 4 || gi[1] != -3+size) { ierr = PetscPrintf(PETSC_COMM_SELF,"1) Error MIN/MAX %" PetscInt_FMT " %" PetscInt_FMT "\n",gi[0],gi[1]);CHKERRQ(ierr); }
21   ierr = PetscGlobalMinMaxInt(PETSC_COMM_WORLD,li,li);CHKERRQ(ierr);
22   if (li[0] != gi[0] || li[1] != gi[1]) { ierr = PetscPrintf(PETSC_COMM_SELF,"2) Error MIN/MAX %" PetscInt_FMT " %" PetscInt_FMT "\n",li[0],li[1]);CHKERRQ(ierr); }
23 
24   if (rank == 0) {
25     li[0] = PETSC_MAX_INT;
26     li[1] = PETSC_MIN_INT;
27   } else if (rank == 1) {
28     li[0] = PETSC_MIN_INT;
29     li[1] = PETSC_MAX_INT;
30   }
31 
32   ierr = PetscGlobalMinMaxInt(PETSC_COMM_WORLD,li,gi);CHKERRQ(ierr);
33   if (gi[0] > li[0] || gi[1] < li[1]) { ierr = PetscPrintf(PETSC_COMM_SELF,"3) Error MIN/MAX %" PetscInt_FMT " %" PetscInt_FMT "\n",gi[0],gi[1]);CHKERRQ(ierr); }
34 
35   lr[0] = 4.0 + rank;
36   lr[1] = -3.0 + size - rank;
37   ierr = PetscGlobalMinMaxReal(PETSC_COMM_WORLD,lr,gr);CHKERRQ(ierr);
38   if (gr[0] != 4.0 || gr[1] != -3.0+size) { ierr = PetscPrintf(PETSC_COMM_SELF,"4) Error MIN/MAX %g %g\n",(double)gr[0],(double)gr[1]);CHKERRQ(ierr); }
39   ierr = PetscGlobalMinMaxReal(PETSC_COMM_WORLD,lr,lr);CHKERRQ(ierr);
40   if (lr[0] != gr[0] || lr[1] != gr[1]) { ierr = PetscPrintf(PETSC_COMM_SELF,"5) Error MIN/MAX %g %g\n",(double)lr[0],(double)li[1]);CHKERRQ(ierr); }
41 
42   ierr = PetscFinalize();
43   return ierr;
44 }
45 
46 /*TEST
47 
48    test:
49      output_file: output/ex58_1.out
50 
51    test:
52      suffix: 2
53      output_file: output/ex58_1.out
54      nsize: 2
55 
56 TEST*/
57