As a follow up question, linked to this post, why is locate_dofs_geometrical returning a list of two arrays ? I expected a single 1D arrray of integers pointing to dofs in a direct manner, going from 0 to self.Space.dofmap.index_map.size_global * self.Space.dofmap.index_map_bs-1 and identifying each dof in a unique manner.
Now I’m unsure my above hack works as expected, given freeinds has been computed naively from a set difference.