xref: /petsc/src/sys/tests/ex42.c (revision b122ec5aa1bd4469eb4e0673542fb7de3f411254)
1c4762a1bSJed Brown static char help[] = "Test scalability of PetscHSetI hash set.\n\n";
2c4762a1bSJed Brown 
3c4762a1bSJed Brown #include <petscsys.h>
4c4762a1bSJed Brown #include <petsctime.h>
5c4762a1bSJed Brown #include <petsc/private/hashseti.h>
6c4762a1bSJed Brown 
7c4762a1bSJed Brown int main(int argc, char **argv)
8c4762a1bSJed Brown {
9c4762a1bSJed Brown   PetscHSetI     table;
10c4762a1bSJed Brown   PetscInt       N = 0, i, j, n;
11c4762a1bSJed Brown   PetscBool      flag;
12c4762a1bSJed Brown   PetscLogDouble t_add = 0;
13c4762a1bSJed Brown   PetscLogDouble t_has = 0;
14c4762a1bSJed Brown   PetscLogDouble t_del = 0;
15c4762a1bSJed Brown 
16*b122ec5aSJacob Faibussowitsch   CHKERRQ(PetscInitialize(&argc,&argv,NULL,help));
175f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscOptionsGetInt(NULL,NULL,"-N",&N,NULL));
185f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscHSetICreate(&table));
19c4762a1bSJed Brown 
20c4762a1bSJed Brown   /* The following line silences warnings from Clang Static Analyzer */
215f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscHSetIResize(table,0));
22c4762a1bSJed Brown 
235f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscTimeSubtract(&t_add));
24c4762a1bSJed Brown   for (i = 0; i < N; ++i) {
25c4762a1bSJed Brown     for (j = 0; j < N; ++j) {
26c4762a1bSJed Brown       PetscInt key = i+j*N;
275f80ce2aSJacob Faibussowitsch       CHKERRQ(PetscHSetIQueryAdd(table, key, &flag));
28c4762a1bSJed Brown     }
29c4762a1bSJed Brown   }
305f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscTimeAdd(&t_add));
31c4762a1bSJed Brown 
325f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscHSetIGetSize(table,&n));
33c4762a1bSJed Brown 
345f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscTimeSubtract(&t_has));
35c4762a1bSJed Brown   for (i = 0; i < N; ++i) {
36c4762a1bSJed Brown     for (j = 0; j < N; ++j) {
37c4762a1bSJed Brown       PetscInt key = i+j*N;
385f80ce2aSJacob Faibussowitsch       CHKERRQ(PetscHSetIHas(table, key, &flag));
39c4762a1bSJed Brown     }
40c4762a1bSJed Brown   }
415f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscTimeAdd(&t_has));
42c4762a1bSJed Brown 
435f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscTimeSubtract(&t_del));
44c4762a1bSJed Brown   for (i = 0; i < N; ++i) {
45c4762a1bSJed Brown     for (j = 0; j < N; ++j) {
46c4762a1bSJed Brown       PetscInt key = i+j*N;
475f80ce2aSJacob Faibussowitsch       CHKERRQ(PetscHSetIQueryDel(table, key, &flag));
48c4762a1bSJed Brown     }
49c4762a1bSJed Brown   }
505f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscTimeAdd(&t_del));
51c4762a1bSJed Brown 
525f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscPrintf(PETSC_COMM_WORLD,"N = %" PetscInt_FMT " - table size: %" PetscInt_FMT ", add: %g, has: %g, del: %g\n",N,n,t_add,t_has,t_del));
53c4762a1bSJed Brown 
545f80ce2aSJacob Faibussowitsch   CHKERRQ(PetscHSetIDestroy(&table));
55*b122ec5aSJacob Faibussowitsch   CHKERRQ(PetscFinalize());
56*b122ec5aSJacob Faibussowitsch   return 0;
57c4762a1bSJed Brown }
58c4762a1bSJed Brown 
59c4762a1bSJed Brown /*TEST
60c4762a1bSJed Brown 
61c4762a1bSJed Brown    test:
62c4762a1bSJed Brown      args: -N 32
63c4762a1bSJed Brown 
64c4762a1bSJed Brown TEST*/
65