CGAL 5.4.4 - 2D and Surface Function Interpolation
Has Model Relationships