[math] ODE step size

View: New views
2 Messages — Rating Filter:   Alert me  

[math] ODE step size

by Angelo Port :: Rate this Message:

Reply to Author | View Threaded | Show Only this Message

Hello,

Is it possible to change the step size parameters on an ODE integrator prior to calling/running the integrator? For example, I would like to control the step size according to the size of the integration time. I am using the DormandPrince54Integrator.

Angelo

Re: [math] ODE step size

by Luc Maisonobe :: Rate this Message:

Reply to Author | View Threaded | Show Only this Message

Angelo Port a écrit :
> Hello,
>
> Is it possible to change the step size parameters on an ODE integrator prior
> to calling/running the integrator? For example, I would like to control the
> step size according to the size of the integration time. I am using the
> DormandPrince54Integrator.

The Dormand-Prince 5(4) integrator is an adaptive stepsize integrator.
This means it handles the step size by itself automatically, changing it
as needed during integration according to error control parameters. The
min and max value of the step size as well as the error control tuning
parameters are set at construction time. They cannot be changed
afterwards (and the step size will change within these prescribed limits).

What you can change however is an output step size (for example to print
the current state at regular intervals). This has nothing to do with the
internal varying integration step size. Changing the output step size is
done by registering a StepNormalizer which itself would wrap your own
fixed step handler. The normalizer will take care to convert the varying
internal step size to your desired fixed step size. Is it this kind of
feature you need ? We can explain in more details how to use it if you want.

There are also fixed step size integrators but their step size is also
fixed at build time and cannot be changed afterwards.

Luc

>
> Angelo


---------------------------------------------------------------------
To unsubscribe, e-mail: user-unsubscribe@...
For additional commands, e-mail: user-help@...