emerge 0.4.10__py3-none-any.whl → 0.5.0__py3-none-any.whl

emerge/_emerge/solver.py

@@ -32,32 +32,14 @@ from enum import Enum
 
 _PARDISO_AVAILABLE = False
 _UMFPACK_AVAILABLE = False
-_PARDISO_ERROR_CODES = """
-   0 | No error.
-  -1 | Input inconsistent.
-  -2 | Not enough memory.
-  -3 | Reordering problem.
-  -4 | Zero pivot, numerical fac. or iterative refinement problem.
-  -5 | Unclassified (internal) error.
-  -6 | Preordering failed (matrix types 11(real and nonsymmetric), 13(complex and nonsymmetric) only).
-  -7 | Diagonal Matrix problem.
-  -8 | 32-bit integer overflow problem.
- -10 | No license file pardiso.lic found.
- -11 | License is expired.
- -12 | Wrong username or hostname.
--100 | Reached maximum number of Krylov-subspace iteration in iterative solver.
--101 | No sufficient convergence in Krylov-subspace iteration within 25 iterations.
--102 | Error in Krylov-subspace iteration.
--103 | Bread-Down in Krylov-subspace iteration
-"""
+
 """ Check if the PC runs on a non-ARM architechture
 If so, attempt to import PyPardiso (if its installed)
 """
 
 if 'arm' not in platform.processor():
     try:
-        from pypardiso import spsolve as pardiso_solve
-        from pypardiso.pardiso_wrapper import PyPardisoError
+        from .solve_interfaces.pardiso_interface import PardisoInterface
         _PARDISO_AVAILABLE = True
     except ModuleNotFoundError as e:
         logger.info('Pardiso not found, defaulting to SuperLU')
@@ -67,6 +49,11 @@ try:
 except ModuleNotFoundError as e:
     logger.debug('UMFPACK not found, defaulting to SuperLU')
 
+
+############################################################
+#                 EIGENMODE FILTER ROUTINE                #
+############################################################
+
 def filter_real_modes(eigvals, eigvecs, k0, ermax, urmax, sign):
     """
     Given arrays of eigenvalues `eigvals` and eigenvectors `eigvecs` (cols of shape (N,)),
@@ -103,6 +90,11 @@ def filter_real_modes(eigvals, eigvecs, k0, ermax, urmax, sign):
     filtered_vecs = filtered_vecs[:, order] 
     return filtered_vals, filtered_vecs
 
+
+############################################################
+#               EIGENMODE ORTHOGONALITY CHECK              #
+############################################################
+
 def filter_unique_eigenpairs(eigen_values: list[complex], eigen_vectors: list[np.ndarray], tol=-3) -> tuple[list[complex], list[np.ndarray]]:
     """
     Filters eigenvectors by orthogonality using dot-product tolerance.
@@ -133,6 +125,11 @@ def filter_unique_eigenpairs(eigen_values: list[complex], eigen_vectors: list[np
 
     return unique_values, unique_vectors
 
+
+############################################################
+#         COMPLEX MATRIX TO REAL MATRIX CONVERSION        #
+############################################################
+
 def complex_to_real_block(A, b):
     """Return (Â,  b̂) real-augmented representation of A x = b."""
     A_r = sparse.csr_matrix(A.real)
@@ -153,7 +150,14 @@ def real_to_complex_block(x):
     x_i = x[n:]
     return x_r + 1j * x_i
 
-    
+
+############################################################
+#                  BASE CLASS DEFINITIONS                 #
+############################################################
+
+class SimulationError(Exception):
+    pass
+
 class Sorter:
     """ A Generic class that executes a sort on the indices.
     It must implement a sort and unsort method.
@@ -238,7 +242,11 @@ class EigSolver:
         pass
 
 
-## -----  SORTERS ----------------------------------------------
+
+############################################################
+#                          SORTERS                         #
+############################################################
+
 @dataclass
 class SolveReport:
     solver: str
@@ -270,7 +278,11 @@ class ReverseCuthillMckee(Sorter):
         logger.debug('Reversing Reverse Cuthill-Mckee sorting.')
         return  x[self.inv_perm]
     
-## -----  PRECONS ----------------------------------------------
+
+############################################################
+#                      PRECONDITIONERS                     #
+############################################################
+
 
 class ILUPrecon(Preconditioner):
     """ Implements the incomplete LU preconditioner on matrix A. """
@@ -285,7 +297,11 @@ class ILUPrecon(Preconditioner):
         self.ilu = sparse.linalg.spilu(A, drop_tol=1e-2, fill_factor=self.fill_factor, permc_spec='MMD_AT_PLUS_A', diag_pivot_thresh=0.001, options=self.options)
         self.M = sparse.linalg.LinearOperator(A.shape, self.ilu.solve)
 
-## ----- ITERATIVE SOLVERS -------------------------------------
+
+############################################################
+#                     ITERATIVE SOLVERS                    #
+############################################################
+
 
 class SolverBicgstab(Solver):
     """ Implements the Bi-Conjugate Gradient Stabilized method"""
@@ -359,7 +375,11 @@ class SolverGMRES(Solver):
             x, info = gmres(A, b, atol=self.atol, callback=self.callback, restart=500, callback_type='pr_norm')
         return x, info
 
-## -----  DIRECT SOLVERS ----------------------------------------
+
+############################################################
+#                      DIRECT SOLVERS                     #
+############################################################
+
 
 class SolverSuperLU(Solver):
     """ Implements Scipi's direct SuperLU solver."""
@@ -392,16 +412,15 @@ class SolverUMFPACK(Solver):
         super().__init__()
         self.A: np.ndarray = None
         self.b: np.ndarray = None
-        self.up: um.UmfpackContext = um.UmfpackContext('zl')
-        self.up.control[um.UMFPACK_PRL] = 0  #less terminal printing
-        self.up.control[um.UMFPACK_IRSTEP] = 2
-        self.up.control[um.UMFPACK_STRATEGY] = um.UMFPACK_STRATEGY_SYMMETRIC
-        self.up.control[um.UMFPACK_ORDERING] = 3
-        self.up.control[um.UMFPACK_PIVOT_TOLERANCE] = 0.001
-        self.up.control[um.UMFPACK_SYM_PIVOT_TOLERANCE] = 0.001
-        self.up.control[um.UMFPACK_BLOCK_SIZE] = 64
-        self.up.control[um.UMFPACK_FIXQ] = -1
-        #self.up.control[um.UMFPACK_]
+        self.umfpack: um.UmfpackContext = um.UmfpackContext('zl')
+        self.umfpack.control[um.UMFPACK_PRL] = 0
+        self.umfpack.control[um.UMFPACK_IRSTEP] = 2
+        self.umfpack.control[um.UMFPACK_STRATEGY] = um.UMFPACK_STRATEGY_SYMMETRIC
+        self.umfpack.control[um.UMFPACK_ORDERING] = 3
+        self.umfpack.control[um.UMFPACK_PIVOT_TOLERANCE] = 0.001
+        self.umfpack.control[um.UMFPACK_SYM_PIVOT_TOLERANCE] = 0.001
+        self.umfpack.control[um.UMFPACK_BLOCK_SIZE] = 64
+        self.umfpack.control[um.UMFPACK_FIXQ] = -1
 
         self.fact_symb: bool = False
 
@@ -414,39 +433,49 @@ class SolverUMFPACK(Solver):
         A.indices = A.indices.astype(np.int64)
         if self.fact_symb is False:
             logger.debug('Executing symbollic factorization.')
-            self.up.symbolic(A)
+            self.umfpack.symbolic(A)
             #self.up.report_symbolic()
             self.fact_symb = True
         if not reuse_factorization:
             #logger.debug('Executing numeric factorization.')
-            self.up.numeric(A)
+            self.umfpack.numeric(A)
             self.A = A
-        x = self.up.solve(um.UMFPACK_A, self.A, b, autoTranspose = False )
+        x = self.umfpack.solve(um.UMFPACK_A, self.A, b, autoTranspose = False )
         return x, 0
 
 class SolverPardiso(Solver):
     """ Implements the PARDISO solver through PyPardiso. """
-    real_only: bool = True
+    real_only: bool = False
     req_sorter: bool = False
 
     def __init__(self):
         super().__init__()
-
+        self.solver: PardisoInterface = PardisoInterface()
+        self.fact_symb: bool = False
         self.A: np.ndarray = None
         self.b: np.ndarray = None
     
     def solve(self, A, b, precon, reuse_factorization: bool = False, id: int = -1):
         logger.info(f'Calling Pardiso Solver. ID={id}')
-        self.A = A
-        self.b = b
-        try:
-            x = pardiso_solve(A, b)
-        except PyPardisoError as e:
-            print('Error Codes:')
-            print(_PARDISO_ERROR_CODES)
-        return x, 0
+        if self.fact_symb is False:
+            logger.debug('Executing symbollic factorization.')
+            self.solver.symbolic(A)
+            self.fact_symb = True
+        if not reuse_factorization:
+            self.solver.numeric(A)
+            self.A = A
+        x, error = self.solver.solve(A, b)
+        if error != 0:
+            logger.error(f'Terminated with error code {error}')
+            logger.error(self.solver.get_error(error))
+            raise SimulationError(f'PARDISO Terminated with error code {error}')
+        return x, error
     
-## -----  DIRECT EIG SOLVERS --------------------------------------
+
+############################################################
+#                 DIRECT EIGENMODE SOLVERS                #
+############################################################
+
 class SolverLAPACK(EigSolver):
 
     def __init__(self):
@@ -480,7 +509,11 @@ class SolverLAPACK(EigSolver):
         lam, vecs = filter_real_modes(lam, vecs, target_k0, 2, 2, sign=sign)
         return lam, vecs
     
-## -----  ITER EIG SOLVERS ---------------------------------------
+
+############################################################
+#                  ITERATIVE EIGEN SOLVERS                 #
+############################################################
+
 
 class SolverARPACK(EigSolver):
     """ Implements the Scipy ARPACK iterative eigenmode solver."""
@@ -600,7 +633,11 @@ class SmartARPACK(EigSolver):
         return eigen_values, eigen_modes
 
 
-## ----- SOLVE ENUMS ---------------------------------------------
+
+############################################################
+#                        SOLVER ENUM                       #
+############################################################
+
 
 class EMSolver(Enum):
     SUPERLU = 1
@@ -633,7 +670,11 @@ class EMSolver(Enum):
         elif self==EMSolver.SMART_ARPACK_BMA:
             return SmartARPACK_BMA()
     
-## -----  SOLVE ROUTINES -----------------------------------------
+
+############################################################
+#                       SOLVE ROUTINE                      #
+############################################################
+
 
 class SolveRoutine:
     """ A generic class describing a solve routine.
@@ -874,11 +915,12 @@ class SolveRoutine:
         NS = solve_ids.shape[0]
 
         A = A.tocsc()
-
-        logger.debug(f'    Removing {NF-NS} prescribed DOFs ({NS} left)')
-
+        
         Asel = A[np.ix_(solve_ids, solve_ids)]
         bsel = b[solve_ids]
+        nnz = Asel.nnz
+
+        logger.debug(f'    Removed {NF-NS} prescribed DOFs ({NS:,} left, {nnz:,} non-zero)')
 
         if solver.real_only:
             logger.debug('    Converting to real matrix')
@@ -900,6 +942,7 @@ class SolveRoutine:
                 precon = str(self.precon)
 
         start = time.time()
+        
         x_solved, code = solver.solve(Asorted, bsorted, self.precon, reuse_factorization=reuse, id=id)
         end = time.time()
         simtime = end-start
@@ -1053,4 +1096,10 @@ class AutomaticRoutine(SolveRoutine):
             return self._try_solver(EMSolver.SUPERLU)
         return self._try_solver(EMSolver.SUPERLU)
     
+
+
+############################################################
+#                    DEFAULT DEFINITION                   #
+############################################################
+
 DEFAULT_ROUTINE = AutomaticRoutine()

emerge/lib.py

@@ -10,48 +10,283 @@ MU0 = 1.2566370612720e-6
 
 ### MATERIALS 
 VACUUM = Material(color="#2d8cd5", opacity=0.05)
-FR4 = Material(er=4.4, tand=0.001, color="#3c9747")
-# Conductors
-ALUMINUM  = Material(cond=3.5e7,         color="#aaaaaa")
-GOLD      = Material(cond=4.1e7,         color="#A07130")
-SILVER    = Material(cond=6.3e7,         color="#bfbfbf")
-TIN       = Material(cond=9.17e6,        color="#d6d6d6")
-NICKEL    = Material(cond=1.43e7,        color="#78706c")
-IRON      = Material(cond=1.0e7,  ur=5000, color="#666666")
-STEEL     = Material(cond=1.45e6, ur=100,  color="#808080")
-
-# Semiconductors & Dielectrics
-SILICON   = Material(er=11.68,    cond=0.1, color="#333333")
-SIO2      = Material(er=3.9,                color="#e6e6e6")
-GAAS      = Material(er=12.9,    cond=0.0, color="#404071")
-
-# Dielectrics
-PTFE      = Material(er=2.1,     tand=0.0002,       color="#cccccc")
-POLYIMIDE = Material(er=3.4,     tand=0.02,         color="#b8b8b8")
-CERAMIC   = Material(er=6.0,     tand=0.001,        color="#efead1")
-
-# Liquids
-WATER     = Material(er=80.1,    cond=0.0,  color="#0080ff", opacity=0.3)
-FERRITE   = Material(er=12.0,    ur=2000, tand=0.02, color="#994d4d")
-
-# Specialty RF Substrates
-ROGERS_4350B = Material(er=3.66, tand=0.0037,       color="#3c9747")
-ROGERS_5880  = Material(er=2.2,  tand=0.0009,       color="#3c9747")
-ROGERS_RO3003   = Material(er=3.0,  tand=0.0013,      color="#3c9747")
-ROGERS_RO3010   = Material(er=10.2, tand=0.0023,      color="#3c9747")
-ROGERS_RO4003C  = Material(er=3.55, tand=0.0027,      color="#3c9747")
-ROGERS_DUROID6002 = Material(er=2.94, tand=0.0012,    color="#3c9747")
-ROGERS_DUROID6010 = Material(er=6.15, tand=0.0027,    color="#3c9747")
-ROGERS_RT5880  = Material(er=2.2,  tand=0.0009,       color="#3c9747")
-TACONIC_RF35   = Material(er=3.5,  tand=0.0018,       color="#3c9747")
-TACONIC_TLC30  = Material(er=3.0,  tand=0.0020,       color="#3c9747")
-ISOLA_I_TERA_MT = Material(er=3.45, tand=0.0030,      color="#3c9747")
-ISOLA_NELCO_4000_13 = Material(er=3.77, tand=0.008,  color="#3c9747")
-VENTEC_VERDELAY_400HR = Material(er=4.0,  tand=0.02,  color="#3c9747")
+
+############################################################
+#                          METALS                         #
+############################################################
+
+GREY = "#bfbfbf"
+MET_ALUMINUM         = Material(cond=3.77e7,  color=GREY)
+MET_CARBON           = Material(cond=3.33e4,  color=GREY)
+MET_CHROMIUM         = Material(cond=5.56e6,  color=GREY)
+MET_COPPER           = Material(cond=5.8e7,   color="#62290c")
+MET_GOLD             = Material(cond=4.10e7,  color=GREY)
+MET_INDIUM           = Material(cond=6.44e6,  color=GREY)
+MET_IRIDIUM          = Material(cond=2.13e7,  color=GREY)
+MET_IRON             = Material(cond=1.04e7,  color=GREY)
+MET_LEAD             = Material(cond=4.84e6,  color=GREY)
+MET_MAGNESIUM        = Material(cond=2.38e7,  color=GREY)
+MET_NICKEL           = Material(cond=1.14e7,  color=GREY)
+MET_NICHROME         = Material(cond=9.09e5,  color=GREY)
+MET_PALLADIUM        = Material(cond=9.42e6,  color=GREY)
+MET_PLATINUM         = Material(cond=9.42e6,  color=GREY)
+MET_RHODIUM          = Material(cond=2.22e7,  color=GREY)
+MET_SILVER           = Material(cond=6.29e7,  color=GREY)
+MET_TANTALUM         = Material(cond=6.44e6,  color=GREY)
+MET_TANTALUM_NITRIDE = Material(cond=3.97e5,  color=GREY)
+MET_TIN              = Material(cond=8.66e6,  color=GREY)
+MET_TITANIUM         = Material(cond=1.82e6,  color=GREY)
+MET_TUNGSTEN         = Material(cond=1.79e7,  color=GREY)
+MET_ZINC             = Material(cond=1.76e7,  color=GREY)
+MET_ZIRCONIUM        = Material(cond=2.44e7,  color=GREY)
+
+
+############################################################
+#                      SEMICONDUCTORS                     #
+############################################################
+
+SEMI_SILICON       = Material(er=11.7,  tand=0.005,     color="#b4b4b4")   # Crystalline Si
+SEMI_SILICON_N     = Material(er=7.5,   tand=0.0003,    color="#a0a0a0")   # Silicon Nitride (Si₃N₄)
+SEMI_SILICON_OXIDE = Material(er=3.9,   tand=0.0001,    color="#e0e0e0")   # Silicon Dioxide (SiO₂)
+SEMI_GERMANIUM     = Material(er=16.0,  tand=0.001,     color="#787878")
+SEMI_GAAS          = Material(er=13.1,  tand=0.0016,    color="#aa8888")   # Gallium Arsenide
+SEMI_GA_N          = Material(er=8.9,   tand=0.002,     color="#8888cc")   # Gallium Nitride
+SEMI_INP           = Material(er=12.5,  tand=0.0015,    color="#cc99aa")   # Indium Phosphide
+SEMI_ALN           = Material(er=8.6,   tand=0.0003,    color="#ccccee")   # Aluminum Nitride
+SEMI_AL2O3         = Material(er=9.8,   tand=0.0002,    color="#eaeaea")   # Alumina
+SEMI_SAPPHIRE      = Material(er=9.4,   tand=0.0001,    color="#ddddff")
+SEMI_DIAMOND       = Material(er=5.5,   tand=0.00005,   color="#cceeff")   # Synthetic CVD diamond
+SEMI_HBN           = Material(er=4.0,   tand=0.0001,    color="#eeeeff")   # Hexagonal Boron Nitride
+SEMI_SIOXNY        = Material(er=5.0,   tand=0.002,     color="#ddddee")   # Silicon Oxynitride (SiOxNy)
+
+############################################################
+#                          LIQUIDS                         #
+############################################################
+
+LIQ_WATER     = Material(er=80.1,    cond=0.0,  color="#0080ff", opacity=0.3)
+LIQ_FERRITE   = Material(er=12.0,    ur=2000, tand=0.02, color="#994d4d")
+
+
+############################################################
+#                        DIELECTRICS                       #
+############################################################
+DIEL_PTFE      = Material(er=2.1,     tand=0.0002,       color="#21912b")
+DIEL_POLYIMIDE = Material(er=3.4,     tand=0.02,         color="#b8b8b8")
+DIEL_CERAMIC   = Material(er=6.0,     tand=0.001,        color="#efead1")
+DIEL_AD10 = Material(er=10.2, tand=0.0078, color="#21912b")
+DIEL_AD1000 = Material(er=10.2, tand=0.0023, color="#21912b")
+DIEL_AD250 = Material(er=2.5, tand=0.0018, color="#21912b")
+DIEL_AD250_PIM = Material(er=2.5, tand=0.0018, color="#21912b")
+DIEL_AD250A = Material(er=2.50, tand=0.0015, color="#21912b")
+DIEL_AD250C = Material(er=2.50, tand=0.0014, color="#21912b")
+DIEL_AD255 = Material(er=2.55, tand=0.0018, color="#21912b")
+DIEL_AD255A = Material(er=2.55, tand=0.0015, color="#21912b")
+DIEL_AD255C = Material(er=2.55, tand=0.0014, color="#21912b")
+DIEL_AD260A = Material(er=2.60, tand=0.0017, color="#21912b")
+DIEL_AD270 = Material(er=2.7, tand=0.0023, color="#21912b")
+DIEL_AD300 = Material(er=3, tand=0.003, color="#21912b")
+DIEL_AD300_PIM = Material(er=3, tand=0.003, color="#21912b")
+DIEL_AD300A = Material(er=3.00, tand=0.002, color="#21912b")
+DIEL_AD300C = Material(er=2.97, tand=0.002, color="#21912b")
+DIEL_AD320 = Material(er=3.2, tand=0.0038, color="#21912b")
+DIEL_AD320_PIM = Material(er=3.2, tand=0.003, color="#21912b")
+DIEL_AD320A = Material(er=3.20, tand=0.0032, color="#21912b")
+DIEL_AD350 = Material(er=3.5, tand=0.003, color="#21912b")
+DIEL_AD350_PIM = Material(er=3.5, tand=0.003, color="#21912b")
+DIEL_AD350A = Material(er=3.50, tand=0.003, color="#21912b")
+DIEL_AD410 = Material(er=4.1, tand=0.003, color="#21912b")
+DIEL_AD430 = Material(er=4.3, tand=0.003, color="#21912b")
+DIEL_AD450 = Material(er=4.5, tand=0.0035, color="#21912b")
+DIEL_AD450A = Material(er=4.5, tand=0.0035, color="#21912b")
+DIEL_AD5 = Material(er=5.1, tand=0.003, color="#21912b")
+DIEL_AD600 = Material(er=5.90, tand=0.003, color="#21912b")
+DIEL_AR1000 = Material(er=9.8, tand=0.003, color="#21912b")
+DIEL_CER_10 = Material(er=10.00, tand=0.0035, color="#21912b")
+DIEL_CLTE = Material(er=2.96, tand=0.0023, color="#21912b")
+DIEL_CLTE_AT = Material(er=3.00, tand=0.0013, color="#21912b")
+DIEL_CLTE_LC = Material(er=2.94, tand=0.0025, color="#21912b")
+DIEL_CLTE_XT = Material(er=2.94, tand=0.0012, color="#21912b")
+DIEL_COMCLAD_HF_ER2 = Material(er=2, tand=0.0025, color="#21912b")
+DIEL_COMCLAD_HF_ER3 = Material(er=3, tand=0.0025, color="#21912b")
+DIEL_COMCLAD_HF_ER4 = Material(er=4, tand=0.0025, color="#21912b")
+DIEL_COMCLAD_HF_ER5 = Material(er=5, tand=0.0025, color="#21912b")
+DIEL_COMCLAD_HF_ER6 = Material(er=6, tand=0.0025, color="#21912b")
+DIEL_COPPER_CLAD_ULTEM = Material(er=3.05, tand=0.003, color="#21912b")
+DIEL_CUCLAD_217LX = Material(er=2.17, tand=0.0009, color="#21912b")
+DIEL_CUCLAD_233LX = Material(er=2.33, tand=0.0013, color="#21912b")
+DIEL_CUCLAD_250GT = Material(er=2.5, tand=0.0018, color="#21912b")
+DIEL_CUCLAD_250GX = Material(er=2.4, tand=0.0018, color="#21912b")
+DIEL_CUFLON = Material(er=2.05, tand=0.00045, color="#21912b")
+DIEL_DICLAD_522 = Material(er=2.4, tand=0.0018, color="#21912b")
+DIEL_DICLAD_527 = Material(er=2.4, tand=0.0018, color="#21912b")
+DIEL_DICLAD_870 = Material(er=2.33, tand=0.0013, color="#21912b")
+DIEL_DICLAD_880 = Material(er=2.17, tand=0.0009, color="#21912b")
+DIEL_DICLAD_880_PIM = Material(er=2.17, tand=0.0009, color="#21912b")
+DIEL_GETEK = Material(er=3.5, tand=0.01, color="#21912b")
+DIEL_GETEK = Material(er=3.8, tand=0.01, color="#21912b")
+DIEL_IS6802_80 = Material(er=2.80, tand=0.003, color="#21912b")
+DIEL_IS6803_00 = Material(er=3.00, tand=0.003, color="#21912b")
+DIEL_IS6803_20 = Material(er=3.20, tand=0.003, color="#21912b")
+DIEL_IS6803_33 = Material(er=3.33, tand=0.003, color="#21912b")
+DIEL_IS6803_38 = Material(er=3.38, tand=0.0032, color="#21912b")
+DIEL_IS6803_45 = Material(er=3.45, tand=0.0035, color="#21912b")
+DIEL_ISOCLAD_917 = Material(er=2.17, tand=0.0013, color="#21912b")
+DIEL_ISOCLAD_933 = Material(er=2.33, tand=0.0016, color="#21912b")
+DIEL_ISOLA_I_TERA_MT = Material(er=3.45, tand=0.0030,      color="#3c9747")
+DIEL_ISOLA_NELCO_4000_13 = Material(er=3.77, tand=0.008,  color="#3c9747")
+DIEL_MAT_25N = Material(er=3.38, tand=0.0025, color="#21912b")
+DIEL_MAT25FR = Material(er=3.58, tand=0.0035, color="#21912b")
+DIEL_MEGTRON6R5775 = Material(er=3.61, tand=0.004, color="#21912b")
+DIEL_MERCURYWAVE_9350 = Material(er=3.5, tand=0.004, color="#21912b")
+DIEL_MULTICLAD_HF = Material(er=3.7, tand=0.0045, color="#21912b")
+DIEL_N_8000 = Material(er=3.5, tand=0.011, color="#21912b")
+DIEL_N4350_13RF = Material(er=3.5, tand=0.0065, color="#21912b")
+DIEL_N4380_13RF = Material(er=3.8, tand=0.007, color="#21912b")
+DIEL_N8000Q = Material(er=3.2, tand=0.006, color="#21912b")
+DIEL_N9300_13RF = Material(er=3, tand=0.004, color="#21912b")
+DIEL_N9320_13RF = Material(er=3.2, tand=0.0045, color="#21912b")
+DIEL_N9338_13RF = Material(er=3.38, tand=0.0046, color="#21912b")
+DIEL_N9350_13RF = Material(er=3.48, tand=0.0055, color="#21912b")
+DIEL_NH9294 = Material(er=2.94, tand=0.0022, color="#21912b")
+DIEL_NH9300 = Material(er=3.00, tand=0.0023, color="#21912b")
+DIEL_NH9320 = Material(er=3.20, tand=0.0024, color="#21912b")
+DIEL_NH9338 = Material(er=3.38, tand=0.0025, color="#21912b")
+DIEL_NH9348 = Material(er=3.48, tand=0.003, color="#21912b")
+DIEL_NH9350 = Material(er=3.50, tand=0.003, color="#21912b")
+DIEL_NH9410 = Material(er=4.10, tand=0.003, color="#21912b")
+DIEL_NH9450 = Material(er=4.50, tand=0.003, color="#21912b")
+DIEL_NORCLAD = Material(er=2.55, tand=0.0011, color="#21912b")
+DIEL_NX9240 = Material(er=2.40, tand=0.0016, color="#21912b")
+DIEL_NX9245 = Material(er=2.45, tand=0.0016, color="#21912b")
+DIEL_NX9250 = Material(er=2.50, tand=0.0017, color="#21912b")
+DIEL_NX9255 = Material(er=2.55, tand=0.0018, color="#21912b")
+DIEL_NX9260 = Material(er=2.60, tand=0.0019, color="#21912b")
+DIEL_NX9270 = Material(er=2.70, tand=0.002, color="#21912b")
+DIEL_NX9294 = Material(er=2.94, tand=0.0022, color="#21912b")
+DIEL_NX9300 = Material(er=3.00, tand=0.0023, color="#21912b")
+DIEL_NX9320 = Material(er=3.20, tand=0.0024, color="#21912b")
+DIEL_NY9208 = Material(er=2.08, tand=0.0006, color="#21912b")
+DIEL_NY9217 = Material(er=2.17, tand=0.0008, color="#21912b")
+DIEL_NY9220 = Material(er=2.20, tand=0.0009, color="#21912b")
+DIEL_NY9233 = Material(er=2.33, tand=0.0011, color="#21912b")
+DIEL_POLYGUIDE = Material(er=2.320, tand=0.0005, color="#21912b")
+DIEL_RF_30 = Material(er=3.00, tand=0.0019, color="#21912b")
+DIEL_RF_301 = Material(er=2.97, tand=0.0018, color="#21912b")
+DIEL_RF_35 = Material(er=3.50, tand=0.0025, color="#21912b")
+DIEL_RF_35A2 = Material(er=3.50, tand=0.0015, color="#21912b")
+DIEL_RF_35P = Material(er=3.50, tand=0.0034, color="#21912b")
+DIEL_RF_35TC = Material(er=3.50, tand=0.0011, color="#21912b")
+DIEL_RF_41 = Material(er=4.10, tand=0.0038, color="#21912b")
+DIEL_RF_43 = Material(er=4.30, tand=0.0033, color="#21912b")
+DIEL_RF_45 = Material(er=4.50, tand=0.0037, color="#21912b")
+DIEL_RF_60A = Material(er=6.15, tand=0.0038, color="#21912b")
+DIEL_RO3003 = Material(er=3.00, tand=0.0011, color="#21912b")
+DIEL_RO3006 = Material(er=6.15, tand=0.002, color="#21912b")
+DIEL_RO3010 = Material(er=10.2, tand=0.0022, color="#21912b")
+DIEL_RO3035 = Material(er=3.50, tand=0.0017, color="#21912b")
+DIEL_RO3203 = Material(er=3.02, tand=0.0016, color="#21912b")
+DIEL_RO3206 = Material(er=6.15, tand=0.0027, color="#21912b")
+DIEL_RO3210 = Material(er=10.2, tand=0.0027, color="#21912b")
+DIEL_RO3730 = Material(er=3.00, tand=0.0016, color="#21912b")
+DIEL_RO4003C = Material(er=3.38, tand=0.0029, color="#21912b")
+DIEL_RO4350B = Material(er=3.48, tand=0.0037, color="#21912b")
+DIEL_RO4350B_TX = Material(er=3.48, tand=0.0034, color="#21912b")
+DIEL_RO4360 = Material(er=6.15, tand=0.0038, color="#21912b")
+DIEL_RO4533 = Material(er=3.30, tand=0.0025, color="#21912b")
+DIEL_RO4534 = Material(er=3.40, tand=0.0027, color="#21912b")
+DIEL_RO4535 = Material(er=3.50, tand=0.0037, color="#21912b")
+DIEL_RO4730 = Material(er=3.00, tand=0.0033, color="#21912b")
+DIEL_RT_Duroid_5870 = Material(er=2.33, tand=0.0012, color="#21912b")
+DIEL_RT_Duroid_5880 = Material(er=2.20, tand=0.0009, color="#21912b")
+DIEL_RT_Duroid_5880LZ = Material(er=1.96, tand=0.0019, color="#21912b")
+DIEL_RT_Duroid_6002 = Material(er=2.94, tand=0.0012, color="#21912b")
+DIEL_RT_Duroid_6006 = Material(er=6.15, tand=0.0027, color="#21912b")
+DIEL_RT_Duroid_6010_2LM = Material(er=10.2, tand=0.0023, color="#21912b")
+DIEL_RT_Duroid_6035HTC = Material(er=3.50, tand=0.0013, color="#21912b")
+DIEL_RT_Duroid_6202 = Material(er=2.94, tand=0.0015, color="#21912b")
+DIEL_RT_Duroid_6202PR = Material(er=2.90, tand=0.002, color="#21912b")
+DIEL_SYRON_70000_002IN_Thick = Material(er=3.4, tand=0.0045, color="#21912b")
+DIEL_SYRON_71000_004IN_THICK = Material(er=3.39, tand=0.005, color="#21912b")
+DIEL_SYRON_71000INCH = Material(er=3.61, tand=0.006, color="#21912b")
+DIEL_TACLAMPLUS= Material(er=2.10, tand=0.0004, color="#21912b")
+DIEL_TC350 = Material(er=3.50, tand=0.002, color="#21912b")
+DIEL_TC600 = Material(er=6.15, tand=0.002, color="#21912b")
+DIEL_THETA = Material(er=3.85, tand=0.0123, color="#21912b")
+DIEL_TLA_6 = Material(er=2.62, tand=0.0017, color="#21912b")
+DIEL_TLC_27 = Material(er=2.75, tand=0.003, color="#21912b")
+DIEL_TLC_30 = Material(er=3.00, tand=0.003, color="#21912b")
+DIEL_TLC_32 = Material(er=3.20, tand=0.003, color="#21912b")
+DIEL_TLC_338 = Material(er=3.38, tand=0.0034, color="#21912b")
+DIEL_TLC_35 = Material(er=3.50, tand=0.0037, color="#21912b")
+DIEL_TLE_95 = Material(er=2.95, tand=0.0028, color="#21912b")
+DIEL_TLF_34 = Material(er=3.40, tand=0.002, color="#21912b")
+DIEL_TLF_35 = Material(er=3.50, tand=0.002, color="#21912b")
+DIEL_TLG_29 = Material(er=2.87, tand=0.0027, color="#21912b")
+DIEL_TLG_30 = Material(er=3, tand=0.0038, color="#21912b")
+DIEL_TLP_3 = Material(er=2.33, tand=0.0009, color="#21912b")
+DIEL_TLP_5 = Material(er=2.20, tand=0.0009, color="#21912b")
+DIEL_TLP_5A = Material(er=2.17, tand=0.0009, color="#21912b")
+DIEL_TLT_0 = Material(er=2.45, tand=0.0019, color="#21912b")
+DIEL_TLT_6 = Material(er=2.65, tand=0.0019, color="#21912b")
+DIEL_TLT_7 = Material(er=2.60, tand=0.0019, color="#21912b")
+DIEL_TLT_8 = Material(er=2.55, tand=0.0019, color="#21912b")
+DIEL_TLT_9 = Material(er=2.50, tand=0.0019, color="#21912b")
+DIEL_TLX_0 = Material(er=2.45, tand=0.0019, color="#21912b")
+DIEL_TLX_6 = Material(er=2.65, tand=0.0019, color="#21912b")
+DIEL_TLX_7 = Material(er=2.60, tand=0.0019, color="#21912b")
+DIEL_TLX_8 = Material(er=2.55, tand=0.0019, color="#21912b")
+DIEL_TLX_9 = Material(er=2.50, tand=0.0019, color="#21912b")
+DIEL_TLY_3 = Material(er=2.33, tand=0.0012, color="#21912b")
+DIEL_TLY_3F = Material(er=2.33, tand=0.0012, color="#21912b")
+DIEL_TLY_5 = Material(er=2.20, tand=0.0009, color="#21912b")
+DIEL_TLY_5_L = Material(er=2.20, tand=0.0009, color="#21912b")
+DIEL_TLY_5A = Material(er=2.17, tand=0.0009, color="#21912b")
+DIEL_TLY_5AL = Material(er=2.17, tand=0.0009, color="#21912b")
+DIEL_TMM_10 = Material(er=9.20, tand=0.0022, color="#21912b")
+DIEL_TMM_10i = Material(er=9.80, tand=0.002, color="#21912b")
+DIEL_TMM_3 = Material(er=3.27, tand=0.002, color="#21912b")
+DIEL_TMM_4 = Material(er=4.50, tand=0.002, color="#21912b")
+DIEL_TMM_6 = Material(er=6.00, tand=0.0023, color="#21912b")
+DIEL_TRF_41 = Material(er=4.10, tand=0.0035, color="#21912b")
+DIEL_TRF_43 = Material(er=4.30, tand=0.0035, color="#21912b")
+DIEL_TRF_45 = Material(er=4.50, tand=0.0035, color="#21912b")
+DIEL_TSM_26 = Material(er=2.60, tand=0.0014, color="#21912b")
+DIEL_TSM_29 = Material(er=2.94, tand=0.0013, color="#21912b")
+DIEL_TSM_30 = Material(er=3.00, tand=0.0013, color="#21912b")
+DIEL_TSM_DS = Material(er=2.85, tand=0.001, color="#21912b")
+DIEL_TSM_DS3 = Material(er=3.00, tand=0.0011, color="#21912b")
+DIEL_ULTRALAM_2000 = Material(er=2.4, tand=0.0019, color="#21912b")
+DIEL_ULTRALAM_3850 = Material(er=2.9, tand=0.0025, color="#21912b")
+DIEL_XT_Duroid_80000_002IN_Thick = Material(er=3.23, tand=0.0035, color="#21912b")
+DIEL_XT_Duroid_8100 = Material(er=3.54, tand=0.0049, color="#21912b")
+DIEL_XT_Duroid_81000_004IN_Thick = Material(er=3.32, tand=0.0038, color="#21912b")
 
 # Legacy FR Materials
-FR1            = Material(er=4.8,  tand=0.025,        color="#3c9747")
-FR2            = Material(er=4.8,  tand=0.02,        color="#3c9747")
+DIEL_FR1 = Material(er=4.8,  tand=0.025, color="#3c9747")  # Paper + phenolic resin
+DIEL_FR2 = Material(er=4.8,  tand=0.02,  color="#3c9747")  # Paper + phenolic resin
+DIEL_FR3 = Material(er=4.5,  tand=0.02,  color="#2b7a4b")  # Paper + epoxy resin
+DIEL_FR4 = Material(er=4.4,  tand=0.015, color="#1e8449")  # Woven glass + epoxy resin (industry standard)
+DIEL_FR5 = Material(er=4.2,  tand=0.012, color="#156e38")  # Woven glass + high-temp epoxy resin
+DIEL_FR6 = Material(er=5.2,  tand=0.030, color="#145a32")  # Paper + unknown resin, poor thermal performance
 
 # Magnetic Materials
-MU_METAL  = Material(cond=1.0e6, ur=200000, color="#666680")
+MU_METAL  = Material(cond=1.0e6, ur=200000, color="#666680")
+
+
+############################################################
+#                           FOAMS                          #
+############################################################
+
+FOAM_ROHACELL_31    = Material(er=1.05,  tand=0.0005,  color="#f0e1a1")  # PMI-based structural foam
+FOAM_ROHACELL_51    = Material(er=1.07,  tand=0.0006,  color="#f0dea0")  # denser version
+FOAM_ROHACELL_71    = Material(er=1.10,  tand=0.0007,  color="#e5d199")
+FOAM_PEI       = Material(er=1.15,  tand=0.0035,  color="#e0b56f")  # polyetherimide-based foam
+FOAM_PMI       = Material(er=1.10,  tand=0.0008,  color="#d9c690")  # polymethacrylimide
+FOAM_PVC       = Material(er=1.20,  tand=0.0040,  color="#cccccc")
+FOAM_EPS            = Material(er=1.03,  tand=0.0050,  color="#f7f7f7")  # expanded polystyrene
+FOAM_XPS            = Material(er=1.05,  tand=0.0030,  color="#e0e0e0")  # extruded polystyrene
+FOAM_PU        = Material(er=1.10,  tand=0.0080,  color="#d0d0d0")  # polyurethane foam
+FOAM_GLAS       = Material(er=3.10,  tand=0.0050,  color="#888888")  # cellular glass, denser
+FOAM_AIREX_C70      = Material(er=1.10,  tand=0.0010,  color="#f7e7a3")  # PET closed cell
+FOAM_AIREX_T92      = Material(er=1.10,  tand=0.0020,  color="#f6d08a")  # higher strength PET
+FOAM_PVC_CORECELL   = Material(er=1.56,  tand=0.0025,  color="#aaaaaa")  # structural core PVC

{emerge-0.4.10.dist-info → emerge-0.5.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: emerge
-Version: 0.4.10
+Version: 0.5.0
 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
@@ -10,13 +10,12 @@ Requires-Dist: gmsh>=4.13.1
 Requires-Dist: joblib>=1.5.1
 Requires-Dist: loguru>=0.7.3
 Requires-Dist: matplotlib>=3.8.0
+Requires-Dist: mkl!=2024.0; platform_machine == 'x86_64' or platform_machine == 'AMD64'
 Requires-Dist: numba-progress>=1.1.3
 Requires-Dist: numba>=0.57.0
 Requires-Dist: numpy<2.3,>=1.24
 Requires-Dist: pyvista>=0.45.2
 Requires-Dist: scipy>=1.14.0
-Provides-Extra: pypardiso
-Requires-Dist: pypardiso; (platform_machine == 'x86_64' or platform_machine == 'AMD64') and extra == 'pypardiso'
 Provides-Extra: umfpack
 Requires-Dist: scikit-umfpack; (sys_platform != 'win32') and extra == 'umfpack'
 Description-Content-Type: text/markdown