1 #ifdef PETSC_RCS_HEADER 2 static char vcid[] = "$Id: PetscMemcmp.c,v 1.9 1998/03/31 23:33:57 balay Exp bsmith $"; 3 #endif 4 5 #include "petsc.h" 6 7 #undef __FUNC__ 8 #define __FUNC__ "main" 9 int main( int argc, char **argv) 10 { 11 PLogDouble x, y, z; 12 Scalar A[10000], B[10000]; 13 int i,ierr; 14 15 PetscInitialize(&argc, &argv,0,0); 16 17 for (i=0; i<10000; i++) { 18 A[i] = i%61897; 19 B[i] = i%61897; 20 } 21 /* To take care of paging effects */ 22 PetscMemcmp(A,B,sizeof(Scalar)*0); 23 ierr = PetscGetTime(&x); CHKERRA(ierr); 24 25 ierr = PetscGetTime(&x); CHKERRA(ierr); 26 PetscMemcmp(A,B,sizeof(Scalar)*10000); 27 PetscMemcmp(A,B,sizeof(Scalar)*10000); 28 PetscMemcmp(A,B,sizeof(Scalar)*10000); 29 PetscMemcmp(A,B,sizeof(Scalar)*10000); 30 PetscMemcmp(A,B,sizeof(Scalar)*10000); 31 PetscMemcmp(A,B,sizeof(Scalar)*10000); 32 PetscMemcmp(A,B,sizeof(Scalar)*10000); 33 PetscMemcmp(A,B,sizeof(Scalar)*10000); 34 PetscMemcmp(A,B,sizeof(Scalar)*10000); 35 PetscMemcmp(A,B,sizeof(Scalar)*10000); 36 ierr = PetscGetTime(&y); CHKERRA(ierr); 37 PetscMemcmp(A,B,sizeof(Scalar)*0); 38 PetscMemcmp(A,B,sizeof(Scalar)*0); 39 PetscMemcmp(A,B,sizeof(Scalar)*0); 40 PetscMemcmp(A,B,sizeof(Scalar)*0); 41 PetscMemcmp(A,B,sizeof(Scalar)*0); 42 PetscMemcmp(A,B,sizeof(Scalar)*0); 43 PetscMemcmp(A,B,sizeof(Scalar)*0); 44 PetscMemcmp(A,B,sizeof(Scalar)*0); 45 PetscMemcmp(A,B,sizeof(Scalar)*0); 46 PetscMemcmp(A,B,sizeof(Scalar)*0); 47 ierr = PetscGetTime(&z); CHKERRA(ierr); 48 49 fprintf(stderr,"%s : \n","PetscMemcmp"); 50 fprintf(stderr," %-11s : %e sec\n","Latency",(z-y)/10.0); 51 fprintf(stderr," %-11s : %e sec\n","Per Scalar",(2*y-x-z)/100000); 52 53 PetscFinalize(); 54 PetscFunctionReturn(0); 55 } 56