Suppose we have a domain $\Omega$ such that
RegionDimension[[\Omega]] == RegionEmbeddingDimension[\[Omega]] == 2 is
We also have a point $P$ such that
RegionMember[RegionBoundary[\[Omega]], P] is
There is a way to find the direction of the outward normal to $\Omega$ in $P$?
I know how to mathematically obtain the outward normal if the region is described by a parametric boundary for example, but I don't know how to compute the outward normal for a generic
The method should be reasonably efficient because I need to apply to many points lying on the boundary.