emerge 0.6.5__tar.gz → 0.6.6__tar.gz

{emerge-0.6.5 → emerge-0.6.6}/.bumpversion.toml

@@ -1,5 +1,5 @@
 [tool.bumpversion]
-current_version = "0.6.5"
+current_version = "0.6.6"
 parse = "(?P<major>\\d+)\\.(?P<minor>\\d+)\\.(?P<patch>\\d+)"
 serialize = ["{major}.{minor}.{patch}"]
 search = "{current_version}"

{emerge-0.6.5 → emerge-0.6.6}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: emerge
-Version: 0.6.5
+Version: 0.6.6
 Summary: An open source EM FEM simulator in Python
 Project-URL: Homepage, https://github.com/FennisRobert/EMerge
 Project-URL: Issues, https://github.com/FennisRobert/EMerge/issues

{emerge-0.6.5 → emerge-0.6.6}/emerge/__init__.py

@@ -18,7 +18,7 @@ along with this program; if not, see
 """
 import os
 
-__version__ = "0.6.5"
+__version__ = "0.6.6"
 
 ############################################################
 #               HANDLE ENVIRONMENT VARIABLES              #
@@ -52,6 +52,7 @@ from ._emerge.cs import CoordinateSystem, CS, GCS, Plane, Axis, XAX, YAX, ZAX, X
 from ._emerge.coord import Line
 from ._emerge import geo
 from ._emerge.selection import Selection, FaceSelection, DomainSelection, EdgeSelection
+from ._emerge.geometry import select
 from ._emerge.mth.common_functions import norm, coax_rout, coax_rin
 from ._emerge.physics.microwave.sc import stratton_chu
 from ._emerge.periodic import RectCell, HexCell

{emerge-0.6.5 → emerge-0.6.6}/emerge/_emerge/geo/__init__.py

@@ -19,6 +19,6 @@ from .pcb import PCB
 from .pmlbox import pmlbox
 from .horn import Horn
 from .shapes import Cylinder, CoaxCylinder, Box, XYPlate, HalfSphere, Sphere, Plate, OldBox, Alignment, Cone
-from .operations import subtract, add, embed, remove, rotate, mirror, change_coordinate_system, translate, intersect
+from .operations import subtract, add, embed, remove, rotate, mirror, change_coordinate_system, translate, intersect, unite
 from .polybased import XYPolygon, GeoPrism, Disc, Curve
 from .step import STEPItems

{emerge-0.6.5 → emerge-0.6.6}/emerge/_emerge/geo/operations.py

@@ -15,8 +15,8 @@
 # along with this program; if not, see
 # <https://www.gnu.org/licenses/>.
 
-from typing import TypeVar
-from ..geometry import GeoSurface, GeoVolume, GeoObject, GeoPoint, GeoEdge
+from typing import TypeVar, overload
+from ..geometry import GeoSurface, GeoVolume, GeoObject, GeoPoint, GeoEdge, GeoPolygon
 from ..cs import CoordinateSystem, GCS
 import gmsh
 import numpy as np
@@ -144,7 +144,8 @@ def rotate(main: GeoVolume,
            c0: tuple[float, float, float],
            ax: tuple[float, float, float],
            angle: float,
-           degree=True) -> GeoVolume:
+           make_copy: bool = False,
+           degree=True) -> GeoObject:
     """Rotates a GeoVolume object around an axist defined at a coordinate.
 
     Args:
@@ -159,17 +160,23 @@ def rotate(main: GeoVolume,
     """
     if degree:
         angle = angle * np.pi/180
-    gmsh.model.occ.rotate(main.dimtags, *c0, *ax, -angle)
-
+    
+    if make_copy:
+        rotate_obj = main.make_copy()
+    else:
+        rotate_obj = main
+     
+    gmsh.model.occ.rotate(rotate_obj.dimtags, *c0, *ax, -angle)
     # Rotate the facepointers
-    for fp in main._all_pointers:
+    for fp in rotate_obj._all_pointers:
         fp.rotate(c0, ax, angle)
-    return main
+    return rotate_obj
 
 def translate(main: GeoVolume,
               dx: float = 0,
               dy: float = 0,
-              dz: float = 0) -> GeoVolume:
+              dz: float = 0,
+              make_copy: bool = False) -> GeoObject:
     """Translates the GeoVolume object along a given displacement
 
     Args:
@@ -177,17 +184,23 @@ def translate(main: GeoVolume,
         dx (float, optional): The X-displacement in meters. Defaults to 0.
         dy (float, optional): The Y-displacement in meters. Defaults to 0.
         dz (float, optional): The Z-displacement in meters. Defaults to 0.
+        make_copy (bool, optional): Whether to make a copy first before translating.
 
     Returns:
-        GeoVolume: The translated object
+        GeoObject: The translated object
     """
-    gmsh.model.occ.translate(main.dimtags, dx, dy, dz)
+    
+    if make_copy:
+        trans_obj = main.make_copy()
+    else:
+        trans_obj = main
+    gmsh.model.occ.translate(trans_obj.dimtags, dx, dy, dz)
     
      # Rotate the facepointers
-    for fp in main._all_pointers:
+    for fp in trans_obj._all_pointers:
         fp.translate(dx, dy, dz)
 
-    return main
+    return trans_obj
 
 def mirror(main: GeoObject,
            origin: tuple[float, float, float] = (0.0, 0.0, 0.0),
@@ -199,6 +212,7 @@ def mirror(main: GeoObject,
         main (GeoVolume): The object to mirror
         origin (tuple[float, float, float], optional): The point of origin in meters. Defaults to (0.0, 0.0, 0.0).
         direction (tuple[float, float, float], optional): The normal axis defining the plane of reflection. Defaults to (1.0, 0.0, 0.0).
+        make_copy (bool, optional): Whether to make a copy first before mirroring.
 
     Returns:
         GeoVolume: The mirrored GeoVolume object
@@ -211,14 +225,12 @@ def mirror(main: GeoObject,
     d = -(a*x0 + b*y0 + c*z0)
     if (a==0) and (b==0) and (c==0):
         return main
+    
     mirror_obj = main
     if make_copy:
-        new_obj = main.make_copy()
-        gmsh.model.occ.mirror(new_obj.dimtags, a,b,c,d)
-        mirror_obj = new_obj
-    else:
-        gmsh.model.occ.mirror(main.dimtags, a,b,c,d)
-
+        mirror_obj = main.make_copy()
+    gmsh.model.occ.mirror(mirror_obj.dimtags, a,b,c,d)
+    
     for fp in mirror_obj._all_pointers:
         fp.mirror(origin, direction)
     return mirror_obj
@@ -237,7 +249,7 @@ def change_coordinate_system(main: GeoObject,
         old_cs (CoordinateSystem, optional): The old coordinate system. Defaults to GCS.
 
     Returns:
-        _type_: _description_
+        GeoObject: The output object
     """
     if new_cs._is_global and old_cs._is_global:
         return main
@@ -254,4 +266,39 @@ def change_coordinate_system(main: GeoObject,
     for fp in main._all_pointers:
         fp.affine_transform(M1)
         fp.affine_transform(M2)
-    return main
+    return main
+
+@overload
+def unite(*objects: GeoVolume) -> GeoVolume: ...
+
+@overload
+def unite(*objects: GeoSurface) -> GeoSurface: ...
+
+@overload
+def unite(*objects: GeoEdge) -> GeoEdge: ...
+
+@overload
+def unite(*objects: GeoPolygon) -> GeoSurface: ...
+
+def unite(*objects: GeoObject) -> GeoObject:
+    """Applies a fusion consisting of all geometries in the argument.
+
+    Returns:
+        GeoObject: The resultant object
+    """
+    main, *rest = objects
+    if not rest:
+        return main
+    main._exists = False
+    dts = []
+    for other in rest:
+        dts.extend(other.dimtags)
+        other._exists = False
+    
+    new_dimtags, mapping = gmsh.model.occ.fuse(main.dimtags, dts)
+    
+    new_obj = GeoObject.from_dimtags(new_dimtags)._take_tools(*objects)
+    new_obj.set_material(main.material)
+    new_obj.prio_set(main._priority)
+    
+    return new_obj

{emerge-0.6.5 → emerge-0.6.6}/emerge/_emerge/geo/pcb.py

@@ -22,7 +22,7 @@ from ..geometry import GeoPolygon, GeoVolume, GeoSurface
 from ..material import Material, AIR, COPPER
 from .shapes import Box, Plate, Cylinder
 from .polybased import XYPolygon
-from .operations import change_coordinate_system
+from .operations import change_coordinate_system, unite
 from .pcb_tools.macro import parse_macro
 from .pcb_tools.calculator import PCBCalculator
 
@@ -310,7 +310,6 @@ class StripCurve(StripTurn):
         points: list[tuple[float, float]] = []
         Npts = int(np.ceil(abs(self.angle/self.dang)))
         R = self.radius-np.sign(self.angle)*self.width/2
-        print(R, self.circ_origin, self._xhat, self._yhat)
         for i in range(Npts):
             ang = abs((i+1)/Npts * self.angle * np.pi/180)
             pnew = self.circ_origin + R*(self._xhat*np.cos(ang)+self._yhat*np.sin(ang))
@@ -324,7 +323,6 @@ class StripCurve(StripTurn):
 
         Npts = int(np.ceil(abs(self.angle/self.dang)))
         R = self.radius+np.sign(self.angle)*self.width/2
-        print(R, self.circ_origin, self._xhat, self._yhat)
         for i in range(Npts):
             ang = abs((i+1)/Npts * self.angle * np.pi/180)
             pnew = self.circ_origin + R*(self._xhat*np.cos(ang)+self._yhat*np.sin(ang))
@@ -1247,7 +1245,13 @@ class PCB:
         plate = change_coordinate_system(plate, self.cs)
         return plate # type: ignore
 
-    def generate_vias(self, merge=False) -> list[Cylinder] | Cylinder:
+    @overload
+    def generate_vias(self, merge=Literal[True]) -> GeoVolume: ...
+    
+    @overload
+    def generate_vias(self, merge=Literal[False]) -> list[Cylinder]: ...
+        
+    def generate_vias(self, merge=False) -> list[Cylinder] | GeoVolume:
         """Generates the via objects.
 
         Args:
@@ -1321,7 +1325,7 @@ class PCB:
         Returns:
             list[Polygon] | GeoSurface: The output stripline polygons possibly merged if merge = True.
         """
-        polys = []
+        polys: list[GeoSurface] = []
         allx = []
         ally = []
 
@@ -1359,13 +1363,14 @@ class PCB:
 
         self.traces = polys
         if merge:
-            tags = []
-            for p in polys:
-                tags.extend(p.tags)
-                if p.material != COPPER:
-                    logger.warning(f'Merging a polygon with material {p.material} into a single polygon that will be COPPER.')
-            polys = GeoSurface(tags)
-            polys.material = COPPER
+            polys = unite(*polys)
+            # tags = []
+            # for p in polys:
+            #     tags.extend(p.tags)
+            #     if p.material != COPPER:
+            #         logger.warning(f'Merging a polygon with material {p.material} into a single polygon that will be COPPER.')
+            # polys = GeoSurface(tags)
+            # polys.material = COPPER
         return polys
                 
 ############################################################

{emerge-0.6.5 → emerge-0.6.6}/emerge/_emerge/geo/shapes.py

@@ -31,9 +31,18 @@ class Alignment(Enum):
     CORNER = 2
 
 class Box(GeoVolume):
-    """ A class that represents a box shaped volume
+    """Creates a box volume object.
+        Specify the alignment of the box with the provided position. The options are CORNER (default)
+        for the front-left-bottom node of the box or CENTER for the center of the box.
 
-    """
+        Args:
+            width (float): The x-size
+            depth (float): The y-size
+            height (float): The z-size
+            position (tuple, optional): The position of the box. Defaults to (0,0,0).
+            alignment (Alignment, optional): Which point of the box is placed at the position. 
+                Defaults to Alignment.CORNER.
+        """
 
     def __init__(self, 
                  width: float, 
@@ -94,7 +103,12 @@ class Box(GeoVolume):
 
     
 class Sphere(GeoVolume):
+    """Generates a sphere objected centered ont he position with the given radius
 
+    Args:
+        radius (float): The sphere radius
+        position (tuple, optional): The center position. Defaults to (0,0,0).
+    """
     def __init__(self, 
                  radius: float,
                  position: tuple = (0,0,0)):
@@ -109,6 +123,17 @@ class Sphere(GeoVolume):
         self.tags: list[int] = [gmsh.model.occ.addSphere(x,y,z,radius),]
 
 class XYPlate(GeoSurface):
+    """Generates and XY-plane oriented plate
+        
+        Specify the alignment of the plate with the provided position. The options are CORNER (default)
+        for the front-left node of the plate or CENTER for the center of the plate.
+
+        Args:
+            width (float): The x-size of the plate
+            depth (float): The y-size of the plate
+            position (tuple, optional): The position of the alignment node. Defaults to (0,0,0).
+            alignment (Alignment, optional): Which node to align to. Defaults to Alignment.CORNER.
+        """
     def __init__(self, 
                  width: float, 
                  depth: float, 
@@ -134,7 +159,19 @@ class XYPlate(GeoSurface):
 
 
 class Plate(GeoSurface):
-        
+    """A generalized 2D rectangular plate in XYZ-space.
+
+        The plate is specified by an origin (o) in meters coordinate plus two vectors (u,v) in meters
+        that span two of the sides such that all points of the plate are defined by:
+            p1 = o
+            p2 = o+u
+            p3 = o+v
+            p4 = o+u+v
+        Args:
+            origin (tuple[float, float, float]): The origin of the plate in meters
+            u (tuple[float, float, float]): The u-axis of the plate
+            v (tuple[float, float, float]): The v-axis of the plate
+        """
     def __init__(self,
                 origin: tuple[float, float, float],
                 u: tuple[float, float, float],
@@ -297,7 +334,6 @@ class CoaxCylinder(GeoVolume):
 
         xo, yo, zo = self.cs.in_global_cs(x.flatten(), y.flatten(), z.flatten())
         return xo, yo, zo
-        return super().boundary()
     
 class HalfSphere(GeoVolume):
 

{emerge-0.6.5 → emerge-0.6.6}/emerge/_emerge/geometry.py

@@ -411,8 +411,8 @@ class GeoObject:
         self._priority = _GEOMANAGER.highest_priority()+10
         return self
     
-    def outside(self, *exclude: FaceNames, tags: list[int] | None = None) -> FaceSelection:
-        """Returns the complete set of outside faces.
+    def boundary(self, exclude: tuple[FaceNames,...] | None = None, tags: list[int] | None = None) -> FaceSelection:
+        """Returns the complete set of boundary faces.
 
         If implemented, it is possible to exclude a set of faces based on their name
         or a list of tags (integers)
@@ -420,8 +420,13 @@ class GeoObject:
         Returns:
             FaceSelection: The selected faces
         """
+        if exclude is None:
+            exclude = tuple()
+            
         if tags is None:
             tags = []
+        
+        
         for name in exclude:
             tags.extend(self.face(name).tags)
         dimtags = gmsh.model.get_boundary(self.dimtags, True, False)
@@ -441,6 +446,22 @@ class GeoObject:
         
         return FaceSelection(self._face_tags(name, tool))
 
+    def faces(self, *names: FaceNames, tool: GeoObject | None = None) -> FaceSelection:
+        """Returns the FaceSelection for a given face names.
+        
+        The face name must be defined for the type of geometry.
+
+        Args:
+            name (FaceNames): The name of the face to select.
+
+        Returns:
+            FaceSelection: The selected face
+        """
+        tags = []
+        for name in names:
+            tags.extend(self._face_tags(name, tool))
+        return FaceSelection(tags)
+        
     @property
     def dimtags(self) -> list[tuple[int, int]]:
         return [(self.dim, tag) for tag in self.tags]
@@ -448,15 +469,6 @@ class GeoObject:
     @property
     def embeddings(self) -> list[tuple[int,int]]:
         return []
-    
-    def boundary(self) -> FaceSelection:
-        if self.dim == 3:
-            tags = gmsh.model.get_boundary(self.dimtags, oriented=False)
-            return FaceSelection([t[1] for t in tags])
-        if self.dim == 2:
-            return FaceSelection(self.tags)
-        else:
-            raise ValueError('Can only generate faces for objects of dimension 2 or higher.')
 
     @staticmethod
     def from_dimtags(dimtags: list[tuple[int,int]]) -> GeoVolume | GeoSurface | GeoObject:
@@ -549,6 +561,25 @@ class GeoPolygon(GeoSurface):
         self.points: list[int] = []
         self.lines: list[int] = []
 
-    
-T = TypeVar('T', GeoVolume, GeoEdge, GeoPoint, GeoSurface)
-    
+
+############################################################
+#                      SHORT FUNCTIONS                     #
+############################################################
+
+def select(*items: Selection | GeoObject) -> Selection:
+    """Generate a selection from a series of selections and/or objects that share the same dimension.
+
+    Raises:
+        ValueError: Raised if the dimensions provided are not consistent
+
+    Returns:
+        Selection: An output selection object.
+    """
+    dim = items[0].dim
+    tags = []
+    for item in items:
+        if item.dim!=dim:
+            raise ValueError(f'Cannot group of objects with a dissimilar dimensions. Trying to include {item} in a list of dimension {dim}.')
+        tags.extend(item.tags)
+    return Selection.from_dim_tags(dim, tags)
+            

{emerge-0.6.5 → emerge-0.6.6}/emerge/_emerge/howto.py

@@ -69,7 +69,7 @@ class _HowtoClass:
         
         The naming convention is left(-X), right(+X), front(-Y), back(+Y), bottom(-Z), top(+Z)
         All outside faces can be selected using
-        >>> outside = object.outside()
+        >>> outside = object.boundary()
         
         If objects are the results from operations, you can access the faces from the
         source objects using the optional tool argument

{emerge-0.6.5 → emerge-0.6.6}/emerge/_emerge/selection.py

@@ -590,3 +590,4 @@ class Selector:
         return FaceSelection([lowest_key,])
 
 SELECTOR_OBJ = Selector()
+

{emerge-0.6.5 → emerge-0.6.6}/examples/demo10_sgh.py

@@ -29,7 +29,7 @@ dx = 2 * mm              # distance from horn exit to PML start
 
 # Create simulation object
 m = em.Simulation('HornAntenna')
-m.check_version("0.6.5") # Checks version compatibility.
+m.check_version("0.6.6") # Checks version compatibility.
 # --- Coordinate system for horn geometry -------------------------------
 hornCS = em.CS(em.YAX, em.ZAX, em.XAX)
 
@@ -81,7 +81,7 @@ m.generate_mesh()
 # --- Boundary conditions ------------------------------------------------
 p1 = m.mw.bc.ModalPort(feed.face('left'), 1)     # excite TE10 in feed
 PMC = m.mw.bc.PMC(m.select.face.inplane(0, 0, 0, 0, 1, 0))  # perfect magnetic on symmetry
-radiation_boundary = air.outside('front', 'left', 'bottom')  # open faces
+radiation_boundary = air2.faces('back','top','right', tool=air)  # open faces
 abc = m.mw.bc.AbsorbingBoundary(m.select.face.inplane(Lhorn-dx,0,0,1,0,0))
 # View mesh and BC selections
 m.view(selections=[p1.selection, PMC.selection, radiation_boundary])

{emerge-0.6.5 → emerge-0.6.6}/examples/demo11_lumped_element_filter.py

@@ -29,7 +29,7 @@ def Cf(C):
 pack = '0603'         # package footprint for lumped components
 # Create simulation and PCB layouter with substrate thickness and material
 m = em.Simulation('LumpedFilter')
-m.check_version("0.6.5") # Checks version compatibility.
+m.check_version("0.6.6") # Checks version compatibility.
 
 th = 0.5         # substrate thickness (meters)
 Hair = 2.0
@@ -99,7 +99,7 @@ for le in LEs:
     m.mw.bc.LumpedElement(le)
 # Perfect conductor on copper traces and vias
 m.mw.bc.PEC(traces)
-m.mw.bc.PEC(vias.outside())
+m.mw.bc.PEC(vias.boundary())
 
 # --- Run frequency-domain simulation ------------------------------------
 data = m.mw.run_sweep(parallel=True, njobs=4, frequency_groups=8)

{emerge-0.6.5 → emerge-0.6.6}/examples/demo12_mode_alignment.py

@@ -12,7 +12,7 @@ direction (+Z vs -Z). In EMerge we can alignm modes using the .align_mode() meth
 
 # First we define our simulation
 m = em.Simulation('aligntest')
-m.check_version("0.6.5") # Checks version compatibility.
+m.check_version("0.6.6") # Checks version compatibility.
 
 # We create a cyllindrical waveguide in the Y-axis.
 cyl = em.geo.Cylinder(0.012, 0.05, em.CS(em.XAX, em.ZAX, em.YAX))

{emerge-0.6.5 → emerge-0.6.6}/examples/demo13_helix_antenna.py

@@ -25,7 +25,7 @@ porth = 2*mm                    # vertical height of the feed extrusion
 
 # --- Simulation object -------------------------------------------------------
 model = em.Simulation('helix')
-model.check_version("0.6.5") # Checks version compatibility.
+model.check_version("0.6.6") # Checks version compatibility.
 
 dfeed = 3*mm                    # straight feed length before the helix starts
 
@@ -78,8 +78,8 @@ model.generate_mesh()
 model.view()
 
 # --- Boundary selections for BCs & port -------------------------------------
-abc_sel = airbox.outside('bottom')     # absorbing boundary applied at bottom face of airbox
-port_sel = feed.outside('front','back')# port faces on the feed (two opposite faces)
+abc_sel = airbox.boundary(exclude=('bottom',))     # absorbing boundary applied at bottom face of airbox
+port_sel = feed.boundary(exclude=('front','back'))# port faces on the feed (two opposite faces)
 
 # --- Boundary conditions -----------------------------------------------------
 abc = model.mw.bc.AbsorbingBoundary(abc_sel)                     # open-space termination

{emerge-0.6.5 → emerge-0.6.6}/examples/demo14_boundary_selection.py

@@ -37,7 +37,7 @@ wgb = 10.16*mm
 L = 50*mm
 
 model = em.Simulation('Test Mode')
-model.check_version("0.6.5") # Checks version compatibility.
+model.check_version("0.6.6") # Checks version compatibility.
 
 # first lets define a WR90 waveguide
 wg_box = em.geo.Box(L, wga, wgb, position=(-L, -wga/2, -wgb/2))
@@ -71,7 +71,7 @@ feed_port = wg_box_new.face('left', tool=wg_box)
 
 # We can also select the outside and exclude a given face. Because our airbox
 # is not modified, we don't have to work with tools.
-radiation_boundary = airbox.outside('left')
+radiation_boundary = airbox.boundary(exclude=('left',))
 
 # Lets view our result
 model.view(selections=[feed_port, radiation_boundary])

{emerge-0.6.5 → emerge-0.6.6}/examples/demo1_stepped_imp_filter.py

@@ -28,7 +28,7 @@ pcbmat = em.Material(er=2.2, tand=0.00, color="#217627", opacity=0.2)
 # We start by creating our simulation object.
 
 m = em.Simulation('Demo1_SIF', loglevel='INFO')
-m.check_version("0.6.5") # Checks version compatibility.
+m.check_version("0.6.6") # Checks version compatibility.
 # To accomodate PCB routing we make use of the PCBLayouter class. To use it we need to 
 # supply it with a thickness, the desired air-box height, the units at which we supply
 # the dimensions and the PCB material.

{emerge-0.6.5 → emerge-0.6.6}/examples/demo2_combline_filter.py

@@ -43,7 +43,7 @@ lfeed = 100*mil
 
 # A usual we start our simulation file
 model = em.Simulation('Combline_DEMO')
-model.check_version("0.6.5") # Checks version compatibility.
+model.check_version("0.6.6") # Checks version compatibility.
 
 # The filter consists of quarter lamba cylindrical pins inside an airbox.
 # First we create the airbox

{emerge-0.6.5 → emerge-0.6.6}/examples/demo3_coupled_line_filter.py

@@ -38,7 +38,7 @@ extra = 100             # extra margin (mil)
 
 # --- Simulation setup ----------------------------------------------------
 model = em.Simulation('Demo3')
-model.check_version("0.6.5") # Checks version compatibility.
+model.check_version("0.6.6") # Checks version compatibility.
 # --- Material and layouter -----------------------------------------------
 mat = em.Material(er=3.55, color="#488343", opacity=0.1)
 # Create PCB layouter with given substrate thickness and units

{emerge-0.6.5 → emerge-0.6.6}/examples/demo4_patch_antenna.py

@@ -36,7 +36,7 @@ f2 = 1.6e9              # stop frequency
 # --- Create simulation object -------------------------------------------
 # Using PVDisplay backend for 3D visualization
 model = em.Simulation('MyPatchAntenna')
-model.check_version("0.6.5") # Checks version compatibility.
+model.check_version("0.6.6") # Checks version compatibility.
 # --- Define geometry primitives -----------------------------------------
 # Substrate block centered at origin in XY, thickness in Z (negative down)
 dielectric = em.geo.Box(wsub, hsub, th,
@@ -103,7 +103,7 @@ port = model.mw.bc.LumpedPort(
     active=True, Z0=50
 )
 # Apply absorbing boundary on underside of airbox to simulate open space
-boundary_selection = air.outside('bottom')
+boundary_selection = air.boundary(exclude=('bottom',))
 
 model.view(selections=[boundary_selection,])
 abc = model.mw.bc.AbsorbingBoundary(boundary_selection)

{emerge-0.6.5 → emerge-0.6.6}/examples/demo5_revolve.py

@@ -4,7 +4,7 @@ import emerge as em
 For now it just shows you how to work with the revolve system.
 """
 model = em.Simulation('Revolve test')
-model.check_version("0.6.5") # Checks version compatibility.
+model.check_version("0.6.6") # Checks version compatibility.
 
 mm = 0.001
 rad_feed = 30*mm

{emerge-0.6.5 → emerge-0.6.6}/examples/demo6_striplines_with_vias.py

@@ -12,7 +12,7 @@ mm = 0.001
 th = 1
 
 model = em.Simulation('Stripline_test')
-model.check_version("0.6.5") # Checks version compatibility.
+model.check_version("0.6.6") # Checks version compatibility.
 
 # As usual we start by creating our layouter
 ly = em.geo.PCB(th, mm, em.GCS, em.lib.DIEL_RO4350B)
@@ -59,9 +59,9 @@ model.mesher.set_boundary_size(trace, 0.001)
 
 model.generate_mesh()
 
-# We display the geometry with extra attention to the vias. With the vias.outside() method we can
+# We display the geometry with extra attention to the vias. With the vias.boundary() method we can
 # specifically show the outside faces of the via.
-model.view(selections=[vias.outside()])
+model.view(selections=[vias.boundary()])
 
 # We setup the lumped port boundary conditions. Because of an added functionality in the PCBLayouter 
 # class, you don't have to specify the width, height and direction of the lumped port, this information
@@ -74,7 +74,7 @@ p2 = model.mw.bc.LumpedPort(lp2, 2)
 pec = model.mw.bc.PEC(trace)
 
 #We also add a PEC for the outsides of our via.
-pecvia = model.mw.bc.PEC(vias.outside())
+pecvia = model.mw.bc.PEC(vias.boundary())
 
 # Finally we run the simulation!
 data = model.mw.run_sweep(True, 4, frequency_groups=8)

{emerge-0.6.5 → emerge-0.6.6}/examples/demo7_periodic_cells.py

@@ -29,7 +29,7 @@ fl = 25*mm
 
 # We start again by defining our simulation model
 model = em.Simulation('Periodic')
-model.check_version("0.6.5") # Checks version compatibility.
+model.check_version("0.6.6") # Checks version compatibility.
 
 # Next we will create a PeriodicCell class (in our case a hexagonal cell). This class
 # is simply meant to simplify our lives and improve the simulation setup flow.

{emerge-0.6.5 → emerge-0.6.6}/examples/demo8_waveguide_bpf_synthesis.py

@@ -80,7 +80,7 @@ wgaps = np.linspace(1*mm, 20*mm, 21)
 Ks = []
 hphis = []
 with em.Simulation('IrisSim') as sim:
-    sim.check_version("0.6.5") # Checks version compatibility.
+    sim.check_version("0.6.6") # Checks version compatibility.
     for (wgap,) in sim.parameter_sweep(True, wgap=wgaps):
         # Define two short waveguide sections separated by iris plate
         wg1 = em.geo.Box(wga, Lfeed, wgb, (-wga/2, -Lfeed - t_thickness/2, 0))

{emerge-0.6.5 → emerge-0.6.6}/examples/demo9_dielectric_resonator.py

@@ -33,7 +33,7 @@ Nmodes = 5
 
 # --- Create simulation ---------------------------------------------------
 model = em.Simulation('DielectricResonatorFilter')
-model.check_version("0.6.5") # Checks version compatibility.
+model.check_version("0.6.6") # Checks version compatibility.
 
 # --- Build geometry ------------------------------------------------------
 # Metal enclosure box (PEC by default)

{emerge-0.6.5 → emerge-0.6.6}/pyproject.toml

@@ -3,7 +3,7 @@ allow-direct-references = true
 
 [project]
 name = "emerge"
-version = "0.6.5"
+version = "0.6.6"
 description = "An open source EM FEM simulator in Python"
 readme = "README.md"
 requires-python = ">=3.10, <4.0"

{emerge-0.6.5 → emerge-0.6.6}/uv.lock

@@ -236,7 +236,7 @@ wheels = [
 
 [[package]]
 name = "emerge"
-version = "0.6.5"
+version = "0.6.6"
 source = { editable = "." }
 dependencies = [
     { name = "gmsh" },