Actual source code: PetscMemcmp.c

  1: #include <petscsys.h>
  2: #include <petsctime.h>

  4: int main(int argc,char **argv)
  5: {
  6:   PetscLogDouble x,y,z;
  7:   PetscScalar    A[10000],B[10000];
  8:   PetscInt       i;
  9:   PetscBool      flg;

 11:   PetscCall(PetscInitialize(&argc,&argv,0,0));
 12:   for (i=0; i<10000; i++) {
 13:     A[i] = i%61897;
 14:     B[i] = i%61897;
 15:   }
 16:   /* To take care of paging effects */
 17:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 18:   PetscCall(PetscTime(&x));

 20:   PetscCall(PetscTime(&x));
 21:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 22:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 23:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 24:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 25:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 26:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 27:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 28:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 29:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 30:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 31:   PetscCall(PetscTime(&y));
 32:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 33:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 34:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 35:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 36:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 37:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 38:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 39:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 40:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 41:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 42:   PetscCall(PetscTime(&z));

 44:   fprintf(stdout,"%s : \n","PetscMemcmp");
 45:   fprintf(stdout,"    %-15s : %e sec\n","Latency",(z-y)/10.0);
 46:   fprintf(stdout,"    %-15s : %e sec\n","Per PetscScalar",(2*y-x-z)/100000);

 48:   PetscCall(PetscFinalize());
 49:   return 0;
 50: }