from __future__ import annotations
import gdsfactory as gf
from gdsfactory.component import Component
from gdsfactory.components.mzit import mzit
from gdsfactory.typings import ComponentSpec
[docs]
@gf.cell
def mzit_lattice(
coupler_lengths: tuple[float, ...] = (10.0, 20.0),
coupler_gaps: tuple[float, ...] = (0.2, 0.3),
delta_lengths: tuple[float, ...] = (10.0,),
mzi: ComponentSpec = mzit,
) -> Component:
r"""Mzi fab tolerant lattice filter.
.. code::
cp1
o4 o2 __ __ o3___w0_t2 _w2___
\ / \
\ length1 / |
============== gap1 |
/ \ |
__/ \_____w0___t1 _w1 |
o3 o1 o4 \ | .
... | | .
o2 o2 o3 | | .
__ _____w0___t1___w1__/ |
\ / |
\ lengthN / |
============== gapN |
/ \ |
__/ \_ |
o1 o1 \___w0___t2___w1_____/
cpN o4
"""
if len(coupler_lengths) != len(coupler_gaps):
raise ValueError(
f"Got {len(coupler_lengths)} coupler_lengths and "
f"{len(coupler_gaps)} coupler_gaps"
)
if len(coupler_lengths) != len(delta_lengths) + 1:
raise ValueError(
f"Got {len(coupler_lengths)} coupler_lengths and "
f"{len(delta_lengths)} delta_lengths. "
"You need one more coupler_length than delta_lengths "
)
assert len(coupler_lengths) >= 2
c = Component()
cp1 = coupler0 = c << gf.get_component(
mzi,
coupler_gap1=coupler_gaps[0],
coupler_gap2=coupler_gaps[1],
coupler_length1=coupler_lengths[0],
coupler_length2=coupler_lengths[1],
delta_length=delta_lengths[0],
)
couplers = [
c
<< gf.get_component(
mzi,
coupler_gap2=coupler_gap,
coupler_length2=coupler_length,
coupler1=None,
delta_length=delta_length,
)
for coupler_length, coupler_gap, delta_length in zip(
coupler_lengths[2:], coupler_gaps[2:], delta_lengths[1:]
)
]
for i, coupler in enumerate(couplers):
if i % 2 == 0:
coupler.mirror()
coupler.connect("o3", coupler0.ports["o1"])
coupler.connect("o4", coupler0.ports["o2"])
coupler0 = coupler
c.add_port("o1", port=coupler0.ports["o1"])
c.add_port("o2", port=coupler0.ports["o2"])
c.add_port("o3", port=cp1.ports["o3"])
c.add_port("o4", port=cp1.ports["o4"])
return c
if __name__ == "__main__":
# cpl = [10, 20, 30]
# cpg = [0.2, 0.3, 0.5]
# dl0 = [10, 20]
cpl = [10, 20, 30, 40]
cpg = [0.2, 0.3, 0.5, 0.5]
dl0 = [10, 20, 30]
c = mzit_lattice(coupler_lengths=cpl, coupler_gaps=cpg, delta_lengths=dl0)
# c = mzit_lattice()
c.show(show_ports=True)
