Subdivision schemes; Vector valued; Interpolatory; Data preservation

In surface restoration usually a damaged region of a surface has to be replaced by a surface patch which restores the region in a suitable way. In particular one aims for C1-continuity at the patch boundary. The Willmore energy is considered to measure fairness and to allow appropriate boundary conditions to ensure continuity of the normal field. The corresponding L2-gradient flow as the actual restoration process leads to a system of fourth order partial differential equations, which can also be written as a system of two coupled second order equations. As it is well known, fourth order problems require an implicit time discretization. Here a semi-implicit approach is presented which allows large time steps. For the discretization of the boundary condition, two different numerical methods are introduced. Finally, we show applications to different surface restoration problems.