def test_spiral_creation(self):
top = gdspy.Cell("t4")
wgt = WaveguideTemplate(bend_radius=50, resist='-')
sp1 = Spiral(wgt, 1000.0, 1000.0, 10000.0)
tk.add(top, sp1)
# print("Spiral area = "+str(top.area()))
self.assertTrue(len(top.elements) == 1)
self.assertTrue(abs(top.area() - 241199.988149) <= 1e-6)
def test_metal_creation(self):
top = gdspy.Cell("t-electrical")
mt = MetalTemplate(bend_radius=0, resist="+", fab="ETCH")
mt1 = MetalRoute([(0, 0), (0, 250), (100, 250), (100, 500),
(400, 500)], mt)
tk.add(top, mt1)
print("MetalRoute area = " + str(top.area()))
self.assertTrue(len(top.references) == 1)
self.assertTrue(abs(top.area() - 88800.0) <= AREA_TOL)
def test_waveguide_creation(self):
top = gdspy.Cell("t1")
wgt = WaveguideTemplate(bend_radius=50, resist='+', fab="LIFTOFF")
wg1 = Waveguide([(50, 0), (250, 0), (250, 500), (500, 500)], wgt)
wg2 = Waveguide([(0, 0), (0, 100), (-250, 100), (-250, -100)], wgt)
tk.add(top, wg1)
tk.add(top, wg2)
# print("Waveguide area = "+str(top.area()))
self.assertTrue(len(top.elements) == 2)
self.assertTrue(abs(top.area() - 33939.818854) <= 1e-6)
def test_waveguide_creation(self):
top = gdspy.Cell("t1")
wgt = WaveguideTemplate(bend_radius=50, resist="+", fab="LIFTOFF")
wg1 = Waveguide([(50, 0), (250, 0), (250, 500), (500, 500)], wgt)
wg2 = Waveguide([(0, 0), (0, 100), (-250, 100), (-250, -100)], wgt)
tk.add(top, wg1)
tk.add(top, wg2)
print("Waveguide area = " + str(top.area()))
self.assertTrue(len(top.references) == 2)
self.assertTrue(abs(top.area() - 33939.719956968205) <= AREA_TOL)
def test_disk_creation(self):
top = gdspy.Cell("t8")
wgt = WaveguideTemplate(bend_radius=50, wg_width=1.0, resist='+')
wg1 = Waveguide([(0, 0), (100, 0)], wgt)
tk.add(top, wg1)
d1 = Disk(wgt, 60.0, 1.0, parity=1, **wg1.portlist["output"])
tk.add(top, d1)
# print("Disk area = "+str(top.area()))
print(len(top.elements))
self.assertTrue(len(top.elements) == 2)
self.assertTrue(abs(top.area() - 31953.5046652) <= 1e-6)
def test_ring_creation(self):
top = gdspy.Cell("t7")
wgt = WaveguideTemplate(bend_radius=50, wg_width=1.0, resist='+')
wg1 = Waveguide([(0, 0), (100, 0)], wgt)
tk.add(top, wg1)
r1 = Ring(wgt, 60.0, 1.0, parity=1, **wg1.portlist["output"])
tk.add(top, r1)
print("Ring area = " + str(top.area()))
print(len(top.elements))
self.assertTrue(len(top.elements) == 2)
self.assertTrue(abs(top.area() - 13133.8016946) <= 1e-6)
def test_disk_creation(self):
top = gdspy.Cell("t8")
wgt = WaveguideTemplate(bend_radius=50, wg_width=1.0, resist="+")
wg1 = Waveguide([(0, 0), (100, 0)], wgt)
tk.add(top, wg1)
d1 = Disk(wgt, 60.0, 1.0, parity=1, **wg1.portlist["output"])
tk.add(top, d1)
print("Disk area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 2)
self.assertTrue(abs(top.area() - 31543.258086158858) <= AREA_TOL)
def test_ring_creation(self):
top = gdspy.Cell("t7")
wgt = WaveguideTemplate(bend_radius=50, wg_width=1.0, resist="+")
wg1 = Waveguide([(0, 0), (100, 0)], wgt)
tk.add(top, wg1)
r1 = Ring(wgt, 60.0, 1.0, parity=1, **wg1.portlist["output"])
tk.add(top, r1)
print("Ring area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 2)
self.assertTrue(abs(top.area() - 13133.710126840513) <= AREA_TOL)
def test_stripslotcoupler_creation(self):
top = gdspy.Cell("t-stripslotycoupler")
wgt_strip = WaveguideTemplate(bend_radius=50,
wg_type='strip',
wg_width=0.7)
wgt_slot = WaveguideTemplate(bend_radius=50,
wg_type='slot',
wg_width=0.7,
slot=0.2)
wg1 = Waveguide([(0, 0), (100, 0)], wgt_strip)
tk.add(top, wg1)
ycoup = StripSlotYCoupler(wgt_strip,
wgt_slot,
10.0,
0.2,
end_slot_width=0,
**wg1.portlist["output"])
tk.add(top, ycoup)
(x1, y1) = ycoup.portlist["output"]["port"]
wg2 = Waveguide([(x1, y1), (x1 + 100, y1)], wgt_slot)
tk.add(top, wg2)
ycoup2 = StripSlotMMICoupler(wgt_strip, wgt_slot, 2.5, 6.0, 20.0,
**wg2.portlist["output"])
tk.add(top, ycoup2)
# print("StripSlotCoupler area = "+str(top.area()))
print(len(top.elements))
self.assertTrue(len(top.elements) == 4)
self.assertTrue(abs(top.area() - 4921.6) <= 1e-6)
def test_mmi1x2_creation(self):
top = gdspy.Cell("t-mmi1x2")
wgt = WaveguideTemplate(bend_radius=50, resist='+')
wg1 = Waveguide([(0, 0), (250, 0)], wgt)
tk.add(top, wg1)
mmi = MMI1x2(wgt,
length=50,
width=10,
taper_width=2.0,
wg_sep=3,
**wg1.portlist["output"])
tk.add(top, mmi)
# print("MMI1x2 area = "+str(top.area()))
print(len(top.elements))
self.assertTrue(len(top.elements) == 2)
self.assertTrue(abs(top.area() - 11623.2729314) <= 1e-6)
def test_spiral_creation(self):
top = gdspy.Cell("t4")
wgt = WaveguideTemplate(bend_radius=50, resist='-')
sp1 = Spiral(wgt,
width=2000.0,
length=20000.0,
spacing=50.0,
parity=1,
port=(0, -10000),
direction='WEST')
tk.add(top, sp1)
# print("Spiral area = "+str(top.area()))
self.assertTrue(len(top.elements) == 1)
self.assertTrue(abs(top.area() - 479999.98235) <= 1e-6)
def test_mmi2x2_creation(self):
top = gdspy.Cell("t6")
wgt = WaveguideTemplate(bend_radius=50, wg_width=1.0, resist='+')
wg1 = Waveguide([(0, 0), (0, -100)], wgt)
tk.add(top, wg1)
mmi = MMI2x2(wgt,
length=50,
width=10,
taper_width=2.0,
wg_sep=3.0,
port=(0, 0),
direction='EAST')
tk.add(top, mmi)
# print("MMI2x2 area = "+str(top.area()))
print(len(top.elements))
self.assertTrue(len(top.elements) == 2)
self.assertTrue(abs(top.area() - 9482.38678993) <= 1e-6)
def test_spiral_creation(self):
top = gdspy.Cell("t4")
wgt = WaveguideTemplate(bend_radius=50, resist="-")
sp1 = Spiral(
wgt,
width=2000.0,
length=20000.0,
spacing=50.0,
parity=1,
port=(0, -10000),
direction="WEST",
)
tk.add(top, sp1)
print("Spiral area = " + str(top.area()))
self.assertTrue(len(top.references) == 1)
self.assertTrue(abs(top.area() - 479999.4878648407) <= AREA_TOL)
def test_broadbanddc_creation(self):
top = gdspy.Cell("t-bdc")
wgt = WaveguideTemplate(wg_width=2.0, bend_radius=100, resist='+')
wg1 = Waveguide([(0, 0), (100, 0)], wgt)
tk.add(top, wg1)
bdc = BroadbandDirectionalCoupler(wgt,
20.0,
0.5,
1.0,
angle=np.pi / 12.0,
parity=1,
**wg1.portlist["output"])
tk.add(top, bdc)
# print("BroadbandDirectoinalCoupler area = "+str(top.area()))
print(len(top.elements))
self.assertTrue(len(top.elements) == 2)
self.assertTrue(abs(top.area() - 8386.54754147) <= 1e-6)
def test_fullcoupler_creation(self):
top = gdspy.Cell("t-fc")
wgt = WaveguideTemplate(wg_width=2.0, bend_radius=100, resist="+")
wg1 = Waveguide([(0, 0), (100, 0)], wgt)
tk.add(top, wg1)
fc = FullCoupler(wgt,
40.0,
0.5,
1.0,
angle=np.pi / 12.0,
parity=1,
**wg1.portlist["output"])
tk.add(top, fc)
print("FullCoupler area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 2)
self.assertTrue(abs(top.area() - 9346.514088150665) <= AREA_TOL)
def test_stripslotconverter_creation(self):
top = gdspy.Cell("t-stripslotyconverter")
wgt_strip = WaveguideTemplate(bend_radius=50,
wg_type="strip",
wg_width=0.7)
wgt_slot = WaveguideTemplate(bend_radius=50,
wg_type="slot",
wg_width=0.7,
slot=0.2)
wg1 = Waveguide([(0, 0), (100, 0)], wgt_strip)
tk.add(top, wg1)
ycoup = StripSlotYConverter(wgt_strip,
wgt_slot,
10.0,
0.2,
end_slot_width=0,
**wg1.portlist["output"])
tk.add(top, ycoup)
(x1, y1) = ycoup.portlist["output"]["port"]
wg2 = Waveguide([(x1, y1), (x1 + 100, y1)], wgt_slot)
tk.add(top, wg2)
coup2 = StripSlotMMIConverter(wgt_strip, wgt_slot, 2.5, 6.0, 20.0,
**wg2.portlist["output"])
tk.add(top, coup2)
coup3 = StripSlotConverter(wgt_strip,
wgt_slot,
length1=15.0,
length2=15.0,
start_rail_width=0.1,
end_strip_width=0.4,
d=1.0,
**coup2.portlist["output"])
tk.add(top, coup3)
print("StripSlotConverter area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 5)
self.assertTrue(abs(top.area() - 5559.85) <= AREA_TOL)
def test_mmi1x2_creation(self):
top = gdspy.Cell("t-mmi1x2")
wgt = WaveguideTemplate(bend_radius=50, resist="+")
wg1 = Waveguide([(0, 0), (250, 0)], wgt)
tk.add(top, wg1)
mmi = MMI1x2(wgt,
length=20,
width=7,
taper_length=10.0,
taper_width=2.0,
wg_sep=7 / 2,
output_wg_sep=10.0,
output_length=20.0,
output_width=2.0,
**wg1.portlist["output"])
tk.add(top, mmi)
print("MMI1x2 area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 2)
self.assertTrue(abs(top.area() - 8406.848925208546) <= AREA_TOL)
def test_mmi2x2_creation(self):
top = gdspy.Cell("t6")
wgt = WaveguideTemplate(bend_radius=50, wg_width=1.0, resist="+")
wg1 = Waveguide([(0, 0), (0, -100)], wgt)
tk.add(top, wg1)
mmi = MMI2x2(
wgt,
length=50,
width=10,
taper_width=2.0,
wg_sep=3.0,
port=(0, 0),
direction="EAST",
)
tk.add(top, mmi)
print("MMI2x2 area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 2)
self.assertTrue(abs(top.area() - 9482.324973105131) <= AREA_TOL)
def test_adiabaticcoupler_creation(self):
top = gdspy.Cell("t-ac")
wgt = WaveguideTemplate(wg_width=2.0, bend_radius=100, resist="+")
wg1 = Waveguide([(0, 0), (100, 0)], wgt)
tk.add(top, wg1)
ac = AdiabaticCoupler(wgt,
length1=30.0,
length2=50.0,
length3=20.0,
wg_sep=1.0,
input_wg_sep=3.0,
output_wg_sep=3.0,
dw=0.1,
**wg1.portlist["output"])
tk.add(top, ac)
print("AdiabaticCoupler area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 2)
self.assertTrue(abs(top.area() - 7202.844941894362) <= AREA_TOL)
def test_contradc_creation(self):
top = gdspy.Cell("t-contradc")
wgt = WaveguideTemplate(wg_width=1.0, bend_radius=50, resist="+")
wg1 = Waveguide([(0, 0), (20, 0)], wgt)
tk.add(top, wg1)
cdc = ContraDirectionalCoupler(wgt,
length=30.0,
gap=0.5,
period=0.220,
dc=0.5,
angle=np.pi / 12.0,
width_top=3.0,
width_bot=0.75,
input_bot=True,
**wg1.portlist["output"])
tk.add(top, cdc)
print("Contra DC area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 2)
self.assertTrue(abs(top.area() - 4222.097590756624) <= AREA_TOL)
def test_taper_creation(self):
top = gdspy.Cell("t2")
wgt = WaveguideTemplate(bend_radius=50, resist="+")
wg1 = Waveguide([(50, 0), (250, 0), (250, 500), (500, 500)], wgt)
tk.add(top, wg1)
tp1 = Taper(wgt, 100.0, 0.3, **wg1.portlist["input"])
tp2 = Taper(wgt, 100.0, 0.0, **wg1.portlist["output"])
tk.add(top, tp1)
tk.add(top, tp2)
print("Taper area = " + str(top.area()))
self.assertTrue(len(top.references) == 3)
self.assertTrue(abs(top.area() - 26199.859978484103) <= AREA_TOL)
def test_taper_creation(self):
top = gdspy.Cell("t2")
wgt = WaveguideTemplate(bend_radius=50, resist='+')
wg1 = Waveguide([(50, 0), (250, 0), (250, 500), (500, 500)], wgt)
tk.add(top, wg1)
tp1 = Taper(wgt, 100.0, 0.3, **wg1.portlist["input"])
tp2 = Taper(wgt, 100.0, 0.0, **wg1.portlist["output"])
tk.add(top, tp1)
tk.add(top, tp2)
# print("Taper area = "+str(top.area()))
self.assertTrue(len(top.elements) == 3)
self.assertTrue(abs(top.area() - 27005.909427) <= 1e-6)
def test_dbr_creation(self):
top = gdspy.Cell("t-dbr")
wgt = WaveguideTemplate(bend_radius=50, resist="+")
wg1 = Waveguide([(0, 0), (100, 0)], wgt)
tk.add(top, wg1)
dbr1 = DBR(wgt, 10.0, 0.85, 0.5, 0.4, **wg1.portlist["output"])
tk.add(top, dbr1)
(x1, y1) = dbr1.portlist["output"]["port"]
wg2 = Waveguide([(x1, y1), (x1 + 100, y1), (x1 + 100, y1 + 100)], wgt)
tk.add(top, wg2)
dbr2 = DBR(wgt, 10.0, 0.85, 0.5, 0.6, **wg2.portlist["output"])
tk.add(top, dbr2)
print("DBR area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 4)
self.assertTrue(abs(top.area() - 9093.453181758501) <= AREA_TOL)
def test_grating_coupler_creation(self):
top = gdspy.Cell("t3")
wgt = WaveguideTemplate(bend_radius=50, resist="+", fab="ETCH")
wg1 = Waveguide([(0, 0), (250, 0), (250, 500), (500, 500)], wgt)
tk.add(top, wg1)
gc1 = GratingCouplerStraight(wgt,
width=20,
length=50,
taper_length=20,
period=1.0,
dutycycle=0.7,
**wg1.portlist["input"])
tk.add(top, gc1)
gc2 = GratingCouplerFocusing(
wgt,
focus_distance=20.0,
width=20,
length=50,
period=1.0,
dutycycle=0.7,
port=(500, 600),
direction="NORTH",
)
tk.add(top, gc2)
gc3 = GratingCoupler(wgt,
theta=np.pi / 5.0,
length=40.0,
taper_length=20.0,
period=1.0,
dutycycle=0.6,
ridge=True,
**wg1.portlist["output"])
tk.add(top, gc3)
print("Grating coupler area = " + str(top.area()))
self.assertTrue(len(top.references) == 4)
self.assertTrue(abs(top.area() - 32705.633256135465) <= AREA_TOL)
def test_bbend_creation(self):
top = gdspy.Cell("t-bbend")
wgt = WaveguideTemplate(bend_radius=25, resist="+")
wg1 = Waveguide([(0, 0), (100, 0)], wgt)
tk.add(top, wg1)
bb1 = BBend(wgt, [(100, 0), (200, 0), (100, 100), (200, 100)])
tk.add(top, bb1)
x, y = bb1.portlist["output"]["port"]
wg2 = Waveguide([(x, y), (x + 100, y)], wgt)
tk.add(top, wg2)
print("BBend area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 3)
self.assertTrue(abs(top.area() - 8823.676248339165) <= AREA_TOL)
def test_eulersbend_creation(self):
top = gdspy.Cell("t-esbend")
wgt = WaveguideTemplate(bend_radius=25, resist="+")
wg1 = Waveguide([(0, 0), (25, 0)], wgt)
tk.add(top, wg1)
esb = EulerSBend(wgt, 200.0, 100.0, **wg1.portlist["output"])
tk.add(top, esb)
x, y = esb.portlist["output"]["port"]
wg2 = Waveguide([(x, y), (x + 25, y)], wgt)
tk.add(top, wg2)
print("EulerSBend area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 3)
self.assertTrue(abs(top.area() - 6888.198875373508) <= AREA_TOL)
def test_sbend_creation(self):
top = gdspy.Cell("t-sbend")
wgt = WaveguideTemplate(wg_width=2.0, bend_radius=100, resist="+")
wg1 = Waveguide([(0, 0), (100, 0)], wgt)
tk.add(top, wg1)
sb1 = SBend(wgt, 200.0, 100.0, **wg1.portlist["output"])
tk.add(top, sb1)
x, y = sb1.portlist["output"]["port"]
wg2 = Waveguide([(x, y), (x + 100, y)], wgt)
tk.add(top, wg2)
print("SBend area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 3)
self.assertTrue(abs(top.area() - 10271.21589480641) <= AREA_TOL)
def test_mzi_creation(self):
top = gdspy.Cell("t-mzi")
wgt = WaveguideTemplate(bend_radius=50, wg_width=1.0, resist='+')
htr_mt = MetalTemplate(width=25,
clad_width=25,
bend_radius=wgt.bend_radius,
resist='+',
fab="ETCH",
metal_layer=13,
metal_datatype=0,
clad_layer=14,
clad_datatype=0)
mt = MetalTemplate(width=25,
clad_width=25,
resist='+',
fab="ETCH",
metal_layer=11,
metal_datatype=0,
clad_layer=12,
clad_datatype=0)
wg_in = Waveguide([(0, 0), (300, 0)], wgt)
tk.add(top, wg_in)
mzi = MachZehnder(wgt,
MMIlength=50,
MMIwidth=10,
MMItaper_width=2.0,
MMIwg_sep=3,
arm1=0,
arm2=100,
heater=True,
heater_length=400,
mt=htr_mt,
**wg_in.portlist["output"])
tk.add(top, mzi)
wg_out = Waveguide([
mzi.portlist["output"]["port"],
(mzi.portlist["output"]["port"][0] + 300,
mzi.portlist["output"]["port"][1])
], wgt)
tk.add(top, wg_out)
# print("MZI area = "+str(top.area()))
print(len(top.elements))
self.assertTrue(len(top.elements) == 3)
self.assertTrue(abs(top.area() - 184111.305754) <= 1e-6)
def test_ebend_creation(self):
top = gdspy.Cell("t-ebend")
wgt = WaveguideTemplate(bend_radius=25, resist="+")
wg1 = Waveguide([(0, 0), (25, 0)], wgt)
tk.add(top, wg1)
eb1 = EBend(wgt, turnby=np.pi / 4.0, **wg1.portlist["output"])
tk.add(top, eb1)
x, y = eb1.portlist["output"]["port"]
wg2 = Waveguide([(x, y), (x + 25 / np.sqrt(2), y + 25 / np.sqrt(2))],
wgt)
tk.add(top, wg2)
print("EBend area = " + str(top.area()))
print(len(top.references))
self.assertTrue(len(top.references) == 3)
self.assertTrue(abs(top.area() - 2142.4695237593724) <= AREA_TOL)
def test_grating_coupler_creation(self):
top = gdspy.Cell("t3")
wgt = WaveguideTemplate(bend_radius=50, resist='+', fab='ETCH')
wg1 = Waveguide([(0, 0), (250, 0), (250, 500), (500, 500)], wgt)
tk.add(top, wg1)
gc1 = GratingCouplerStraight(wgt,
width=20,
length=50,
taper_length=20,
period=1.0,
dutycycle=0.7,
**wg1.portlist["input"])
tk.add(top, gc1)
gc2 = GratingCouplerFocusing(wgt,
focus_distance=20.0,
width=20,
length=50,
period=1.0,
dutycycle=0.7,
**wg1.portlist["output"])
tk.add(top, gc2)
# print("Grating coupler area = "+str(top.area()))
self.assertTrue(len(top.elements) == 3)
self.assertTrue(abs(top.area() - 30061.6918115) <= 1e-6)
