Tag All Mesh Cells Sufficiently-Close to a Given Point

Hi everyone,

I’d like to find a relatively quick way to tag every cell in a mesh that is within a distance \delta of a given point \mathbf{x}_0 \in \Omega.

Does anyone know of a way to do this?

You can have a look at Bitbucket
focusing in particular on obstacle subdomain (the remaing left, right etcetera are related to assigning boundary labels).

Please make sure in future post to say if you are willing to do this for dolfin or dolfinx. Having read your previous posts, I assumed dolfin.