gplugins 1.3.1__tar.gz → 1.3.3__tar.gz

{gplugins-1.3.1 → gplugins-1.3.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: gplugins
-Version: 1.3.1
+Version: 1.3.3
 Summary: gdsfactory plugins
 Keywords: python
 Author-email: gdsfactory <contact@gdsfactory.com>
@@ -11,24 +11,25 @@ Classifier: Programming Language :: Python :: 3.11
 Classifier: Programming Language :: Python :: 3.12
 Classifier: Operating System :: OS Independent
 License-File: LICENSE
-Requires-Dist: gdsfactory>=9.0.0
+Requires-Dist: gdsfactory>=9.3.5
 Requires-Dist: pint
 Requires-Dist: gdstk
 Requires-Dist: tqdm
 Requires-Dist: numpy
+Requires-Dist: xarray==2025.1.2
 Requires-Dist: pre-commit ; extra == "dev"
 Requires-Dist: pytest ; extra == "dev"
 Requires-Dist: pytest-cov ; extra == "dev"
 Requires-Dist: pytest_regressions ; extra == "dev"
 Requires-Dist: pyswarms ; extra == "dev"
-Requires-Dist: autograd ; extra == "dev"
 Requires-Dist: devsim ; extra == "devsim"
 Requires-Dist: pyvista<=0.43.8 ; extra == "devsim"
-Requires-Dist: tidy3d>=2.7.8 ; extra == "devsim"
+Requires-Dist: tidy3d>=2.8.2,<2.9 ; extra == "devsim"
 Requires-Dist: jupytext ; extra == "docs"
 Requires-Dist: matplotlib ; extra == "docs"
 Requires-Dist: jupyter-book~=1.0 ; extra == "docs"
 Requires-Dist: pyvista[all]<=0.43.8 ; extra == "docs"
+Requires-Dist: bokeh ; extra == "docs"
 Requires-Dist: femwell~=0.1.11 ; extra == "femwell"
 Requires-Dist: meshwell~=1.0.7 ; extra == "femwell"
 Requires-Dist: jinja2 ; extra == "gfviz"
@@ -56,9 +57,8 @@ Requires-Dist: towncrier ; extra == "maintainer"
 Requires-Dist: jax>=0.4.26 ; extra == "meow"
 Requires-Dist: jaxlib>=0.4.26 ; extra == "meow"
 Requires-Dist: flax>=0.8.2 ; extra == "meow"
-Requires-Dist: meow-sim>=0.11,<0.13 ; extra == "meow"
-Requires-Dist: tidy3d>=2.7.8 ; extra == "meow"
-Requires-Dist: scipy<=1.14.1 ; extra == "meow"
+Requires-Dist: meow-sim>=0.13,<0.14 ; extra == "meow"
+Requires-Dist: tidy3d>=2.8.2,<2.9 ; extra == "meow"
 Requires-Dist: jax>=0.4.26 ; extra == "sax"
 Requires-Dist: jaxlib>=0.4.26 ; extra == "sax"
 Requires-Dist: flax>=0.8.2 ; extra == "sax"
@@ -69,7 +69,7 @@ Requires-Dist: ray ; extra == "sax"
 Requires-Dist: bokeh ; extra == "schematic"
 Requires-Dist: ipywidgets ; extra == "schematic"
 Requires-Dist: natsort ; extra == "schematic"
-Requires-Dist: tidy3d>=2.7.8 ; extra == "tidy3d"
+Requires-Dist: tidy3d>=2.8.2,<2.9 ; extra == "tidy3d"
 Requires-Dist: meshio ; extra == "tidy3d"
 Requires-Dist: numpy ; extra == "tidy3d"
 Requires-Dist: meshwell~=1.0.7 ; extra == "tidy3d"
@@ -90,7 +90,7 @@ Provides-Extra: schematic
 Provides-Extra: tidy3d
 Provides-Extra: vlsir
 
-# gplugins 1.3.1
+# gplugins 1.3.3
 
 [![docs](https://github.com/gdsfactory/gplugins/actions/workflows/pages.yml/badge.svg)](https://gdsfactory.github.io/gplugins/)
 [![PyPI](https://img.shields.io/pypi/v/gplugins)](https://pypi.org/project/gplugins/)

{gplugins-1.3.1 → gplugins-1.3.3}/README.md

@@ -1,4 +1,4 @@
-# gplugins 1.3.1
+# gplugins 1.3.3
 
 [![docs](https://github.com/gdsfactory/gplugins/actions/workflows/pages.yml/badge.svg)](https://gdsfactory.github.io/gplugins/)
 [![PyPI](https://img.shields.io/pypi/v/gplugins)](https://pypi.org/project/gplugins/)

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/__init__.py

@@ -1,6 +1,6 @@
 """gplugins - gdsfactory plugins."""
 
-__version__ = "1.3.1"
+__version__ = "1.3.3"
 
 import pathlib
 

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/common/base_models/component.py

@@ -248,14 +248,3 @@ class LayeredComponentBase(BaseModel):
     def get_layer_center(self, layername: str) -> tuple[float, float, float]:
         bbox = self.get_layer_bbox(layername)
         return tuple(np.mean(bbox, axis=0))
-
-    def get_vertices(self, layer_name: str, buffer: float = 0.0):
-        poly = self.polygons[layer_name].buffer(buffer, join_style="mitre")
-        match poly:
-            case MultiPolygon():
-                verts = tuple(tuple(p.exterior.coords) for p in poly.geoms)
-            case Polygon():
-                verts = (tuple(poly.exterior.coords),)
-            case _:
-                raise TypeError(f"Invalid polygon type: {type(poly)}")
-        return verts

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/common/utils/get_component_with_net_layers.py

@@ -1,6 +1,7 @@
 import copy
 
 import gdsfactory as gf
+import klayout.db as kdb
 from gdsfactory import Component
 from gdsfactory.technology import LayerStack, LogicalLayer
 
@@ -35,8 +36,8 @@ def get_component_layer_stack(
 
 
 def get_component_with_net_layers(
-    component,
-    layer_stack,
+    component: Component,
+    layer_stack: LayerStack,
     port_names: list[str],
     delimiter: str = "#",
     new_layers_init: tuple[int, int] = (10010, 0),
@@ -74,7 +75,9 @@ def get_component_with_net_layers(
         net_component = net_component.remove_layers(layers=(port.layer,))
         for polygon in polygons:
             # If polygon belongs to port, create a unique new layer, and add the polygon to it
-            if polygon.sized(3 * gf.kcl.dbu).inside(port.center):
+            if polygon.sized(int(3 * gf.kcl.dbu)).inside(
+                kdb.Point(*port.to_itype().center)
+            ):
                 # if gdstk.inside(
                 #     [port.center],
                 #     gdstk.offset(gdstk.Polygon(polygon), gf.get_active_pdk().grid_size),
@@ -111,8 +114,3 @@ def get_component_with_net_layers(
 
     net_component.name = f"{component.name}_net_layers"
     return net_component
-
-
-if __name__ == "__main__":
-    c = get_component_with_net_layers()
-    c.show()

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/common/utils/get_effective_indices.py

@@ -2,9 +2,10 @@
 
 from __future__ import annotations
 
-from typing import Literal
+from typing import Any, Literal, cast
 
 import numpy as np
+import numpy.typing as npt
 from scipy.optimize import fsolve
 
 
@@ -19,9 +20,9 @@ def get_effective_indices(
     """Returns the effective refractive indices for a 1D mode.
 
     Args:
-        epsilon_core: Relative permittivity of the film.
-        epsilon_substrate: Relative permittivity of the substrate.
-        epsilon_cladding: Relative permittivity of the cladding.
+        core_material: Refractive index of the core material.
+        nsubstrate: Refractive index of the substrate.
+        clad_materialding: Refractive index of the cladding.
         thickness: Thickness of the film in um.
         wavelength: Wavelength in um.
         polarization: Either "te" or "tm".
@@ -66,18 +67,20 @@ def get_effective_indices(
 
     k_0 = 2 * np.pi / wavelength
 
-    def k_f(e_eff):
+    def k_f(e_eff: npt.NDArray[np.floating[Any]]) -> npt.NDArray[np.floating[Any]]:
         return k_0 * np.sqrt(epsilon_core - e_eff) / (epsilon_core if tm else 1)
 
-    def k_s(e_eff):
+    def k_s(e_eff: npt.NDArray[np.floating[Any]]) -> npt.NDArray[np.floating[Any]]:
         return (
             k_0 * np.sqrt(e_eff - epsilon_substrate) / (epsilon_substrate if tm else 1)
         )
 
-    def k_c(e_eff):
+    def k_c(e_eff: npt.NDArray[np.floating[Any]]) -> npt.NDArray[np.floating[Any]]:
         return k_0 * np.sqrt(e_eff - epsilon_cladding) / (epsilon_cladding if tm else 1)
 
-    def objective(e_eff):
+    def objective(
+        e_eff: npt.NDArray[np.floating[Any]],
+    ) -> npt.NDArray[np.floating[Any]]:
         return 1 / np.tan(k_f(e_eff) * thickness) - (
             k_f(e_eff) ** 2 - k_s(e_eff) * k_c(e_eff)
         ) / (k_f(e_eff) * (k_s(e_eff) + k_c(e_eff)))
@@ -92,14 +95,14 @@ def get_effective_indices(
         return []
 
     # and then use fsolve to get exact indices
-    indices_temp = fsolve(objective, indices_temp)
+    indices_temp = cast(npt.NDArray[np.floating[Any]], fsolve(objective, indices_temp))
 
-    indices = []
+    indices: list[float] = []
     for index in indices_temp:
         if not any(np.isclose(index, i, atol=1e-5) for i in indices):
             indices.append(index)
 
-    return np.sqrt(indices).tolist()
+    return cast(list[float], np.sqrt(indices).tolist())
 
 
 def test_effective_index() -> None:

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/common/utils/parse_layer_stack.py

@@ -26,7 +26,7 @@ def list_unique_layer_stack_z(
 def map_unique_layer_stack_z(
     layer_stack: LayerStack,
     include_zmax: bool = True,
-):
+) -> dict[str, set[float]]:
     """Map unique LayerStack z coordinates to various layers.
 
     Args:
@@ -57,7 +57,7 @@ def map_unique_layer_stack_z(
 def get_layer_overlaps_z(
     layer_stack: LayerStack,
     include_zmax: bool = True,
-):
+) -> dict[float, set[str]]:
     """Maps layers to unique LayerStack z coordinates.
 
     Args:
@@ -67,7 +67,7 @@ def get_layer_overlaps_z(
     """
     z_grid = list_unique_layer_stack_z(layer_stack)
     unique_z_dict = map_unique_layer_stack_z(layer_stack, include_zmax)
-    intersection_z_dict = {}
+    intersection_z_dict: dict[float, set[str]] = {}
     for z in z_grid:
         current_layers = {
             layername for layername, layer_zs in unique_z_dict.items() if z in layer_zs
@@ -77,11 +77,12 @@ def get_layer_overlaps_z(
     return intersection_z_dict
 
 
-def get_layers_at_z(layer_stack: LayerStack, z: float):
+def get_layers_at_z(layer_stack: LayerStack, z: float) -> list[str]:
     """Returns layers present at a given z-position.
 
     Args:
         layer_stack: LayerStack
+        z: float
     Returns:
         List of layers
     """
@@ -93,11 +94,11 @@ def get_layers_at_z(layer_stack: LayerStack, z: float):
         raise ValueError("Requested z-value is above the minimum layer_stack z")
     for z_unique in intersection_z_dict.keys():
         if z <= z_unique:
-            return intersection_z_dict[z_unique]
+            return list(intersection_z_dict[z_unique])
     raise AssertionError("Could not find z-value in layer_stack z-range.")
 
 
-def order_layer_stack(layer_stack: LayerStack):
+def order_layer_stack(layer_stack: LayerStack) -> list[str]:
     """Orders layer_stack according to mesh_order.
 
     Args:

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/common/utils/plot.py

@@ -4,13 +4,17 @@ import re
 from collections.abc import Sequence
 from functools import partial
 from itertools import combinations
+from typing import Any
 
-import gdsfactory as gf
 import matplotlib.pyplot as plt
 import numpy as np
+import numpy.typing as npt
+from matplotlib.axes import Axes
 
 
-def _check_ports(sp: dict[str, np.ndarray], ports: Sequence[str]) -> None:
+def _check_ports(
+    sp: dict[str, npt.NDArray[np.floating[Any]]], ports: Sequence[str]
+) -> None:
     """Ensure ports exist in Sparameters."""
     for port in ports:
         if port not in sp:
@@ -18,7 +22,7 @@ def _check_ports(sp: dict[str, np.ndarray], ports: Sequence[str]) -> None:
 
 
 def plot_sparameters(
-    sp: dict[str, np.ndarray],
+    sp: dict[str, npt.NDArray[np.floating[Any]]],
     logscale: bool = True,
     plot_phase: bool = False,
     keys: tuple[str, ...] | None = None,
@@ -39,7 +43,7 @@ def plot_sparameters(
 
     """
     w = sp["wavelengths"] * units
-    keys = keys or [key for key in sp if not key.lower().startswith("wav")]
+    keys = keys or tuple(key for key in sp if not key.lower().startswith("wav"))
 
     for key in keys:
         if with_simpler_input_keys:
@@ -74,9 +78,12 @@ def plot_sparameters(
 
 
 def plot_imbalance(
-    sp: dict[str, np.ndarray], ports: Sequence[str], ax: plt.Axes | None = None
+    sp: dict[str, npt.NDArray[np.floating[Any]]],
+    ports: Sequence[str],
+    ax: Axes | None = None,
 ) -> None:
     """Plots imbalance in dB for coupler.
+
     The imbalance is always defined between two ports, so this function plots the
     imbalance between all unique port combinations.
 
@@ -107,7 +114,9 @@ def plot_imbalance(
 
 
 def plot_loss(
-    sp: dict[str, np.ndarray], ports: Sequence[str], ax: plt.Axes | None = None
+    sp: dict[str, npt.NDArray[np.floating[Any]]],
+    ports: Sequence[str],
+    ax: Axes | None = None,
 ) -> None:
     """Plots loss dB for coupler.
 
@@ -137,7 +146,9 @@ def plot_loss(
 
 
 def plot_reflection(
-    sp: dict[str, np.ndarray], ports: Sequence[str], ax: plt.Axes | None = None
+    sp: dict[str, npt.NDArray[np.floating[Any]]],
+    ports: Sequence[str],
+    ax: Axes | None = None,
 ) -> None:
     """Plots reflection in dB for coupler.
 
@@ -172,11 +183,3 @@ plot_imbalance1x2 = partial(plot_imbalance, ports=["o1@0,o2@0", "o1@0,o3@0"])
 plot_imbalance2x2 = partial(plot_imbalance, ports=["o1@0,o3@0", "o1@0,o4@0"])
 plot_reflection1x2 = partial(plot_reflection, ports=["o1@0,o1@0"])
 plot_reflection2x2 = partial(plot_reflection, ports=["o1@0,o1@0", "o2@0,o1@0"])
-
-if __name__ == "__main__":
-    import gplugins as sim
-
-    sp = sim.get_sparameters_data_tidy3d(component=gf.components.mmi1x2)
-    # plot_sparameters(sp, logscale=False, keys=["o1@0,o2@0"])
-    # plot_sparameters(sp, logscale=False, keys=["S21"])
-    # plt.show()

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/fdtdz/get_ports_fdtdz.py

@@ -224,8 +224,8 @@ if __name__ == "__main__":
         omega=omega,
         port=c.ports["o1"],
         epsilon=epsilon,
-        xmin=c.dxmin,
-        ymin=c.dymin,
+        xmin=c.xmin,
+        ymin=c.ymin,
         nm_per_pixel=nm_per_pixel,
         port_extent_xy=1,
     )
@@ -234,8 +234,8 @@ if __name__ == "__main__":
         port=c.ports["o1"],
         epsilon_port=epsilon_port,
         excitation=excitation,
-        xmin=c.dxmin,
-        ymin=c.dymin,
+        xmin=c.xmin,
+        ymin=c.ymin,
         zmin=zmin,
         nm_per_pixel=1000,
         figsize=(11, 4),

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/fdtdz/get_sparameters_fdtdz.py

@@ -93,8 +93,8 @@ def get_sparameters_fdtdz(
     """
     # Checks from gmeep
     component_ref = component.ref()
-    component_ref.dx = 0
-    component_ref.dy = 0
+    component_ref.x = 0
+    component_ref.y = 0
 
     optical_port_names = list(component_ref.get_ports_dict(port_type="optical").keys())
 
@@ -129,8 +129,8 @@ def get_sparameters_fdtdz(
             omega=omega,
             port=component.ports[portname],
             epsilon=epsilon,
-            xmin=component_extended.dxmin,
-            ymin=component_extended.dymin,
+            xmin=component_extended.xmin,
+            ymin=component_extended.ymin,
             nm_per_pixel=nm_per_pixel,
             port_extent_xy=port_margin,
         )

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/femwell/mode_solver.py

@@ -68,8 +68,8 @@ def compute_cross_section_modes(
     """
     # Get meshable component from cross-section
     c = gf.components.straight(length=10, cross_section=cross_section)
-    dx = c.dxsize
-    dy = c.dysize
+    dx = c.xsize
+    dy = c.ysize
 
     xsection_bounds = [
         [dx / 2, dy - wafer_padding],

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/gfviz/netlist.py

@@ -167,10 +167,10 @@ def get_ports(child, parent=None):
     else:
         ports = ports1
 
-    x0 = c.dxmin
-    x1 = c.dxmax
-    y0 = c.dymin
-    y1 = c.dymax
+    x0 = c.xmin
+    x1 = c.xmax
+    y0 = c.ymin
+    y1 = c.ymax
     x0, x1 = min(x0, x1), max(x0, x1)
     y0, y1 = min(y0, y1), max(y0, y1)
 
@@ -271,7 +271,7 @@ def _line_dists(points, start, end):
 
 
 def sort_ports(ports) -> list[gf.Port]:
-    return natsorted(ports, key=lambda port: port.dx)
+    return natsorted(ports, key=lambda port: port.x)
 
 
 def wrap_component_in_netlist(name):

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/gmeep/get_port_eigenmode.py

@@ -19,10 +19,10 @@ def get_domain_measurements(sim, output_plane, frequency, resolution=0):
     # Get domain measurements
     sim_center, sim_size = get_2D_dimensions(sim, output_plane)
 
-    xmin = sim_center.dx- sim_size.dx/ 2
-    xmax = sim_center.dx+ sim_size.dx/ 2
-    ymin = sim_center.dy- sim_size.dy/ 2
-    ymax = sim_center.dy+ sim_size.dy/ 2
+    xmin = sim_center.x- sim_size.x/ 2
+    xmax = sim_center.x+ sim_size.x/ 2
+    ymin = sim_center.y- sim_size.y/ 2
+    ymax = sim_center.y+ sim_size.y/ 2
     zmin = sim_center.z - sim_size.z / 2
     zmax = sim_center.z + sim_size.z / 2
 
@@ -31,12 +31,12 @@ def get_domain_measurements(sim, output_plane, frequency, resolution=0):
     Ny = int((ymax - ymin) * grid_resolution + 1)
     Nz = int((zmax - zmin) * grid_resolution + 1)
 
-    if sim_size.dx== 0:
+    if sim_size.x== 0:
         # Plot y on x axis, z on y axis (YZ plane)
         xtics = np.array([sim_center.x])
         ytics = np.linspace(ymin, ymax, Ny)
         ztics = np.linspace(zmin, zmax, Nz)
-    elif sim_size.dy== 0:
+    elif sim_size.y== 0:
         # Plot x on x axis, z on y axis (XZ plane)
         xtics = np.linspace(xmin, xmax, Nx)
         ytics = np.array([sim_center.y])
@@ -117,10 +117,10 @@ def get_port_2Dx_eigenmode(
     # The output of this function is slightly different then MPB (there is no mode_solver object)
     # Format like the Mode objects in gdsfactory/simulation/modes to reuse modes' functions
     if not choose_yz:
-        ny = int(size.dy * sim.resolution)
+        ny = int(size.y * sim.resolution)
         nz = int(size.z * sim.resolution)
         y = np.linspace(
-            center.dy - size.dy / 2, center.dy + size.dy / 2, ny
+            center.y - size.y / 2, center.y + size.y / 2, ny
         )  # eigenmode solver and sim res are technically different
         z = np.linspace(center.z - size.z / 2, center.z + size.z / 2, nz)
     yy, zz = np.meshgrid(y, z, indexing="ij")

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/gmeep/get_simulation.py

@@ -144,8 +144,8 @@ def get_simulation(
 
     dummy_component = gf.Component()
     component_ref = dummy_component << component
-    component_ref.dx = 0
-    component_ref.dy = 0
+    component_ref.x = 0
+    component_ref.y = 0
 
     wavelength = (wavelength_start + wavelength_stop) / 2
 
@@ -195,8 +195,8 @@ def get_simulation(
     cell_thickness = tpml + zmargin_bot + t_core + zmargin_top + tpml if is_3d else 0
 
     cell_size = mp.Vector3(
-        component.dxsize + 2 * tpml,
-        component.dysize + 2 * tpml,
+        component.xsize + 2 * tpml,
+        component.ysize + 2 * tpml,
         cell_thickness,
     )
 

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/gmeep/get_simulation_grating_farfield.py

@@ -172,7 +172,7 @@ def get_simulation_grating_farfield(
         0.5 * sz - pml_thickness + y_offset - 1,
     )
     fiber_port_size = mp.Vector3(3 * fiber_core_diameter, 0, 0)
-    # fiber_port_direction = mp.Vector3(y=-1).drotate(mp.Vector3(z=1), -1 * fiber_angle)
+    # fiber_port_direction = mp.Vector3(y=-1).rotate(mp.Vector3(z=1), -1 * fiber_angle)
 
     waveguide_port_center = mp.Vector3(-sxy / 4)
     waveguide_port_size = mp.Vector3(0, 2 * clad_thickness - 0.2)
@@ -193,8 +193,8 @@ def get_simulation_grating_farfield(
                 x=grating_start + fiber_xposition - fiber_offset_from_angle
             ),
             size=mp.Vector3(fiber_clad, hfiber_geom),
-            e1=mp.Vector3(x=1).drotate(mp.Vector3(z=1), -1 * fiber_angle),
-            e2=mp.Vector3(y=1).drotate(mp.Vector3(z=1), -1 * fiber_angle),
+            e1=mp.Vector3(x=1).rotate(mp.Vector3(z=1), -1 * fiber_angle),
+            e2=mp.Vector3(y=1).rotate(mp.Vector3(z=1), -1 * fiber_angle),
         )
     ]
 
@@ -205,8 +205,8 @@ def get_simulation_grating_farfield(
                 x=grating_start + fiber_xposition - fiber_offset_from_angle
             ),
             size=mp.Vector3(fiber_core_diameter, hfiber_geom),
-            e1=mp.Vector3(x=1).drotate(mp.Vector3(z=1), -1 * fiber_angle),
-            e2=mp.Vector3(y=1).drotate(mp.Vector3(z=1), -1 * fiber_angle),
+            e1=mp.Vector3(x=1).rotate(mp.Vector3(z=1), -1 * fiber_angle),
+            e2=mp.Vector3(y=1).rotate(mp.Vector3(z=1), -1 * fiber_angle),
         )
     )
 
@@ -415,8 +415,8 @@ def get_port_1D_eigenmode(
         frequency=fsrc,
     )
     ys_waveguide = np.linspace(
-        center_waveguide.dy - size_waveguide.dy / 2,
-        center_waveguide.dy + size_waveguide.dy / 2,
+        center_waveguide.y - size_waveguide.y / 2,
+        center_waveguide.y + size_waveguide.y / 2,
         int(sim.resolution * size_waveguide.y),
     )
     x_waveguide = center_waveguide.x
@@ -426,14 +426,14 @@ def get_port_1D_eigenmode(
         direction=mp.NO_DIRECTION,
         where=mp.Volume(center=center_fiber, size=size_fiber),
         band_num=band_num,
-        kpoint=mp.Vector3(0, fsrc * 1.45, 0).drotate(
+        kpoint=mp.Vector3(0, fsrc * 1.45, 0).rotate(
             mp.Vector3(z=1), -1 * np.radians(fiber_angle_deg)
         ),  # Hardcoded index for now, pull from simulation eventually
         frequency=fsrc,
     )
     xs_fiber = np.linspace(
-        center_fiber.dx - size_fiber.dx / 2,
-        center_fiber.dx + size_fiber.dx / 2,
+        center_fiber.x - size_fiber.x / 2,
+        center_fiber.x + size_fiber.x / 2,
         int(sim.resolution * size_fiber.x),
     )
     y_fiber = center_fiber.y

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/gmeep/get_simulation_grating_fiber.py

@@ -501,8 +501,8 @@ def get_port_1D_eigenmode(
         frequency=fsrc,
     )
     ys_waveguide = np.linspace(
-        center_waveguide.dy - size_waveguide.dy / 2,
-        center_waveguide.dy + size_waveguide.dy / 2,
+        center_waveguide.y - size_waveguide.y / 2,
+        center_waveguide.y + size_waveguide.y / 2,
         int(sim.resolution * size_waveguide.y),
     )
     x_waveguide = center_waveguide.x
@@ -518,8 +518,8 @@ def get_port_1D_eigenmode(
         frequency=fsrc,
     )
     xs_fiber = np.linspace(
-        center_fiber.dx - size_fiber.dx / 2,
-        center_fiber.dx + size_fiber.dx / 2,
+        center_fiber.x - size_fiber.x / 2,
+        center_fiber.x + size_fiber.x / 2,
         int(sim.resolution * size_fiber.x),
     )
     y_fiber = center_fiber.y

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/gmeep/meep_adjoint_optimization.py

@@ -151,8 +151,8 @@ def get_meep_adjoint_optimizer(
         Vector3(*cell_size)
         if cell_size
         else Vector3(
-            c.dxsize + 2 * sim.boundary_layers[0].thickness,
-            c.dysize + 2 * sim.boundary_layers[0].thickness,
+            c.xsize + 2 * sim.boundary_layers[0].thickness,
+            c.ysize + 2 * sim.boundary_layers[0].thickness,
             cell_thickness,
         )
     )
@@ -262,10 +262,10 @@ def get_component_from_sim(
     """
     grid_resolution = upscale_factor * sim.resolution
     sim_center, sim_size = get_2D_dimensions(sim, output_plane=None)
-    xmin = sim_center.dx - sim_size.dx / 2
-    xmax = sim_center.dx + sim_size.dx / 2
-    ymin = sim_center.dy - sim_size.dy / 2
-    ymax = sim_center.dy + sim_size.dy / 2
+    xmin = sim_center.x - sim_size.x / 2
+    xmax = sim_center.x + sim_size.x / 2
+    ymin = sim_center.y - sim_size.y / 2
+    ymax = sim_center.y + sim_size.y / 2
     Nx = int((xmax - xmin) * grid_resolution + 1)
     Ny = int((ymax - ymin) * grid_resolution + 1)
     xtics = np.linspace(xmin, xmax, Nx)

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/gmsh/define_polysurfaces.py

@@ -6,16 +6,16 @@ from shapely.affinity import scale
 
 
 def define_polysurfaces(
-    polygons_dict: dict,
+    polygons_dict: dict[str, Any],
     layer_stack: LayerStack,
-    layer_physical_map: dict,
-    layer_meshbool_map: dict,
+    layer_physical_map: dict[str, Any],
+    layer_meshbool_map: dict[str, Any],
     model: Any,
-    resolutions: dict,
+    resolutions: dict[str, Any] | None = None,
     scale_factor: float = 1,
-):
+) -> list[PolySurface]:
     """Define meshwell polysurfaces dimtags from gdsfactory information."""
-    polysurfaces_list = []
+    polysurfaces_list: list[PolySurface] = []
 
     if resolutions is None:
         resolutions = {}
@@ -24,6 +24,11 @@ def define_polysurfaces(
         if polygons_dict[layername].is_empty:
             continue
 
+        layer_stack_ = layer_stack.layers.get(layername)
+
+        if layer_stack_ is None:
+            continue
+
         polysurfaces_list.append(
             PolySurface(
                 polygons=scale(
@@ -33,7 +38,7 @@ def define_polysurfaces(
                 ),
                 model=model,
                 resolution=resolutions.get(layername, None),
-                mesh_order=layer_stack.layers.get(layername).mesh_order,
+                mesh_order=layer_stack_.mesh_order,
                 physical_name=layer_physical_map[layername]
                 if layername in layer_physical_map
                 else layername,

{gplugins-1.3.1 → gplugins-1.3.3}/gplugins/gmsh/get_mesh.py

@@ -98,9 +98,7 @@ def get_mesh(
 
     # Default meshing flags (all True)
     if layer_meshbool_map is None:
-        layer_meshbool_map = {
-            layer_name: True for layer_name in layer_stack.layers.keys()
-        }
+        layer_meshbool_map = dict.fromkeys(layer_stack.layers.keys(), True)
     else:
         for layer_name in layer_stack.layers.keys():
             if layer_name not in layer_physical_map.keys():