Vertex integrals

.github/workflows/ccpp.yml

@@ -111,15 +111,15 @@ jobs:
       - name: Install FEniCS Python components (default branches/tags)
         if: github.event_name != 'workflow_dispatch'
         run: |
-          pip install git+https://github.com/fenics/ufl.git@${{ env.ufl_ref }}
-          pip install git+https://github.com/fenics/basix.git@${{ env.basix_ref }}
-          pip install git+https://github.com/fenics/ffcx.git@${{ env.ffcx_ref }}
+          pip install git+https://github.com/${{ env.ufl_repository }}.git@${{ env.ufl_ref }}
+          pip install git+https://github.com/${{ env.basix_repository }}.git@${{ env.basix_ref }}
+          pip install git+https://github.com/${{ env.ffcx_repository }}.git@${{ env.ffcx_ref }}
       - name: Install FEniCS Python components
         if: github.event_name == 'workflow_dispatch'
         run: |
-          pip install git+https://github.com/FEniCS/ufl.git@${{ env.ufl_ref }}
-          pip install git+https://github.com/FEniCS/basix.git@${{ env.basix_ref }}
-          pip install git+https://github.com/FEniCS/ffcx.git@${{ env.ffcx_ref }}
+          pip install git+https://github.com/${{ env.ufl_repository }}.git@${{ env.ufl_ref }}
+          pip install git+https://github.com/${{ env.basix_repository }}.git@${{ env.basix_ref }}
+          pip install git+https://github.com/${{ env.ffcx_repository }}.git@${{ env.ffcx_ref }}
 
       - name: Configure and install C++
         run: |
@@ -184,16 +184,15 @@ jobs:
       - name: Install FEniCS Python components (default branches/tags)
         if: github.event_name != 'workflow_dispatch'
         run: |
-          pip install git+https://github.com/FEniCS/ufl.git@${{ env.ufl_ref }}
-          pip install git+https://github.com/FEniCS/basix.git@${{ env.basix_ref }}
-          pip install git+https://github.com/FEniCS/ffcx.git@${{ env.ffcx_ref }}
-
+          pip install git+https://github.com/${{ env.ufl_repository }}.git@${{ env.ufl_ref }}
+          pip install git+https://github.com/${{ env.basix_repository }}.git@${{ env.basix_ref }}
+          pip install git+https://github.com/${{ env.ffcx_repository }}.git@${{ env.ffcx_ref }}
       - name: Install FEniCS Python components
         if: github.event_name == 'workflow_dispatch'
         run: |
-          pip install git+https://github.com/FEniCS/ufl.git@${{ env.ufl_ref }}
-          pip install git+https://github.com/FEniCS/basix.git@${{ env.basix_ref }}
-          pip install git+https://github.com/FEniCS/ffcx.git@${{ env.ffcx_ref }}
+          pip install git+https://github.com/${{ env.ufl_repository }}.git@${{ env.ufl_ref }}
+          pip install git+https://github.com/${{ env.basix_repository }}.git@${{ env.basix_ref }}
+          pip install git+https://github.com/${{ env.ffcx_repository }}.git@${{ env.ffcx_ref }}
 
       - name: Configure, build and install C++ library
         run: |
@@ -266,15 +265,15 @@ jobs:
       - name: Install FEniCS Python components (default branches/tags)
         if: github.event_name != 'workflow_dispatch'
         run: |
-          pip install git+https://github.com/FEniCS/ufl.git@${{ env.ufl_ref }}
-          pip install git+https://github.com/FEniCS/basix.git@${{ env.basix_ref }}
-          pip install git+https://github.com/FEniCS/ffcx.git@${{ env.ffcx_ref }}
+          pip install git+https://github.com/${{ env.ufl_repository }}.git@${{ env.ufl_ref }}
+          pip install git+https://github.com/${{ env.basix_repository }}.git@${{ env.basix_ref }}
+          pip install git+https://github.com/${{ env.ffcx_repository }}.git@${{ env.ffcx_ref }}
       - name: Install FEniCS Python components
         if: github.event_name == 'workflow_dispatch'
         run: |
-          pip install git+https://github.com/FEniCS/ufl.git@${{ env.ufl_ref }}
-          pip install git+https://github.com/FEniCS/basix.git@${{ env.basix_ref }}
-          pip install git+https://github.com/FEniCS/ffcx.git@${{ env.ffcx_ref }}
+          pip install git+https://github.com/${{ env.ufl_repository }}.git@${{ env.ufl_ref }}
+          pip install git+https://github.com/${{ env.basix_repository }}.git@${{ env.basix_ref }}
+          pip install git+https://github.com/${{ env.ffcx_repository }}.git@${{ env.ffcx_ref }}
 
       - name: Configure C++
         run: |

.github/workflows/fenicsx-refs.env

@@ -1,3 +1,6 @@
+basix_repository=FEniCS/basix
 basix_ref=main
+ufl_repository=FEniCS/ufl
 ufl_ref=main
-ffcx_ref=main
+ffcx_repository=schnellerhase/fenics-ffcx
+ffcx_ref=point_measure

.github/workflows/macos.yml

@@ -73,9 +73,9 @@ jobs:
 
       - name: Install FEniCSx dependencies
         run: |
-          python -m pip install git+https://github.com/fenics/ufl.git@${{ env.ufl_ref }}
-          python -m pip install git+https://github.com/fenics/basix.git@${{ env.basix_ref }}
-          python -m pip install git+https://github.com/fenics/ffcx.git@${{ env.ffcx_ref }}
+          python -m pip install git+https://github.com/${{ env.ufl_repository }}.git@${{ env.ufl_ref }}
+          python -m pip install git+https://github.com/${{ env.basix_repository }}.git@${{ env.basix_ref }}
+          python -m pip install git+https://github.com/${{ env.ffcx_repository }}.git@${{ env.ffcx_ref }}
 
       - name: Build and install DOLFINx C++ library
         run: |

.github/workflows/oneapi.yml

@@ -56,13 +56,13 @@ jobs:
           conda list
 
       - name: Install Basix
-        run: uv pip install --no-build-isolation git+https://github.com/FEniCS/basix.git@${{ env.basix_ref }}
+        run: uv pip install --no-build-isolation git+https://github.com/${{ env.basix_repository }}.git@${{ env.basix_ref }}
 
       - name: Clone FFCx
         uses: actions/checkout@v4
         with:
           path: ./ffcx
-          repository: FEniCS/ffcx
+          repository: ${{ env.ffcx_repository }}
           ref: ${{ env.ffcx_ref }}
 
       - name: Install UFCx C interface
@@ -82,7 +82,7 @@ jobs:
 
       - name: Install UFL and FFCx modules
         run: |
-          uv pip install --no-build-isolation git+https://github.com/FEniCS/ufl.git@${{ env.ufl_ref }}
+          uv pip install --no-build-isolation git+https://github.com/${{ env.ufl_repository }}.git@${{ env.ufl_ref }}
           uv pip install --no-build-isolation ffcx/
 
       - name: Build and run DOLFINx C++ unit tests (serial and MPI)

.github/workflows/pyvista.yml

@@ -39,9 +39,9 @@ jobs:
 
       - name: Install FEniCS Python components
         run: |
-          python -m pip install git+https://github.com/fenics/ufl.git@${{ env.ufl_ref }}
-          python -m pip install git+https://github.com/fenics/basix.git@${{ env.basix_ref }}
-          python -m pip install git+https://github.com/fenics/ffcx.git@${{ env.ffcx_ref }}
+          pip install git+https://github.com/${{ env.ufl_repository }}.git@${{ env.ufl_ref }}
+          pip install git+https://github.com/${{ env.basix_repository }}.git@${{ env.basix_ref }}
+          pip install git+https://github.com/${{ env.ffcx_repository }}.git@${{ env.ffcx_ref }}
           apt-get update
           apt-get install -y --no-install-recommends libgl1-mesa-dev xvfb  # pyvista
           apt-get install -y --no-install-recommends libqt5gui5t64 libgl1 # pyvistaqt

.github/workflows/redhat.yml

@@ -33,9 +33,9 @@ jobs:
 
       - name: Install FEniCS Python components
         run: |
-          python3 -m pip install git+https://github.com/fenics/ufl.git@${{ env.ufl_ref }}
-          python3 -m pip install git+https://github.com/fenics/basix.git@${{ env.basix_ref }}
-          python3 -m pip install git+https://github.com/fenics/ffcx.git@${{ env.ffcx_ref }}
+          pip install git+https://github.com/${{ env.ufl_repository }}.git@${{ env.ufl_ref }}
+          pip install git+https://github.com/${{ env.basix_repository }}.git@${{ env.basix_ref }}
+          pip install git+https://github.com/${{ env.ffcx_repository }}.git@${{ env.ffcx_ref }}
 
       - name: Configure C++
         run: cmake -G Ninja -DCMAKE_BUILD_TYPE=Developer -B build -S cpp/

.github/workflows/sonarcloud.yml

@@ -62,9 +62,9 @@ jobs:
           echo "$HOME/.sonar/build-wrapper-linux-x86" >> $GITHUB_PATH
       - name: Install FEniCS Python components
         run: |
-          python -m pip install git+https://github.com/fenics/ufl.git@${{ env.ufl_ref }}
-          python -m pip install git+https://github.com/fenics/basix.git@${{ env.basix_ref }}
-          python -m pip install git+https://github.com/fenics/ffcx.git@${{ env.ffcx_ref }}
+          pip install git+https://github.com/${{ env.ufl_repository }}.git@${{ env.ufl_ref }}
+          pip install git+https://github.com/${{ env.basix_repository }}.git@${{ env.basix_ref }}
+          pip install git+https://github.com/${{ env.ffcx_repository }}.git@${{ env.ffcx_ref }}
       - name: Run build-wrapper
         run: |
           mkdir build

cpp/dolfinx/fem/assemble_scalar_impl.h

@@ -1,4 +1,4 @@
-// Copyright (C) 2019-2025 Garth N. Wells
+// Copyright (C) 2019-2025 Garth N. Wells and Paul T. Kühner
 //
 // This file is part of DOLFINx (https://www.fenicsproject.org)
 //
@@ -150,6 +150,43 @@ T assemble_interior_facets(
   return value;
 }
 
+/// Assemble functional over vertices
+template <dolfinx::scalar T>
+T assemble_vertices(mdspan2_t x_dofmap,
+                    md::mdspan<const scalar_value_t<T>,
+                               md::extents<std::size_t, md::dynamic_extent, 3>>
+                        x,
+                    md::mdspan<const std::int32_t,
+                               md::extents<std::size_t, md::dynamic_extent, 2>>
+                        vertices,
+                    FEkernel<T> auto fn, std::span<const T> constants,
+                    md::mdspan<const T, md::dextents<std::size_t, 2>> coeffs)
+{
+  T value(0);
+  if (vertices.empty())
+    return value;
+
+  // Create data structures used in assembly
+  std::vector<scalar_value_t<T>> cdofs(3 * x_dofmap.extent(1));
+
+  // Iterate over all cells
+  for (std::size_t index = 0; index < vertices.extent(0); ++index)
+  {
+    std::int32_t cell = vertices(index, 0);
+    std::int32_t local_vertex_index = vertices(index, 1);
+
+    // Get cell coordinates/geometry
+    auto x_dofs = md::submdspan(x_dofmap, cell, md::full_extent);
+    for (std::size_t i = 0; i < x_dofs.size(); ++i)
+      std::copy_n(&x(x_dofs[i], 0), 3, std::next(cdofs.begin(), 3 * i));
+
+    fn(&value, &coeffs(index, 0), constants.data(), cdofs.data(),
+       &local_vertex_index, nullptr, nullptr);
+  }
+
+  return value;
+}
+
 /// Assemble functional into an scalar with provided mesh geometry.
 template <dolfinx::scalar T, std::floating_point U>
 T assemble_scalar(
@@ -228,6 +265,25 @@ T assemble_scalar(
         perms);
   }
 
+  for (int i : M.integral_ids(IntegralType::vertex))
+  {
+    auto fn = M.kernel(IntegralType::vertex, i, 0);
+    assert(fn);
+
+    auto& [coeffs, cstride] = coefficients.at({IntegralType::vertex, i});
+
+    std::span<const std::int32_t> data = M.domain(IntegralType::vertex, i, 0);
+    assert(data.size() * cstride == coeffs.size());
+
+    md::mdspan<const std::int32_t,
+               md::extents<std::size_t, md::dynamic_extent, 2>>
+        cell_and_vertex(data.data(), data.size() / 2, 2);
+
+    value += impl::assemble_vertices(
+        x_dofmap, x, cell_and_vertex, fn, constants,
+        md::mdspan(coeffs.data(), data.size() / 2, cstride));
+  }
+
   return value;
 }
 

cpp/dolfinx/fem/assemble_vector_impl.h

@@ -1,4 +1,4 @@
-// Copyright (C) 2018-2021 Garth N. Wells
+// Copyright (C) 2018-2025 Garth N. Wells and Paul T. Kühner
 //
 // This file is part of DOLFINx (https://www.fenicsproject.org)
 //
@@ -907,6 +907,95 @@ void assemble_interior_facets(
   }
 }
 
+/// @brief Execute kernel over a set of vertices and accumulate result in
+/// vector.
+///
+/// @tparam T  Scalar type
+/// @tparam _bs Block size of the form test function dof map. If less
+/// than zero the block size is determined at runtime. If `_bs` is
+/// positive the block size is used as a compile-time constant, which
+/// has performance benefits.
+/// @param[in] P0 Function that applies transformation `P0.b` in-place
+/// to `b` to transform test degrees-of-freedom.
+/// @param[in,out] b Array to accumulate into.
+/// @param[in] x_dofmap Dofmap for the mesh geometry.
+/// @param[in] x Mesh geometry (coordinates).
+/// @param[in] vertices Vertex indices `(vertices.size(), 2)` - first entry
+/// holds the cell index over which the vertex is to evaluated, and the second
+/// stores the local index of the vertex within the cell.
+/// @param[in] dofmap Test function (row) degree-of-freedom data holding
+/// the (0) dofmap, (1) dofmap block size and (2) dofmap cell indices.
+/// @param[in] kernel Kernel function to execute over each cell.
+/// @param[in] constants Constant coefficient data in the kernel.
+/// @param[in] coeffs Coefficient data in the kernel. It has shape
+/// `(vertices.size(), num_cell_coeffs)`. `coeffs(i, j)` is the `j`th
+/// coefficient for cell `i`.
+/// @param[in] cell_info0 Cell permutation information for the test
+/// function mesh.
+template <dolfinx::scalar T, int _bs = -1>
+void assemble_vertices(
+    fem::DofTransformKernel<T> auto P0, std::span<T> b, mdspan2_t x_dofmap,
+    md::mdspan<const scalar_value_t<T>,
+               md::extents<std::size_t, md::dynamic_extent, 3>>
+        x,
+    md::mdspan<const std::int32_t,
+               md::extents<std::size_t, md::dynamic_extent, 2>>
+        vertices,
+    std::tuple<mdspan2_t, int,
+               md::mdspan<const std::int32_t,
+                          md::extents<std::size_t, md::dynamic_extent, 2>>>
+        dofmap,
+    FEkernel<T> auto kernel, std::span<const T> constants,
+    md::mdspan<const T, md::dextents<std::size_t, 2>> coeffs,
+    std::span<const std::uint32_t> cell_info0)
+{
+  if (vertices.empty())
+    return;
+
+  const auto [dmap, bs, vertices0] = dofmap;
+  assert(_bs < 0 or _bs == bs);
+
+  // Create data structures used in assembly
+  std::vector<scalar_value_t<T>> cdofs(3 * x_dofmap.extent(1));
+  std::vector<T> be(bs * dmap.extent(1));
+  std::span<T> _be(be);
+
+  // Iterate over active cells
+  for (std::size_t index = 0; index < vertices.extent(0); ++index)
+  {
+    // Integration domain cell, local index, and test function cell
+    std::int32_t cell = vertices(index, 0);
+    std::int32_t local_index = vertices(index, 1);
+    std::int32_t c0 = vertices0(index, 0);
+
+    // Get cell coordinates/geometry
+    auto x_dofs = md::submdspan(x_dofmap, cell, md::full_extent);
+    for (std::size_t i = 0; i < x_dofs.size(); ++i)
+      std::copy_n(&x(x_dofs[i], 0), 3, std::next(cdofs.begin(), 3 * i));
+
+    // Tabulate vector for cell
+    std::ranges::fill(be, 0);
+    kernel(be.data(), &coeffs(index, 0), constants.data(), cdofs.data(),
+           &local_index, nullptr, nullptr);
+    P0(_be, cell_info0, c0, 1);
+
+    // Scatter cell vector to 'global' vector array
+    auto dofs = md::submdspan(dmap, c0, md::full_extent);
+    if constexpr (_bs > 0)
+    {
+      for (std::size_t i = 0; i < dofs.size(); ++i)
+        for (int k = 0; k < _bs; ++k)
+          b[_bs * dofs[i] + k] += be[_bs * i + k];
+    }
+    else
+    {
+      for (std::size_t i = 0; i < dofs.size(); ++i)
+        for (int k = 0; k < bs; ++k)
+          b[bs * dofs[i] + k] += be[bs * i + k];
+    }
+  }
+}
+
 /// Modify RHS vector to account for boundary condition such that:
 ///
 /// b <- b - alpha * A.(x_bc - x0)
@@ -1340,6 +1429,32 @@ void assemble_vector(
             cell_info0, perms);
       }
     }
+
+    for (int i : L.integral_ids(IntegralType::vertex))
+    {
+      auto fn = L.kernel(IntegralType::vertex, i, 0);
+      assert(fn);
+
+      std::span vertices = L.domain(IntegralType::vertex, i, cell_type_idx);
+      std::span vertices0
+          = L.domain_arg(IntegralType::vertex, 0, i, cell_type_idx);
+
+      auto& [coeffs, cstride] = coefficients.at({IntegralType::vertex, i});
+
+      assert(vertices.size() * cstride == coeffs.size());
+
+      impl::assemble_vertices<T>(
+          P0, b, x_dofmap, x,
+          md::mdspan<const std::int32_t,
+                     md::extents<std::size_t, md::dynamic_extent, 2>>(
+              vertices.data(), vertices.size() / 2, 2),
+          {dofs, bs,
+           md::mdspan<const std::int32_t,
+                      md::extents<std::size_t, md::dynamic_extent, 2>>(
+               vertices0.data(), vertices0.size() / 2, 2)},
+          fn, constants,
+          md::mdspan(coeffs.data(), vertices.size() / 2, cstride), cell_info0);
+    }
   }
 }