Let M ⊂ R d be a compact, smooth and boundaryless manifold with dimension m and unit reach. We show how to construct a function ϕ : R d → R d−m from a uniform (ε, κ)sample P of M that offers several guarantees. Let Z ϕ denote the zero set of ϕ. Let M denote the set of points at distance ε or less from M. There exists ε 0 ∈ (0, 1) that decreases as d increases such that if ε ≤ ε 0 , the following guarantees hold. First, Z ϕ ∩ M is a faithful approximation of M in the sense that Z ϕ ∩ M is homeomorphic to M, the Hausdorff distance between Z ϕ ∩ M and M is O(m 5/2 ε 2 ), and the normal spaces at nearby points in Z ϕ ∩ M and M make an angle O(m 2 √ κε). Second, ϕ has local support; in particular, the value of ϕ at a point is affected only by sample points in P that lie within a distance of O(mε). Third, we give a projection operator that only uses sample points in P at distance O(mε) from the initial point. The projection operator maps any initial point near P onto Z ϕ ∩ M in the limit by repeated applications. * A preliminary version appears in where e g / R · ∆x converges to the zero vector as R · ∆x → 0. Since R is fixed, it means that e g / ∆x tends to the zero vector as ∆x → 0. We multiply both sides of (27) by R t and then subtract the resulting equation from (26). Some terms cancel each other because f (x + ∆x) = R t · g (R · (x + ∆x)) and f (x) = R t · g(x ) = R t · g (R · x). We obtainTherefore,We are free to choose the direction of ∆x. We choose it such that· R · ∆x , i.e., ∆x is an eigenvector corresponding to the largest eigenvalue of J f (x) − R t · J g (x ) · R. Then,