1d8e9fea7SSatish Balay 2c6db04a5SJed Brown #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; 38*6f2b61bcSKarl Rupp for (j = 0; j<10; j++) 39*6f2b61bcSKarl Rupp for (i=0; i<8000000; i++) B[i] = A[i];} 40e24b481bSBarry Smith 41ac355199SBarry Smith ierr = PetscGetTime(&y);CHKERRQ(ierr); 4287828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 4387828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 4487828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);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); 52ac355199SBarry Smith ierr = PetscGetTime(&z);CHKERRQ(ierr); 53b0878937SSatish Balay 5435d8aa7fSBarry Smith fprintf(stdout,"%s : \n","PetscMemcpy"); 55b4d8b9abSSatish Balay fprintf(stdout," %-15s : %e MB/s\n","Bandwidth",10.0*8*8/(y-x)); 56b4d8b9abSSatish Balay fprintf(stdout," %-15s : %e sec\n","Latency",(z-y)/10.0); 57b4d8b9abSSatish Balay fprintf(stdout," %-15s : %e sec\n","Per PetscScalar",(2*y-x-z)/8000000.0); 58b0878937SSatish Balay 59b4d8b9abSSatish Balay ierr = PetscFree(A);CHKERRQ(ierr); 60b4d8b9abSSatish Balay ierr = PetscFree(B);CHKERRQ(ierr); 61f3fe499bSBarry Smith ierr = PetscFinalize(); 623a40ed3dSBarry Smith PetscFunctionReturn(0); 63b0878937SSatish Balay } 64