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