30 #include "component.h"
39 using namespace device;
60 void diac::calcTheModel (
bool last) {
74 isOn = Ud > (Ibo / gi);
76 nr_double_t Ut, Ieq, Vd;
81 Ut = Ubo /
log (Ibo / Is);
88 Id *=
exp (80.0) * (1.0 + Ud - 80.0) - 1.0;
94 gd = Is / Ut *
exp (Ud);
143 matrix diac::calcMatrixY (nr_double_t frequency) {
169 #define qState 0 // charge state
170 #define cState 1 // current state
181 if (time_prev < time) {