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