Preserving discontinuity
We want to preserve discontinuity because they often form noticeable features
sample an additional set of points Xdisc from sharp edges of initial mesh. Compute Edisc by projecting Xdisc onto the corresponding sharp edges
disallow/penalize collapse of boundary and discontinuity edges