diff --git a/src/ts/examples/tutorials/ex2.c b/src/ts/examples/tutorials/ex2.c index ba3a007..92ff84c 100644 --- a/src/ts/examples/tutorials/ex2.c +++ b/src/ts/examples/tutorials/ex2.c @@ -175,10 +175,16 @@ int main(int argc,char **argv) - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */ ierr = TSSetType(ts,TSBEULER);CHKERRQ(ierr); + PetscInt i; + for (i=0;i<30;i++) { + double tt0 = MPI_Wtime(); + ierr = TSSetTimeStep(ts,dt);CHKERRQ(ierr); + ierr = TSSetTime(ts,0);CHKERRQ(ierr); + ierr = TSSetStepNumber(ts,0);CHKERRQ(ierr); ierr = TSSetMaxSteps(ts,time_steps_max);CHKERRQ(ierr); ierr = TSSetMaxTime(ts,time_total_max);CHKERRQ(ierr); ierr = TSSetExactFinalTime(ts,TS_EXACTFINALTIME_STEPOVER);CHKERRQ(ierr); - ierr = TSSetFromOptions(ts);CHKERRQ(ierr); + //ierr = TSSetFromOptions(ts);CHKERRQ(ierr); /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Solve the problem @@ -193,7 +199,8 @@ int main(int argc,char **argv) Run the timestepping solver */ ierr = TSSolve(ts,u);CHKERRQ(ierr); - + PetscPrintf(PETSC_COMM_WORLD,"%g\n",MPI_Wtime()-tt0); + } /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Free work space. All PETSc objects should be destroyed when they are no longer needed.