Elements of Bi-cubic Polynomial Natural Spline Interpolation for Scattered Data: Boundary Conditions Meet Partition of Unity Technique. Stat., optim. inf. comput. 2020;8(4):994-1010. doi:10.19139/soic-2310-5070-1083