Dean Huang
106 lines
2.9 KiB
Python
106 lines
2.9 KiB
Python
# Copyright 2013 DEVSIM LLC
|
|
#
|
|
# SPDX-License-Identifier: Apache-2.0
|
|
|
|
#### circ4.tcl
|
|
####
|
|
#### copy of circ1.tcl using a voltage source
|
|
|
|
import devsim
|
|
import test_common
|
|
|
|
devices = ("MyDevice1", "MyDevice2")
|
|
region = "MyRegion"
|
|
|
|
####
|
|
#### Meshing
|
|
####
|
|
test_common.CreateSimpleMesh(device=devices[0], region=region)
|
|
devsim.create_device(mesh="dog", device=devices[1])
|
|
|
|
for device in devices:
|
|
test_common.SetupResistorConstants(device=device, region="")
|
|
test_common.SetupInitialResistorSystem(device, region)
|
|
|
|
devsim.add_circuit_node(name="cnode0", variable_update="default")
|
|
devsim.add_circuit_node(name="cnode1", variable_update="default")
|
|
devsim.circuit_element(name="R1", n1="cnode1", n2=0, value=1e15)
|
|
devsim.circuit_element(name="V1", n1="cnode0", n2=0, value=0.0)
|
|
|
|
test_common.SetupInitialResistorContact(
|
|
device="MyDevice1",
|
|
contact="top",
|
|
use_circuit_bias=True,
|
|
circuit_node="MyDevice1_top",
|
|
)
|
|
test_common.SetupInitialResistorContact(
|
|
device="MyDevice1",
|
|
contact="bot",
|
|
use_circuit_bias=True,
|
|
circuit_node="MyDevice1_bot",
|
|
)
|
|
test_common.SetupInitialResistorContact(
|
|
device="MyDevice2",
|
|
contact="top",
|
|
use_circuit_bias=True,
|
|
circuit_node="MyDevice2_top",
|
|
)
|
|
test_common.SetupInitialResistorContact(
|
|
device="MyDevice2",
|
|
contact="bot",
|
|
use_circuit_bias=True,
|
|
circuit_node="MyDevice2_bot",
|
|
)
|
|
|
|
|
|
devsim.circuit_node_alias(node="cnode0", alias="MyDevice1_top")
|
|
devsim.circuit_node_alias(node="cnode1", alias="MyDevice2_top")
|
|
devsim.circuit_node_alias(node="cnode1", alias="MyDevice1_bot")
|
|
devsim.circuit_node_alias(node="GND", alias="MyDevice2_bot")
|
|
|
|
devsim.solve(type="dc", absolute_error=1.0, relative_error=1e-14, maximum_iterations=30)
|
|
|
|
for device in devices:
|
|
print(device)
|
|
for name in ("Potential", "IntrinsicElectrons"):
|
|
devsim.print_node_values(device=device, region=region, name=name)
|
|
|
|
for device in devices:
|
|
test_common.SetupCarrierResistorSystem(device=device, region=region)
|
|
|
|
test_common.SetupCarrierResistorContact(
|
|
device="MyDevice1",
|
|
contact="top",
|
|
use_circuit_bias=True,
|
|
circuit_node="MyDevice1_top",
|
|
)
|
|
test_common.SetupCarrierResistorContact(
|
|
device="MyDevice1",
|
|
contact="bot",
|
|
use_circuit_bias=True,
|
|
circuit_node="MyDevice1_bot",
|
|
)
|
|
test_common.SetupCarrierResistorContact(
|
|
device="MyDevice2",
|
|
contact="top",
|
|
use_circuit_bias=True,
|
|
circuit_node="MyDevice2_top",
|
|
)
|
|
test_common.SetupCarrierResistorContact(
|
|
device="MyDevice2",
|
|
contact="bot",
|
|
use_circuit_bias=True,
|
|
circuit_node="MyDevice2_bot",
|
|
)
|
|
|
|
|
|
for v in (0.0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.10):
|
|
devsim.circuit_alter(name="V1", value=v)
|
|
devsim.solve(
|
|
type="dc", absolute_error=1.0, relative_error=1e-9, maximum_iterations=30
|
|
)
|
|
test_common.print_circuit_solution()
|
|
for device in devices:
|
|
test_common.printResistorCurrent(device=device, contact="top")
|
|
test_common.printResistorCurrent(device=device, contact="bot")
|