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*9566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc,&argv,NULL,help)); 17*9566063dSJacob Faibussowitsch PetscCall(PetscOptionsGetInt(NULL,NULL,"-N",&N,NULL)); 18*9566063dSJacob Faibussowitsch PetscCall(PetscHSetICreate(&table)); 19c4762a1bSJed Brown 20c4762a1bSJed Brown /* The following line silences warnings from Clang Static Analyzer */ 21*9566063dSJacob Faibussowitsch PetscCall(PetscHSetIResize(table,0)); 22c4762a1bSJed Brown 23*9566063dSJacob Faibussowitsch PetscCall(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; 27*9566063dSJacob Faibussowitsch PetscCall(PetscHSetIQueryAdd(table, key, &flag)); 28c4762a1bSJed Brown } 29c4762a1bSJed Brown } 30*9566063dSJacob Faibussowitsch PetscCall(PetscTimeAdd(&t_add)); 31c4762a1bSJed Brown 32*9566063dSJacob Faibussowitsch PetscCall(PetscHSetIGetSize(table,&n)); 33c4762a1bSJed Brown 34*9566063dSJacob Faibussowitsch PetscCall(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; 38*9566063dSJacob Faibussowitsch PetscCall(PetscHSetIHas(table, key, &flag)); 39c4762a1bSJed Brown } 40c4762a1bSJed Brown } 41*9566063dSJacob Faibussowitsch PetscCall(PetscTimeAdd(&t_has)); 42c4762a1bSJed Brown 43*9566063dSJacob Faibussowitsch PetscCall(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; 47*9566063dSJacob Faibussowitsch PetscCall(PetscHSetIQueryDel(table, key, &flag)); 48c4762a1bSJed Brown } 49c4762a1bSJed Brown } 50*9566063dSJacob Faibussowitsch PetscCall(PetscTimeAdd(&t_del)); 51c4762a1bSJed Brown 52*9566063dSJacob Faibussowitsch PetscCall(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 54*9566063dSJacob Faibussowitsch PetscCall(PetscHSetIDestroy(&table)); 55*9566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 56b122ec5aSJacob 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