[FEATURE] - add fold constraints to finite difference interpolator

[lachlangrose]

trafficstars

Is your feature request related to a problem? Please describe. Add the folding constraints to the finite different interpolator following a similar structure to the discrete fold interpolator.

Need to work out how to regularise the gaussian curvature wrt the fold frame coordinates