By the way, what is the status of derivative discontinuity detection these days in MTK/Dyad. Would it be now possible to implement a piecewise linear diode model like I posted back in late 2023?
If not, perhaps it’s possible to blend the smooth-but-ill-conditioned Shockley with a small-but-not-too-small resistance Ron in series (to limit the di/dv derivative of the I-V curve when conducting) and a large-but-not-too-large resistance in parallel (to limit reverse voltage when blocked). Or as put in Modelica Standard Library, a small conductance Goff.
I feel that using resistive “numerical snubber” is safer than capacitors (unless of course the goal is to model a genuine physical snubber in the circuit), because I find it easier to reason about (how much voltage or current error does it bring for a given forward current or reverse voltage).