[petsc-users] [TimeStepping] Eventhandler

Zhang, Hong hongzhang at anl.gov
Sat Feb 2 18:59:42 CST 2019


Hi Moritz,

Are you using adaptive time stepping or fixed time stepping? When a fixed stepsize is used, the next time step after an event will be set to the remaining part of the time step automatically. In your example, presumably the stepsize is dt = 0.2, the trajectory should look like
t=0.9
t=1.0 (event time)
t=1.1
t=1.3
...

So a very small stepsize after an event is expected if the event occurs near the end of the original time step. 

Hong (Mr.)

> On Feb 2, 2019, at 12:08 PM, Moritz Huck via petsc-users <petsc-users at mcs.anl.gov> wrote:
> 
> Hi,
> I am trying to use SetEventHandler.
> I hav eencountered two "problems":
> 
> 1. after an event is encountered the next time step is needlessly small (h=1e-6), even if a don't change anything in the function (for testing).
> 
> 
> 2. the integration does not continue at the event encounter
> e.g. my eventfunction is g(t,u)=t-1,
> the event is found in the step from t= 0.9 to 1.1 at tevent=1
> but the next step is not 1+dt but it is 1.1 + dt
> I would expect the next step to be 1.+dt
> 
> Best Regards,
> Moritz
> 	



More information about the petsc-users mailing list