1a5eb4965SSatish Balay #ifdef PETSC_RCS_HEADER 2*549d3d68SSatish Balay static char vcid[] = "$Id: PetscMemcpy.c,v 1.12 1999/03/19 21:24:35 bsmith Exp balay $"; 3d8e9fea7SSatish Balay #endif 4d8e9fea7SSatish Balay 5b0878937SSatish Balay #include "petsc.h" 6b0878937SSatish Balay 75a655dc6SBarry Smith #undef __FUNC__ 85a655dc6SBarry 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 */ 25*549d3d68SSatish Balay ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr); 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); 48*549d3d68SSatish Balay ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr); 49*549d3d68SSatish Balay ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr); 50*549d3d68SSatish Balay ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr); 51*549d3d68SSatish Balay ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr); 52*549d3d68SSatish Balay ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr); 53*549d3d68SSatish Balay ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr); 54*549d3d68SSatish Balay ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr); 55*549d3d68SSatish Balay ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr); 56*549d3d68SSatish Balay ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr); 57*549d3d68SSatish Balay ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr); 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