1 2 #include <petscsys.h> 3 #include <petsctime.h> 4 5 #undef __FUNCT__ 6 #define __FUNCT__ "main" 7 int main(int argc,char **argv) 8 { 9 PetscLogDouble x,y,z; 10 PetscScalar A[10000],B[10000]; 11 PetscInt i; 12 PetscErrorCode ierr; 13 PetscBool flg; 14 15 ierr = PetscInitialize(&argc,&argv,0,0);if (ierr) return ierr; 16 for (i=0; i<10000; i++) { 17 A[i] = i%61897; 18 B[i] = i%61897; 19 } 20 /* To take care of paging effects */ 21 PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg); 22 ierr = PetscTime(&x);CHKERRQ(ierr); 23 24 ierr = PetscTime(&x);CHKERRQ(ierr); 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 PetscMemcmp(A,B,sizeof(PetscScalar)*10000,&flg); 35 ierr = PetscTime(&y);CHKERRQ(ierr); 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 PetscMemcmp(A,B,sizeof(PetscScalar)*0,&flg); 46 ierr = PetscTime(&z);CHKERRQ(ierr); 47 48 fprintf(stdout,"%s : \n","PetscMemcmp"); 49 fprintf(stdout," %-15s : %e sec\n","Latency",(z-y)/10.0); 50 fprintf(stdout," %-15s : %e sec\n","Per PetscScalar",(2*y-x-z)/100000); 51 52 ierr = PetscFinalize(); 53 return ierr; 54 } 55