Broken FEEC projections on polar domains

.github/actions/install_petsc4py/action.yml

@@ -6,7 +6,7 @@ runs:
   steps:
 
     - name: Cache PETSc
-      uses: actions/cache@v4
+      uses: actions/cache@v6
       id: cache-petsc
       env:
         cache-name: cache-PETSc

.github/workflows/documentation.yml

@@ -82,10 +82,10 @@ jobs:
         python docs/update_links.py
 
     - name: Setup Pages
-      uses: actions/configure-pages@v5
+      uses: actions/configure-pages@v6
 
     - name: Upload artifact
-      uses: actions/upload-pages-artifact@v3
+      uses: actions/upload-pages-artifact@v5
       with:
         path: 'docs/build/html'
 
@@ -99,4 +99,4 @@ jobs:
     steps:
     - name: Deploy to GitHub Pages
       id: deployment
-      uses: actions/deploy-pages@v4
+      uses: actions/deploy-pages@v5

CHANGELOG.md

@@ -6,11 +6,16 @@ All notable changes to this project will be documented in this file.
 
 ### Added
 
+-   #576 : Add module `psydac.feec.polar.conga_projections` with broken FEEC polar projections in 2D
+-   #576 : Add 2D Poisson and TE Maxwell examples on polar mapped domains in `psydac/feec/polar/examples`
+-   #576 : Add module `psydac.utilities.parallel_utils` for parallel execution and gathering variable-length arrays
+-   #576 : Add module `psydac.utilities.operators` with class `Laplacian` used in some old examples
 -   #577 : Add an installation configuration option to choose the backend language
 -   [DEVELOPER] Create action `install_petsc4py` to install PETSc & `petsc4py` w/ complex support
 
 ### Fixed
 
+-   #576 : Require `sympde==0.19.3` which fixes a bug in the linearity checks
 -   #579 : Require `h5py>=3.16` which installs correctly with `setuptools>=81.0`
 -   #579 : Don't run postprocessing unit tests with `pytest-xdist` because `h5py` is not thread-safe
 -   #579 : Return error code on failure of the `psydac test` and `psydac compile` commands
@@ -26,12 +31,14 @@ All notable changes to this project will be documented in this file.
 
 ### Changed
 
+-   #576 : Use latest version of Igakit (commit dalcinl/igakit@92ee097 of 2026/07/24) which supports NumPy >= 2.4
 -   #595 : Use PETSc 3.25.5 whose Python bindings `petsc4py` are built correctly with `cython>=3`
 -   #580 : Use PETSc 3.25.0 whose Python bindings `petsc4py` install correctly with `setuptools>=81.0`
 -   #579 : Require `pyccel>=2.2.3` which can compile all kernels with C
 -   #579 : Require `numpy>=2.1` to support Python >= 3.10
 -   #579 : Require `pytest>=9.0` and use `pytest.toml` instead of `pytest.ini` for Pytest configuration
 -   #579 : Move coverage configuration from `pyproject.toml` to `psydac/pytest.toml`
+-   [DEVELOPER] Update GitHub Actions for repository checkout and Python setup
 -   [DEVELOPER] Rename actions: `macos/ubuntu_install` -> `macos/ubuntu_installations`
 -   [DEVELOPER] Do not check file changes to trigger testing workflow on PRs
 -   [DEVELOPER] Run documentation workflow on pushes to `devel` whenever `README.md` is modified

docs/source/modules/feec.polar.rst

@@ -0,0 +1,10 @@
+feec.polar
+==========
+
+.. currentmodule:: psydac.feec
+.. autosummary::
+    :nosignatures:
+    :toctree: STUBDIR
+    :template: autosummary/module.rst
+
+    polar.conga_projections

docs/source/modules/feec.rst

@@ -14,3 +14,10 @@ feec
     feec.multipatch_domain_utilities
     feec.pull_push
     feec.pushforward
+
+feec submodules
+--------------
+.. toctree::
+    :maxdepth: 1
+
+    feec.polar

docs/source/modules/utilities.rst

@@ -7,6 +7,8 @@ utilities
     :toctree: STUBDIR
     :template: autosummary/module.rst
 
+    utilities.operators
+    utilities.parallel_utils
     utilities.quadratures
     utilities.utils
     utilities.vtk

examples/old_examples/poisson_2d_mapping.py

@@ -21,43 +21,12 @@
 from psydac.fem.tensor           import TensorFemSpace
 from psydac.fem.basic            import FemField
 from psydac.mapping.discrete     import SplineMapping
+from psydac.utilities.operators  import Laplacian
 from psydac.utilities.utils      import refine_array_1d
 from psydac.cad.geometry         import Geometry
 from psydac.ddm.cart             import DomainDecomposition
 from psydac.polar.c1_projections import C1Projector
 
-#==============================================================================
-class Laplacian:
-
-    def __init__(self, mapping):
-
-        assert isinstance(mapping, CallableMapping)
-
-        sym = mapping.symbolic_mapping
-
-        self._eta        = sym.logical_coordinates
-        self._metric     = sym.metric_expr
-        self._metric_det = sym.metric_det_expr
-
-    # ...
-    def __call__(self, phi):
-
-        from sympy import sqrt, Matrix
-
-        u      = self._eta
-        G      = self._metric
-        sqrt_g = sqrt(self._metric_det)
-
-        # Store column vector of partial derivatives of phi w.r.t. uj
-        dphi_du = Matrix([phi.diff(uj) for uj in u])
-
-        # Compute gradient of phi in tangent basis: A = G^(-1) dphi_du
-        A = G.LUsolve(dphi_du)
-
-        # Compute Laplacian of phi using formula for divergence of vector A
-        lapl = sum((sqrt_g * Ai).diff(ui) for ui, Ai in zip(u, A)) / sqrt_g
-
-        return lapl
 
 #==============================================================================
 class Poisson2D:

examples/old_examples/poisson_3d_target_torus.py

@@ -12,45 +12,12 @@
 
 from psydac.api.postprocessing import OutputManager, PostProcessManager
 from psydac.fem.basic import FemField
+from psydac.utilities.operators import Laplacian
 from psydac.utilities.utils import refine_array_1d
 
 NUM_DIMS_LOGICAL  = 3
 NUM_DIMS_PHYSICAL = 3
 
-#==============================================================================
-class Laplacian:
-
-    def __init__(self, mapping):
-
-        assert isinstance(mapping, Mapping)
-
-        sym = mapping
-
-        self._eta        = sym.logical_coordinates
-        self._metric     = sym.metric_expr
-        self._metric_det = sym.metric_det_expr
-
-    # ...
-    def __call__(self, phi):
-
-        from sympy import sqrt, Matrix
-
-        u      = self._eta
-        G      = self._metric
-        sqrt_g = sqrt( self._metric_det )
-
-        # Store column vector of partial derivatives of phi w.r.t. uj
-        dphi_du = Matrix( [phi.diff( uj ) for uj in u] )
-
-        # Compute gradient of phi in tangent basis: A = G^(-1) dphi_du
-        A = G.LUsolve( dphi_du )
-
-        # Compute Laplacian of phi using formula for divergence of vector A
-        lapl = sum( (sqrt_g*Ai).diff( ui ) for ui,Ai in zip( u,A ) ) / sqrt_g
-
-        return lapl
-
-
 #==============================================================================
 # Define the Spherical coordinate system
 class TargetTorusMapping(Mapping):

psydac/feec/polar/__init__.py

@@ -0,0 +1,5 @@
+#---------------------------------------------------------------------------#
+# This file is part of PSYDAC which is released under MIT License. See the  #
+# LICENSE file or go to https://github.com/pyccel/psydac/blob/devel/LICENSE #
+# for full license details.                                                 #
+#---------------------------------------------------------------------------#