[Paper Review] Regularity of Non-Stationary Multivariate Subdivision

Maria Charina, Costanza Conti+2 more|arXiv (Cornell University)|Jun 27, 2014

Advanced Numerical Analysis Techniques59 references3 citations

TL;DR

This paper presents a matrix-based approach to analyze the convergence and Hölder regularity of multivariate non-stationary subdivision schemes with integer dilation matrices, leveraging asymptotic similarity and approximate sum rules to extend joint spectral radius techniques to non-stationary settings. It proves a conjecture on the Hölder regularity of generalized Daubechies wavelets by establishing a sharp formula for the Hölder exponent in the case $ M = mI $, $ m \geq 2 $.

ABSTRACT

In this paper, we study scalar multivariate non-stationary subdivision schemes with integer dilation matrix M=mI, m >=2, and present a general approach for checking their convergence and for determining their Hölder regularity. The combination of the concepts of asymptotic similarity and approximate sum rules allows us to link stationary and non-stationary settings and to employ recent advances in methods for exact computation of the joint spectral radius. As an application, we prove a recent conjecture on the Hölder regularity of the generalized Daubechies wavelets. We illustrate our results with several examples.

Motivation & Objective

• Address the lack of general convergence and regularity analysis tools for non-stationary multivariate subdivision schemes, which are crucial for modeling complex geometric and exponential functions.
• Overcome the long-standing belief that joint spectral radius methods are inapplicable in non-stationary settings by establishing conditions under which they can be effectively used.
• Unify the analysis of stationary and non-stationary schemes through the concepts of asymptotic similarity and approximate sum rules.
• Provide a general framework for determining Hölder regularity of limit functions, especially in the case $ M = mI $, $ m \geq 2 $, where a closed-form formula for the Hölder exponent is derived.
• Prove a recent conjecture by N. Dyn et al. regarding the Hölder regularity of generalized Daubechies wavelets using the proposed theoretical framework.

Proposed method

• The paper introduces a matrix-based approach that reduces convergence and regularity analysis of non-stationary subdivision schemes to the computation of the joint spectral radius of a sequence of transition matrices.
• By assuming that level-dependent masks have bounded support and satisfy approximate sum rules, the method ensures the existence of a limiting subspace for the sequence of linear subspaces associated with the transition matrices.
• Approximate sum rules are shown to be 'almost necessary' for convergence and regularity, providing a bridge between stationary and non-stationary settings.
• The framework uses asymptotic similarity to relate the behavior of non-stationary schemes to their stationary counterparts, enabling the application of advanced joint spectral radius computation techniques.
• Key equations include the infinite product representation of the Fourier transform of the limit function: $ \widehat{\phi}(\omega) = \prod_{k=1}^{\infty} p_k(2^{-k}\omega) $, and the asymptotic condition $ p_k(1/2 + 2^{-k}t) = o(2^{-\ell k}) $ as $ k \to \infty $, which links smoothness to decay of derivatives at $ \omega = 1/2 $.
• Smoothness analysis relies on Taylor expansion of $ p_k(1/2 + 2^{-k}t) $ and the resulting decay condition $ \max_{j=0,\dots,\ell} 2^{-jk} |D^j p_k(1/2)| = o(2^{-\ell k}) $, which implies the Hölder regularity of the limit function.

Experimental results

Research questions

• RQ1 What conditions allow the joint spectral radius approach to be extended from stationary to non-stationary multivariate subdivision schemes?
• RQ2 How can the concepts of asymptotic similarity and approximate sum rules be used to unify the analysis of stationary and non-stationary schemes?
• RQ3 What is the precise Hölder regularity of the limit functions generated by non-stationary subdivision schemes when the dilation matrix is $ M = mI $, $ m \geq 2 $?
• RQ4 Can the proposed framework be used to verify the Hölder regularity of generalized Daubechies wavelets, as conjectured by N. Dyn et al.?
• RQ5 How do the decay properties of the symbols $ p_k $ at $ \omega = 1/2 $ relate to the smoothness of the limit function?

Key findings

• The paper establishes that the joint spectral radius technique is applicable to non-stationary subdivision schemes under mild assumptions: bounded support of masks and satisfaction of approximate sum rules.
• A sharp formula for the Hölder exponent of limit functions is derived in the case $ M = mI $, $ m \geq 2 $, enabling precise regularity estimation.
• The proposed framework proves the conjecture by N. Dyn et al. on the Hölder regularity of generalized Daubechies wavelets, confirming their smoothness properties.
• An asymptotic condition on the derivatives of the symbols $ p_k $ at $ \omega = 1/2 $, namely $ \max_{j=0,\dots,\ell} 2^{-jk} |D^j p_k(1/2)| = o(2^{-\ell k}) $, is shown to be necessary and sufficient for $ C^\ell $ regularity of the limit function.
• The existence of a limiting subspace for the sequence of linear subspaces associated with the transition matrices is rigorously proven under convergence assumptions, enabling spectral analysis.
• The method successfully links the Fourier domain behavior of the limit function to the decay of derivatives of the mask symbols, providing a practical criterion for regularity assessment.