xref: /petsc/src/benchmarks/PetscMemcpy.c (revision b0a32e0c6855ee6a6cd3495fa7da12ea9885bc5d)
1*b0a32e0cSBarry Smith /*$Id: PetscMemcpy.c,v 1.16 2000/11/28 17:32:38 bsmith Exp bsmith $*/
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 {
9*b0a32e0cSBarry 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 
15*b0a32e0cSBarry Smith ierr = PetscMalloc(8000000*sizeof(Scalar),&(  A ));
16*b0a32e0cSBarry Smith ierr = PetscMalloc(8000000*sizeof(Scalar),&(  B ));
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 */
23549d3d68SSatish Balay   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr);
24d3093643SSatish Balay   ierr = PetscGetTime(&x);CHKERRA(ierr);
25b0878937SSatish Balay 
26d3093643SSatish Balay   ierr = PetscGetTime(&x);CHKERRA(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 
45d3093643SSatish Balay   ierr = PetscGetTime(&y);CHKERRA(ierr);
46549d3d68SSatish Balay   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr);
47549d3d68SSatish Balay   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr);
48549d3d68SSatish Balay   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr);
49549d3d68SSatish Balay   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr);
50549d3d68SSatish Balay   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr);
51549d3d68SSatish Balay   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr);
52549d3d68SSatish Balay   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr);
53549d3d68SSatish Balay   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr);
54549d3d68SSatish Balay   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr);
55549d3d68SSatish Balay   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRA(ierr);
56d3093643SSatish Balay   ierr = PetscGetTime(&z);CHKERRA(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 
63b0878937SSatish Balay   PetscFinalize();
643a40ed3dSBarry Smith   PetscFunctionReturn(0);
65b0878937SSatish Balay }
66