def loss_deembedding_ch14_23(
pitch: float = 127.0,
grating_coupler: ComponentFactory = grating_coupler_te,
input_port_indexes: Tuple[int, ...] = (0, 1),
**kwargs) -> Component:
"""Grating coupler test structure for fiber array.
Connects channel 1->4, 2->3
Args:
pitch:
grating_coupler:
input_port_indexes:
Keyword Args:
cross_section settings
"""
gc = grating_coupler()
c = gf.Component()
dx = pitch
gcs = [
gc.ref(position=(i * dx, 0), port_id="o1", rotation=-90)
for i in range(4)
]
gc_ports = [g.ports["o1"] for g in gcs]
c.add(gcs)
c.add(
get_route(gc_ports[0],
gc_ports[3],
start_straight_length=40.0,
taper=None,
**kwargs).references)
c.add(
get_route(gc_ports[1],
gc_ports[2],
start_straight_length=30.0,
taper=None,
**kwargs).references)
for i, index in enumerate(input_port_indexes):
label = get_input_label(gc_ports[index],
gc,
i,
component_name=inspect.stack()[0][3])
label.position = gc_ports[index].position
c.add(label)
return c
def loss_deembedding_ch13_24(
pitch: float = 127.0,
R: float = 10.0,
grating_coupler_factory: ComponentFactory = grating_coupler_te,
input_port_indexes: Tuple[int, ...] = (0, 1),
cross_section: CrossSectionFactory = strip,
**kwargs) -> Component:
gc = grating_coupler_factory()
c = gf.Component()
dx = pitch
gcs = [
gc.ref(position=(i * dx, 0), port_id="o1", rotation=-90)
for i in range(4)
]
gc_ports = [g.ports["o1"] for g in gcs]
c.add(gcs)
c.add(
get_route(gc_ports[0],
gc_ports[2],
start_straight=40.0,
taper_factory=None,
cross_section=cross_section,
**kwargs).references)
gsi = gc.size_info
p1 = gc_ports[1]
p3 = gc_ports[3]
a = R + 5.0 # 0.5
b = max(2 * a, pitch / 2)
y_bot_align_route = -gsi.width - 5.0
c.add(
connect_loopback(p1,
p3,
a,
b,
R,
y_bot_align_route,
cross_section=cross_section,
**kwargs))
for i, index in enumerate(input_port_indexes):
label = get_input_label(gc_ports[index],
gc,
i,
component_name=inspect.stack()[0][3])
label.position = gc_ports[index].position
c.add(label)
return c
def loss_deembedding_ch14_23(
pitch: float = 127.0,
R: float = 10.0,
grating_coupler_factory: ComponentFactory = grating_coupler_te,
input_port_indexes: Tuple[int, ...] = (0, 1),
**kwargs) -> Component:
gc = grating_coupler_factory()
c = gf.Component()
dx = pitch
gcs = [
gc.ref(position=(i * dx, 0), port_id="o1", rotation=-90)
for i in range(4)
]
gc_ports = [g.ports["o1"] for g in gcs]
c.add(gcs)
c.add(
get_route(gc_ports[0],
gc_ports[3],
start_straight=40.0,
taper_factory=None,
**kwargs).references)
c.add(
get_route(gc_ports[1],
gc_ports[2],
start_straight=30.0,
taper_factory=None,
**kwargs).references)
for i, index in enumerate(input_port_indexes):
label = get_input_label(gc_ports[index],
gc,
i,
component_name=inspect.stack()[0][3])
label.position = gc_ports[index].position
c.add(label)
return c
def grating_coupler_loss_fiber_array(
pitch: float = 127.0,
grating_coupler: ComponentFactory = grating_coupler_te,
input_port_indexes: Tuple[int, ...] = (0, 1),
**kwargs) -> Component:
"""Returns Grating coupler fiber array loopback.
Args:
pitch:
grating_coupler:
cross_section:
input_port_indexes: adds test labels
kwargs: cross_section settings
"""
gc = grating_coupler()
c = gf.Component()
dx = pitch
gcs = [
gc.ref(position=(i * dx, 0), port_id="o1", rotation=-90)
for i in range(2)
]
gc_ports = [g.ports["o1"] for g in gcs]
c.add(gcs)
c.add(
get_route(gc_ports[0],
gc_ports[1],
start_straight_length=40.0,
taper=None,
**kwargs).references)
for i, index in enumerate(input_port_indexes):
label = get_input_label(gc_ports[index],
gc,
i,
component_name=inspect.stack()[0][3])
label.position = gc_ports[index].position
c.add(label)
return c
def loss_deembedding_ch13_24(
pitch: float = 127.0,
grating_coupler: ComponentFactory = grating_coupler_te,
input_port_indexes: Tuple[int, ...] = (0, 1),
cross_section: CrossSectionFactory = strip,
**kwargs) -> Component:
"""Grating coupler test structure for fiber array.
Connects channel 1->3, 2->4
Args:
pitch:
grating_coupler:
input_port_indexes: adds test labels
cross_section:
kwargs: cross_section settings
"""
gc = grating_coupler()
c = gf.Component()
dx = pitch
gcs = [
gc.ref(position=(i * dx, 0), port_id="o1", rotation=-90)
for i in range(4)
]
gc_ports = [g.ports["o1"] for g in gcs]
c.add(gcs)
c.add(
get_route(gc_ports[0],
gc_ports[2],
start_straight_length=40.0,
taper=None,
cross_section=cross_section,
**kwargs).references)
x = cross_section(**kwargs)
radius = x.info["radius"]
gsi = gc.size_info
p1 = gc_ports[1]
p3 = gc_ports[3]
a = radius + 5.0 # 0.5
b = max(2 * a, pitch / 2)
y_bot_align_route = -gsi.width - 5.0
c.add(
connect_loopback(p1,
p3,
a,
b,
y_bot_align_route,
cross_section=cross_section,
**kwargs))
for i, index in enumerate(input_port_indexes):
label = get_input_label(gc_ports[index],
gc,
i,
component_name=inspect.stack()[0][3])
label.position = gc_ports[index].position
c.add(label)
return c
def route_south(
component: Component,
optical_routing_type: int = 1,
excluded_ports: Optional[Tuple[str, ...]] = None,
straight_separation: Number = 4.0,
io_gratings_lines: Optional[List[List[ComponentReference]]] = None,
gc_port_name: str = 1,
bend: ComponentFactory = bend_euler,
straight: ComponentFactory = straight_function,
taper: Optional[ComponentFactory] = taper_function,
auto_widen: bool = True,
select_ports: Callable = select_ports_optical,
cross_section: CrossSectionFactory = strip,
**kwargs,
) -> Routes:
"""Returns Routes
Args:
component: component to route
optical_routing_type: routing heuristic `1` or `2`
`1` uses the component size info to estimate the box size.
`2` only looks at the optical port positions to estimate the size
excluded_ports=[]: list of port names to NOT route
straight_separation
io_gratings_lines: list of ports to which the ports produced by this
function will be connected. Supplying this information helps
avoiding straight collisions
gc_port_name: grating port name
Returns:
list of references, list of ports
Works well if the component looks rougly like a rectangular box with
north ports on the north of the box
south ports on the south of the box
east ports on the east of the box
west ports on the west of the box
"""
excluded_ports = excluded_ports or []
assert optical_routing_type in [
1,
2,
], f"optical_routing_type = {optical_routing_type}, not supported "
optical_ports = list(select_ports(component.ports).values())
optical_ports = [p for p in optical_ports if p.name not in excluded_ports]
csi = component.size_info
references = []
lengths = []
bend90 = bend(cross_section=cross_section, **kwargs) if callable(bend) else bend
dy = abs(bend90.info.dy)
# Handle empty list gracefully
if not optical_ports:
return [], []
conn_params = dict(
bend=bend,
straight=straight,
taper=taper,
auto_widen=auto_widen,
cross_section=cross_section,
**kwargs,
)
# Used to avoid crossing between straights in special cases
# This could happen when abs(x_port - x_grating) <= 2 * dy
delta_gr_min = 2 * dy + 1
sep = straight_separation
# Get lists of optical ports by orientation
direction_ports = direction_ports_from_list_ports(optical_ports)
north_ports = direction_ports["N"]
north_start = north_ports[0 : len(north_ports) // 2]
north_finish = north_ports[len(north_ports) // 2 :]
west_ports = direction_ports["W"]
west_ports.reverse()
east_ports = direction_ports["E"]
south_ports = direction_ports["S"]
north_finish.reverse() # Sort right to left
north_start.reverse() # Sort right to left
ordered_ports = north_start + west_ports + south_ports + east_ports + north_finish
def get_index_port_closest_to_x(x, list_ports):
return np.array([abs(x - p.ports[gc_port_name].x) for p in list_ports]).argmin()
def gen_port_from_port(x, y, p):
new_p = pd.Port(name=p.name, midpoint=(x, y), orientation=90.0, width=p.width)
return new_p
west_ports.reverse()
y0 = min([p.y for p in ordered_ports]) - dy - 0.5
ports_to_route = []
i = 0
optical_xs_tmp = [p.x for p in ordered_ports]
x_optical_min = min(optical_xs_tmp)
x_optical_max = max(optical_xs_tmp)
# Set starting ``x`` on the west side
# ``x`` is the x-coord of the waypoint where the current component port is connected.
# x starts as close as possible to the component.
# For each new port, the distance is increased by the separation.
# The starting x depends on the heuristic chosen : ``1`` or ``2``
if optical_routing_type == 1:
# use component size to know how far to route
x = csi.west - dy - 1
elif optical_routing_type == 2:
# use optical port to know how far to route
x = x_optical_min - dy - 1
else:
raise ValueError("Invalid optical routing type")
# First route the ports facing west
# In case we have to connect these ports to a line of gratings,
# Ensure that the port is aligned with the grating port or
# has enough space for manhattan routing (at least two bend radius)
for p in west_ports:
if io_gratings_lines:
i_grating = get_index_port_closest_to_x(x, io_gratings_lines[-1])
x_gr = io_gratings_lines[-1][i_grating].ports[gc_port_name].x
if abs(x - x_gr) < delta_gr_min:
if x > x_gr:
x = x_gr
elif x < x_gr:
x = x_gr - delta_gr_min
tmp_port = gen_port_from_port(x, y0, p)
ports_to_route.append(tmp_port)
route = get_route(input_port=p, output_port=tmp_port, **conn_params)
references.extend(route.references)
lengths.append(route.length)
x -= sep
i += 1
start_straight_length = 0.5
# First-half of north ports
# This ensures that north ports are routed above the top west one
north_start.reverse() # We need them from left to right
if len(north_start) > 0:
y_max = max([p.y for p in west_ports + north_start])
for p in north_start:
tmp_port = gen_port_from_port(x, y0, p)
route = get_route(
input_port=p,
output_port=tmp_port,
start_straight_length=start_straight_length + y_max - p.y,
**conn_params,
)
references.extend(route.references)
lengths.append(route.length)
ports_to_route.append(tmp_port)
x -= sep
start_straight_length += sep
# Set starting ``x`` on the east side
if optical_routing_type == 1:
# use component size to know how far to route
x = csi.east + dy + 1
elif optical_routing_type == 2:
# use optical port to know how far to route
x = x_optical_max + dy + 1
else:
raise ValueError(
f"Invalid optical routing type. Got {optical_routing_type}, only (1, 2 supported) "
)
i = 0
# Route the east ports
# In case we have to connect these ports to a line of gratings,
# Ensure that the port is aligned with the grating port or
# has enough space for manhattan routing (at least two bend radius)
start_straight_length = 0.5
for p in east_ports:
if io_gratings_lines:
i_grating = get_index_port_closest_to_x(x, io_gratings_lines[-1])
x_gr = io_gratings_lines[-1][i_grating].ports[gc_port_name].x
if abs(x - x_gr) < delta_gr_min:
if x < x_gr:
x = x_gr
elif x > x_gr:
x = x_gr + delta_gr_min
tmp_port = gen_port_from_port(x, y0, p)
route = get_route(
p, tmp_port, start_straight_length=start_straight_length, **conn_params
)
references.extend(route.references)
lengths.append(route.length)
ports_to_route.append(tmp_port)
x += sep
i += 1
# Route the remaining north ports
start_straight_length = 0.5
if len(north_finish) > 0:
y_max = max([p.y for p in east_ports + north_finish])
for p in north_finish:
tmp_port = gen_port_from_port(x, y0, p)
ports_to_route.append(tmp_port)
route = get_route(
input_port=p,
output_port=tmp_port,
start_straight_length=start_straight_length + y_max - p.y,
**conn_params,
)
references.extend(route.references)
lengths.append(route.length)
x += sep
start_straight_length += sep
# Add south ports
ports = [flip(p) for p in ports_to_route] + south_ports
return Routes(references=references, ports=ports, lengths=lengths)