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