1 /*$Id: PetscMemcmp.c,v 1.17 2001/01/23 20:58:37 balay Exp balay $*/ 2 3 #include "petsc.h" 4 5 #undef __FUNCT__ 6 #define __FUNCT__ "main" 7 int main(int argc,char **argv) 8 { 9 PetscLogDouble x,y,z; 10 Scalar A[10000],B[10000]; 11 int i,ierr; 12 13 PetscInitialize(&argc,&argv,0,0); 14 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(Scalar)*0); 21 ierr = PetscGetTime(&x);CHKERRQ(ierr); 22 23 ierr = PetscGetTime(&x);CHKERRQ(ierr); 24 PetscMemcmp(A,B,sizeof(Scalar)*10000); 25 PetscMemcmp(A,B,sizeof(Scalar)*10000); 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 ierr = PetscGetTime(&y);CHKERRQ(ierr); 35 PetscMemcmp(A,B,sizeof(Scalar)*0); 36 PetscMemcmp(A,B,sizeof(Scalar)*0); 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 ierr = PetscGetTime(&z);CHKERRQ(ierr); 46 47 fprintf(stdout,"%s : \n","PetscMemcmp"); 48 fprintf(stdout," %-11s : %e sec\n","Latency",(z-y)/10.0); 49 fprintf(stdout," %-11s : %e sec\n","Per Scalar",(2*y-x-z)/100000); 50 51 ierr = PetscFinalize();CHKERRQ(ierr); 52 PetscFunctionReturn(0); 53 } 54