Migrate project to use uv (dependency manager)

.flake8

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

.github/workflows/publish.yml

@@ -0,0 +1,35 @@
+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

.github/workflows/python-app.yml

@@ -1,6 +1,3 @@
-# This workflow will install Python dependencies, run tests, and lint with a single version of Python
-# For more information, see: https://help.github.com/actions/language-and-framework-guides/using-python-with-github-actions
-
 name: Python application
 
 on:
@@ -17,29 +14,27 @@ jobs:
     strategy:
       fail-fast: false
       matrix:
-        python-version: ["3.10", "3.11", "3.12", "3.13"]
+        python-version: ["3.10", "3.11", "3.12"]
 
     steps:
     - uses: actions/checkout@v6
-    - name: Set up Python ${{ matrix.python-version }}
-      uses: actions/setup-python@v6
+
+    - name: Install uv and set the Python version
+      uses: astral-sh/setup-uv@v7
       with:
-        python-version: "${{ matrix.python-version }}"
-    - name: Install dependencies
-      run: |
-        python -m pip install --upgrade pip
-        pip install flake8 mypy pytest
-        pip install -r requirements.txt
-        pip install --editable .
+        python-version: ${{ matrix.python-version }}
+
+    - name: Install the project
+      run: uv sync --locked --all-extras --dev
+
     - name: Lint with flake8
-      run: |
-        flake8 EMpy tests examples scripts
+      run: uv run flake8 EMpy tests examples scripts
+
     - name: Lint with pyflakes
-      run: |
-        pyflakes EMpy tests examples scripts
+      run: uv run pyflakes EMpy tests examples scripts
+
     - name: Lint with mypy
-      run: |
-        mypy EMpy tests examples scripts
+      run: uv run mypy EMpy tests examples scripts
+
     - name: Run tests
-      run: |
-        pytest
+      run: uv run pytest

.gitignore

@@ -1,10 +1,7 @@
-*.pyc
-*.egg-info/
-*.egg
-.eggs/
-.tox/
-.vscode/
-dist
+*.egg-info
+.mypy_cache
+.venv
+.vscode
 build
-.idea
-.DS_Store
+dist
+__pycache__

EMpy/RCWA.py

@@ -12,6 +12,7 @@
 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 @@ def get_pitch(self):
 
 
 class IsotropicRCWA(RCWA):
-
     """Isotropic RCWA solver."""
 
     def solve(self, wls):
@@ -215,7 +214,11 @@ def solve(self, wls):
         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 @@ def solve(self, wls):
                     # 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 @@ def __str__(self):
 
 
 class AnisotropicRCWA(RCWA):
-
     """Anisotropic RCWA solver."""
 
     def solve(self, wls):
@@ -487,14 +491,20 @@ def solve(self, wls):
         # 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 @@ def solve(self, wls):
             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 @@ def solve(self, wls):
             )
             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 @@ def solve(self, wls):
             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

EMpy/__init__.py

@@ -10,6 +10,7 @@
 It is based on U{numpy<http://numpy.scipy.org>} and U{scipy<http://www.scipy.org>}.
 
 """
+
 __all__ = [
     "constants",
     "devices",