How to apply a point based load by modifying the RHS?

[dokken]

I would suggest considering:

Point sources in DOLFINx — scifem

or

Point Sources (Redux)

For dolfinx 0.8.0, simply use cmap = mesh.geometry.cmap. For dolfinx main branch (future 0.9.0), a few further changes in the call to dolfinx.cpp.geometry.determine_point_ownership are required due to Make standalone function for interpolation over non-matching meshes (… · FEniCS/dolfinx@ab79530 · GitHub. I copy the updated code below # Create a point source for Poisson problem # Author: Jørgen S. Dokken # SPDX-License-Identifier: MIT from mpi4py import MPI from petsc4py import PETSc import …