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_linear():
var('R1 R2 V0')
c = SubCircuit(toolkit=symbolic)
c['V0'] = VS(1, gnd, v=V0, vac=1, toolkit=symbolic)
c['L'] = L(1,2, L=1e-3)
c['R1'] = R(2, 3, r = Symbol('R1'))
c['R2'] = R(3, gnd, r = Symbol('R2'))
dc = SymbolicDC(c)
res = dc.solve()
assert_equal(sympy.simplify(res.v(3, gnd) - V0*R2/(R1+R2)), 0)
def test_nonlinear():
var('k qelectron I0 Isat qelectron T', positive=True, real=True)
c = SubCircuit(toolkit=symbolic)
c['I0'] = IS(gnd, 'net1', i=I0, toolkit=symbolic)
c['D'] = Diode('net1', gnd, IS=Isat, toolkit=symbolic)
dc = SymbolicDC(c)
dc.epar.T = T
res = dc.solve()
assert_equal(sympy.simplify(res.v('net1') - k * T / qelectron * log(I0/Isat+1)), 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])