1 #ifdef PETSC_RCS_HEADER 2 static char vcid[] = "$Id: PetscMemcpy.c,v 1.9 1997/10/19 03:30:47 bsmith Exp balay $"; 3 #endif 4 5 #include "petsc.h" 6 7 int main( int argc, char **argv) 8 { 9 PLogDouble x, y, z; 10 int i,ierr; 11 Scalar A[10000], B[10000]; 12 13 PetscInitialize(&argc, &argv,0,0); 14 for (i=0; i<10000; i++) { 15 A[i] = i%61897; 16 B[i] = i%61897; 17 } 18 /* To take care of paging effects */ 19 PetscMemcpy(A,B,sizeof(Scalar)*0); 20 ierr = PetscGetTime(&x); CHKERRA(ierr); 21 22 ierr = PetscGetTime(&x); CHKERRA(ierr); 23 PetscMemcpy(A,B,sizeof(Scalar)*10000); 24 PetscMemcpy(A,B,sizeof(Scalar)*10000); 25 PetscMemcpy(A,B,sizeof(Scalar)*10000); 26 PetscMemcpy(A,B,sizeof(Scalar)*10000); 27 PetscMemcpy(A,B,sizeof(Scalar)*10000); 28 PetscMemcpy(A,B,sizeof(Scalar)*10000); 29 PetscMemcpy(A,B,sizeof(Scalar)*10000); 30 PetscMemcpy(A,B,sizeof(Scalar)*10000); 31 PetscMemcpy(A,B,sizeof(Scalar)*10000); 32 PetscMemcpy(A,B,sizeof(Scalar)*10000); 33 ierr = PetscGetTime(&y); CHKERRA(ierr); 34 PetscMemcpy(A,B,sizeof(Scalar)*0); 35 PetscMemcpy(A,B,sizeof(Scalar)*0); 36 PetscMemcpy(A,B,sizeof(Scalar)*0); 37 PetscMemcpy(A,B,sizeof(Scalar)*0); 38 PetscMemcpy(A,B,sizeof(Scalar)*0); 39 PetscMemcpy(A,B,sizeof(Scalar)*0); 40 PetscMemcpy(A,B,sizeof(Scalar)*0); 41 PetscMemcpy(A,B,sizeof(Scalar)*0); 42 PetscMemcpy(A,B,sizeof(Scalar)*0); 43 PetscMemcpy(A,B,sizeof(Scalar)*0); 44 ierr = PetscGetTime(&z); CHKERRA(ierr); 45 46 fprintf(stderr,"%s : \n","PetscMemcpy"); 47 fprintf(stderr," %-11s : %e sec\n","Latency",(z-y)/10.0); 48 fprintf(stderr," %-11s : %e sec\n","Per Scalar",(2*y-x-z)/100000.0); 49 50 PetscFinalize(); 51 PetscFunctionReturn(0); 52 } 53