1 #ifdef PETSC_RCS_HEADER 2 static char vcid[] = "$Id: PetscVecNorm.c,v 1.5 1998/04/21 18:34:13 balay Exp bsmith $"; 3 #endif 4 5 #include "vec.h" 6 7 #undef __FUNC__ 8 #define __FUNC__ "main" 9 int main( int argc, char **argv) 10 { 11 Vec x; 12 double norm; 13 PLogDouble t1,t2; 14 int ierr,n = 10000,flg; 15 16 PetscInitialize(&argc, &argv,0,0); 17 ierr = OptionsGetInt(PETSC_NULL,"-n",&n,&flg); CHKERRA(ierr); 18 19 ierr = VecCreateSeq(PETSC_COMM_SELF,n,&x); CHKERRA(ierr); 20 21 /* To take care of paging effects */ 22 ierr = PetscGetTime(&t1); CHKERRA(ierr); 23 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 24 25 ierr = PetscGetTime(&t1); 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 = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 35 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 36 ierr = PetscGetTime(&t2); CHKERRA(ierr); 37 ierr = VecNorm(x,NORM_2,&norm); CHKERRA(ierr); 38 39 fprintf(stderr,"%s : \n","PetscMemcpy"); 40 fprintf(stderr," Time %g\n",t2-t1); 41 42 PetscFinalize(); 43 PetscFunctionReturn(0); 44 } 45