1d8e9fea7SSatish Balay 2d382aafbSBarry Smith #include "petscsys.h" 3b0878937SSatish Balay 44a2ae208SSatish Balay #undef __FUNCT__ 54a2ae208SSatish Balay #define __FUNCT__ "main" 6b0878937SSatish Balay int main(int argc,char **argv) 7b0878937SSatish Balay { 8b0a32e0cSBarry Smith PetscLogDouble x,y,z; 9d3093643SSatish Balay int i,ierr; 10ea709b57SSatish Balay PetscScalar *A,*B; 11b0878937SSatish Balay 1277c4ece6SBarry Smith PetscInitialize(&argc,&argv,0,0); 13e24b481bSBarry Smith 1487828ca2SBarry Smith ierr = PetscMalloc(8000000*sizeof(PetscScalar),&A);CHKERRQ(ierr); 1587828ca2SBarry Smith ierr = PetscMalloc(8000000*sizeof(PetscScalar),&B);CHKERRQ(ierr); 16e24b481bSBarry Smith 17e24b481bSBarry Smith for (i=0; i<8000000; i++) { 18b0878937SSatish Balay A[i] = i%61897; 19b0878937SSatish Balay B[i] = i%61897; 20b0878937SSatish Balay } 21b0878937SSatish Balay /* To take care of paging effects */ 2287828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 23ac355199SBarry Smith ierr = PetscGetTime(&x);CHKERRQ(ierr); 24b0878937SSatish Balay 25ac355199SBarry Smith ierr = PetscGetTime(&x);CHKERRQ(ierr); 26e24b481bSBarry Smith /* 2787828ca2SBarry Smith PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 2887828ca2SBarry Smith PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 2987828ca2SBarry Smith PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 3087828ca2SBarry Smith PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 3187828ca2SBarry Smith PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 3287828ca2SBarry Smith PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 3387828ca2SBarry Smith PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 3487828ca2SBarry Smith PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 3587828ca2SBarry Smith PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 3687828ca2SBarry Smith PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); */ 37e24b481bSBarry Smith { int j; 38e24b481bSBarry Smith for (j = 0; j<10; j++) { 39e24b481bSBarry Smith for (i=0; i<8000000; i++) { 40e24b481bSBarry Smith B[i] = A[i]; 41e24b481bSBarry Smith } 42e24b481bSBarry Smith }} 43e24b481bSBarry Smith 44ac355199SBarry Smith ierr = PetscGetTime(&y);CHKERRQ(ierr); 4587828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 4687828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 4787828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 4887828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 4987828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 5087828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 5187828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 5287828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 5387828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 5487828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 55ac355199SBarry Smith ierr = PetscGetTime(&z);CHKERRQ(ierr); 56b0878937SSatish Balay 5735d8aa7fSBarry Smith fprintf(stdout,"%s : \n","PetscMemcpy"); 58b4d8b9abSSatish Balay fprintf(stdout," %-15s : %e MB/s\n","Bandwidth",10.0*8*8/(y-x)); 59b4d8b9abSSatish Balay fprintf(stdout," %-15s : %e sec\n","Latency",(z-y)/10.0); 60b4d8b9abSSatish Balay fprintf(stdout," %-15s : %e sec\n","Per PetscScalar",(2*y-x-z)/8000000.0); 61b0878937SSatish Balay 62b4d8b9abSSatish Balay ierr = PetscFree(A);CHKERRQ(ierr); 63b4d8b9abSSatish Balay ierr = PetscFree(B);CHKERRQ(ierr); 64*f3fe499bSBarry Smith ierr = PetscFinalize(); 653a40ed3dSBarry Smith PetscFunctionReturn(0); 66b0878937SSatish Balay } 67