1d8e9fea7SSatish Balay 2c6db04a5SJed Brown #include <petscsys.h> 3*8563dfccSBarry Smith #include <petsctime.h> 4b0878937SSatish Balay 54a2ae208SSatish Balay #undef __FUNCT__ 64a2ae208SSatish Balay #define __FUNCT__ "main" 7b0878937SSatish Balay int main(int argc,char **argv) 8b0878937SSatish Balay { 9b0a32e0cSBarry Smith PetscLogDouble x,y,z; 10d3093643SSatish Balay int i,ierr; 11ea709b57SSatish Balay PetscScalar *A,*B; 12b0878937SSatish Balay 1377c4ece6SBarry Smith PetscInitialize(&argc,&argv,0,0); 14e24b481bSBarry Smith 1587828ca2SBarry Smith ierr = PetscMalloc(8000000*sizeof(PetscScalar),&A);CHKERRQ(ierr); 1687828ca2SBarry Smith ierr = PetscMalloc(8000000*sizeof(PetscScalar),&B);CHKERRQ(ierr); 17e24b481bSBarry Smith 18e24b481bSBarry Smith for (i=0; i<8000000; i++) { 19b0878937SSatish Balay A[i] = i%61897; 20b0878937SSatish Balay B[i] = i%61897; 21b0878937SSatish Balay } 22b0878937SSatish Balay /* To take care of paging effects */ 2387828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 24*8563dfccSBarry Smith ierr = PetscTime(&x);CHKERRQ(ierr); 25b0878937SSatish Balay 26*8563dfccSBarry Smith ierr = PetscTime(&x);CHKERRQ(ierr); 27e24b481bSBarry Smith /* 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); 3787828ca2SBarry Smith PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); */ 38e24b481bSBarry Smith { int j; 396f2b61bcSKarl Rupp for (j = 0; j<10; j++) 406f2b61bcSKarl Rupp for (i=0; i<8000000; i++) B[i] = A[i];} 41e24b481bSBarry Smith 42*8563dfccSBarry Smith ierr = PetscTime(&y);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); 5287828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 53*8563dfccSBarry Smith ierr = PetscTime(&z);CHKERRQ(ierr); 54b0878937SSatish Balay 5535d8aa7fSBarry Smith fprintf(stdout,"%s : \n","PetscMemcpy"); 56b4d8b9abSSatish Balay fprintf(stdout," %-15s : %e MB/s\n","Bandwidth",10.0*8*8/(y-x)); 57b4d8b9abSSatish Balay fprintf(stdout," %-15s : %e sec\n","Latency",(z-y)/10.0); 58b4d8b9abSSatish Balay fprintf(stdout," %-15s : %e sec\n","Per PetscScalar",(2*y-x-z)/8000000.0); 59b0878937SSatish Balay 60b4d8b9abSSatish Balay ierr = PetscFree(A);CHKERRQ(ierr); 61b4d8b9abSSatish Balay ierr = PetscFree(B);CHKERRQ(ierr); 62f3fe499bSBarry Smith ierr = PetscFinalize(); 633a40ed3dSBarry Smith PetscFunctionReturn(0); 64b0878937SSatish Balay } 65