Modeling piecewise linear diodes (discontinuous behaviors) with Dyad?

OK, so there is ongoing work on a switched system solver then. Good to know.

Some months ago, a nearby researcher @benoit_caillaud pointed me to the work of one of his colleague: the nonsmooth ODE solver Siconos. Very nice contact mechanics simulation videos. I’ve never tried it though. From what I got, the downside is that it is a first order ODE solver.