def test_geteqsys():
var('R1 V0 Isat T')
k = symbolic.kboltzmann
qelectron = symbolic.qelectron
c = SubCircuit(toolkit=symbolic)
c['V0'] = VS('net1', gnd, v=V0, toolkit=symbolic)
c['R1'] = R('net1', 'net2', r=R1)
c['D1'] = Diode('net2', gnd, IS=Isat, toolkit=symbolic)
dc = SymbolicDC(c)
dc.epar.T = T
eqsys, x = dc.get_eqsys()
x0, x2, x3 = x
eqsys_ref = np.array([x3 + x0/R1 - x2/R1,
-Isat*(1 - sympy.exp(qelectron*x2/(T*k))) + x2/R1 - x0/R1,
x0 - V0])
assert sympy.simplify(eqsys_ref[0] - eqsys_ref[0]) == 0
assert sympy.simplify(eqsys_ref[0] - eqsys_ref[0]) == 0
assert sympy.simplify(eqsys_ref[0] - eqsys_ref[0]) == 0
def test_geteqsys():
var('R1 V0 Isat T')
k = symbolic.kboltzmann
qelectron = symbolic.qelectron
c = SubCircuit(toolkit=symbolic)
c['V0'] = VS('net1', gnd, v=V0, toolkit=symbolic)
c['R1'] = R('net1', 'net2', r=R1)
c['D1'] = Diode('net2', gnd, IS=Isat, toolkit=symbolic)
dc = SymbolicDC(c)
dc.epar.T = T
eqsys, x = dc.get_eqsys()
x0, x2, x3 = x
eqsys_ref = np.array([x3 + x0/R1 - x2/R1,
-Isat*(1 - sympy.exp(qelectron*x2/(T*k))) + x2/R1 - x0/R1,
x0 - V0])
assert sympy.simplify(eqsys_ref[0] - eqsys_ref[0]) == 0
assert sympy.simplify(eqsys_ref[0] - eqsys_ref[0]) == 0
assert sympy.simplify(eqsys_ref[0] - eqsys_ref[0]) == 0
def test_geteqsys():
var('k qelectron R1 V0 Isat q T qelectron')
c = SubCircuit(toolkit=symbolic)
c['V0'] = VS('net1', gnd, v=V0, toolkit=symbolic)
c['R1'] = R('net1', 'net2', r=R1)
c['D1'] = Diode('net2', gnd, IS=Isat, toolkit=symbolic)
dc = SymbolicDC(c)
dc.epar.T = Symbol('T')
eqsys, x = dc.get_eqsys()
x0, x2, x3 = x
assert_array_equal(eqsys, [x3 + x0/R1 - x2/R1,
-Isat*(1 - sympy.exp(qelectron*x2/(T*k))) + x2/R1 - x0/R1,
x0 - V0])
