ElectroMagneticPython 2.2.1__tar.gz → 2.2.4__tar.gz

electromagneticpython-2.2.4/.bumpversion.cfg

@@ -0,0 +1,6 @@
+[bumpversion]
+commit = True
+sign_tags = True
+tag = True
+tag_name = {new_version}
+current_version = 2.2.4

electromagneticpython-2.2.4/.flake8

@@ -0,0 +1,3 @@
+[flake8]
+max-line-length = 130
+ignore = E203,E741,E501,W503,W292

electromagneticpython-2.2.4/.github/dependabot.yml

@@ -0,0 +1,13 @@
+# Keep GitHub Actions up to date with GitHub's Dependabot...
+# https://docs.github.com/en/code-security/dependabot/working-with-dependabot/keeping-your-actions-up-to-date-with-dependabot
+# https://docs.github.com/en/code-security/dependabot/dependabot-version-updates/configuration-options-for-the-dependabot.yml-file#package-ecosystem
+# Working with uv
+# https://docs.astral.sh/uv/guides/integration/dependency-bots/#dependabot
+
+version: 2
+
+updates:
+  - package-ecosystem: uv
+    directory: /
+    schedule:
+      interval: weekly

electromagneticpython-2.2.4/.github/workflows/codeql-analysis.yml

@@ -0,0 +1,70 @@
+# For most projects, this workflow file will not need changing; you simply need
+# to commit it to your repository.
+#
+# You may wish to alter this file to override the set of languages analyzed,
+# or to provide custom queries or build logic.
+#
+# ******** NOTE ********
+# We have attempted to detect the languages in your repository. Please check
+# the `language` matrix defined below to confirm you have the correct set of
+# supported CodeQL languages.
+#
+name: "CodeQL"
+
+on:
+  push:
+    branches: [ master ]
+  pull_request:
+    # The branches below must be a subset of the branches above
+    branches: [ master ]
+  schedule:
+    - cron: '43 12 * * 4'
+
+jobs:
+  analyze:
+    name: Analyze
+    runs-on: ubuntu-latest
+    permissions:
+      actions: read
+      contents: read
+      security-events: write
+
+    strategy:
+      fail-fast: false
+      matrix:
+        language: [ 'python' ]
+        # CodeQL supports [ 'cpp', 'csharp', 'go', 'java', 'javascript', 'python', 'ruby' ]
+        # Learn more about CodeQL language support at https://git.io/codeql-language-support
+
+    steps:
+    - name: Checkout repository
+      uses: actions/checkout@v6
+
+    # Initializes the CodeQL tools for scanning.
+    - name: Initialize CodeQL
+      uses: github/codeql-action/init@v4
+      with:
+        languages: ${{ matrix.language }}
+        # If you wish to specify custom queries, you can do so here or in a config file.
+        # By default, queries listed here will override any specified in a config file.
+        # Prefix the list here with "+" to use these queries and those in the config file.
+        # queries: ./path/to/local/query, your-org/your-repo/queries@main
+
+    # Autobuild attempts to build any compiled languages  (C/C++, C#, or Java).
+    # If this step fails, then you should remove it and run the build manually (see below)
+    - name: Autobuild
+      uses: github/codeql-action/autobuild@v4
+
+    # ℹ️ Command-line programs to run using the OS shell.
+    # 📚 https://git.io/JvXDl
+
+    # ✏️ If the Autobuild fails above, remove it and uncomment the following three lines
+    #    and modify them (or add more) to build your code if your project
+    #    uses a compiled language
+
+    #- run: |
+    #   make bootstrap
+    #   make release
+
+    - name: Perform CodeQL Analysis
+      uses: github/codeql-action/analyze@v4

electromagneticpython-2.2.4/.github/workflows/publish.yml

@@ -0,0 +1,36 @@
+name: Release and Publish
+
+
+on:
+  push:
+    tags:
+      - '[0-9]+.[0-9]+.[0-9]+'
+
+jobs:
+  release:
+    runs-on: ubuntu-latest
+    steps:
+      - name: Checkout
+        uses: actions/checkout@v6
+      - name: Release on GitHub
+        uses: softprops/action-gh-release@v2
+  publish:
+    runs-on: ubuntu-latest
+    environment:
+      name: pypi
+    permissions:
+      id-token: write
+      contents: read
+    steps:
+      - uses: actions/checkout@v6
+
+      - name: Install uv and set the Python version
+        uses: astral-sh/setup-uv@v7
+        with:
+          python-version: '3.12'
+
+      - name: Build
+        run: uv build
+
+      - name: Publish
+        run: uv publish

electromagneticpython-2.2.4/.github/workflows/python-app.yml

@@ -0,0 +1,40 @@
+name: Python application
+
+on:
+  push:
+    branches: [ master ]
+  pull_request:
+    branches: [ master ]
+
+jobs:
+  build:
+
+    runs-on: ubuntu-latest
+
+    strategy:
+      fail-fast: false
+      matrix:
+        python-version: ["3.10", "3.11", "3.12", "3.13", "3.14"]
+
+    steps:
+    - uses: actions/checkout@v6
+
+    - name: Install uv and set the Python version
+      uses: astral-sh/setup-uv@v7
+      with:
+        python-version: ${{ matrix.python-version }}
+
+    - name: Install the project
+      run: uv sync --locked --all-extras --dev
+
+    - name: Lint with flake8
+      run: uv run flake8 EMpy tests examples scripts
+
+    - name: Lint with pyflakes
+      run: uv run pyflakes EMpy tests examples scripts
+
+    - name: Lint with mypy
+      run: uv run mypy EMpy tests examples scripts
+
+    - name: Run tests
+      run: uv run pytest

electromagneticpython-2.2.4/.gitignore

@@ -0,0 +1,8 @@
+*.egg-info
+.mypy_cache
+.venv
+.vscode
+build
+dist
+__pycache__
+EMpy/_version.py

{electromagneticpython-2.2.1 → electromagneticpython-2.2.4}/EMpy/RCWA.py

@@ -12,6 +12,7 @@ for every diffraction order and layer thickness, and an anisotropic
 one, only stable for low diffraction orders and thin layers.
 
 """
+
 __author__ = "Lorenzo Bolla"
 
 import numpy as np
@@ -95,7 +96,6 @@ def dispersion_relation_extraordinary(kx, ky, k, nO, nE, c):
 
 
 class RCWA:
-
     """Class to handle the RCWA solvers.
 
     NOTE
@@ -168,7 +168,6 @@ class RCWA:
 
 
 class IsotropicRCWA(RCWA):
-
     """Isotropic RCWA solver."""
 
     def solve(self, wls):
@@ -215,7 +214,11 @@ class IsotropicRCWA(RCWA):
         DE3 = np.zeros_like(DE1)
 
         dirk1 = np.array(
-            [np.sin(alpha) * np.cos(delta), np.sin(alpha) * np.sin(delta), np.cos(alpha)]
+            [
+                np.sin(alpha) * np.cos(delta),
+                np.sin(alpha) * np.sin(delta),
+                np.cos(alpha),
+            ]
         )
 
         # usefull matrices
@@ -323,7 +326,9 @@ class IsotropicRCWA(RCWA):
                     # H_S_reduced = np.dot(Ky, Ky) + np.dot(Kx, linsolve(E, np.dot(Kx, E11))) - E11
                     H_S_reduced = (
                         (ky / k) ** 2 * I
-                        + kx[:, np.newaxis] / k * linsolve(E, kx[:, np.newaxis] / k * E11)
+                        + kx[:, np.newaxis]
+                        / k
+                        * linsolve(E, kx[:, np.newaxis] / k * E11)
                         - E11
                     )
 
@@ -448,7 +453,6 @@ class IsotropicRCWA(RCWA):
 
 
 class AnisotropicRCWA(RCWA):
-
     """Anisotropic RCWA solver."""
 
     def solve(self, wls):
@@ -487,14 +491,20 @@ class AnisotropicRCWA(RCWA):
         # grating on the xy plane
         K = 2 * pi / LAMBDA * np.array([np.sin(phi), 0.0, np.cos(phi)], dtype=complex)
         dirk1 = np.array(
-            [np.sin(alpha) * np.cos(delta), np.sin(alpha) * np.sin(delta), np.cos(alpha)]
+            [
+                np.sin(alpha) * np.cos(delta),
+                np.sin(alpha) * np.sin(delta),
+                np.cos(alpha),
+            ]
         )
 
         # D polarization vector
         u = np.array(
             [
-                np.cos(psi) * np.cos(alpha) * np.cos(delta) - np.sin(psi) * np.sin(delta),
-                np.cos(psi) * np.cos(alpha) * np.sin(delta) + np.sin(psi) * np.cos(delta),
+                np.cos(psi) * np.cos(alpha) * np.cos(delta)
+                - np.sin(psi) * np.sin(delta),
+                np.cos(psi) * np.cos(alpha) * np.sin(delta)
+                + np.sin(psi) * np.cos(delta),
                 -np.cos(psi) * np.sin(alpha),
             ]
         )
@@ -814,7 +824,11 @@ class AnisotropicRCWA(RCWA):
             denom = (k1[2] - np.dot(u, k1) * u[2]).real
             DEO1[:, iwl] = (
                 -(
-                    (np.absolute(ROx) ** 2 + np.absolute(ROy) ** 2 + np.absolute(ROz) ** 2)
+                    (
+                        np.absolute(ROx) ** 2
+                        + np.absolute(ROy) ** 2
+                        + np.absolute(ROz) ** 2
+                    )
                     * np.conj(kO1i[2, :])
                     - (ROx * kO1i[0, :] + ROy * kO1i[1, :] + ROz * kO1i[2, :])
                     * np.conj(ROz)
@@ -823,7 +837,11 @@ class AnisotropicRCWA(RCWA):
             )
             DEE1[:, iwl] = (
                 -(
-                    (np.absolute(REx) ** 2 + np.absolute(REy) ** 2 + np.absolute(REz) ** 2)
+                    (
+                        np.absolute(REx) ** 2
+                        + np.absolute(REy) ** 2
+                        + np.absolute(REz) ** 2
+                    )
                     * np.conj(kE1i[2, :])
                     - (REx * kE1i[0, :] + REy * kE1i[1, :] + REz * kE1i[2, :])
                     * np.conj(REz)
@@ -833,12 +851,14 @@ class AnisotropicRCWA(RCWA):
             DEO3[:, iwl] = (
                 (np.absolute(TOx) ** 2 + np.absolute(TOy) ** 2 + np.absolute(TOz) ** 2)
                 * np.conj(kO3i[2, :])
-                - (TOx * kO3i[0, :] + TOy * kO3i[1, :] + TOz * kO3i[2, :]) * np.conj(TOz)
+                - (TOx * kO3i[0, :] + TOy * kO3i[1, :] + TOz * kO3i[2, :])
+                * np.conj(TOz)
             ).real / denom
             DEE3[:, iwl] = (
                 (np.absolute(TEx) ** 2 + np.absolute(TEy) ** 2 + np.absolute(TEz) ** 2)
                 * np.conj(kE3i[2, :])
-                - (TEx * kE3i[0, :] + TEy * kE3i[1, :] + TEz * kE3i[2, :]) * np.conj(TEz)
+                - (TEx * kE3i[0, :] + TEy * kE3i[1, :] + TEz * kE3i[2, :])
+                * np.conj(TEz)
             ).real / denom
 
         # save the results

{electromagneticpython-2.2.1 → electromagneticpython-2.2.4}/EMpy/__init__.py

@@ -10,6 +10,7 @@ The package contains some useful routines to study electromagnetic problems with
 It is based on U{numpy<http://numpy.scipy.org>} and U{scipy<http://www.scipy.org>}.
 
 """
+
 __all__ = [
     "constants",
     "devices",
@@ -30,6 +31,13 @@ from . import RCWA
 from . import scattering
 from . import transfer_matrix
 from . import utils
-from .version import version
 
-__version__ = version
+try:
+    from ._version import version as __version__
+except ImportError:
+    # Fallback for development environments or edge cases
+    try:
+        from importlib.metadata import version, PackageNotFoundError
+        __version__ = version("electromagneticpython")
+    except PackageNotFoundError:
+        __version__ = "0.0.0.dev0"

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.1 → 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.1 → 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.
 
@@ -31,6 +29,7 @@ devices, new ones can be studied.
     isnumber=8367>}
 
 """
+
 __author__ = "Lorenzo Bolla"
 
 import numpy
@@ -39,7 +38,6 @@ from functools import reduce
 
 
 class DeviceMatrix:
-
     """Device Matrix.
 
     Base class for TransferMatrix and ChainMatrix.
@@ -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")
@@ -83,7 +81,6 @@ class DeviceMatrix:
 
 
 class TransferMatrix(DeviceMatrix):
-
     """Transfer Matrix.
 
     Notation::
@@ -138,6 +135,7 @@ class TransferMatrix(DeviceMatrix):
 
         Notation::
 
+
             b12 = H1 a12 and d12 = H2 c12, but b12 == c12,
 
         therefore::
@@ -154,11 +152,10 @@ class TransferMatrix(DeviceMatrix):
 
         """
 
-        return TransferMatrix(M.to_transfer().data * self.data)
+        return TransferMatrix(M.to_transfer().data @ self.data)
 
 
 class ChainMatrix(DeviceMatrix):
-
     """Chain Matrix.
 
     Notation::
@@ -229,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):
@@ -289,7 +286,6 @@ def composeCMlist(CMlist):
 
 
 class Device:
-
     """Generic device.
 
     Notes
@@ -307,7 +303,6 @@ class Device:
 
 
 class Coupler(Device):
-
     """Coupler.
 
     Examples
@@ -436,7 +431,6 @@ class Coupler(Device):
 
 
 class Line(Device):
-
     """Line.
 
     Notes
@@ -570,7 +564,6 @@ class Line(Device):
 
 
 class MZ(Device):
-
     """Mach-Zehnder Interferometer (MZI).
 
     A MZI is made of two L{Coupler}s connected by a L{Line}.
@@ -679,7 +672,6 @@ class MZ(Device):
 
 
 class APRR(Device):
-
     """All-Pass Ring Resonator.
 
     An APRR is made of a L{Coupler} with a feedback line (the ring).
@@ -746,7 +738,6 @@ class APRR(Device):
 
 
 class SRR(Device):
-
     """Single Ring Resonator.
 
     A SRR is made of the parallel of a L{Coupler} a feedback line (the
@@ -872,7 +863,6 @@ class SRR(Device):
 
 
 class NRR(Device):
-
     """N-Rings Resonator.
 
     A NRR is made of the parallel of two lines connected by N ring resonators.
@@ -996,7 +986,6 @@ class NRR(Device):
 
 
 class T_TCRR(Device):
-
     """Tunable Triple-Coupled Ring Resonator.
 
     A T_TCRR is made of a L{MZ}, with a feedback line.
@@ -1056,8 +1045,7 @@ class T_TCRR(Device):
         t4 = numpy.exp(-1j * beta * l[3])
 
         denom = (
-            q**3 * Kbar[2] * t3 * t4 * (Kbar[0] * Kbar[1] * t1 - K[0] * K[1] * t2)
-            - 1.0
+            q**3 * Kbar[2] * t3 * t4 * (Kbar[0] * Kbar[1] * t1 - K[0] * K[1] * t2) - 1.0
         )
         self.DROP = (
             q**3 * K[2] * t3 * (K[0] * Kbar[1] * t1 + Kbar[0] * K[1] * t2) / denom
@@ -1071,7 +1059,6 @@ class T_TCRR(Device):
 
 
 class T_CRTCRR(Device):
-
     """Tunable Coupled Ring Triple-Coupled Ring Resonator.
 
     NOTE
@@ -1129,16 +1116,7 @@ class T_CRTCRR(Device):
             - q**2 * Kbar[1] * t2**2 * (Kbar[2] - q**3 * t1**2 * N1)
         )
         self.DROP = (
-            1j
-            * q**4
-            * K[0]
-            * K[1]
-            * K[2]
-            * Kbar[0]
-            * t1**1.5
-            * t2
-            * (1 + t31)
-            / denom
+            1j * q**4 * K[0] * K[1] * K[2] * Kbar[0] * t1**1.5 * t2 * (1 + t31) / denom
         )
         self.THRU = (
             q**2
@@ -1156,7 +1134,6 @@ class T_CRTCRR(Device):
 
 
 class T_CTCRR(Device):
-
     """Tunable Coupled Triple Coupled Ring Resonator.
 
     NOTE
@@ -1363,7 +1340,6 @@ class SWG(Device):
 
 
 class Etalon(Device):
-
     """Etalon.
 
     @ivar R: Reflectivity.
@@ -1402,7 +1378,7 @@ class Etalon(Device):
         self.wl0 = numpy.mean(wls)
         self.FSRwl = self.wl0**2 / (2 * n * l * numpy.cos(self.theta) + self.wl0)
 
-        (self.f0, self.FSR) = EMpy.utils.wl2f(self.wl0, self.FSRwl)
+        self.f0, self.FSR = EMpy.utils.wl2f(self.wl0, self.FSRwl)
 
         self.FINESSE = numpy.pi / (2 * numpy.arcsin(1.0 / numpy.sqrt(self.F)))
         # self.FINESSE = numpy.pi * numpy.sqrt(self.F) / 2.

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

@@ -1,10 +1,9 @@
-# pylint: disable=W0622,R0903,R0902,R0913,W0633
-
 """Functions and objects to manipulate materials.
 
 A material is an object with a refractive index function.
 
 """
+
 from functools import partial
 
 import numpy
@@ -15,7 +14,6 @@ __author__ = "Lorenzo Bolla"
 
 
 class Material:
-
     """Generic class to handle materials.
 
     This class is intended to be subclassed to obtain isotropic and
@@ -28,7 +26,6 @@ class Material:
 
 
 class RefractiveIndex:
-
     """Refractive Index.
 
     Unaware of temperature.
@@ -146,7 +143,6 @@ class RefractiveIndex:
 
 
 class ThermalOpticCoefficient:
-
     """Thermal Optic Coefficient."""
 
     def __init__(self, data=None, T0=300):
@@ -168,7 +164,6 @@ class ThermalOpticCoefficient:
 
 
 class IsotropicMaterial(Material):
-
     """Subclasses Material to describe isotropic materials.
 
     Frequency dispersion and thermic aware.
@@ -232,7 +227,6 @@ class EpsilonTensor:
 
 
 class AnisotropicMaterial(Material):
-
     """Subclass Material to describe anisotropic materials.
 
     No frequency dispersion nor thermic aware.
@@ -292,7 +286,6 @@ BK7 = IsotropicMaterial(
 
 
 class LiquidCrystal(Material):
-
     """Liquid Crystal.
 
     A liquid crystal is determined by it ordinary and extraordinary