Scribd
Upload
19 views

backup.py

Uploaded by

adnansami987606
Copyright
© © All Rights Reserved
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
19 views

backup.py

Uploaded by

adnansami987606
Copyright
© © All Rights Reserved
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
You are on page 1/ 2

import tidy3d as td

import numpy as np
import gdstk
from tidy3d_lambda import entrypoint

@entrypoint
def generate_object(param):
# parameters
x0=0 # x coordinate of directional coupler center (um)
y0=0 # y coordinate of directional coupler center (um)
z0=0 # z coordinate of directional coupler center (um)
wg_width=0.5 # width of the waveguide (um)
wg_thickness=0.22 # thickness of the waveguide (um)
wg_spacing_in=0 # spacing between the input straight waveguides (um)
wg_spacing_coup=1.0 # gap size of the coupling region (um)
wg_length=4 # length of the input straight waveguides (um)
coup_length=4 #length of the coupling region (um)
bend_length=4 # length of the bending (um)
sidewall_angle=0 # side wall angle of the waveguide (rad)
# medium of the waveguide
medium=td.material_library["cSi"]["Li1993_293K"]
cell = gdstk.Cell("dc") # define a gds cell

# define a path
coup = gdstk.RobustPath(
(-wg_length - bend_length - coup_length / 2 + x0, wg_spacing_in / 2 + y0),
wg_width,
layer=1,
datatype=0,
)

# add a curved segment instead of a straight one


coup.arc(radius=.85*bend_length, initial_angle=np.pi, final_angle=0)

A = (wg_spacing_in - wg_spacing_coup - wg_width) / 4

coup.segment(
(-coup_length /2 + x0, wg_spacing_in / 2 + y0), offset=lambda u: A/100 *
2*np.cos(np.pi * u) - A
)

# Modify this part to create a full circular loop replacing the blue-marked
portion.
coup.arc(radius=coup_length/1.5, initial_angle=-np.pi/2, final_angle=
0.05*np.pi/2)

# add path to the cell


cell.add(coup)

# add the mirrored path to the cell


cell.add(coup.copy().mirror((x0 + 1, y0), (x0, y0)))

# define geometry from the gds cell


dc_geo = td.PolySlab.from_gds(
cell,
gds_layer=1,
axis=2,
slab_bounds=(z0 - wg_thickness / 2, z0 + wg_thickness / 2),
sidewall_angle=sidewall_angle,
)

# define tidy3d structure of the bend


dc = td.Structure(geometry=td.GeometryGroup(geometries=dc_geo), medium=medium)

return [dc]

You might also like