xref: /petsc/src/benchmarks/PetscMemcpy.c (revision e24b481b4a15e258c09dabebcc19da3246e723e1)
1a5eb4965SSatish Balay #ifdef PETSC_RCS_HEADER
2*e24b481bSBarry Smith static char vcid[] = "$Id: PetscMemcpy.c,v 1.10 1998/03/31 23:34:02 balay Exp $";
3d8e9fea7SSatish Balay #endif
4d8e9fea7SSatish Balay 
5b0878937SSatish Balay #include "petsc.h"
6b0878937SSatish Balay 
7b0878937SSatish Balay int main( int argc, char **argv)
8b0878937SSatish Balay {
947794344SBarry Smith   PLogDouble x, y, z;
10d3093643SSatish Balay   int        i,ierr;
11*e24b481bSBarry Smith   Scalar     *A, *B;
12b0878937SSatish Balay 
1377c4ece6SBarry Smith   PetscInitialize(&argc, &argv,0,0);
14*e24b481bSBarry Smith 
15*e24b481bSBarry Smith   A = (Scalar *) PetscMalloc(8000000*sizeof(Scalar));
16*e24b481bSBarry Smith   B = (Scalar *) PetscMalloc(8000000*sizeof(Scalar));
17*e24b481bSBarry Smith 
18*e24b481bSBarry 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 */
237fb9f68cSSatish Balay   PetscMemcpy(A,B,sizeof(Scalar)*0);
24d3093643SSatish Balay   ierr = PetscGetTime(&x); CHKERRA(ierr);
25b0878937SSatish Balay 
26d3093643SSatish Balay   ierr = PetscGetTime(&x); CHKERRA(ierr);
27*e24b481bSBarry Smith   /*
28*e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000);
29*e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000);
30*e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000);
31*e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000);
32*e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000);
33*e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000);
34*e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000);
35*e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000);
36*e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000);
37*e24b481bSBarry Smith   PetscMemcpy(A,B,sizeof(Scalar)*8000000); */
38*e24b481bSBarry Smith   { int j;
39*e24b481bSBarry Smith   for (j = 0; j<10; j++ ) {
40*e24b481bSBarry Smith     for ( i=0; i<8000000; i++) {
41*e24b481bSBarry Smith       B[i] = A[i];
42*e24b481bSBarry Smith     }
43*e24b481bSBarry Smith   }}
44*e24b481bSBarry Smith 
45d3093643SSatish Balay   ierr = PetscGetTime(&y); CHKERRA(ierr);
467fb9f68cSSatish Balay   PetscMemcpy(A,B,sizeof(Scalar)*0);
477fb9f68cSSatish Balay   PetscMemcpy(A,B,sizeof(Scalar)*0);
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);
56d3093643SSatish Balay   ierr = PetscGetTime(&z); CHKERRA(ierr);
57b0878937SSatish Balay 
58b0878937SSatish Balay   fprintf(stderr,"%s : \n","PetscMemcpy");
59*e24b481bSBarry Smith   fprintf(stderr,"    %-11s : %e MB/s\n","Bandwidth",10.0*8*8/(y-x));
60b0878937SSatish Balay   fprintf(stderr,"    %-11s : %e sec\n","Latency",(z-y)/10.0);
61*e24b481bSBarry Smith   fprintf(stderr,"    %-11s : %e sec\n","Per Scalar",(2*y-x-z)/8000000.0);
62b0878937SSatish Balay 
63b0878937SSatish Balay   PetscFinalize();
643a40ed3dSBarry Smith   PetscFunctionReturn(0);
65b0878937SSatish Balay }
66