I fixed most of the problems with Scene.plot(). But Polyslabs are not transposed correctly. It looks like I need to modify the "polyslab.py" file. (I haven't touched this file yet.)

Author: jewettaijfc

tidy3d/components/geometry/base.py

@@ -2181,7 +2181,7 @@ def intersections_with(self, other, transpose: bool = False):
         dim = "xyz"[normal_ind]
         pos = self.center[normal_ind]
         xyz_kwargs = {dim: pos}
-        shapes_plane = other.intersections_plane(**xyz_kwargs)
+        shapes_plane = other.intersections_plane(**xyz_kwargs, transpose=transpose)
 
         # intersect all shapes with the input self
         bs_min, bs_max = (self.pop_axis(bounds, axis=normal_ind, transpose=transpose)[1] for bounds in self.bounds)

tidy3d/components/scene.py

@@ -319,7 +319,7 @@ def intersecting_media(
 
     @staticmethod
     def intersecting_structures(
-        test_object: Box, structures: tuple[Structure, ...]
+        test_object: Box, structures: tuple[Structure, ...], transpose: bool = False,
     ) -> tuple[Structure, ...]:
         """From a given list of structures, returns a list of :class:`.Structure` that intersect
         with the ``test_object``, if it is a surface, or its surfaces, if it is a volume.
@@ -346,7 +346,7 @@ def intersecting_structures(
 
             structures_merged = []
             for structure in structures:
-                intersections = structure.geometry.intersections_plane(**xyz_kwargs)
+                intersections = structure.geometry.intersections_plane(**xyz_kwargs, transpose=transpose)
                 if len(intersections) > 0:
                     structures_merged.append(structure)
             return structures_merged
@@ -648,8 +648,8 @@ def _get_structures_2dbox(
         """
         # if no hlim and/or vlim given, the bounds will then be the usual pml bounds
         axis, _ = Box.parse_xyz_kwargs(x=x, y=y, z=z)
-        _, (hmin, vmin) = Box.pop_axis(self.bounds[0], axis=axis)
-        _, (hmax, vmax) = Box.pop_axis(self.bounds[1], axis=axis)
+        _, (hmin, vmin) = Box.pop_axis(self.bounds[0], axis=axis, transpose=transpose)
+        _, (hmax, vmax) = Box.pop_axis(self.bounds[1], axis=axis, transpose=transpose)
 
         if hlim is not None:
             (hmin, hmax) = hlim