1 2 #include <petscsys.h> 3 #include <petsctime.h> 4 5 int main(int argc,char **argv) 6 { 7 PetscLogDouble x,y,z; 8 PetscInt i; 9 PetscErrorCode ierr; 10 PetscScalar *A,*B; 11 12 ierr = PetscInitialize(&argc,&argv,0,0);if (ierr) return ierr; 13 ierr = PetscMalloc1(8000000,&A);CHKERRQ(ierr); 14 ierr = PetscMalloc1(8000000,&B);CHKERRQ(ierr); 15 16 for (i=0; i<8000000; i++) { 17 A[i] = i%61897; 18 B[i] = i%61897; 19 } 20 /* To take care of paging effects */ 21 ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 22 ierr = PetscTime(&x);CHKERRQ(ierr); 23 24 ierr = PetscTime(&x);CHKERRQ(ierr); 25 /* 26 PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 27 PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 28 PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 29 PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 30 PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 31 PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 32 PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 33 PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 34 PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); 35 PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); */ 36 { int j; 37 for (j = 0; j<10; j++) 38 for (i=0; i<8000000; i++) B[i] = A[i];} 39 40 ierr = PetscTime(&y);CHKERRQ(ierr); 41 ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 42 ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 43 ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 44 ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 45 ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 46 ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 47 ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 48 ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 49 ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 50 ierr = PetscMemcpy(A,B,sizeof(PetscScalar)*0);CHKERRQ(ierr); 51 ierr = PetscTime(&z);CHKERRQ(ierr); 52 53 fprintf(stdout,"%s : \n","PetscMemcpy"); 54 fprintf(stdout," %-15s : %e MB/s\n","Bandwidth",10.0*8*8/(y-x)); 55 fprintf(stdout," %-15s : %e sec\n","Latency",(z-y)/10.0); 56 fprintf(stdout," %-15s : %e sec\n","Per PetscScalar",(2*y-x-z)/8000000.0); 57 58 ierr = PetscFree(A);CHKERRQ(ierr); 59 ierr = PetscFree(B);CHKERRQ(ierr); 60 ierr = PetscFinalize(); 61 return ierr; 62 } 63