xref: /petsc/src/benchmarks/PetscMemcpy.c (revision ac355199758bcb0816d12ec02ec1ad655f49044a)
1*ac355199SBarry Smith /*$Id: PetscMemcpy.c,v 1.18 2001/01/16 18:21:52 balay 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 {
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 
15*ac355199SBarry Smith   ierr = PetscMalloc(8000000*sizeof(Scalar),&A);CHKERRQ(ierr);
16*ac355199SBarry 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 */
23*ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
24*ac355199SBarry Smith   ierr = PetscGetTime(&x);CHKERRQ(ierr);
25b0878937SSatish Balay 
26*ac355199SBarry 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 
45*ac355199SBarry Smith   ierr = PetscGetTime(&y);CHKERRQ(ierr);
46*ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
47*ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
48*ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
49*ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
50*ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
51*ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
52*ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
53*ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
54*ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
55*ac355199SBarry Smith   ierr = PetscMemcpy(A,B,sizeof(Scalar)*0);CHKERRQ(ierr);
56*ac355199SBarry 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 
63b0878937SSatish Balay   PetscFinalize();
643a40ed3dSBarry Smith   PetscFunctionReturn(0);
65b0878937SSatish Balay }
66