Actual source code: PetscMemcpy.c

  1: /*$Id: PetscMemcpy.c,v 1.24 2001/08/29 20:59:33 balay Exp $*/

 3:  #include petsc.h

  7: int main(int argc,char **argv)
  8: {
  9:   PetscLogDouble  x,y,z;
 10:   int             i,ierr;
 11:   PetscScalar     *A,*B;

 13:   PetscInitialize(&argc,&argv,0,0);

 15:   PetscMalloc(8000000*sizeof(PetscScalar),&A);
 16:   PetscMalloc(8000000*sizeof(PetscScalar),&B);

 18:   for (i=0; i<8000000; i++) {
 19:     A[i] = i%61897;
 20:     B[i] = i%61897;
 21:   }
 22:   /* To take care of paging effects */
 23:   PetscMemcpy(A,B,sizeof(PetscScalar)*0);
 24:   PetscGetTime(&x);

 26:   PetscGetTime(&x);
 27:   /*
 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:   PetscMemcpy(A,B,sizeof(PetscScalar)*8000000);
 37:   PetscMemcpy(A,B,sizeof(PetscScalar)*8000000); */
 38:   { int j;
 39:   for (j = 0; j<10; j++) {
 40:     for (i=0; i<8000000; i++) {
 41:       B[i] = A[i];
 42:     }
 43:   }}

 45:   PetscGetTime(&y);
 46:   PetscMemcpy(A,B,sizeof(PetscScalar)*0);
 47:   PetscMemcpy(A,B,sizeof(PetscScalar)*0);
 48:   PetscMemcpy(A,B,sizeof(PetscScalar)*0);
 49:   PetscMemcpy(A,B,sizeof(PetscScalar)*0);
 50:   PetscMemcpy(A,B,sizeof(PetscScalar)*0);
 51:   PetscMemcpy(A,B,sizeof(PetscScalar)*0);
 52:   PetscMemcpy(A,B,sizeof(PetscScalar)*0);
 53:   PetscMemcpy(A,B,sizeof(PetscScalar)*0);
 54:   PetscMemcpy(A,B,sizeof(PetscScalar)*0);
 55:   PetscMemcpy(A,B,sizeof(PetscScalar)*0);
 56:   PetscGetTime(&z);

 58:   fprintf(stdout,"%s : \n","PetscMemcpy");
 59:   fprintf(stdout,"    %-15s : %e MB/s\n","Bandwidth",10.0*8*8/(y-x));
 60:   fprintf(stdout,"    %-15s : %e sec\n","Latency",(z-y)/10.0);
 61:   fprintf(stdout,"    %-15s : %e sec\n","Per PetscScalar",(2*y-x-z)/8000000.0);

 63:   PetscFree(A);
 64:   PetscFree(B);
 65:   PetscFinalize();
 66:   return(0);
 67: }