Actual source code: PetscMemcmp.c

petsc-3.3-p7 2013-05-11
  2: #include <petscsys.h>

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

 13:   PetscInitialize(&argc,&argv,0,0);

 15:   for (i=0; i<10000; i++) {
 16:     A[i] = i%61897;
 17:     B[i] = i%61897;
 18:   }
 19:   /* To take care of paging effects */
 20:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 21:   PetscGetTime(&x);

 23:   PetscGetTime(&x);
 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:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 32:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 33:   PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg);
 34:   PetscGetTime(&y);
 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:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 43:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 44:   PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg);
 45:   PetscGetTime(&z);

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

 51:   PetscFinalize();
 52:   return(0);
 53: }