1 #ifdef PETSC_RCS_HEADER 2 static char vcid[] = "$Id: PetscVecNorm.c,v 1.3 1997/10/19 03:30:47 bsmith Exp bsmith $"; 3 #endif 4 5 #include "vec.h" 6 7 int main( int argc, char **argv) 8 { 9 Vec x; 10 double norm; 11 PLogDouble t1,t2; 12 int ierr,n = 10000,flg; 13 14 PetscInitialize(&argc, &argv,0,0); 15 ierr = OptionsGetInt(PETSC_NULL,"-n",&n,&flg); CHKERRA(ierr); 16 17 ierr = VecCreateSeq(PETSC_COMM_SELF,n,&x); CHKERRA(ierr); 18 19 /* To take care of paging effects */ 20 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 21 22 t1 = PetscGetTime(); 23 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 24 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 25 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 26 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 27 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 28 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 29 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 30 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 31 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 32 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 33 t2 = PetscGetTime(); 34 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 35 36 fprintf(stderr,"%s : \n","PetscMemcpy"); 37 fprintf(stderr," Time %g\n",t2-t1); 38 39 PetscFinalize(); 40 PetscFunctionReturn(0); 41 } 42