xref: /petsc/src/benchmarks/PetscMemcpy.c (revision 2b7fea2a9991bf5a64900aaa80f2645ded1d7a90)
1*2b7fea2aSSatish Balay /*$Id: PetscMemcpy.c,v 1.19 2001/01/17 22:28:38 bsmith Exp balay $*/
2d8e9fea7SSatish Balay 
3b0878937SSatish Balay #include "petsc.h"
4b0878937SSatish Balay 
55a655dc6SBarry Smith #undef __FUNC__
65a655dc6SBarry Smith #define __FUNC__ "main"
7b0878937SSatish Balay int main(int argc,char **argv)
8b0878937SSatish Balay {
9b0a32e0cSBarry Smith   PetscLogDouble x,y,z;
10d3093643SSatish Balay   int        i,ierr;
11e24b481bSBarry Smith   Scalar     *A,*B;
12b0878937SSatish Balay 
1377c4ece6SBarry Smith   PetscInitialize(&argc,&argv,0,0);
14e24b481bSBarry Smith 
15ac355199SBarry Smith   ierr = PetscMalloc(8000000*sizeof(Scalar),&A);CHKERRQ(ierr);
16ac355199SBarry Smith   ierr = PetscMalloc(8000000*sizeof(Scalar),&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 */
23ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
24ac355199SBarry Smith   ierr = PetscGetTime(&x);CHKERRQ(ierr);
25b0878937SSatish Balay 
26ac355199SBarry Smith   ierr = PetscGetTime(&x);CHKERRQ(ierr);
27e24b481bSBarry Smith   /*
28e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000);
29e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000);
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   { 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);
46ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
47ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
48ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
49ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
50ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
51ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
52ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
53ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
54ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
55ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*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);
6135d8aa7fSBarry Smith   fprintf(stdout,"    %-11s : %e sec\n","Per Scalar",(2*y-x-z)/8000000.0);
62b0878937SSatish Balay 
63*2b7fea2aSSatish Balay   ierr = PetscFinalize();CHKERRQ(ierr);
643a40ed3dSBarry Smith   PetscFunctionReturn(0);
65b0878937SSatish Balay }
66