1 Introduction
Subdivision schemes are iterative methods for the construction of curves and surfaces exploited in various applications, ranging from computer-aided geometric design, computer graphics and animation (see, e.g., MR2415757 ; Warren:2001:SMG:580358 ) to the construction of wavelets and frames (see, e.g., MR1968118 ; MR1971300 ). Here we focus on univariate stationary subdivision MR1079033 . Convergent subdivision schemes are characterized by a compactly supported function , called basic limit function, which satisfies a refinement equation of the type
(1) |
where is the arity, , is a compactly supported sequence of real values called mask and is a shift parameter. If a compactly supported function satisfies (1), we can always consider it to have , , changing the indexing and the shift parameter properly. In particular, denoting
it is easy to prove that has a symmetric support if and only if one of the following is true:
(2) |
Because of (2), univariate subdivision schemes are divided in two major families: primal schemes, satisfying , and dual schemes, satisfying .
By constructing the subdivision symbol
(3) |
associated to the mask , it is well-known (see, e.g., MR2775138 ; MR2843037 ; MR3071114 ) that the shift parameter satisfies
The target of this work are interpolatory schemes, i.e., schemes with a basic limit function satisfying
(4) |
Within the class of interpolatory schemes, primal ones are characterized by a simple polynomial equation which involves the sub-symbols of the scheme, i.e., the Laurent polynomials
(5) |
In particular, a convergent subdivision scheme is primal interpolatory if and only if MR1397613
(6) |
Due to (6), primal interpolatory schemes have the so called stepwise interpolation property, i.e., at each subdivision step they maintain the data of the previous one.
Beside the use of primal interpolatory schemes, there exist other approaches to interpolate points via subdivision, such as those that apply an approximating scheme after suitably preprocessing the data to be interpolated (see, e.g., Deng2010137 ; MR3921224 ; Zheng2006301 ).
However, before LUCIA , none of the existing approaches took into account the possibility of constructing a native dual interpolatory scheme which does not have the property of retaining the initial data at each iteration,
but achieves the interpolation in the sense that the initial data are still preserved in the limit function.
All dual subdivision schemes we can find in literature are indeed not interpolatory MR3071114 ,
except for the family of dual quaternary schemes introduced in LUCIA
and for the class of dual -point subdivision schemes proposed in DENG2019344 , which is shown to possess the interpolation property only when tends to infinity, i.e., when the subdivision mask has infinite length.
To fill this theoretical gap in the literature, in this paper we investigate dual interpolatory schemes with finite masks and present a complete algebraic characterization of their symbols, which covers every arity.
Additionally, for any arbitrary arity greater than
, this characterization is used as a general constructive method to produce a great amount of new interpolatory schemes with different features. In fact, the method we propose allows the user to tune a good amount of degrees of freedom: the arity
, the length of the mask , the desired degree of polynomial reproduction and some samples of the resulting basic limit function .The remainder of the paper is organized as follows. In Section 2 we state a necessary condition for a basic limit function to have prescribed values at lattices of the form , for and satisfying . In Section 3, we show an explicit algebraic characterization of dual interpolatory schemes in terms of polynomial equalities involving sub-symbols and other polynomials related to the evaluation of the basic limit function at . Exploiting this fact we are able to construct, in Section 4, new and interesting interpolatory schemes never considered in the literature, which possess comparable or even superior properties, in terms of regularity, length of the support and/or polynomial reproduction, with respect to their primal counterparts. Conclusions are drawn in Section 5.
2 Necessary Algebraic Condition for Refinability
We start proving a general necessary condition for a function to be the basic limit function of a convergent subdivision scheme, having prescribed values over the lattice , . The proof exploits the notation
(7) |
to refer to the Fourier transform of a function
, and is based on a fundamental result of harmonic analysis, i.e. the Poisson summation formula in its generalized form MR1420504 , which states:if a function , satisfies
(8) |
for some , then,
(9) |
Remark 2.1.
Condition (8) is easily satisfied by continuous compactly supported functions such as the basic limit functions of convergent subdivision schemes.
Theorem 2.2.
Let be the basic limit function of a convergent subdivision scheme of arity , compactly supported mask and sub-symbols , , satisfying (1) with . Then, for every such that , the following polynomial identity holds
(10) |
where
(11) |
with .
Proof.
It is well known that, on the Fourier side, the refinement equation (1) reads as
(12) |
where
(13) |
is the symbol associated to the refinement equation (1). Then, for every , we get
Summing over , we obtain
where we set . Since , we have
Thus,
(14) |
At this point we are ready to apply Poisson summation formula to both sides of (14) obtaining, for the left-hand-side,
(15) |
and, for the right-hand-side,
(16) |
Combining (14) with (15) and (16) we obtain
Since
we finally arrive at (10). ∎
Remark 2.3.
Remark 2.4.
In Theorem 2.2, once and are fixed, the bigger is , the more evaluations of we need to take into account. In particular, since the aim is to construct a refinable function given some of its values, it makes sense to consider
(17) |
Indeed, if there are no trigonometric polynomials satisfying (10) with as in (17) for a fixed set of values , there is no hope about the existence of solutions considering , , instead of , since .
3 Algebraic Characterization
Let us consider the dual case, i.e., in (2). Since , according to (17) we choose . In Lemma 3.5 we start specializing the necessary condition of Theorem 2.2 to the dual interpolatory case. Then we split the characterization based on the arity
being odd (Theorem
3.8) or even (Theorem 3.9).Lemma 3.5.
Consider a convergent -ary dual interpolatory subdivision scheme with compactly supported mask , sub-symbols , , and basic limit function with compact symmetric support, that satisfies , . Then
(18) |
Proof.
As a consequence of Lemma 3.5, it is easy to prove that for the arity there are no continuous refinable functions. This could be the reason why dual interpolatory schemes have not been investigated before.
Corollary 3.6.
For there are no convergent dual interpolatory subdivision schemes.
Proof.
Consider the necessary condition (18) for :
In particular, the left-hand side can be rewritten as
Then, recalling (5), we rewrite the right-hand side as
Thus, since
we get
(20) |
Now both and (and so ) are compactly supported. Then, because of (20), the first and the last non-zero elements of must be equal to . There follows that the associated subdivision scheme can not be convergent, since its difference scheme cannot be contractive CHOI2006351 ; MR1172120 ; MR2008967 . ∎
Remark 3.7.
Actually a dual interpolatory scheme was almost known since 1884. The scheme is the one with arity and mask , whose basic limit function is related to the Cantor function (see, e.g., MR2195181 ). Indeed, the resulting basic limit function , Figure 1, has and is divided into three parts: over it is exactly the well-known ascending Cantor function, over it is constant equal to and over it is equal to the descending Cantor function. This scheme however reproduces only constants and . So it is not really useful for applications. This is the shortest possible basic limit function obtained by a converging dual interpolatory scheme.

At this point, to proceed with the algebraic characterization of dual interpolatory schemes we have to split computations into two cases: the case with odd (Theorem 3.8) and the case with even (Theorem 3.9). The changes are due to the equivalences that the indices of the sums in the right-hand-side of (18) must satisfy, which depend on . We thus treat the two cases separately.
Theorem 3.8.
Let . A convergent -ary subdivision scheme is a dual interpolatory scheme if and only if
(21) |
Proof.
We start from (18) and we need to further study the right-hand-side. We observe that the only satisfying is . On the other hand, the only possibility to satisfy for is that , i.e., if , and . Thus, the expressions in the right-hand side of (18) can be rewritten as
where the last equality holds since . Thus, equation (21) follows from (18).
Conversely, if (21) holds, we can subtract it from the necessary condition for refinability (10) with and , getting, for the left-hand-side,
(22) |
and, for the right-hand-side,
(23) |
where we used the fact that . Combining (22) and (23), we obtain
(24) |
Recalling (5) and (11), one realizes that all the powers of of the second term of the right-hand-side of (24) are odd, and thus, since the left-hand side presents a unique power of which is even, it must hold that
(25) |
which implies , , and this concludes the proof. ∎
Theorem 3.9.
Let . A convergent -ary subdivision scheme is a dual interpolatory scheme if and only if
(26) |
Proof.
As for the proof of Theorem 3.8, we start from (18) and we need to further study the right-hand-side. We first observe that no satisfies , being even. On the other hand, the only possibility to satisfy for is that , namely or , , . Thus,
where the last equality holds since and . Equation (26) follows then from (18).
4 Constructive Examples
In this section we show a linear algebra approach exploiting (21) and (26) for the construction of dual interpolatory schemes with arity , and , pointing out the pros and cons with respect to known primal interpolatory schemes. Exploiting the various degrees of freedom, this method can be used to search for dual interpolatory schemes with given wanted properties such as given arity, values of the basic limit function , specific support length of the mask/basic limit function and/or degree of polynomial reproduction. The cases when the resulting linear system has no solution, are equivalent to the non-existence of a dual interpolatory scheme that satisfies the required properties.
In what follows matrices and vectors are considered to be bi-infinite where not specified and are labeled by bold uppercase and lowercase letters respectively. When useful, subsets of those matrices and vectors will be denoted using the convenient MatLab notation. Furthermore, we let the reader know that the given estimates of the regularity of the proposed schemes are obtained via joint spectral radius techniques. For the approximation of the joint spectral radius, the MatLab package
t-toolboxes, with parameter , has been used. This package implements the modified invariant polytope method by Guglielmi, Mejstrik and Protasov (see MR3886713 ; MR3009529 ; THOMAS2 ). The Hölder regularities are given with a precision of .4.1 General Strategy
The linear system to be solved arises from matching the coefficients of the Laurent polynomial in the left- and right-hand side of (21) (when is odd) or (26) (when is even). The result is described in the following Proposition.
Proposition 4.10.
Let , . Necessary condition for a bi-infinite vector to be the mask of an -ary dual interpolatory scheme with basic limit function with prescribed values
Comments
There are no comments yet.