ElectroMagneticPython 2.2.3__tar.gz → 2.2.4__tar.gz

{electromagneticpython-2.2.3 → electromagneticpython-2.2.4}/.bumpversion.cfg

@@ -3,4 +3,4 @@ commit = True
 sign_tags = True
 tag = True
 tag_name = {new_version}
-current_version = 2.2.3
+current_version = 2.2.4

electromagneticpython-2.2.4/EMpy/_version.py

@@ -0,0 +1,24 @@
+# file generated by vcs-versioning
+# don't change, don't track in version control
+from __future__ import annotations
+
+__all__ = [
+    "__version__",
+    "__version_tuple__",
+    "version",
+    "version_tuple",
+    "__commit_id__",
+    "commit_id",
+]
+
+version: str
+__version__: str
+__version_tuple__: tuple[int | str, ...]
+version_tuple: tuple[int | str, ...]
+commit_id: str | None
+__commit_id__: str | None
+
+__version__ = version = '2.2.4'
+__version_tuple__ = version_tuple = (2, 2, 4)
+
+__commit_id__ = commit_id = 'gd8c276251'

{electromagneticpython-2.2.3 → electromagneticpython-2.2.4}/EMpy/constants.py

@@ -1,5 +1,3 @@
-# pylint: disable=C0301
-
 """Useful constants.
 
 Constants used in mathematics and electromagnetism.

{electromagneticpython-2.2.3 → electromagneticpython-2.2.4}/EMpy/devices.py

@@ -1,5 +1,3 @@
-# pylint: disable=C0302,W0622,W1001,R0902,R0903,W0201,R0913,R0914,C0325,W0221
-
 """Set of objects and functions to describe the behaviour of some
 widely known devices and build new ones.
 
@@ -62,7 +60,7 @@ class DeviceMatrix:
         return self.__data
 
     def setdata(self, value):
-        self.__data = numpy.asmatrix(value)
+        self.__data = numpy.array(value)
 
     __data = None
     data = property(fget=getdata, fset=setdata, doc="data")
@@ -137,6 +135,7 @@ class TransferMatrix(DeviceMatrix):
 
         Notation::
 
+
             b12 = H1 a12 and d12 = H2 c12, but b12 == c12,
 
         therefore::
@@ -153,7 +152,7 @@ class TransferMatrix(DeviceMatrix):
 
         """
 
-        return TransferMatrix(M.to_transfer().data * self.data)
+        return TransferMatrix(M.to_transfer().data @ self.data)
 
 
 class ChainMatrix(DeviceMatrix):
@@ -227,7 +226,7 @@ class ChainMatrix(DeviceMatrix):
 
         """
 
-        return ChainMatrix(self.data * M.to_chain().data)
+        return ChainMatrix(self.data @ M.to_chain().data)
 
 
 def composeTM(A, B):

{electromagneticpython-2.2.3 → electromagneticpython-2.2.4}/EMpy/materials.py

@@ -1,5 +1,3 @@
-# pylint: disable=W0622,R0903,R0902,R0913,W0633
-
 """Functions and objects to manipulate materials.
 
 A material is an object with a refractive index function.

{electromagneticpython-2.2.3 → electromagneticpython-2.2.4}/EMpy/modesolvers/FD.py

@@ -1,7 +1,3 @@
-# pylint: disable=line-too-long,too-many-locals,too-many-statements,too-many-branches
-# pylint: disable=redefined-builtin,wildcard-import,unused-wildcard-import
-# pylint: disable=attribute-defined-outside-init,too-many-instance-attributes
-# pylint: disable=arguments-differ,too-many-arguments
 """Finite Difference Modesolver.
 
 @see: Fallahkhair, "Vector Finite Difference Modesolver for Anisotropic Dielectric Waveguides",

{electromagneticpython-2.2.3 → electromagneticpython-2.2.4}/EMpy/transfer_matrix.py

@@ -39,8 +39,8 @@ class IsotropicTransferMatrix(TransferMatrix):
 
         INPUT
         multilayer = Multilayer obj describing the sequence of layers.
-        theta_inc = angle of the incident wave (in radiant) wrt the
-        normal.
+        theta_inc = angle of the incident wave (in radians) wrt the
+        normal, measured in the first layer (multilayer[0]).
         """
         if not multilayer.isIsotropic():
             raise ValueError(
@@ -118,12 +118,14 @@ class IsotropicTransferMatrix(TransferMatrix):
 
             Rs[iwl] = np.absolute(rs) ** 2
             Ts[iwl] = (
-                np.absolute((n[-1] * np.cos(theta[-1])) / (n[0] * np.cos(theta[0])))
+                np.real(n[-1] * np.cos(theta[-1]))
+                / np.real(n[0] * np.cos(theta[0]))
                 * np.absolute(ts) ** 2
             )
             Rp[iwl] = np.absolute(rp) ** 2
             Tp[iwl] = (
-                np.absolute((n[-1] * np.cos(theta[-1])) / (n[0] * np.cos(theta[0])))
+                np.real(np.conj(n[-1]) * np.cos(theta[-1]))
+                / np.real(np.conj(n[0]) * np.cos(theta[0]))
                 * np.absolute(tp) ** 2
             )
 
@@ -160,7 +162,9 @@ class AnisotropicTransferMatrix(TransferMatrix):
         INPUT
         multilayer = Multilayer obj describing the sequence of layers.
         theta_inc_x and theta_inc_y = angles of the incident wave (in
-        radiant) wrt the normal.
+        radians) wrt the normal in vacuum, i.e. kx = k0 * sin(theta_inc_x)
+        and ky = k0 * sin(theta_inc_y) (unlike IsotropicTransferMatrix,
+        where theta_inc is measured in multilayer[0]).
         """
         TransferMatrix.__init__(self, multilayer)
         self.theta_inc_x = theta_inc_x

{electromagneticpython-2.2.3 → electromagneticpython-2.2.4}/EMpy/utils.py

@@ -1,5 +1,3 @@
-# pylint: disable=R0913,R0914,W0201,W0622,C0302,R0902,R0903,W1001,W0612,W0613
-
 """Useful functions and objects used more or less everywhere."""
 
 __author__ = "Lorenzo Bolla"

{electromagneticpython-2.2.3 → electromagneticpython-2.2.4/ElectroMagneticPython.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ElectroMagneticPython
-Version: 2.2.3
+Version: 2.2.4
 Summary: EMpy - ElectroMagnetic Python
 Author-email: Lorenzo Bolla <code@lbolla.info>
 License-Expression: MIT

{electromagneticpython-2.2.3 → electromagneticpython-2.2.4}/ElectroMagneticPython.egg-info/SOURCES.txt

@@ -29,6 +29,8 @@ ElectroMagneticPython.egg-info/PKG-INFO
 ElectroMagneticPython.egg-info/SOURCES.txt
 ElectroMagneticPython.egg-info/dependency_links.txt
 ElectroMagneticPython.egg-info/requires.txt
+ElectroMagneticPython.egg-info/scm_file_list.json
+ElectroMagneticPython.egg-info/scm_version.json
 ElectroMagneticPython.egg-info/top_level.txt
 examples/ex_APRR.py
 examples/ex_NRR.py
@@ -54,4 +56,5 @@ tests/__init__.py
 tests/devices_test.py
 tests/materials_test.py
 tests/modesolvers_test.py
+tests/transfer_matrix_test.py
 tests/utils_test.py

electromagneticpython-2.2.4/ElectroMagneticPython.egg-info/scm_file_list.json

@@ -0,0 +1,56 @@
+{
+  "files": [
+    ".bumpversion.cfg",
+    ".flake8",
+    ".github/dependabot.yml",
+    ".github/workflows/codeql-analysis.yml",
+    ".github/workflows/publish.yml",
+    ".github/workflows/python-app.yml",
+    ".gitignore",
+    "EMpy/RCWA.py",
+    "EMpy/__init__.py",
+    "EMpy/constants.py",
+    "EMpy/devices.py",
+    "EMpy/materials.py",
+    "EMpy/modesolvers/FD.py",
+    "EMpy/modesolvers/FMM.py",
+    "EMpy/modesolvers/__init__.py",
+    "EMpy/modesolvers/geometries.py",
+    "EMpy/modesolvers/interface.py",
+    "EMpy/scattering.py",
+    "EMpy/transfer_matrix.py",
+    "EMpy/utils.py",
+    "LICENSE",
+    "Makefile",
+    "README.md",
+    "examples/ex_APRR.py",
+    "examples/ex_NRR.py",
+    "examples/ex_RCWA.py",
+    "examples/ex_RCWA_2.py",
+    "examples/ex_RCWA_3.py",
+    "examples/ex_SRR.py",
+    "examples/ex_laser_etch_monitor.py",
+    "examples/ex_modesolver.py",
+    "examples/ex_modesolver_2.py",
+    "examples/ex_modesolver_3.py",
+    "examples/ex_modesolver_4.py",
+    "examples/ex_modesolver_dch.py",
+    "examples/ex_modesolver_dch_anis.py",
+    "examples/ex_transfer_matrix.py",
+    "examples/ex_transfer_matrix_2.py",
+    "examples/ex_transfer_matrix_2.py.png",
+    "examples/ex_transfer_matrix_3.py",
+    "examples/lc_data.dat",
+    "examples/nk.py",
+    "mypy.ini",
+    "pyproject.toml",
+    "scripts/FDTD.py",
+    "tests/__init__.py",
+    "tests/devices_test.py",
+    "tests/materials_test.py",
+    "tests/modesolvers_test.py",
+    "tests/transfer_matrix_test.py",
+    "tests/utils_test.py",
+    "uv.lock"
+  ]
+}

electromagneticpython-2.2.4/ElectroMagneticPython.egg-info/scm_version.json

@@ -0,0 +1,8 @@
+{
+  "tag": "2.2.4",
+  "distance": 0,
+  "node": "gd8c276251899f9a6cf5bdd352a4fe68a53c60e61",
+  "dirty": false,
+  "branch": "HEAD",
+  "node_date": "2026-11-14"
+}

{electromagneticpython-2.2.3/ElectroMagneticPython.egg-info → electromagneticpython-2.2.4}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: ElectroMagneticPython
-Version: 2.2.3
+Version: 2.2.4
 Summary: EMpy - ElectroMagnetic Python
 Author-email: Lorenzo Bolla <code@lbolla.info>
 License-Expression: MIT

{electromagneticpython-2.2.3 → electromagneticpython-2.2.4}/examples/ex_transfer_matrix_2.py

@@ -40,5 +40,5 @@ plt.xlabel("wavelength /m")
 plt.ylabel("Power /dB")
 plt.grid()
 plt.xlim(wls.min(), wls.max())
-plt.savefig(__file__ + ".png")
+# plt.savefig(__file__ + ".png")
 plt.show()

electromagneticpython-2.2.4/tests/transfer_matrix_test.py

@@ -0,0 +1,134 @@
+from unittest import TestCase
+
+import numpy as np
+from numpy.testing import assert_allclose
+
+import EMpy
+
+
+def _make_iso_multilayer(n_list, d_list):
+    layers = EMpy.utils.Multilayer()
+    for n, d in zip(n_list, d_list):
+        mat = EMpy.materials.IsotropicMaterial(
+            "mat", n0=EMpy.materials.RefractiveIndex(n0_const=n)
+        )
+        layers.append(EMpy.utils.Layer(mat, d))
+    return layers
+
+
+def _make_aniso_multilayer(n_list, d_list):
+    layers = EMpy.utils.Multilayer()
+    for n, d in zip(n_list, d_list):
+        eps = EMpy.materials.EpsilonTensor(
+            epsilon_tensor_const=(n**2) * EMpy.constants.eps0 * np.eye(3)
+        )
+        mat = EMpy.materials.AnisotropicMaterial("mat", epsilon_tensor=eps)
+        layers.append(EMpy.utils.Layer(mat, d))
+    return layers
+
+
+class TransferMatrixTest(TestCase):
+
+    def test_fresnel_brewster(self):
+        n1, n2 = 1.0, 1.5
+        theta_b = np.arctan2(n2, n1)
+        layers = _make_iso_multilayer([n1, n2], [np.inf, np.inf])
+        sol = EMpy.transfer_matrix.IsotropicTransferMatrix(layers, theta_b).solve(
+            [1.55e-6]
+        )
+        assert_allclose(sol.Rp, [0.0], atol=1e-15)
+        assert_allclose(sol.Tp, [1.0], atol=1e-15)
+        assert_allclose(sol.Rs + sol.Ts, [1.0], atol=1e-15)
+
+    def test_total_internal_reflection(self):
+        layers = _make_iso_multilayer([1.5, 1.0], [np.inf, np.inf])
+        theta_inc = EMpy.utils.deg2rad(60.0)
+        sol = EMpy.transfer_matrix.IsotropicTransferMatrix(layers, theta_inc).solve(
+            [1.55e-6]
+        )
+        assert_allclose(sol.Rs, [1.0], atol=1e-15)
+        assert_allclose(sol.Ts, [0.0], atol=1e-15)
+        assert_allclose(sol.Rp, [1.0], atol=1e-15)
+        assert_allclose(sol.Tp, [0.0], atol=1e-15)
+
+    def test_frustrated_total_internal_reflection(self):
+        layers = _make_iso_multilayer([1.5, 1.0, 1.5], [np.inf, 0.2e-6, np.inf])
+        theta_inc = EMpy.utils.deg2rad(60.0)
+        sol = EMpy.transfer_matrix.IsotropicTransferMatrix(layers, theta_inc).solve(
+            [1.55e-6]
+        )
+        self.assertTrue((0.0 < sol.Ts[0] < 1.0) and (0.0 < sol.Tp[0] < 1.0))
+        assert_allclose(sol.Rs + sol.Ts, [1.0], atol=1e-15)
+        assert_allclose(sol.Rp + sol.Tp, [1.0], atol=1e-15)
+
+    def test_absorbing_exit_medium(self):
+        wls = [1.55e-6]
+        theta_inc = EMpy.utils.deg2rad(40.0)
+        n_list = [1.0, 1.5 - 0.5j]
+        d_list = [np.inf, np.inf]
+
+        iso_layers = _make_iso_multilayer(n_list, d_list)
+        aniso_layers = _make_aniso_multilayer(n_list, d_list)
+
+        sol_iso = EMpy.transfer_matrix.IsotropicTransferMatrix(
+            iso_layers, theta_inc
+        ).solve(wls)
+        sol_aniso = EMpy.transfer_matrix.AnisotropicTransferMatrix(
+            aniso_layers, theta_inc, 0.0
+        ).solve(wls)
+
+        # Single interface from lossless medium: net flux across z=0 is conserved
+        assert_allclose(sol_iso.Rs + sol_iso.Ts, [1.0], atol=1e-15)
+        assert_allclose(sol_iso.Rp + sol_iso.Tp, [1.0], atol=1e-15)
+
+        assert_allclose(sol_iso.Rs, sol_aniso.R[0, 0, :], rtol=1e-6)
+        assert_allclose(sol_iso.Ts, sol_aniso.T[0, 0, :], rtol=1e-6)
+        assert_allclose(sol_iso.Rp, sol_aniso.R[1, 1, :], rtol=1e-6)
+        assert_allclose(sol_iso.Tp, sol_aniso.T[1, 1, :], rtol=1e-6)
+
+    def test_absorbing_slab(self):
+        wls = [1.55e-6]
+        theta_inc = EMpy.utils.deg2rad(20.0)
+        n_list = [1.0, 2.0 - 0.1j, 1.0]
+        d_list = [np.inf, 0.5e-6, np.inf]
+
+        iso_layers = _make_iso_multilayer(n_list, d_list)
+        aniso_layers = _make_aniso_multilayer(n_list, d_list)
+
+        sol_iso = EMpy.transfer_matrix.IsotropicTransferMatrix(
+            iso_layers, theta_inc
+        ).solve(wls)
+        sol_aniso = EMpy.transfer_matrix.AnisotropicTransferMatrix(
+            aniso_layers, theta_inc, 0.0
+        ).solve(wls)
+
+        self.assertLess(sol_iso.Rs[0] + sol_iso.Ts[0], 1.0)
+        self.assertLess(sol_iso.Rp[0] + sol_iso.Tp[0], 1.0)
+        assert_allclose(sol_iso.Rs, sol_aniso.R[0, 0, :], rtol=1e-6)
+        assert_allclose(sol_iso.Ts, sol_aniso.T[0, 0, :], rtol=1e-6)
+        assert_allclose(sol_iso.Rp, sol_aniso.R[1, 1, :], rtol=1e-6)
+        assert_allclose(sol_iso.Tp, sol_aniso.T[1, 1, :], rtol=1e-6)
+
+    def test_iso_vs_aniso_non_unit_incident_index(self):
+        wls = np.linspace(1.4e-6, 1.7e-6, 5)
+        n_list = [1.5, 2.1, 1.0]
+        d_list = [np.inf, 0.3e-6, np.inf]
+        theta_inc = EMpy.utils.deg2rad(30.0)
+        theta_vac = np.arcsin(n_list[0] * np.sin(theta_inc))
+
+        iso_layers = _make_iso_multilayer(n_list, d_list)
+        aniso_layers = _make_aniso_multilayer(n_list, d_list)
+
+        sol_iso = EMpy.transfer_matrix.IsotropicTransferMatrix(
+            iso_layers, theta_inc
+        ).solve(wls)
+        sol_aniso = EMpy.transfer_matrix.AnisotropicTransferMatrix(
+            aniso_layers, theta_vac, 0.0
+        ).solve(wls)
+
+        assert_allclose(sol_iso.Rs + sol_iso.Ts, 1.0, atol=1e-14)
+        assert_allclose(sol_iso.Rp + sol_iso.Tp, 1.0, atol=1e-14)
+        assert_allclose(sol_iso.Rs, sol_aniso.R[0, 0, :], rtol=1e-6)
+        assert_allclose(sol_iso.Ts, sol_aniso.T[0, 0, :], rtol=1e-6)
+        assert_allclose(sol_iso.Rp, sol_aniso.R[1, 1, :], rtol=1e-6)
+        assert_allclose(sol_iso.Tp, sol_aniso.T[1, 1, :], rtol=1e-6)