Resolve #3167 (#3389)

python/dolfinx/fem/bcs.py

@@ -132,6 +132,21 @@ def function_space(self) -> dolfinx.fem.FunctionSpace:
         """The function space on which the boundary condition is defined"""
         return self._cpp_object.function_space
 
+    def dof_indices(self) -> tuple[np.ndarray, int]:
+        """Access dof indices `(local indices, unrolled)`, including ghosts, to
+        which a Dirichlet condition is applied, and the index to the first
+        non-owned (ghost) index. The array of indices is sorted.
+
+        Note:
+            The returned array is read-only.
+
+        Returns:
+            Sorted array of dof indices (unrolled) and index to the
+            first entry in the dof index array that is not owned. Entries
+            `dofs[:pos]` are owned and entries `dofs[pos:]` are ghosts.
+        """
+        return self._cpp_object.dof_indices()
+
 
 def dirichletbc(
     value: typing.Union[Function, Constant, np.ndarray],