1 #ifdef PETSC_RCS_HEADER 2 static char vcid[] = "$Id: PetscVecNorm.c,v 1.4 1997/11/28 16:22:21 bsmith Exp balay $"; 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 = PetscGetTime(&t1); CHKERRA(ierr); 21 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 22 23 ierr = PetscGetTime(&t1); 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 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 34 ierr = PetscGetTime(&t2); CHKERRA(ierr); 35 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 36 37 fprintf(stderr,"%s : \n","PetscMemcpy"); 38 fprintf(stderr," Time %g\n",t2-t1); 39 40 PetscFinalize(); 41 PetscFunctionReturn(0); 42 } 43