Facet-Based Smoothing for Mesh Generation

Starting from the triangular, linear mesh that approximates the glacier surface, a smooth surface is constructed, using the Walton-Meeks Patch, as in ​Mesh-Based Geometry and ​Smooth Facets Evaluations

In this method, surface normals are approximated at each vertex, and Bezier-type curves and surfaces are build on top of triangulated surface. Operations like closest point and normal evaluations for the surface allow using a quad mesh generation algorithm.

An iGeom type interface is implemented on the smooth evaluator. The evaluator is used by the quad mesh generator and by the solver in the coupling algorithm.