Local lagrange interpolation with cubic splines on tetrahedral partitions

We describe an algorithm for constructing a Lagrange interpolation pair based on C¹ cubic splines defined on tetrahedral partitions. In particular, given a set of points , we construct a set P containing and a spline space based on a tetrahedral partition whose set of vertices include such that interpolation at the points of P is well-defined and unique. Earlier results are extended in two ways: (1) here we allow arbitrary sets , and (2) the method provides optimal approximation order of smooth functions.