DAE with continuous callback and DImplicitEuler
results in incorrect integrator.u
in affect!()
#1029
Labels
DImplicitEuler
results in incorrect integrator.u
in affect!()
#1029
When using continuous callbacks with
DAEProblem
andDImplicitEuler
, the values inintegrator.u
don't match the values computed by the root finder in thecondition()
function. The root finder successfully computes the correct root but then this root is not passed on to theaffect!()
function.With
IDA
, however, everything works as expected.Expected behavior
When a
ContinuousCallback
is triggered, the values inintegrator.u
(in the call toaffect!(integrator)
) should be identical to the root found in thecondition()
function.With
IDA
this is the case. This is demonstrated in the example output below.Also note that sometimes
DImplicitEuler
does in fact pass on the correct root.Minimal Reproducible Example 👇
Output (abbreviated)
Output of
using Pkg; Pkg.status()
Output of
using Pkg; Pkg.status(; mode = PKGMODE_MANIFEST)
Output of
versioninfo()
The text was updated successfully, but these errors were encountered: