Polynomial-based non-uniform interpolatory subdivision with features control

Starting from a well-known construction of polynomial-based interpolatory 4-point schemes, in this paper we present an original affine combination of quadratic polynomial samples that leads to a non-uniform 4-point scheme with edge parameters. This blending-type formulation is then further generalized to provide a powerful subdivision algorithm that combines the fairing curve of a non-uniform refinement with the advantages of a shape-controlled interpolation method and an arbitrary point insertion rule. The result is a non-uniform interpolatory 4-point scheme that is unique in combining a number of distinctive properties. In fact it generates visually-pleasing limit curves where special features ranging from cusps and flat edges to point/edge tension effects may be included without creating undesired undulations. Moreover such a scheme is capable of inserting new points at any positions of existing intervals, so that the most convenient parameter values may be chosen as well as the intervals for insertion. Such a fully flexible curve scheme is a fundamental step towards the construction of high-quality interpolatory subdivision surfaces with features control

Multigrid methods and stationary subdivisions

Multigridmethods are fast iterative solvers for sparse large ill-conditioned linear systems of equations derived, for instance, via discretization of PDEs in fluid dynamics, electrostatics and continuummechanics problems. Subdivision schemes are simple iterative algorithms for generation of smooth curves and surfaces with applications in 3D computer graphics and animation industry. This thesis presents the first definition and analysis of subdivision based multigrid methods. The main goal is to improve the convergence rate and the computational cost of multigrid taking advantage of the reproduction and regularity properties of underlying subdivision. The analysis focuses on the grid transfer operators appearing at the coarse grid correction step in the multigrid procedure. The convergence of multigrid is expressed in terms of algebraic properties of the trigonometric polynomial associated to the grid transfer operator. We interpreter the coarse-to-fine grid transfer operator as one step of subdivision. We reformulate the algebraic properties ensuring multigrid convergence in terms of regularity and generation properties of subdivision. The theoretical analysis is supported by numerical experiments for both algebraic and geometric multigrid. The numerical tests with the bivariate anisotropic Laplacian ask for subdivision schemes with anisotropic dilation. We construct a family of interpolatory subdivision schemes with such dilation which are optimal in terms of the size of the support versus their polynomial generation properties. The numerical tests confirmthe validity of our theoretical analysis