1d8e9fea7SSatish Balay 2c6db04a5SJed Brown #include <petscsys.h> 38563dfccSBarry Smith #include <petsctime.h> 4b0878937SSatish Balay 5b0878937SSatish Balay int main(int argc,char **argv) 6b0878937SSatish Balay { 7580bdb30SBarry Smith PetscLogDouble x,y; 8ea709b57SSatish Balay PetscScalar *A,*B; 9b0878937SSatish Balay 10*9566063dSJacob Faibussowitsch PetscCall(PetscInitialize(&argc,&argv,0,0)); 11*9566063dSJacob Faibussowitsch PetscCall(PetscCalloc1(8000000,&A)); 12*9566063dSJacob Faibussowitsch PetscCall(PetscMalloc1(8000000,&B)); 13e24b481bSBarry Smith 14e24b481bSBarry Smith for (i=0; i<8000000; i++) { 15b0878937SSatish Balay A[i] = i%61897; 16b0878937SSatish Balay B[i] = i%61897; 17b0878937SSatish Balay } 18b0878937SSatish Balay /* To take care of paging effects */ 19*9566063dSJacob Faibussowitsch PetscCall(PetscArraycpy(A,B,8000000)); 20*9566063dSJacob Faibussowitsch PetscCall(PetscTime(&x)); 21*9566063dSJacob Faibussowitsch PetscCall(PetscArraycpy(A,B,8000000)); 22*9566063dSJacob Faibussowitsch PetscCall(PetscTime(&x)); 23b0878937SSatish Balay 2435d8aa7fSBarry Smith fprintf(stdout,"%s : \n","PetscMemcpy"); 25b4d8b9abSSatish Balay fprintf(stdout," %-15s : %e MB/s\n","Bandwidth",10.0*8*8/(y-x)); 26b0878937SSatish Balay 27*9566063dSJacob Faibussowitsch PetscCall(PetscFree(A)); 28*9566063dSJacob Faibussowitsch PetscCall(PetscFree(B)); 29*9566063dSJacob Faibussowitsch PetscCall(PetscFinalize()); 30b122ec5aSJacob Faibussowitsch return 0; 31b0878937SSatish Balay } 32