emerge 1.0.2__tar.gz → 1.0.3__tar.gz

{emerge-1.0.2 → emerge-1.0.3}/.bumpversion.toml

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

{emerge-1.0.2 → emerge-1.0.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: emerge
-Version: 1.0.2
+Version: 1.0.3
 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-1.0.2 → emerge-1.0.3}/emerge/__init__.py

@@ -18,7 +18,7 @@ along with this program; if not, see
 """
 import os
 
-__version__ = "1.0.2"
+__version__ = "1.0.3"
 
 ############################################################
 #               HANDLE ENVIRONMENT VARIABLES              #

{emerge-1.0.2 → emerge-1.0.3}/emerge/_emerge/physics/microwave/microwave_bc.py

@@ -96,6 +96,15 @@ class MWBoundaryConditionSet(BoundaryConditionSet):
         self._cell._ports.append(port)
         return port
 
+    # Checks
+    def _is_excited(self) -> bool:
+        for bc in self.boundary_conditions:
+            if not isinstance(bc, RobinBC):
+                continue
+            if bc._include_force:
+                return True
+            
+        return False
 
 ############################################################
 #                    BOUNDARY CONDITIONS                   #

{emerge-1.0.2 → emerge-1.0.3}/emerge/_emerge/plot/pyvista/display.py

@@ -213,11 +213,13 @@ class _AnimObject:
                  field: np.ndarray,
                  T: Callable,
                  grid: pv.Grid,
+                 filtered_grid: pv.Grid,
                  actor: pv.Actor,
                  on_update: Callable):
         self.field: np.ndarray = field
         self.T: Callable = T
         self.grid: pv.Grid = grid
+        self.fgrid: pv.Grid = filtered_grid
         self.actor: pv.Actor = actor
         self.on_update: Callable = on_update
 
@@ -611,7 +613,8 @@ class PVDisplay(BaseDisplay):
         
         grid = pv.StructuredGrid(x,y,z)
         field_flat = field.flatten(order='F')
-
+        
+        
 
         if scale=='log':
             T = lambda x: np.log10(np.abs(x+1e-12))
@@ -626,10 +629,12 @@ class PVDisplay(BaseDisplay):
         else:
             name = _fieldname
         self._ctr += 1
+        
         grid[name] = static_field
 
         grid_no_nan = grid.threshold(scalars=name)
-
+        
+        
         # Determine color limits
         if clim is None:
             fmin = np.nanmin(static_field)
@@ -647,7 +652,9 @@ class PVDisplay(BaseDisplay):
             def on_update(obj: _AnimObject, phi: complex):
                 field_anim = obj.T(np.real(obj.field * phi))
                 obj.grid[name] = field_anim
-            self._objs.append(_AnimObject(field_flat, T, grid_no_nan, actor, on_update))
+                obj.fgrid[name] = obj.grid.threshold(scalars=name)[name]
+                #obj.fgrid replace with thresholded scalar data.
+            self._objs.append(_AnimObject(field_flat, T, grid, grid_no_nan, actor, on_update))
 
 
     def add_title(self, title: str) -> None:

{emerge-1.0.2 → emerge-1.0.3}/emerge/_emerge/solver.py

@@ -29,10 +29,6 @@ import platform
 import time
 from typing import Literal, Callable
 from enum import Enum
-import importlib.util
-
-def module_exists(name: str) -> bool:
-    return importlib.util.find_spec(name) is not None
 
 _PARDISO_AVAILABLE = False
 _UMFPACK_AVAILABLE = False
@@ -56,26 +52,26 @@ if 'arm' not in platform.processor():
 #                          UMFPACK                         #
 ############################################################
 
-if module_exists('scikits.umfpack'):
-    try:
-        import scikits.umfpack as um # type: ignore
-        _UMFPACK_AVAILABLE = True
-    except ModuleNotFoundError:
-        logger.debug('UMFPACK not found, defaulting to SuperLU')
+
+try:
+    import scikits.umfpack as um # type: ignore
+    _UMFPACK_AVAILABLE = True
+except ModuleNotFoundError:
+    logger.debug('UMFPACK not found, defaulting to SuperLU')
 
 ############################################################
 #                           CUDSS                          #
 ############################################################
 
-if module_exists('cupy'):
-    try:
-        from .solve_interfaces.cudss_interface import CuDSSInterface
-        _CUDSS_AVAILABLE = True
-    except ModuleNotFoundError:
-        pass
-    except ImportError as e:
-        logger.error('Error while importing CuDSS dependencies:')
-        logger.exception(e)
+
+try:
+    from .solve_interfaces.cudss_interface import CuDSSInterface
+    _CUDSS_AVAILABLE = True
+except ModuleNotFoundError:
+    pass
+except ImportError as e:
+    logger.error('Error while importing CuDSS dependencies:')
+    logger.exception(e)
     
 ############################################################
 #                       SOLVE REPORT                       #

{emerge-1.0.2 → emerge-1.0.3}/examples/demo10_sgh.py

@@ -33,7 +33,7 @@ dx = 2 * mm              # distance from horn exit to PML start
 
 # Create simulation object
 m = em.Simulation('HornAntenna')
-m.check_version("1.0.2") # Checks version compatibility.
+m.check_version("1.0.3") # Checks version compatibility.
 # --- Coordinate system for horn geometry -------------------------------
 hornCS = em.CS(em.YAX, em.ZAX, em.XAX)
 

{emerge-1.0.2 → emerge-1.0.3}/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("1.0.2") # Checks version compatibility.
+m.check_version("1.0.3") # Checks version compatibility.
 
 th = 0.5         # substrate thickness (meters)
 Hair = 2.0

{emerge-1.0.2 → emerge-1.0.3}/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("1.0.2") # Checks version compatibility.
+m.check_version("1.0.3") # Checks version compatibility.
 
 # We create a cyllindrical waveguide in the Y-axis.
 cyl = em.geo.Cylinder(0.012, 0.05, em.YAX.construct_cs())

{emerge-1.0.2 → emerge-1.0.3}/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("1.0.2") # Checks version compatibility.
+model.check_version("1.0.3") # Checks version compatibility.
 
 dfeed = 3*mm                    # straight feed length before the helix starts
 

{emerge-1.0.2 → emerge-1.0.3}/examples/demo14_boundary_selection.py

@@ -37,7 +37,7 @@ wgb = 10.16*mm
 L = 50*mm
 
 model = em.Simulation('Test Mode')
-model.check_version("1.0.2") # Checks version compatibility.
+model.check_version("1.0.3") # Checks version compatibility.
 
 # first lets define a WR90 waveguide
 wg_box = em.geo.Box(L, wga, wgb, position=(-L, -wga/2, -wgb/2))

{emerge-1.0.2 → emerge-1.0.3}/examples/demo1_stepped_imp_filter.py

@@ -28,7 +28,7 @@ pcbmat = em.Material(er=er, color="#217627", opacity=0.2)
 # We start by creating our simulation object.
 
 m = em.Simulation('Demo1_SIF')
-m.check_version("1.0.2") # Checks version compatibility.
+m.check_version("1.0.3") # 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-1.0.2 → emerge-1.0.3}/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("1.0.2") # Checks version compatibility.
+model.check_version("1.0.3") # Checks version compatibility.
 
 # The filter consists of quarter lamba cylindrical pins inside an airbox.
 # First we create the airbox

{emerge-1.0.2 → emerge-1.0.3}/examples/demo3_coupled_line_filter.py

@@ -38,7 +38,7 @@ extra = 100             # extra margin (mil)
 
 # --- Simulation setup ----------------------------------------------------
 model = em.Simulation('Demo3')
-model.check_version("1.0.2") # Checks version compatibility.
+model.check_version("1.0.3") # Checks version compatibility.
 # --- Material and layouter -----------------------------------------------
 mat = em.Material(er=3.55, color="#488343", opacity=0.4)
 

{emerge-1.0.2 → emerge-1.0.3}/examples/demo4_patch_antenna.py

@@ -34,7 +34,7 @@ f2 = 1.60e9             # stop frequency
 # --- Create simulation object -------------------------------------------
 model = em.Simulation('MyPatchAntenna')
 
-model.check_version("1.0.2") # Checks version compatibility.
+model.check_version("1.0.3") # Checks version compatibility.
 
 # --- Define geometry primitives -----------------------------------------
 # Substrate block centered at origin in XY, thickness in Z (negative down)
@@ -81,7 +81,7 @@ dielectric.material = em.Material(3.38, color="#207020", opacity=0.9)
 model.mw.resolution = 0.15
 
 # Frequency sweep across the resonance
-model.mw.set_frequency_range(f1, f2, 11)
+model.mw.set_frequency_range(f1, f2, 3)
 
 # --- Combine geometry into simulation -----------------------------------
 model.commit_geometry()
@@ -142,4 +142,5 @@ ff3d = data.field.find(freq=1.575e9).farfield_3d(boundary_selection)
 surf = ff3d.surfplot('normE', rmax=60 * mm,
                       offset=(0, 0, 20 * mm))
 model.display.add_surf(*surf, cmap='viridis', symmetrize=False)
+model.display.animate().add_surf(*data.field.find(freq=1.575e9).cutplane(2*mm, x=0).scalar('Ey','real'), scale='symlog', cmap='coolwarm')
 model.display.show()

{emerge-1.0.2 → emerge-1.0.3}/examples/demo5_revolve.py

@@ -29,7 +29,7 @@ airbox_width = 15*cm                # airbox width
 
 # --- Create simulation object -------------------------------------------
 model = em.Simulation('ConicalHornAntenna')
-model.check_version("1.0.2") # Checks version compatibility
+model.check_version("1.0.3") # Checks version compatibility
 
 # --- Feed geometry -------------------------------------------------------
 feed = em.geo.Cylinder(

{emerge-1.0.2 → emerge-1.0.3}/examples/demo6_striplines_with_vias.py

@@ -14,7 +14,7 @@ mm = 0.001    # Define a millimeter
 th = 1.0      # mm
 
 model = em.Simulation('Stripline_test')
-model.check_version("1.0.2") # Checks version compatibility.
+model.check_version("1.0.3") # Checks version compatibility.
 
 # As usual we start by creating our layouter
 ly = em.geo.PCB(th, mm, em.GCS, em.lib.DIEL_RO4350B, trace_material=em.lib.MET_COPPER)

{emerge-1.0.2 → emerge-1.0.3}/examples/demo7_periodic_cells.py

@@ -32,7 +32,7 @@ fl      = 25*mm  # Feed length
 
 # We start again by defining our simulation model
 model = em.Simulation('Periodic')
-model.check_version("1.0.2") # Checks version compatibility.
+model.check_version("1.0.3") # 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-1.0.2 → emerge-1.0.3}/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("1.0.2") # Checks version compatibility.
+    sim.check_version("1.0.3") # 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-1.0.2 → emerge-1.0.3}/examples/demo9_dielectric_resonator.py

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

{emerge-1.0.2 → emerge-1.0.3}/pyproject.toml

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

{emerge-1.0.2 → emerge-1.0.3}/uv.lock

@@ -245,7 +245,7 @@ wheels = [
 
 [[package]]
 name = "emerge"
-version = "1.0.2"
+version = "1.0.3"
 source = { editable = "." }
 dependencies = [
     { name = "cloudpickle" },