def test_Simulation_cascade(self):
simu = sim.Simulation(self.nl)
simu._cascade()
assert simu.combined.nets == [0, 1]
assert len(simu.combined.f) == 2000
assert len(simu.combined.s) == 2000
assert len(simu.combined.nets) == 2
assert simu.combined.s.shape[1] == 2
def test_mzi(self):
y1 = core.ComponentInstance(dev.ebeam_y_1550)
y2 = core.ComponentInstance(dev.ebeam_y_1550)
wg1 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':50e-6})
wg2 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':150e-6})
c1 = [y1, y1, y2, y2]
p1 = [1, 2, 2, 1]
c2 = [wg1, wg2, wg1, wg2]
p2 = [0, 0, 1, 1]
con = zip(c1, p1, c2, p2)
nl = core.Netlist()
nl.load(con, formatter='ll')
simu = sim.Simulation(nl)
freq = simu.freq_array
zero2zero = abs(simu.s_parameters()[:, 0, 0])**2
zero2one = abs(simu.s_parameters()[:, 0, 1])**2
one2zero = abs(simu.s_parameters()[:, 1, 0])**2
one2one = abs(simu.s_parameters()[:, 1, 1])**2
# np.savez('test_simphony_test_mzi', freq=freq, zero2zero=zero2zero, zero2one=zero2one, one2zero=one2zero, one2one=one2one)
expected = np.load(os.path.join(os.path.dirname(os.path.realpath(__file__)), 'benchmarks', 'test_simphony_test_mzi.npz'))
assert np.all(freq == expected['freq'])
assert np.all(zero2zero == expected['zero2zero'])
assert np.all(zero2one == expected['zero2one'])
assert np.all(one2zero == expected['one2zero'])
assert np.all(one2one == expected['one2one'])
# import matplotlib.pyplot as plt
# plt.subplot(221)
# plt.plot(freq, zero2zero)
# plt.subplot(222)
# plt.plot(freq, zero2one)
# plt.subplot(223)
# plt.plot(freq, one2zero)
# plt.subplot(224)
# plt.plot(freq, one2one)
# plt.suptitle("MZI")
# plt.show()
wg2 = core.ComponentInstance(dev.ebeam_wg_integral_1550,
extras={
'length': 150e-6,
'ne': 12.15,
'ng': -3.7
})
c1 = [y1, y1, y2, y2, bg_in, bg_out]
p1 = [1, 2, 2, 1, 0, 0]
c2 = [wg1, wg2, wg1, wg2, y1, y2]
p2 = [0, 0, 1, 1, 0, 0]
con = zip(c1, p1, c2, p2)
nl = core.Netlist()
nl.load(con, formatter='ll')
simu = sim.Simulation(nl)
freq = simu.freq_array
zero2zero = np.log10(abs(simu.s_parameters()[:, 0, 0])**2)
zero2one = np.log10(abs(simu.s_parameters()[:, 0, 1])**2)
one2zero = abs(simu.s_parameters()[:, 1, 0])**2
one2one = abs(simu.s_parameters()[:, 1, 1])**2
# assert np.all(freq == expected['freq'])
# assert np.all(zero2zero == expected['zero2zero'])
# assert np.all(zero2one == expected['zero2one'])
# assert np.all(one2zero == expected['one2zero'])
# assert np.all(one2one == expected['one2one'])
import matplotlib.pyplot as plt
plt.subplot(221)
def test_Simulation_rearrange(self):
simu = sim.Simulation(self.nl)
simu._cascade()
def test_caching(self):
simu = sim.Simulation(self.nl)
assert np.array_equal(simu.freq_array,
np.linspace(1.88e+14, 1.99e+14, 2000))
assert len(simu.cache) == len(
[item for item in self.components if item.model.cachable])
def test_Simulation_initialization(self):
length = 1000
s = sim.Simulation(self.nl, num=length)
assert len(s.freq_array == length)
def test_4Port_Circuit(self):
gc1 = core.ComponentInstance(dev.ebeam_gc_te1550)
gc2 = core.ComponentInstance(dev.ebeam_gc_te1550)
gc3 = core.ComponentInstance(dev.ebeam_gc_te1550)
gc4 = core.ComponentInstance(dev.ebeam_gc_te1550)
y1 = core.ComponentInstance(dev.ebeam_y_1550)
y2 = core.ComponentInstance(dev.ebeam_y_1550)
y3 = core.ComponentInstance(dev.ebeam_y_1550)
bdc1 = core.ComponentInstance(dev.ebeam_bdc_te1550)
bdc2 = core.ComponentInstance(dev.ebeam_bdc_te1550)
term1 = core.ComponentInstance(dev.ebeam_terminator_te1550)
wg1 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':165.51e-6})
wg2 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':247.73e-6})
wg3 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':642.91e-6})
wg4 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':391.06e-6})
wg5 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':10.45e-6})
wg6 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':10.45e-6})
wg7 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':10.45e-6})
wg8 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':10.45e-6})
wg9 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':162.29e-6})
wg10 = core.ComponentInstance(dev.ebeam_wg_integral_1550, extras={'length':205.47e-6})
connections = []
connections.append([gc1, 0, wg1, 1])
connections.append([gc3, 0, wg2, 1])
connections.append([bdc1, 3, wg1, 0])
connections.append([bdc1, 2, wg2, 0])
connections.append([gc2, 0, y1, 0])
connections.append([y1, 1, wg3, 0])
connections.append([y1, 2, wg4, 0])
connections.append([y2, 0, wg4, 1])
connections.append([y3, 0, wg3, 1])
connections.append([y2, 1, wg5, 1])
connections.append([bdc1, 0, wg5, 0])
connections.append([bdc1, 1, wg6, 1])
connections.append([y3, 2, wg6, 0])
connections.append([y2, 2, wg7, 0])
connections.append([y3, 1, wg8, 1])
connections.append([bdc2, 2, wg7, 1])
connections.append([bdc2, 3, wg8, 0])
connections.append([bdc2, 0, wg9, 0])
connections.append([term1, 0, wg9, 1])
connections.append([bdc2, 1, wg10, 0])
connections.append([gc4, 0, wg10, 1])
nl = core.Netlist()
nl.load(connections, formatter='ll')
simu = sim.Simulation(nl)
freq = simu.freq_array
two2zero = abs(simu.s_parameters()[:, 2, 0])**2
two2one = abs(simu.s_parameters()[:, 2, 1])**2
two2two = abs(simu.s_parameters()[:, 2, 2])**2
two2three = abs(simu.s_parameters()[:, 2, 3])**2
# np.savez('test_simphony_test_4Port_Circuit', freq=freq, two2zero=two2zero, two2one=two2one, two2two=two2two, two2three=two2three)
expected = np.load(os.path.join(os.path.dirname(os.path.realpath(__file__)), 'benchmarks', 'test_simphony_test_4Port_Circuit.npz'))
assert np.all(freq == expected['freq'])
assert np.all(two2zero == expected['two2zero'])
assert np.all(two2one == expected['two2one'])
assert np.all(two2two == expected['two2two'])
assert np.all(two2three == expected['two2three'])
for i in range(3):
connections.append([data_line[i], 1, selectors[i][0], 2])
for i in range(1, 3):
connections.append([data_line[i], 0, selectors[i - 1][0], 0])
for i in range(3):
connections.append([selectors[i][0], 3, selectors[i][1], 1])
connections.append([selectors[i][0], 1, selectors[i][1], 3])
for i in range(3):
connections.append([outputs[i], 0, selectors[i][1], 0])
connections.append([terminators[i], 0, selectors[i][1], 2])
connections.append([selectors[2][0], 0, terminators[-1], 0])
nl = core.Netlist()
nl.load(connections, formatter='ll')
simu = sim.Simulation(nl,
start_freq=wl2freq(1551.15e-9),
stop_freq=wl2freq(1524.5e-9))
import matplotlib.pyplot as plt
import matplotlib.gridspec as gridspec
freq, s = simu.freq_array, simu.s_parameters()
fig = plt.figure(tight_layout=True)
gs = gridspec.GridSpec(1, 3)
ax = fig.add_subplot(gs[0, :2])
for inport in range(1):
for outport in range(1, 4):
ax.plot(freq2wl(freq) * 1e9,
np.abs(s[:, outport, inport])**2,
plt.legend()
plt.tight_layout()
plt.show()
# local()
# -----------------------------------------------------------------------------
#
# Run the actual simulation (over some optional frequency range)
#
nl = core.Netlist()
nl.load(connections, formatter='ll')
# simu = sim.Simulation(nl, start_freq=1.925e+14, stop_freq=1.945e+14)
simu = sim.Simulation(nl,
start_freq=wl2freq(1.5501e-6),
stop_freq=wl2freq(1.5499e-6))
# simu = sim.Simulation(nl)
# Get the simulation results
freq, s = simu.freq_array, simu.s_parameters()
# -----------------------------------------------------------------------------
#
# We're interested in investigating behavior at this frequency
#
set_wl = 1550e-9
set_freq = wl2freq(set_wl)
# set_freq = 1.93e+14
# -----------------------------------------------------------------------------