1*ea709b57SSatish Balay /*$Id: PetscMemcpy.c,v 1.22 2001/08/06 21:19:06 bsmith Exp balay $*/ 2d8e9fea7SSatish Balay 3b0878937SSatish Balay #include "petsc.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; 11*ea709b57SSatish 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); 24ac355199SBarry Smith ierr = PetscGetTime(&x);CHKERRQ(ierr); 25b0878937SSatish Balay 26ac355199SBarry Smith ierr = PetscGetTime(&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; 39e24b481bSBarry Smith for (j = 0; j<10; j++) { 40e24b481bSBarry Smith for (i=0; i<8000000; i++) { 41e24b481bSBarry Smith B[i] = A[i]; 42e24b481bSBarry Smith } 43e24b481bSBarry Smith }} 44e24b481bSBarry Smith 45ac355199SBarry Smith ierr = PetscGetTime(&y);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); 5587828ca2SBarry Smith ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 56ac355199SBarry Smith ierr = PetscGetTime(&z);CHKERRQ(ierr); 57b0878937SSatish Balay 5835d8aa7fSBarry Smith fprintf(stdout,"%s : \n","PetscMemcpy"); 5935d8aa7fSBarry Smith fprintf(stdout," %-11s : %e MB/s\n","Bandwidth",10.0*8*8/(y-x)); 6035d8aa7fSBarry Smith fprintf(stdout," %-11s : %e sec\n","Latency",(z-y)/10.0); 61*ea709b57SSatish Balay fprintf(stdout," %-11s : %e sec\n","Per PetscScalar",(2*y-x-z)/8000000.0); 62b0878937SSatish Balay 632b7fea2aSSatish Balay ierr = PetscFinalize();CHKERRQ(ierr); 643a40ed3dSBarry Smith PetscFunctionReturn(0); 65b0878937SSatish Balay } 66