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])