Martin Otter

Zero crossing termination not correct with zero crossing functions crossing closely together

6

I have a simulation where zero crossing gives a wrong result. Its hard to construct a simple example with this effect, but the reason for this issue can be clearly...

FunctionCallingCallback not correctly synchronized with events

3

The FunctionCallingCallback seems to be not correctly synchronized with events as shown with the following test model: ``` module Test_bouncingBall_withFunctionCallBack using DifferentialEquations function f(du,u,p,t) #u[1]: height, u[2]:velocity du[1] = u[2]...

State events: Too many condition calls

1

DifferentialEquations uses **too many "condition" calls** to compute zero crossing functions, when compared to a Modelica tool. For a multi-domain model this is critical, because potentially every variable in a...

CVODE_BDF: Statistics is wrong if zero crossings occur

2

If CVODE_BDF is used with zero crossing functions and state events occur, then the reported statistics is wrong ( sol.destats.naccept, sol.destats.nreject, sol.destats.nf, sol.destats.njacs, sol.destats.nreject ), for example njacs=1, although 20...

saveat and DiscreteCallback may save the some time point twice

5

Assume saveat =[0.0, 0.1, 0.2, 0.3, 0.4, 0.5] and DiscreteCallback triggers a time event at t = 0.1, 0.3, 0.5, then ```julia sol.t = [0.0, 0.1, 0.1, 0.1, 0.2, 0.3,...

SavedValues not stored at events

1

DifferentialEquations supports storing additional values (besides u, du) at a pre-defined time-grid. However, **it seems not possible to store additional values also at event points**. Note, this is important for...

SavingCallback for DAEs: Not clear how to get du

7

The callback `save_func(u, t, integrator)` can be provided to compute additional outputs `y` based on the current time t and the current state u (and the model parameters integrator.p): `y=f(u,t,p)`....

SavedValues seems to be not in solution sol (as needed for ensemble simulations)

1

Assume additional variables are computed with the SavedValues callback: ``` prob = ODEProblem((du,u,p,t) -> du .= u, rand(4,4), (0.0,1.0)) saved_values = SavedValues(Float64, Tuple{Float64,Float64}) cb = SavingCallback((u,t,integrator)->(tr(u),norm(u)), saved_values, saveat=0.0:0.1:1.0) sol =...

CVODE_BDF: Event points are not stored

3

Take the following bouncing ball model ``` using DifferentialEquations, Sundials function f(du,u,p,t) du[1] = u[2] du[2] = -9.81 end function condition(u,t,integrator) z = -u[1] return z end function affect(integrator) integrator.u[2]...

## Description The following asciidoc text ``` [source,C] ---- /* This is a long comment */ ---- Text after the source is in bold ``` Changes the font to "bold":...