Design of Hermite Subdivision Schemes Aided by Spectral Radius Optimization

• Authors: Bin Han;Michael L. Overton;Thomas P.-Y. Yu
• Venue: SIAM Journal on Scientific Computing
• Year: 2003

We present a method for constructing multivariate refinable Hermite interpolants and their associated subdivision algorithms based on a combination of analytical and numerical approaches. As the limit of a linear iterative procedure, the critical L2 Sobolev smoothness of a refinable Hermite interpolant is given by the spectral radius of a matrix dependent upon the refinement mask. The design question is, Given certain constraints (support size, symmetry type, refinement pattern, etc.), how can one choose the refinement mask so that the resulting refinable function has optimal smoothness? This question naturally gives rise to a spectral radius optimization problem. In general, the objective function is not convex and may not be differentiable, or even Lipschitz, at a local minimizer. Nonetheless, a recently developed robust solver for nonsmooth optimization problems may be applied to find local minimizers of the spectral radius objective function. In fact, we find that in specific cases that are of particular interest in the present context, the objective function is smooth at local minimizers and may be accurately minimized by standard techniques. We present two necessary mathematical tricks that make the method practical: (i) compression of matrix size based on symmetry and (ii) efficient computation of gradients of the objective function. We conclude by reporting some computational results.