1 #ifdef PETSC_RCS_HEADER 2 static char vcid[] = "$Id: PetscVecNorm.c,v 1.6 1999/03/19 21:24:35 bsmith Exp balay $"; 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