PointSource in dolfinx

[yucesoya]

Hey @dokken ,

Is there any update regarding the pointsource? I’ve also seen your comment above which you mentioned simple implementation working in serial. Could you share the copy of this? I’d greatly appreciate it.

[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.