1a5eb4965SSatish Balay #ifdef PETSC_RCS_HEADER 2*5a655dc6SBarry Smith static char vcid[] = "$Id: PetscMemcpy.c,v 1.11 1998/11/20 15:31:19 bsmith Exp bsmith $"; 3d8e9fea7SSatish Balay #endif 4d8e9fea7SSatish Balay 5b0878937SSatish Balay #include "petsc.h" 6b0878937SSatish Balay 7*5a655dc6SBarry Smith #undef __FUNC__ 8*5a655dc6SBarry Smith #define __FUNC__ "main" 9b0878937SSatish Balay int main( int argc, char **argv) 10b0878937SSatish Balay { 1147794344SBarry Smith PLogDouble x, y, z; 12d3093643SSatish Balay int i,ierr; 13e24b481bSBarry Smith Scalar *A, *B; 14b0878937SSatish Balay 1577c4ece6SBarry Smith PetscInitialize(&argc, &argv,0,0); 16e24b481bSBarry Smith 17e24b481bSBarry Smith A = (Scalar *) PetscMalloc(8000000*sizeof(Scalar)); 18e24b481bSBarry Smith B = (Scalar *) PetscMalloc(8000000*sizeof(Scalar)); 19e24b481bSBarry Smith 20e24b481bSBarry Smith for (i=0; i<8000000; i++) { 21b0878937SSatish Balay A[i] = i%61897; 22b0878937SSatish Balay B[i] = i%61897; 23b0878937SSatish Balay } 24b0878937SSatish Balay /* To take care of paging effects */ 257fb9f68cSSatish Balay PetscMemcpy(A,B,sizeof(Scalar)*0); 26d3093643SSatish Balay ierr = PetscGetTime(&x); CHKERRA(ierr); 27b0878937SSatish Balay 28d3093643SSatish Balay ierr = PetscGetTime(&x); CHKERRA(ierr); 29e24b481bSBarry Smith /* 30e24b481bSBarry Smith PetscMemcpy(A,B,sizeof(Scalar)*8000000); 31e24b481bSBarry Smith PetscMemcpy(A,B,sizeof(Scalar)*8000000); 32e24b481bSBarry Smith PetscMemcpy(A,B,sizeof(Scalar)*8000000); 33e24b481bSBarry Smith PetscMemcpy(A,B,sizeof(Scalar)*8000000); 34e24b481bSBarry Smith PetscMemcpy(A,B,sizeof(Scalar)*8000000); 35e24b481bSBarry Smith PetscMemcpy(A,B,sizeof(Scalar)*8000000); 36e24b481bSBarry Smith PetscMemcpy(A,B,sizeof(Scalar)*8000000); 37e24b481bSBarry Smith PetscMemcpy(A,B,sizeof(Scalar)*8000000); 38e24b481bSBarry Smith PetscMemcpy(A,B,sizeof(Scalar)*8000000); 39e24b481bSBarry Smith PetscMemcpy(A,B,sizeof(Scalar)*8000000); */ 40e24b481bSBarry Smith { int j; 41e24b481bSBarry Smith for (j = 0; j<10; j++ ) { 42e24b481bSBarry Smith for ( i=0; i<8000000; i++) { 43e24b481bSBarry Smith B[i] = A[i]; 44e24b481bSBarry Smith } 45e24b481bSBarry Smith }} 46e24b481bSBarry Smith 47d3093643SSatish Balay ierr = PetscGetTime(&y); CHKERRA(ierr); 487fb9f68cSSatish Balay PetscMemcpy(A,B,sizeof(Scalar)*0); 497fb9f68cSSatish Balay PetscMemcpy(A,B,sizeof(Scalar)*0); 507fb9f68cSSatish Balay PetscMemcpy(A,B,sizeof(Scalar)*0); 517fb9f68cSSatish Balay PetscMemcpy(A,B,sizeof(Scalar)*0); 527fb9f68cSSatish Balay PetscMemcpy(A,B,sizeof(Scalar)*0); 537fb9f68cSSatish Balay PetscMemcpy(A,B,sizeof(Scalar)*0); 547fb9f68cSSatish Balay PetscMemcpy(A,B,sizeof(Scalar)*0); 557fb9f68cSSatish Balay PetscMemcpy(A,B,sizeof(Scalar)*0); 567fb9f68cSSatish Balay PetscMemcpy(A,B,sizeof(Scalar)*0); 577fb9f68cSSatish Balay PetscMemcpy(A,B,sizeof(Scalar)*0); 58d3093643SSatish Balay ierr = PetscGetTime(&z); CHKERRA(ierr); 59b0878937SSatish Balay 60b0878937SSatish Balay fprintf(stderr,"%s : \n","PetscMemcpy"); 61e24b481bSBarry Smith fprintf(stderr," %-11s : %e MB/s\n","Bandwidth",10.0*8*8/(y-x)); 62b0878937SSatish Balay fprintf(stderr," %-11s : %e sec\n","Latency",(z-y)/10.0); 63e24b481bSBarry Smith fprintf(stderr," %-11s : %e sec\n","Per Scalar",(2*y-x-z)/8000000.0); 64b0878937SSatish Balay 65b0878937SSatish Balay PetscFinalize(); 663a40ed3dSBarry Smith PetscFunctionReturn(0); 67b0878937SSatish Balay } 68