Regularity of multiwavelets

The motivation for this paper is an interesting observation made by Plonka concerning the factorization of the matrix symbol associated with the refinement equation for B-splines with equally spaced multiple knots at integers and subsequent developments which relate this factorization to regularity of refinable vector fields over the real line. Our intention is to contribute to this train of ideas which is partially driven by the importance of refinable vector fields in the construction of multiwavelets. The use of subdivision methods will allow us to consider the problem almost entirely in the spatial domain and leads to exact characterizations of differentiability and Hölder regularity in arbitrary L _p spaces. We first study the close relationship between vector subdivision schemes and a generalized notion of scalar subdivision schemes based on bi-infinite matrices with certain periodicity properties. For the latter type of subdivision scheme we will derive criteria for convergence and Hölder regularity of the limit function, which mainly depend on the spectral radius of a bi-infinite matrix induced by the subdivision operator, and we will show that differentiability of the limit functions can be characterized by factorization properties of the subdivision operator. By switching back to vector subdivision we will transfer these results to refinable vectors fields and obtain characterizations of regularity by factorization and spectral radius properties of the symbol associated to the refinable vector field. Finally, we point out how multiwavelets can be generated from orthonormal refinable bi-infinite vector fields.     

Analysis of Hermite subdivision using piecewise polynomials

In this paper we begin to explore a new method of analyzing the regularity of Hermite subdivision schemes that are defined from local polynomial interpolants. The idea of the method is to view the limit of the scheme as the limit of splines formed by these local interpolants rather than as the limit of polygons. We demonstrate the success of the method by obtaining the precise Hölder regularity of the simple, but non-trivial scheme in which the data are uniformly spaced and the refinement rule is defined by quintic interpolation of four values and two derivatives.     

On refinable functions and subdivision with positive masks

We consider aspects of the analysis of refinement equations with positive mask coefficients. First we derive, explicitly in terms of the mask, estimates for the geometric convergence rate of both the cascade algorithm and the corresponding subdivision scheme, as well as the Hölder continuity exponent of the resulting refinable function. Moreover, we show that the subdivision scheme converges for a class of unbounded initial sequences. Finally, we present a regularity result containing sufficient conditions on the mask for the refinable function to possess continuous derivatives up to a given order.     

Refinable bivariate quartic -splines for multi-level data representation and surface display

In this paper, a second-order Hermite basis of the space of -quartic splines on the six-directional mesh is constructed and the refinable mask of the basis functions is derived. In addition, the extra parameters of this basis are modified to extend the Hermite interpolating property at the integer lattices by including Lagrange interpolation at the half integers as well. We also formulate a compactly supported super function in terms of the basis functions to facilitate the construction of quasi-interpolants to achieve the highest (i.e., fifth) order of approximation in an efficient way. Due to the small (minimum) support of the basis functions, the refinable mask immediately yields (up to) four-point matrix-valued coefficient stencils of a vector subdivision scheme for efficient display of -quartic spline surfaces. Finally, this vector subdivision approach is further modified to reduce the size of the coefficient stencils to two-point templates while maintaining the second-order Hermite interpolating property.

     

Noninterpolatory Hermite subdivision schemes

Bivariate interpolatory Hermite subdivision schemes have recently been applied to build free-form subdivision surfaces. It is well known to geometric modelling practitioners that interpolatory schemes typically lead to ``unfair" surfaces--surfaces with unwanted wiggles or undulations--and noninterpolatory (a.k.a. approximating in the CAGD community) schemes are much preferred in geometric modelling applications. In this article, we introduce, analyze and construct noninterpolatory Hermite subdivision schemes, a class of vector subdivision schemes which can be applied to iteratively refine Hermite data in a not necessarily interpolatory fashion. We also study symmetry properties of such subdivision schemes which are crucial for application in free-form subdivision surfaces.

A key step in our mathematical analysis of Hermite type subdivision schemes is that we make use of the strong convergence theory of refinement equations to convert a prescribed geometric condition on the subdivision scheme--namely, the subdivision scheme is of Hermite type--to an algebraic condition on the subdivision mask. The latter algebraic condition can then be used in a computational framework to construct specific schemes.

     

On the regularity analysis of interpolatory Hermite subdivision schemes

It is well known that the critical Hölder regularity of a subdivision schemes can typically be expressed in terms of the joint-spectral radius (JSR) of two operators restricted to a common finite-dimensional invariant subspace. In this article, we investigate interpolatory Hermite subdivision schemes in dimension one and specifically those with optimal accuracy orders. The latter include as special cases the well-known Lagrange interpolatory subdivision schemes by Deslauriers and Dubuc. We first show how to express the critical Hölder regularity of such a scheme in terms of the joint-spectral radius of a matrix pair {F₀,F₁} given in a very explicit form. While the so-called finiteness conjecture for JSR is known to be not true in general, we conjecture that for such matrix pairs arising from Hermite interpolatory schemes of optimal accuracy orders a “strong finiteness conjecture” holds: ρ(F₀,F₁)=ρ(F₀)=ρ(F₁). We prove that this conjecture is a consequence of another conjectured property of Hermite interpolatory schemes which, in turn, is connected to a kind of positivity property of matrix polynomials. We also prove these conjectures in certain new cases using both time and frequency domain arguments; our study here strongly suggests the existence of a notion of “positive definiteness” for non-Hermitian matrices.     

Interpolatory wavelets for manifold-valued data

Geometric wavelet-like transforms for univariate and multivariate manifold-valued data can be constructed by means of nonlinear stationary subdivision rules which are intrinsic to the geometry under consideration. We show that in an appropriate vector bundle setting for a general class of interpolatory wavelet transforms, which applies to Riemannian geometry, Lie groups and other geometries, Hölder smoothness of functions is characterized by decay rates of their wavelet coefficients.     

Multivariate Compactly Supported Fundamental Refinable Functions, Duals, and Biorthogonal Wavelets

In areas of geometric modeling and wavelets, one often needs to construct a compactly supported refinable function φ which has sufficient regularity and which is fundamental for interpolation [that means, φ(0)=1 and φ(α)=0 for all α∈ Z ^s ∖{0}].
Low regularity examples of such functions have been obtained numerically by several authors, and a more general numerical scheme was given in [1]. This article presents several schemes to construct compactly supported fundamental refinable functions, which have higher regularity, directly from a given, continuous, compactly supported, refinable fundamental function φ. Asymptotic regularity analyses of the functions generated by the constructions are given.The constructions provide the basis for multivariate interpolatory subdivision algorithms that generate highly smooth surfaces.
A very important consequence of the constructions is a natural formation of pairs of dual refinable functions, a necessary element in constructing biorthogonal wavelets. Combined with the biorthogonal wavelet construction algorithm for a pair of dual refinable functions given in [2], we are able to obtain symmetrical compactly supported multivariate biorthogonal wavelets which have arbitrarily high regularity. Several examples are computed.     

Optimal two-dimensional interpolatory ternary subdivision schemes with two-ring stencils

For any interpolatory ternary subdivision scheme with two-ring stencils for a regular triangular or quadrilateral mesh, we show that the critical Hölder smoothness exponent of its basis function cannot exceed , where the critical Hölder smoothness exponent of a function is defined to be

On the other hand, for both regular triangular and quadrilateral meshes, we present several examples of interpolatory ternary subdivision schemes with two-ring stencils such that the critical Hölder smoothness exponents of their basis functions do achieve the optimal smoothness upper bound . Consequently, we obtain optimal smoothest interpolatory ternary subdivision schemes with two-ring stencils for the regular triangular and quadrilateral meshes. Our computation and analysis of optimal multidimensional subdivision schemes are based on the projection method and the -norm joint spectral radius.

     

Existence and regularity of isometries

We use local harmonic coordinates to establish sharp results on the regularity of isometric maps between Riemannian manifolds whose metric tensors have limited regularity (e.g., are Hölder continuous). We also discuss the issue of local flatness and of local isometric embedding with given first and second fundamental form, in the context of limited smoothness.

     

Palindromic Matrices of Order Two and Three-Point Subdivision Schemes

We introduce a family of three-point subdivision schemes related to palindromic pairs of matrices of order 2. We apply the Mößner theorem on palindromic matrices to the C ⁰ convergence of these subdivision schemes. We study the Hölder regularity of their limit functions. The Hölder exponent which is found in the regular case is sharp for most limit functions. In the singular case, the modulus of continuity of the limit functions is of order δlogδ. These results can be used for studying the C ¹ convergence of the Merrien family of Hermite subdivision schemes.     

From Hermite to stationary subdivision schemes in one and several variables

Vector and Hermite subdivision schemes both act on vector data, but since the latter one interprets the vectors as function values and consecutive derivatives they differ by the “renormalization” of the Hermite scheme in any step. In this paper we give an algebraic factorization method in one and several variables to relate any Hermite subdivision scheme that satisfies the so–called spectral condition to a vector subdivision scheme. These factorizations are natural extensions of the “zero at π” condition known for the masks of refinable functions. Moreover, we show how this factorization can be used to investigate different forms of convergence of the Hermite scheme and why the multivariate situation is conceptionally more intricate than the univariate one. Finally, we give some examples of such factorizations.     

Stability Implies Convergence of Cascade Algorithms in Sobolev Space

This article proves that the stability of the shifts of a refinable function vector ensures the convergence of the corresponding cascade algorithm in Sobolev space to which the refinable function vector belongs. An example of Hermite interpolants is presented to illustrate the result.     

Differentiability of multivariate refinable functions and factorization

The paper develops a necessary condition for the regularity of a multivariate refinable function in terms of a factorization property of the associated subdivision mask. The extension to arbitrary isotropic dilation matrices necessitates the introduction of the concepts of restricted and renormalized convergence of a subdivision scheme as well as the notion of subconvergence, i.e., the convergence of only a subsequence of the iterations of the subdivision scheme. Since, in addition, factorization methods pass even from scalar to matrix valued refinable functions, those results have to be formulated in terms of matrix refinable functions or vector subdivision schemes, respectively, in order to be suitable for iterated application. Moreover, it is shown for a particular case that the the condition is not only a necessary but also a sufficient one. Dedicated to Charles A. Micchelli, a unique person, friend, mathematician and collaborator, on the occasion of his sixtieth birthday Mathematics subject classifications (2000) 65T60, 65D99.     

A Construction of Biorthogonal Functions to B-Splines with Multiple Knots

We present a construction of a refinable compactly supported vector of functions which is biorthogonal to the vector of B-splines of a given degree with multiple knots at the integers with prescribed multiplicity. The construction is based on Hermite interpolatory subdivision schemes, and on the relation between B-splines and divided differences. The biorthogonal vector of functions is shown to be refinable, with a mask related to that of the Hermite scheme. For simplicity of presentation the special (scalar) case, corresponding to B-splines with simple knots, is treated separately.     

The Sobolev Regularity of Refinable Functions

Refinable functions underlie the theory and constructions of wavelet systems on the one hand and the theory and convergence analysis of uniform subdivision algorithms on the other. The regularity of such functions dictates, in the context of wavelets, the smoothness of the derived wavelet system and, in the subdivision context, the smoothness of the limiting surface of the iterative process. Since the refinable function is, in many circumstances, not known analytically, the analysis of its regularity must be based on the explicitly known mask. We establish in this paper a formula that computes, for isotropic dilation and in any number of variables, the sharp L₂-regularity of the refinable function φ in terms of the spectral radius of the restriction of the associated transfer operator to a specific invariant subspace. For a compactly supported refinable function φ, the relevant invariant space is proved to be finite dimensional and is completely characterized in terms of the dependence relations among the shifts of φ together with the polynomials that these shifts reproduce. The previously known formula for this compact support case requires the further assumptions that the mask is finitely supported and that the shifts of φ are stable. Adopting a stability assumption (but without assuming the finiteness of the mask), we derive that known formula from our general one. Moreover, we show that in the absence of stability, the lower bound provided by that previously known formula may be abysmal. Our characterization is further extended to the FSI (i.e., vector) case, to the unisotropic dilation matrix case, and to even snore general setups. We also establish corresponding results for refinable distributions.     

A construction of interpolating wavelets on invariant sets

We introduce the concept of a refinable set relative to a family of contractive mappings on a metric space, and demonstrate how such sets are useful to recursively construct interpolants which have a multiscale structure. The notion of a refinable set parallels that of a refinable function, which is the basis of wavelet construction. The interpolation points we recursively generate from a refinable set by a set-theoretic multiresolution are analogous to multiresolution for functions used in wavelet construction. We then use this recursive structure for the points to construct multiscale interpolants. Several concrete examples of refinable sets which can be used for generating interpolatory wavelets are included.

     

INTERPOLATION BY LOOP＇S SUBDIVISION FUNCTIONS

For the problem of constructing smooth functions over arbitrary surfaces from discretedata, we propose to use Loop‘s subdivision functions as the interpolants. Results on theexistence, uniqueness and error bound of the interpolants are established. An efficientprogressive computation algorithm for the interpolants is also presented.     

INTERPOLATION BY LOOP''''S SUBDIVISION FUNCTIONS

For the problem of constructing smooth functions over arbitrary surfaces from discretedata,we propose to use Loop's subdivision functions as the interpolants.Results on theexistence,uniqueness and error bound of the interpolants are established.An efficientprogressive computation algorithm for the interpolants is also presented.     

A note on the engulfing property and the -regularity of convex functions in Carnot groups

We study the engulfing property for convex functions in Carnot groups. As an application we show that the horizontal gradient of functions with this property is Hölder continuous.