Nira Dyn

Nira Dyn, Nir Sharon

Subdivision schemes have become an important tool for approximation of manifold-valued functions. In this paper, we describe a construction of manifold-valued subdivision schemes for geodesically complete manifolds. Our construction is based upon the adaptation of linear subdivision schemes using the notion of repeated binary averaging, where as a repeated binary average we propose to use the geodesic inductive mean. We derive conditions on the adapted schemes which guarantee convergence from any initial manifold-valued sequence. The definition and analysis of convergence are intrinsic to the manifold. The adaptation technique and the convergence analysis are demonstrated by several important examples.

Albert Cohen, Nira Dyn, Frédéric Hecht et al.

A simple greedy refinement procedure for the generation of data-adapted triangulations is proposed and studied. Given a function of two variables, the algorithm produces a hierarchy of triangulations and piecewise polynomial approximations on these triangulations. The refinement procedure consists in bisecting a triangle T in a direction which is chosen so as to minimize the local approximation error in some prescribed norm between the approximated function and its piecewise polynomial approximation after T is bisected. The hierarchical structure allows us to derive various approximation tools such as multiresolution analysis, wavelet bases, adaptive triangulations based either on greedy or optimal CART trees, as well as a simple encoding of the corresponding triangulations. We give a general proof of convergence in the Lp norm of all these approximations. Numerical tests performed in the case of piecewise linear approximation of functions with analytic expressions or of numerical images illustrate the fact that the refinement procedure generates triangles with an optimal aspect ratio (which is dictated by the local Hessian of of the approximated function in case of C2 functions).

Nira Dyn, Nir Sharon

A refinement of manifold data is a computational process, which produces a denser set of discrete data from a given one. Such refinements are closely related to multiresolution representations of manifold data by pyramid transforms, and approximation of manifold-valued functions by repeated refinements schemes. Most refinement methods compute each refined element separately, independently of the computations of the other elements. Here we propose a global method which computes all the refined elements simultaneously, using geodesic averages. We analyse repeated refinements schemes based on this global approach, and derive conditions guaranteeing strong convergence.

Caroline Moosmüller, Nira Dyn

In this paper we suggest a method for transforming a vector subdivision scheme generating $C^{\ell}$ limits to another such scheme of the same dimension, generating $C^{\ell+1}$ limits. In scalar subdivision, it is well known that a scheme generating $C^{\ell}$ limit curves can be transformed to a new scheme producing $C^{\ell+1}$ limit curves by multiplying the scheme's symbol with the smoothing factor $\tfrac{z+1}{2}$. We extend this approach to vector and Hermite subdivision schemes, by manipulating symbols. The algorithms presented in this paper allow to construct vector (Hermite) subdivision schemes of arbitrarily high regularity from a convergent vector scheme (from a Hermite scheme whose Taylor scheme is convergent with limit functions of vanishing first component).

Shay Kels, Nira Dyn

In this work we construct subdivision schemes refining general subsets of R^n and study their applications to the approximation of set-valued functions. Differently from previous works on set-valued approximation, our methods are developed and analyzed in the metric space of Lebesgue measurable sets endowed with the symmetric difference metric. The construction of the set-valued subdivision schemes is based on a new weighted average of two sets, which is defined for positive weights (corresponding to interpolation) and also when one weight is negative (corresponding to extrapolation). Using the new average with positive weights, we adapt to sets spline subdivision schemes computed by the Lane-Riesenfeld algorithm, which requires only averages of pairs of numbers. The averages of numbers are then replaced by the new averages of pairs of sets. Among other features of the resulting set-valued subdivision schemes, we prove their monotonicity preservation property. Using the new weighted average of sets with both positive and negative weights, we adapt to sets the 4-point interpolatory subdivision scheme. Finally we discuss the extension of the results obtained in the metric spaces of sets, to general metric spaces endowed with an averaging operation satisfying certain properties.

Nira Dyn, Xiaosheng Zhuang

Multiscale transforms for real-valued data, based on interpolatory subdivision operators have been studied in recent year. They are easy to define, and can be extended to other types of data, for example to manifold-valued data. In this paper we define linear multiscale transforms, based on certain linear, non-interpolatory subdivision operators, termed "even-reversible". For such operators, we prove, using Wiener's lemma, the existence of an inverse to the linear operator defined by the even part of the subdivision mask, and termed it "even-inverse". We show that the non-interpolatory subdivision operators, with spline or with pseudo-spline masks, are even-reversible, and derive explicitly, for the quadratic and cubic spline subdivision operators, the symbols of the corresponding even-inverse operators. We also analyze properties of the multiscale transforms based on even-reversible subdivision operators, in particular, their stability and the rate of decay of the details.

Costanza Conti, Nira Dyn, Lucia Romani

In this paper we give an elementary proof of the convergence of corner cutting algorithms refining points, in case the corner cutting weights are taken from the rather general class of weights considered by Gregory and Qu (1996). We then use similar ideas, adapted to nets of functions, to prove the convergence of corner cutting algorithms refining nets of functions, in case the corner cutting weights are taken from a stricter class of weights than in the refinement of points.

Shay Kels, Nira Dyn

We study the approximation of univariate and multivariate set-valued functions (SVFs) by the adaptation to SVFs of positive samples-based approximation operators for real-valued functions. To this end, we introduce a new weighted average of several sets and study its properties. The approximation results are obtained in the space of Lebesgue measurable sets with the symmetric difference metric. In particular, we apply the new average of sets to adapt to SVFs the classical Bernstein approximation operators, and obtain a set-valued analog of the Weierstrass approximation theorem. The rate of approximation of Hölder continuous SVFs by Bernstein operators is studied and shown to be asymptotically equal to that for real-valued functions. Finally, the results obtained in the metric space of sets are generalized to metric spaces endowed with an average satisfying certain properties.

Nira Dyn, Uri Itai, Nir Sharon

Given a set of matrices, modeled as samples of a matrix-valued function, we suggest a method to approximate the underline function using a product approximation operator. This operator extends known approximation methods by exploiting the structure of the matrices in the samples set, and based on decomposition theorems. We introduce our approach in detail and discuss its advantages using a few examples. In addition, we provide basic tools for analyzing properties of the matrix functions generated by our approximation operators.

Costanza Conti, Nira Dyn, Carla Manni et al.

A new equivalence notion between non-stationary subdivision schemes, termed asymptotical similarity, which is weaker than asymptotical equivalence, is introduced and studied. It is known that asymptotical equivalence between a non-stationary subdivision scheme and a convergent stationary scheme guarantees the convergence of the non-stationary scheme. We show that for non-stationary schemes reproducing constants, the condition of asymptotical equivalence can be relaxed to asymptotical similarity. This result applies to a wide class of non-stationary schemes of importance in theory and applications.

Nira Dyn, Maria Skopina

We study multivariate trigonometric polynomials, satisfying a set of constraints close to the known Strung-Fix conditions. Based on the polyphase representation of these polynomials relative to a general dilation matrix, we develop a simple constructive method for a special type of decomposition of such polynomials. These decompositions are of interest to the analysis of convergence and smoothness of multivariate subdivision schemes associated with general dilation matrices. We apply these decompositions, by verifying sufficient conditions for the convergence and smoothness of multivariate scalar subdivision schemes, proved here. For the convergence analysis our sufficient conditions apply to arbitrary dilation matrices, while the previously known necessary and sufficient conditions are relevant only in case of dilation matrices with a self similar tiling. For the analysis of smoothness, we state and prove two theorems on multivariate matrix subdivision schemes, which lead to sufficient conditions for C^1 limits of scalar multivariate subdivision schemes associated with isotropic dilation matrices. Although similar results are stated in the literature, we give here detailed proofs of the results, which we could not find elsewhere.