1 #ifndef lint 2 static char vcid[] = "$Id: PetscMemcpy.c,v 1.5 1996/03/08 23:27:35 balay Exp bsmith $"; 3 #endif 4 5 #include "stdio.h" 6 #include "petsc.h" 7 8 int main( int argc, char **argv) 9 { 10 double x, y,z; 11 int i; 12 Scalar A[10000], B[10000]; 13 14 PetscInitialize(&argc, &argv,0,0); 15 for (i=0; i<10000; i++) { 16 A[i] = i%61897; 17 B[i] = i%61897; 18 } 19 /* To take care of paging effects */ 20 PetscMemcpy(A,B,sizeof(Scalar)*0); 21 x = PetscGetTime(); 22 23 x = PetscGetTime(); 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 PetscMemcpy(A,B,sizeof(Scalar)*10000); 34 y = PetscGetTime(); 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 PetscMemcpy(A,B,sizeof(Scalar)*0); 45 z = PetscGetTime(); 46 47 fprintf(stderr,"%s : \n","PetscMemcpy"); 48 fprintf(stderr," %-11s : %e sec\n","Latency",(z-y)/10.0); 49 fprintf(stderr," %-11s : %e sec\n","Per Scalar",(2*y-x-z)/100000.0); 50 51 PetscFinalize(); 52 return 0; 53 } 54