Actual source code: PetscMemcpy.c
2: #include <petscsys.h>
3: #include <petsctime.h>
5: int main(int argc,char **argv)
6: {
7: PetscLogDouble x,y;
8: PetscScalar *A,*B;
10: PetscInitialize(&argc,&argv,0,0);
11: PetscCalloc1(8000000,&A);
12: PetscMalloc1(8000000,&B);
14: for (i=0; i<8000000; i++) {
15: A[i] = i%61897;
16: B[i] = i%61897;
17: }
18: /* To take care of paging effects */
19: PetscArraycpy(A,B,8000000);
20: PetscTime(&x);
21: PetscArraycpy(A,B,8000000);
22: PetscTime(&x);
24: fprintf(stdout,"%s : \n","PetscMemcpy");
25: fprintf(stdout," %-15s : %e MB/s\n","Bandwidth",10.0*8*8/(y-x));
27: PetscFree(A);
28: PetscFree(B);
29: PetscFinalize();
30: return 0;
31: }