DAE integration with events (rootfinding) within Pyomo

[François LESAGE]

Hello Everybody,

I have developed a model including a DAE system (about 20 equations) and relying heavily on events (physical modifications of the system, which requires a new computation of some states). For now, I use Matlab solvers to integrate this system, which works well.

I wish however to migrate the model on Python, and more precisely I wish to use Pyomo (in order not only to simulate but also do some dynamic optimization e.g.).

My question is: can I use the rootfinding features of IDAS (with CasADi) within Pyomo, and if yes, how?

I read Pyomo, CasADi and Sundials docs, but could not find the answer.

An example should be sufficient, but any help is welcome.

Thanks a lot!

François