PointSource in dolfinx

[dokken]

I’ve made a version that works in parallel:

How to construct a sparse matrix used for calculating function values at multiple points

What you are searching for is how to create a point source. Following is code for getting the basis functions evaluated at a point in the cell: # Author: Jørgen S. Dokken # # Create a point source import dolfinx from mpi4py import MPI import numpy as np import ufl def compute_cell_contributions(V, points): # Determine what process owns a point and what cells it lies within _, _, owning_points, cells = dolfinx.cpp.geometry.determine_point_ownership( mesh._cpp_object, points, …

The only remaining task is inserting the contributions from a cell into the matrix/vector.