Hello! I am simulating a DAE model as an ODEProblem using a constant mass matrix. A model parameter p changes from 0 to 1 at time t=t1 using DiscreteCallback. The solver returns dt<=dtmin at t=t1 when p changes. However, running the model with p=1 from the start works fine. Am I missing something…

The issue was due to reinitialization which has now been improved for DAEs.

DiffEq: Discrete callback causes dt<=dtmin

aaliqureshi October 13, 2023, 9:39pm 5

Closing this thread.

DAE mass matrix formulation fails but implicit formulation works