Analytic Functions in Local Shift-Invariant Spaces and Analytic Limits of Level Dependent Subdivision

Journal of Fourier Analysis and Applications - We develop a theory of Sobolev orthogonal polynomials on the Sierpiński gasket ( $$SG$$ ), which is a fractal set that can be viewed as a limit...     

Tight wavelet frames for irregular multiresolution analysis

An important tool for the construction of tight wavelet frames is the Unitary Extension Principle first formulated in the Fourier-domain by Ron and Shen. We show that the time-domain analogue of this principle provides a unified approach to the construction of tight frames based on many variations of multiresolution analyses, e.g., regular refinements of bounded L-shaped domains, refinements of subdivision surfaces around irregular vertices, and nonstationary subdivision. We consider the case of nonnegative refinement coefficients and develop a fully local construction method for tight frames. Especially, in the shift-invariant setting, our construction produces the same tight frame generators as the Unitary Extension Principle.     

Regularity of multivariate vector subdivision schemes

In this paper we discuss methods for investigating the convergence of multivariate vector subdivision schemes and the regularity of the associated limit functions. Specifically, we consider difference vector subdivision schemes whose restricted contractivity determines the convergence of the original scheme and describes the connection between the regularity of the limit functions of the difference subdivision scheme and the original subdivision scheme.     

Adaptive frame methods for nonlinear variational problems

In this paper we develop adaptive numerical solvers for certain nonlinear variational problems. The discretization of the variational problems is done by a suitable frame decomposition of the solution, i.e., a complete, stable, and redundant expansion. The discretization yields an equivalent nonlinear problem on the space of frame coefficients. The discrete problem is then adaptively solved using approximated nested fixed point and Richardson type iterations. We investigate the convergence, stability, and optimal complexity of the scheme. A theoretical advantage, for example, with respect to adaptive finite element schemes is that convergence and complexity results for the latter are usually hard to prove. The use of frames is further motivated by their redundancy, which, at least numerically, has been shown to improve the conditioning of the discretization matrices. Also frames are usually easier to construct than Riesz bases. We present a construction of divergence-free wavelet frames suitable for applications in fluid dynamics and magnetohydrodynamics. M. Fornasier acknowledges the financial support provided through the Intra-European Individual Marie Curie Fellowship Programme, under contract MOIF-CT-2006-039438. All of the authors acknowledge the hospitality of Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate, Università di Roma “La Sapienza”, Italy, during the early preparation of this work. The authors want to thank Daniele Boffi, Dorina Mitrea, and Karsten Urban for the helpful and fruitful discussions on divergence-free function spaces.     

Tight wavelet frames for subdivision

In this paper we construct multivariate tight wavelet frame decompositions for scalar and vector subdivision schemes with nonnegative masks. The constructed frame generators have one vanishing moment and are obtained by factorizing certain positive semi-definite matrices. The construction is local and allows us to obtain framelets even in the vicinity of irregular vertices. Constructing tight frames, instead of wavelet bases, we avoid extra computations of the dual masks. In addition, the frame decomposition algorithm is stable as the discrete frame transform is an isometry on _?2${?}_2$, if the data are properly normalized.     

Tight frames of compactly supported multivariate multi-wavelets

This paper is devoted to the study and construction of compactly supported tight frames of multivariate multi-wavelets. In particular, a necessary condition for their existence is derived to provide some useful guide for constructing such MRA tight frames, by reducing the factorization task of the associated polyphase matrix-valued Laurent polynomial to that of certain scalar-valued non-negative ones. We illustrate our construction method with examples of both multivariate scalar- and vector-valued subdivision schemes. Since our constructions for C¹ and C² piecewise cubic schemes are quite involved, we also include the corresponding Matlab code in the Appendix.     

Reproduction of exponential polynomials by multivariate non-stationary subdivision schemes with a general dilation matrix

We study scalar multivariate non-stationary subdivision schemes with a general integer dilation matrix. We characterize the capability of such schemes to reproduce exponential polynomials in terms of simple algebraic conditions on their symbols. These algebraic conditions provide a useful theoretical tool for checking the reproduction properties of existing schemes and for constructing new schemes with desired reproduction capabilities and other enhanced properties. We illustrate our results with several examples.     

An Algebraic Perspective on Multivariate Tight Wavelet Frames

Recent advances in real algebraic geometry and in the theory of polynomial optimization are applied to answer some open questions in the theory of multivariate tight wavelet frames whose generators have at least one vanishing moment. Namely, several equivalent formulations of the so-called Unitary Extension Principle (UEP) are given in terms of Hermitian sums of squares of certain nonnegative Laurent polynomials and in terms of semidefinite programming. These formulations merge recent advances in real algebraic geometry and wavelet frame theory and lead to an affirmative answer to the long-standing open question of the existence of tight wavelet frames in dimension d=2. They also provide, for every d, efficient numerical methods for checking the existence of tight wavelet frames and for their construction. A class of counterexamples in dimension d=3 show that, in general, the so-called sub-QMF condition is not sufficient for the existence of tight wavelet frames. Stronger sufficient conditions for determining the existence of tight wavelet frames in dimension d≥3 are derived. The results are illustrated on several examples.