Given a compact manifold $M$ in $R^n$, $M = f(x)$, f(x) is infinitely differentiable. $x$ $\in$ $R^n$, I want to find a bunch of samples on the manifold.
Currently, I'm setting up an SQP optimization problem to solve for samples on the manifold, with a set of random samples in $R^n$. The optimization goal is $min |x_r - x|$, the optimization constraint is f(x) = 0, where $x_r$ is a random sample.
Is there a better way to sample the manifold or project random points onto the manifold? I know that computation geometry works have some methods like Coxeter triangulation, which triangulates the manifold and could be used for sampling points on the manifold as well. But I'm not sure whether it will be faster than my current method since it also computes a triangulation, which I don't need.