1 2 #include <petscsys.h> 3 #include <petsctime.h> 4 5 #undef __FUNCT__ 6 #define __FUNCT__ "main" 7 int main(int argc,char **argv) 8 { 9 PetscLogDouble x,y,z; 10 PetscScalar A[10000]; 11 PetscErrorCode ierr; 12 13 ierr = PetscInitialize(&argc,&argv,0,0);if (ierr) return ierr; 14 /* To take care of paging effects */ 15 ierr = PetscMemzero(A,sizeof(PetscScalar)*0);CHKERRQ(ierr); 16 ierr = PetscTime(&x);CHKERRQ(ierr); 17 18 ierr = PetscTime(&x);CHKERRQ(ierr); 19 ierr = PetscMemzero(A,sizeof(PetscScalar)*10000);CHKERRQ(ierr); 20 ierr = PetscMemzero(A,sizeof(PetscScalar)*10000);CHKERRQ(ierr); 21 ierr = PetscMemzero(A,sizeof(PetscScalar)*10000);CHKERRQ(ierr); 22 ierr = PetscMemzero(A,sizeof(PetscScalar)*10000);CHKERRQ(ierr); 23 ierr = PetscMemzero(A,sizeof(PetscScalar)*10000);CHKERRQ(ierr); 24 ierr = PetscMemzero(A,sizeof(PetscScalar)*10000);CHKERRQ(ierr); 25 ierr = PetscMemzero(A,sizeof(PetscScalar)*10000);CHKERRQ(ierr); 26 ierr = PetscMemzero(A,sizeof(PetscScalar)*10000);CHKERRQ(ierr); 27 ierr = PetscMemzero(A,sizeof(PetscScalar)*10000);CHKERRQ(ierr); 28 ierr = PetscMemzero(A,sizeof(PetscScalar)*10000);CHKERRQ(ierr); 29 ierr = PetscTime(&y);CHKERRQ(ierr); 30 ierr = PetscMemzero(A,sizeof(PetscScalar)*0);CHKERRQ(ierr); 31 ierr = PetscMemzero(A,sizeof(PetscScalar)*0);CHKERRQ(ierr); 32 ierr = PetscMemzero(A,sizeof(PetscScalar)*0);CHKERRQ(ierr); 33 ierr = PetscMemzero(A,sizeof(PetscScalar)*0);CHKERRQ(ierr); 34 ierr = PetscMemzero(A,sizeof(PetscScalar)*0);CHKERRQ(ierr); 35 ierr = PetscMemzero(A,sizeof(PetscScalar)*0);CHKERRQ(ierr); 36 ierr = PetscMemzero(A,sizeof(PetscScalar)*0);CHKERRQ(ierr); 37 ierr = PetscMemzero(A,sizeof(PetscScalar)*0);CHKERRQ(ierr); 38 ierr = PetscMemzero(A,sizeof(PetscScalar)*0);CHKERRQ(ierr); 39 ierr = PetscMemzero(A,sizeof(PetscScalar)*0);CHKERRQ(ierr); 40 ierr = PetscTime(&z);CHKERRQ(ierr); 41 42 fprintf(stdout,"%s : \n","PetscMemzero"); 43 fprintf(stdout," %-15s : %e sec\n","Latency",(z-y)/10.0); 44 fprintf(stdout," %-15s : %e sec\n","Per PetscScalar",(2*y-x-z)/100000.0); 45 46 ierr = PetscFinalize(); 47 return ierr; 48 } 49