High Order Cubic-Polynomial Interpolation Schemes on Triangular Meshes. Commun. Comput. Phys. 2012;12(5):1588-1602. doi:10.4208/cicp.201210.250212a