Next Article in Journal
Is Inconsistency in the Association between Frontal Alpha Asymmetry and Depression a Function of Sex, Age, and Peripheral Inflammation?
Previous Article in Journal
Synchronized Cyclograms to Assess Inter-Limb Symmetry during Gait in Women with Anorexia and Bulimia: A Retrospective Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Simple Method for Constructing Symmetric Subdivision Schemes with High-Degree Polynomial Reproduction

1
School of Mathematics, Hefei University of Technology, Hefei 230601, China
2
School of Computer and Information, Hefei University of Technology, Hefei 230601, China
*
Author to whom correspondence should be addressed.
Submission received: 24 October 2023 / Revised: 5 December 2023 / Accepted: 7 December 2023 / Published: 15 December 2023

Abstract

:
In this paper, we present an efficient method for constructing symmetric subdivision schemes reproducing high-degree polynomials, without solving a system of linear equations. Original symmetric subdivision and its deduced subdivisions have similar increasing characteristics to the family of pseudo-splines from the polynomial reproduction point of view. Several examples are given to illustrate the efficiency of the method.

1. Introduction

Subdivision schemes are efficient methods for the generation of free-form curves and surfaces. They are easy to implement and have a wide application in computer graphics and geometric design. Given a control polygon, subdivision can generate smooth curves by recursively refining the polygon according to certain refinement rules. In general, subdivision schemes can be divided into two categories: interpolatory schemes and approximating schemes. Interpolatory schemes obtain better shape control, while approximating schemes have better smoothness. The most well-known interpolatory subdivision is the classical 4-point binary scheme proposed by Dyn et al. [1], and the Lane–Riesenfeld algorithm [2] (LR algorithm for short) provided the most representative approximating subdivision for generating arbitrary-degree, uniform B-spline curves.
An important property of subdivision schemes is polynomial reproduction. A subdivision scheme S a is said to reproduce polynomials of degree d if it is convergent and if S a f 0 = p for any polynomial p π d and initial data f i 0 = p ( t i 0 ) , where π d denotes the space of polynomials of a degree less than or equal to d [3]. This property is desirable because any convergent subdivision scheme that reproduces polynomials of degree d has the approximation order d + 1 . That is to say, if initial data is taken by sampling any function f C d + 1 with f ( d + 1 ) < uniformly with distance h, then the limit function approximates f with an error of the order O ( h d + 1 ) [4]. Using the Taylor expansion of f and the fact that the basic limit function is finitely supported, it is straightforward to show that S a has the approximation order d + 1 if it reproduces polynomials of degree d [5]. That is to say, subdivision with high-degree polynomial reproduction can approximate a given function with high accuracy. Thus, providing a lot of approximating functions, subdivision is an important and valuable method in the field of approximation theory and numerical analysis.
In the literature, Dyn et al. [6] presented a C 2 four-point subdivision scheme that reproduces cubic polynomials, based on local cubic interpolation. Choi et al. [7] constructed a family of subdivision schemes, with its high-order members achieving a higher polynomial reproduction degree. Hormann and Sabin [8] were the first to derive the polynomial reproduction degree for a family of schemes using simple algebraic considerations [5]. Later, Dyn et al. [3] generalized this algebraic method and presented necessary and sufficient conditions for primal and dual binary schemes to have polynomial reproduction of degree d. Conti and Hormann [5] further extended the results in [3] and derived a unified condition for polynomial reproduction that covers symmetric and non-symmetric schemes and naturally applies to m-ary subdivision schemes as well. Charina and Conti [9] extended the results in [5] to the multivariate setting for scalar subdivision with the dilation matrix m I ,   | m | 2 and presented algebraic conditions on the symbol for computing the exact degree of polynomial reproduction. Deng and Ma [10] present an efficient framework for the evaluation of subdivision schemes with polynomial reproduction property. In 2017, Jeong et al. [11] presented a new class of quasi-interpolatory Hermite subdivision schemes of order two, with tension parameters. In 2019, Conti et al. [12] investigated the algebraic conditions that the symbols of a univariate, binary, Hermite subdivision scheme have to fulfill in order to reproduce polynomials. And in 2020, Hüning [13] generalized the result in [12] to Hermite schemes of any order.
Our interest in polynomial reproduction is motivated by pseudo-splines ([3,14]), which are the limits of subdivision schemes with the least possible support among all schemes, with specific degrees of polynomial generation and reproduction that neatly fill the gap between uniform B-splines and interpolatory 2 n -point subdivision schemes [15]. One characteristic of pseudo-splines S a m l or S a ˜ m l is that the polynomial reproduction of its members increases progressively with l increasing. Based on this, we want to know, given a subdivision scheme, if there exists subdivision schemes interconnected with it, and their polynomial reproduction, similar to that of members of a family of pseudo-splines, progressively increases. Thus, different from the available literatures, which mainly focus on deriving the polynomial reproduction degree for given schemes using algebraic considerations on the subdivision symbol and its derivatives (see [5,9,11,12,13,16]), our study provides a direct method for deducing subdivision schemes with higher-order polynomial reproduction from a given symmetric subdivision scheme, with only the subdivision symbol being considered. Furthermore, to construct a high-performance subdivision scheme from approximating subdivision schemes, it is often necessary to solve the linear system, for instance [17], which is not easy. However, using our method, with the factorization of a difference mask for the symbol of a given scheme, the subdivision symbol can be computed directly without solving the linear system. Meanwhile, the characteristic of progressively increasing the property of polynomial reproduction contained between the original and deduced schemes is displayed clearly in the method. The main results of our paper, Propositions 1 and 2, present the computing formulas, and common formulas are listed in Table 1, Table 2, Table 3 and Table 4.
We start by recalling some basic concepts and existing results related to polynomial reproduction. In Section 3, we introduce a method for deriving subdivision schemes with higher-order accuracy from a given subdivision scheme. In Section 4, we apply the proposed method to uniform B-splines, Hormann–Sabin’s family [8], and C-schemes [18] as examples to demonstrate the simplicity and efficiency of our method. Section 5 is devoted to the conclusions and future work.

2. Preliminaries

Following the notation in [3], we consider uniform and stationary subdivision schemes S a that are determined by their masks  a = ( a i ) i Z . Define by
a ( z ) = j Z a j z j , z C 0 ,
the generating function (or Laurent polynomial) whose coefficients are the entries of a. The polynomial in (1) is commonly known as the symbol of the subdivision scheme S a .
Definition 1. 
([3]). A subdivision scheme S a is called odd symmetric if
a i = a i , i Z ,
and even symmetric if
a i = a i 1 , i Z .
In terms of Laurent polynomials, these conditions translate to a ( z ) = a ( 1 / z ) and a ( z ) z = a ( 1 / z ) , respectively.
Definition 2. 
([3]). The primal (or standard) parameterization of a subdivision scheme is based on the parameter values
t i k = i 2 k , i Z , k N 0 ,
while the dual parameterization attaches the data f i k to the parameter values
t i k = i 1 2 2 k , i Z , k N 0 .
Theorem 1. 
([3] Theorem 3.2). For a non-singular subdivision scheme S a , the condition a ( z ) is divisible by ( 1 + z ) d G + 1 , which is equivalent to the property that for any polynomial p of degree d d G , there exists some initial data f 0 such that S a f 0 = p .
Theorem 2. 
([3] Theorem 4.6). If S a is a subdivision scheme that generates polynomials of degree d G , then it reproduces polynomials of degree d R d G with respect to the primal parameterization if and only if
a ( z ) 2 is divisible by ( 1 z ) d R + 1 .
Theorem 3. 
([3] Theorem 4.7). If S a is a subdivision scheme that generates polynomials of degree d G , then it reproduces polynomials of degree d R d G with respect to the dual parameterization if and only if
a ( z 2 ) z 2 is divisible by ( 1 z ) d R + 1 .
This paper investigates symmetric subdivision schemes, and the degree of polynomial reproduction mentioned in the following corresponds to using appropriate parameterizations.

3. A Method for Constructing Symmetric Subdivision Schemes with High-Degree Polynomial Reproduction

Let S a be a symmetric approximating scheme with generation degree d G a and reproduction degree d R a , here d G a > d R a , and denote by S b the new subdivision scheme deduced from S a . We write the symbol of S b as
b ( z ) = a ( z ) T ( z ) ,
where T ( z ) is a Laurent polynomial not divisible by 1 + z . Then, by Theorem 1, S b has the same generation degree as S a . Using a necessary condition for the subdivision convergence: a ( 1 ) = b ( 1 ) = 2 , we have T ( 1 ) = 1 . Denote by P k + 1 = { p j k + 1 } j Z + , k 0 , the sequence of points at level k + 1 produced by S b from initial sequence P 0 = { p j 0 } j Z + . Using the generating function of p j k and p j k + 1 , the formal relation between points at two successive refinement levels is P k + 1 ( z ) = b ( z ) P k ( z 2 ) = a ( z ) P k ( z 2 ) T ( z ) , with P k ( z ) = j Z + p j k z j . According to the Lane–Riesenfeld algorithm [2], this means the ( k + 1 ) -level vertices produced by S b are affine combinations of the refined points P ^ k + 1 produced by S a from P k , and combination coefficients are provided by T ( z ) . Thus, a step of subdivision S b can be regarded as a step of subdivision S a followed by a smoothing stage. Our method is based on a simple observation: if we rewrite T ( z ) as T ( z ) = 1 + ρ ( z ) , then the Laurent polynomial ρ ( z ) satisfies ρ ( 1 ) = 0 , and is divisible by 1 z . At this time, we have
b ( z ) = a ( z ) + a ( z ) ρ ( z ) ,
and P k + 1 ( z ) = b ( z ) P k ( z 2 ) = a ( z ) P k ( z 2 ) + a ( z ) P k ( z 2 ) ρ ( z ) , which illustrates that, given initial sequence P 0 , the ( k + 1 ) -level vertices produced by S b are obtained via the displacement of the ( k + 1 ) -level vertices produced by S a , and the displacement is decided by ρ ( z ) . Now, using Formula (3), we give the simple method for deducing new subdivision with higher-order polynomial reproduction from symmetric approximating subdivision.
Proposition 1. 
Let S a be an odd symmetric subdivision scheme with generation degree d G a and reproduction degree d R a , d G a d R a + 2 . Then, new odd symmetric schemes S b with reproduction degree
d R b = d R a + 2 , , d G a
can be deduced from S a . The symbols of S b take the following form:
b k ( z ) = a ( z ) ( 1 + ρ k ( z ) ) ,
k = 1 , 2 , , d G a d R a 2 , and ρ k ( z ) are Laurent polynomials with the form:
ρ k ( z ) = ρ k 1 ( z ) + l k · ( 1 z ) d R a + 2 k 1 z d R a + 2 k 1 2 ,
where ρ 0 ( z ) = 0 , l k are constants depending on a 0 [ j ] , j = 1 , 2 , , N d R a 1 2 , with
a 0 [ j ] = a j ( 1 ) j [ C N j a 0 C N 2 j 1 a 0 [ 1 ] + C N 4 j 2 a 0 [ 2 ] + + ( 1 ) j 1 C N 2 j + 2 1 a 0 [ j 1 ] ] ,
j = 1 , 2 , , N d R a 1 2 , a 0 [ 0 ] = a 0 .
Specifically, when d R a = 1 ,   l k = l k 1 a 0 [ N 2 1 ] + + l 1 a 0 [ N 2 k + 1 ] + a 0 [ N 2 k ] 2 , with l 0 = 1 .
Proof of Proposition 1. 
Since S a is odd symmetric, by Definition 1, we suppose a ( z ) has the form
a ( z ) = i = 0 N a i z i z N 2 , N Z + ,
with even N, and a j = a N j , j = 0 , 1 , , N 2 1 . Using the Binomial Theorem, we obtain
i = 0 N a i z i a 0 ( 1 z ) N = i = 1 N 1 α i [ 1 ] z i = z i = 0 N 2 a i [ 1 ] z i ,
where
a i 1 [ 1 ] = α i [ 1 ] = a i ( 1 ) i C N i a 0 , i = 1 , 2 , , N 1 .
Since
a j [ 1 ] = a N 2 j [ 1 ] , j = 0 , 1 , , N 2 2 1 ,
conducting similar analysis for i = 0 N 2 a i [ 1 ] z i , we have
i = 0 N 2 a i [ 1 ] z i a 0 [ 1 ] ( 1 z ) N 2 = i = 1 N 3 α i [ 2 ] z i = z i = 0 N 4 a i [ 2 ] z i ,
where
a i 1 [ 2 ] = α i [ 2 ] = a i [ 1 ] ( 1 ) i C N 2 i a 0 [ 1 ] , i = 1 , 2 , , N 3 ,
and
a j [ 2 ] = a N 4 j [ 1 ] , j = 0 , 1 , , N 4 2 1 .
Continue the process and we obtain
a ( z ) = a 0 ( 1 z ) N + z a 0 [ 1 ] ( 1 z ) N 2 + z 2 a 0 [ 2 ] ( 1 z ) N 4 + z N 2 .
Since the reproduction degree of S a is d R a , by Theorem 2, a ( z ) 2 is divisible by ( 1 z ) d R a + 1 , then combined with (7), we know
a ( z ) = a 0 ( 1 z ) N + z a 0 [ 1 ] ( 1 z ) N 2 + z 2 a 0 [ 2 ] ( 1 z ) N 4 + + z N d R a 1 2 a 0 [ N d R a 1 2 ] ( 1 z ) d R a + 1 z N 2 + 2 ,
where
a 0 [ j ] = a j ( 1 ) j [ C N j a 0 C N 2 j 1 a 0 [ 1 ] + C N 4 j 2 a 0 [ 2 ] + + ( 1 ) j 1 C N 2 j + 2 1 a 0 [ j 1 ] ] ,
j = 1 , 2 , , N d R a 1 2 , a 0 [ 0 ] = a 0 .
Then, according to (3), we have
b ( z ) 2 = a ( z ) 2 + a ( z ) · ρ ( z ) = j = 0 N d R a 1 2 z j a 0 [ j ] ( 1 z ) N 2 j z N 2 + j = 0 N d R a 1 2 z j a 0 [ j ] ( 1 z ) N 2 j + 2 z N 2 z N 2 · ρ ( z ) .
If d R b = d R a + 2 , then by Theorem 2, b ( z ) 2 should be divisible by ( 1 z ) d R a + 3 , and ρ ( z ) should satisfy
2 z N 2 · ρ ( z ) = z N d R a 1 2 a 0 [ N d R a 1 2 ] ( 1 z ) d R a + 1 ,
thus, we take ρ ( z ) as
ρ 1 ( z ) = l 1 · ( 1 z ) d R a + 1 z d R a + 1 2 ,
where
l 1 = a 0 [ N d R a 1 2 ] 2 ,
then, corresponding to the symbol
b 1 ( z ) = a ( z ) ( 1 + ρ 1 ( z ) ) = a ( z ) ( 1 + l 1 · ( 1 z ) d R a + 1 z d R a + 1 2 ) ,
we obtain a new scheme S b 1 , which has a higher polynomial reproduction degree d R a + 2 than original scheme S a .
Similarly, to make b ( z ) 2 be divisible by ( 1 z ) d R a + 5 , ρ ( z ) should be taken as
ρ 2 ( z ) = ρ 1 ( z ) + l 2 · ( 1 z ) d R a + 3 z d R a + 3 2 = l 1 · ( 1 z ) d R a + 1 z d R a + 1 2 + l 2 · ( 1 z ) d R a + 3 z d R a + 3 2
such that
[ z N d R a 1 2 a 0 [ N d R a 1 2 ] ( 1 z ) d R a + 1 + 2 z N 2 ] · ρ ( z )
can eliminate
z N d R a 3 2 a 0 [ N d R a 3 2 ] ( 1 z ) d R a + 3 + z N d R a 1 2 a 0 [ N d R a 1 2 ] ( 1 z ) d R a + 1 .
Here, for example, if d R a = 1 , then
l 2 = l 1 a 0 [ N 2 2 ] + a 0 [ N 4 2 ] 2 ;
if d R a = 3 , then
l 2 = a 0 [ N 6 2 ] 2 .
And, corresponding to the symbol
b 2 ( z ) = a ( z ) ( 1 + ρ 2 ( z ) ) = a ( z ) ( 1 + l 1 · ( 1 z ) d R a + 1 z d R a + 1 2 + l 2 · ( 1 z ) d R a + 3 z d R a + 3 2 ) ,
the new scheme S b 2 has polynomial reproduction degree d R a + 4 .
Since the generation degree of S b is also d G a , and the reproduction degree is not more than the generation degree, d G a is the highest reproduction degree the new odd symmetric scheme S b can attain; hence, we can conclude that with the above process continuing, the last new scheme deduced from S a is S b n , n = d G a d R a 2 , and the symbol b k ( z ) , k = 1 , 2 , , d G a d R a 2 , takes the following form:
b k ( z ) = a ( z ) ( 1 + ρ k ( z ) ) ,
where
ρ k ( z ) = ρ k 1 ( z ) + l k · ( 1 z ) d R a + 2 k 1 z d R a + 2 k 1 2 ,
and l k are constants depending on a 0 [ j ] , j = 1 , 2 , , N d R a 1 2 .
Here, if d R a = 1 , then we can obtain via computation that
l k = l k 1 a 0 [ N 2 2 ] + + l 1 a 0 [ N 2 k + 2 2 ] + a 0 [ N 2 k 2 ] 2 .
Corresponding to the symbol b k ( z ) , the new scheme S b k has polynomial reproduction degree d R a + 2 k .
Therefore, we obtain new schemes S b with polynomial reproduction degree
d R b = d R a + 2 , , d G a ,
from the given odd symmetric scheme S a .
Remark 1. 
We conclude from Proposition 1 that, if S a is an odd symmetric subdivision scheme with d G a d R a + 2 , then new odd symmetric schemes S b with d R b = d R a + 2 , , d G a can be deduced. To compute the symbol of S b , b k ( z ) , k = 1 , 2 , , d G a d R a 2 , we first use formula
a 0 [ j ] = a j ( 1 ) j [ C N j a 0 C N 2 j 1 a 0 [ 1 ] + C N 4 j 2 a 0 [ 2 ] + + ( 1 ) j 1 C N 2 j + 2 1 a 0 [ j 1 ] ]
to compute a 0 [ j ] , j = 1 , 2 , , N d R a 1 2 , then l k can be computed using the method in Proposition 1, and then we have
ρ k ( z ) = ρ k 1 ( z ) + l k · ( 1 z ) d R a + 2 k 1 z d R a + 2 k 1 2
where ρ 0 ( z ) = 0 , and finally, we obtain b k ( z ) = a ( z ) ( 1 + ρ k ( z ) ) .
Remark 2. 
In Table 1 and Table 2, we list the specific computational formulas of l k and ρ k ( z ) for the cases d R a =1 and d R a =3, respectively.
Proposition 2. 
Let S a be an even symmetric subdivision scheme with generation degree d G a and reproduction degree d R a , d G a > d R a + 2 . Then, new even symmetric schemes S b with reproduction degree
d R b = d R a + 2 , , d G a 1
can be deduced from S a . The symbols of S b take the following form:
b k ( z ) = a ( z ) ( 1 + ρ k ( z ) ) ,
k = 1 , 2 , , d G a d R a 1 2 , and ρ k ( z ) are Laurent polynomials with the form:
ρ k ( z ) = ρ k 1 ( z ) + t k · ( 1 z ) d R a + 2 k 1 z d R a + 2 k 1 2 ,
where ρ 0 ( z ) = 0 , t k are constants depending on m 0 [ j ] , j = 1 , 2 , , 2 N d R a 1 2 , with
m 0 [ j ] = m j ( 1 ) j [ C 2 N j m 0 C 2 N 2 j 1 m 0 [ 1 ] + C 2 N 4 j 2 m 0 [ 2 ] + + ( 1 ) j 1 C 2 N 2 j + 2 1 m 0 [ j 1 ] ] ,
j = 1 , 2 , , 2 N d R a 1 2 , m 0 [ 0 ] = m 0 . When j is even, m j = a j 2 ; when j is odd, m j = 0 .
Proof of Proposition 2. 
Since S a is even symmetric, by Definition 1, we suppose a ( z ) has the form
a ( z ) = i = 0 N a i z i z N 2 , N Z + ,
with odd N, and a j = a N j , j = 0 , 1 , , N 2 , and
a ( z 2 ) z = i = 0 N a i z 2 i z N .
Let
m ( z ) = i = 0 2 N m i z i z N = a ( z 2 ) z ,
then m 2 i = a i , m 2 i + 1 = 0 , i = 0 , 1 , , N . Since m j = m 2 N j , j = 0 , 1 , , N 1 , m ( z ) is odd symmetric, we can obtain similar results to a ( z ) in (7) for m ( z ) . At the same time, since the reproduction degree of S a is d R a , by Theorem 3, a ( z 2 ) z 2 = m ( z ) 2 is divisible by ( 1 z ) d R a + 1 , then it can be deduced that
b ( z ) 2 = a ( z ) 2 + a ( z ) · ρ ( z ) = ( 1 z ) 4 + 8 ( 1 z ) 2 · z 8 z 2 + ( 1 z ) 4 + 8 ( 1 z ) 2 · z + 16 z 2 8 z 2 · ρ ( z ) .
m ( z ) 2 = m 0 ( 1 z ) 2 N + z m 0 [ 1 ] ( 1 z ) 2 N 2 + z 2 m 0 [ 2 ] ( 1 z ) 2 N 4 + + z 2 N d R a 1 2 m 0 [ 2 N d R a 1 2 ] ( 1 z ) d R a + 1 z N
where
m 0 [ j ] = m j ( 1 ) j [ C 2 N j m 0 C 2 N 2 j 1 m 0 [ 1 ] + C 2 N 4 j 2 m 0 [ 2 ] + + ( 1 ) j 1 C 2 N 2 j + 2 1 m 0 [ j 1 ] ] ,
j = 1 , 2 , , 2 N d R a 1 2 , m 0 [ 0 ] = m 0 .
Using (3), we have
b ( z 2 ) z 2 = a ( z 2 ) z 2 + a ( z 2 ) z · ρ ( z 2 ) = j = 0 2 N d R a 1 2 z j m 0 [ j ] ( 1 z ) 2 N 2 j z N + j = 0 2 N d R a 1 2 z j m 0 [ j ] ( 1 z ) 2 N 2 j + 2 z N z N · ρ ( z 2 ) .
If d R b = d R a + 2 , then by Theorem 3, b ( z 2 ) z 2 should be divisible by ( 1 z ) d R a + 3 , ρ ( z 2 ) should be taken as
ρ 1 ( z 2 ) = t 1 · ( 1 z 2 ) d R a + 1 z d R a + 1
such that 2 z N · ρ ( z 2 ) can eliminate z 2 N d R a 1 2 m 0 [ 2 N d R a 1 2 ] ( 1 z ) d R a + 1 . At this point, if d R a = 1 , then
ρ 1 ( z 2 ) = t 1 · ( 1 z 2 ) 2 z 2 = t 1 · ( 1 z ) 4 + 4 ( 1 z ) 2 · z z 2 ,
thus, 8 t 1 = m 0 [ N 1 ] , and
t 1 = m 0 [ N 1 ] 8 ;
If d R a = 3 , then
ρ 1 ( z 2 ) = t 1 · ( 1 z 2 ) 4 z 4 = t 1 · ( 1 z ) 8 + 8 ( 1 z ) 6 · z + 16 ( 1 z ) 4 · z 2 z 4 ,
thus, 32 t 1 = m 0 [ N 2 ] , and
t 1 = m 0 [ N 2 ] 32 ,
then, corresponding to the symbol
b 1 ( z ) = a ( z ) ( 1 + ρ 1 ( z ) ) = a ( z ) ( 1 + t 1 · ( 1 z ) d R a + 1 z d R a + 1 2 ) ,
the new scheme S b 1 has a higher polynomial reproduction degree d R a + 2 than original scheme S a .
Similarly, to make b ( z 2 ) z 2 be divisible by ( 1 z ) d R a + 5 , ρ ( z 2 ) should be taken as
ρ 2 ( z 2 ) = ρ 1 ( z 2 ) + t 2 · ( 1 z 2 ) d R a + 3 z d R a + 3 = t 1 · ( 1 z 2 ) d R a + 1 z d R a + 1 + t 2 · ( 1 z 2 ) d R a + 3 z d R a + 3
such that
[ z 2 N d R a 1 2 m 0 [ 2 N d R a 1 2 ] ( 1 z ) d R a + 1 + 2 z N ] · ρ ( z 2 )
can eliminate
z 2 N d R a 3 2 m 0 [ 2 N d R a 3 2 ] ( 1 z ) d R a + 3 + z 2 N d R a 1 2 m 0 [ 2 N d R a 1 2 ] ( 1 z ) d R a + 1 .
Here, if d R a = 1 , then
t 2 = 2 t 1 + 4 t 1 m 0 [ N 1 ] + m 0 [ N 2 ] 32 ;
if d R a = 3 , then
t 2 = 16 t 1 + m 0 [ N 3 ] 128 .
Then, corresponding to the symbol
b 2 ( z ) = a ( z ) ( 1 + ρ 2 ( z ) ) = a ( z ) ( 1 + t 1 · ( 1 z ) d R a + 1 z d R a + 1 2 + t 2 · ( 1 z ) d R a + 3 z d R a + 3 2 ) ,
S b 2 has a higher polynomial reproduction degree d R a + 4 than original scheme S a . Since the generation degree of S b is also d G a , and the reproduction degree is not more than the generation degree, d G a 1 is the highest reproduction degree the new even symmetric scheme S b can attain. Hence, we conclude that, with the above process continuing, the last new scheme deduced from S a is S b n , n = d G a d R a 1 2 , and b k ( z ) takes the following form:
b k ( z ) = a ( z ) ( 1 + ρ k ( z ) ) , ( k = 1 , 2 , , d G a d R a 1 2 ) ,
where
ρ k ( z ) = ρ k 1 ( z ) + t k · ( 1 z ) d R a + 2 k 1 z d R a + 2 k 1 2 ,
and t k are constants depending on m 0 [ j ] , j = 1 , 2 , , 2 N d R a 1 2 . Corresponding to the symbol b k ( z ) , new scheme S b k has a higher polynomial reproduction degree d R a + 2 k than original scheme S a .
Therefore, we obtain new schemes S b with the reproduction degree
d R b = d R a + 2 , , d G a 1 ,
from the given even symmetric scheme S a . □
Remark 3. 
We conclude from Proposition 2 that, if S a is an even symmetric subdivision scheme with d G a > d R a + 2 , then new even symmetric schemes S b with d R b = d R a + 2 , , d G a 1 can be deduced. To compute the symbol of S b , b k ( z ) , k = 1 , 2 , , d G a d R a 1 2 , we first use formula
m 0 [ j ] = m j ( 1 ) j [ C 2 N j m 0 C 2 N 2 j 1 m 0 [ 1 ] + C 2 N 4 j 2 m 0 [ 2 ] + + ( 1 ) j 1 C 2 N 2 j + 2 1 m 0 [ j 1 ] ]
to compute m 0 [ j ] , j = 1 , 2 , , 2 N d R a 1 2 , m 0 [ 0 ] = m 0 . When j is even, m j = a j 2 ; when j is odd, m j = 0 ; then t k can be computed using the method in Proposition 2 and then
ρ k ( z ) = ρ k 1 ( z ) + t k · ( 1 z ) d R a + 2 k 1 z d R a + 2 k 1 2 ,
where ρ 0 ( z ) = 0 , and finally, we obtain b k ( z ) = a ( z ) ( 1 + ρ k ( z ) ) .
Remark 4. 
In Table 3 and Table 4, we list the specific computational formulas of t k and ρ k ( z ) for the cases d R a = 1 and d R a = 3, respectively.

4. Applications and Examples

The aim of this section is to exploit the strategy proposed in Section 3 to derive new schemes with higher-order polynomial reproduction from the given symmetric scheme. Let us start by considering odd symmetric subdivision schemes.

4.1. Subdivision Scheme Deduced from Uniform Cubic B-Spline Refinement

The symbol of uniform cubic B-spline refinement is
a ( z ) = ( 1 + z ) 4 8 z 2 = 1 + 4 z + 6 z 2 + 4 z 3 + z 4 8 z 2 .
Here, N = 4 , ( a 0 , a 1 , a 2 , a 3 , a 4 ) = 1 8 ( 1 , 4 , 6 , 4 , 1 ) , d G a = 3, and d R a = 1. By Proposition 1, k = 1 , one scheme S b 1 with d R b 1 = 3 can be deduced from S a .
Using (9) and (8), we have ( a 0 , a 0 [ 1 ] ) = 1 8 ( 1 , 8 ) , and
a ( z ) = ( 1 z ) 4 + 8 ( 1 z ) 2 · z 8 z 2 + 2 .
Thus, if b ( z ) = a ( z ) + a ( z ) · ρ ( z ) , then
b ( z ) 2 = a ( z ) 2 + a ( z ) · ρ ( z ) = ( 1 z ) 4 + 8 ( 1 z ) 2 · z 8 z 2 + ( 1 z ) 4 + 8 ( 1 z ) 2 · z + 16 z 2 8 z 2 · ρ ( z ) .
By (10) and (11), we obtain l 1 = a 0 [ 1 ] 2 = 1 2 , and ρ 1 ( z ) = ( 1 z ) 2 2 z , thus
b 1 ( z ) = a ( z ) · ( 1 + ρ 1 ( z ) ) = ( 1 + z ) 4 8 z 2 · 1 4 z + z 2 2 z = 1 9 z 2 16 z 3 9 z 4 + z 6 16 z 3 .
Therefore, we obtain S b 1 with d R b 1 = 3 from the given scheme S a with d R a = 1 . In fact, S b 1 is precisely the classical 4-point interpolatory scheme proposed in [1].

4.2. New Subdivision Schemes S 8 _ 1 and S 8 _ 2 Deduced from Hormann–Sabin’ Family S 8

The symbol of S 8 [8] is
a ( z ) = ( 1 + z ) 8 128 z 4 · 1 + 3 z z 2 z = ( 1 + z 10 ) 5 ( z + z 9 ) 5 ( z 2 + z 8 ) + 20 ( z 3 + z 7 ) + 70 ( z 4 + z 6 ) + 98 z 5 128 z 5 .
Here, N = 10 , ( a 0 , a 1 , , a 10 ) = 1 128 ( 1 , 5 , 5 , 20 , 70 , 98 , 70 , 20 , 5 , 5 , 1 ) , d G a = 7, and d R a = 3. From Proposition 1, k = 1 , 2 , two schemes S b 1 with d R b 1 = 5 and S b 2 with d R b 2 = 7 can be deduced from S a .
Using (8) and (9), we have ( a 0 , a 0 [ 1 ] , a 0 [ 2 ] , a 0 [ 3 ] ) = 1 128 ( 1 , 15 , 80 , 160 ) , and
a ( z ) = ( 1 z ) 10 15 ( 1 z ) 8 · z 80 ( 1 z ) 6 · z 2 160 ( 1 z ) 4 · z 3 128 z 5 + 2 .
By (10) and (11), we obtain l 1 = a 0 [ 3 ] 2 = 5 8 , and ρ 1 ( z ) = 5 ( 1 z ) 4 8 z 2 , thus
b 1 ( z ) = ( 1 + z ) 8 128 z 4 · 1 + 3 z z 2 z · ( 1 + ρ 1 ( z ) ) = ( 1 + z ) 8 128 z 4 · 1 + 3 z z 2 z · 5 20 z + 38 z 2 20 z 3 + 5 z 4 8 z 2 .
And this is the first new scheme S b 1 , denoted by S 8 _ 1 , with d R b 1 = 5 deduced from S 8 with d R a = 3 .
Further, by (12) and (14), we have l 2 = a 0 [ 2 ] 2 = 5 16 , and ρ 2 ( z ) = 5 ( 1 z ) 4 8 z 2 + 5 ( 1 z ) 6 16 z 3 . Hence,
b 2 ( z ) = a ( z ) · ( 1 + ρ 2 ( z ) ) = ( 1 + z ) 8 128 z 4 · 1 + 3 z z 2 z · 5 20 z + 35 z 2 24 z 3 + 35 z 4 20 z 5 + 5 z 6 16 z 3 ,
and we obtain the second new scheme S b 2 , denoted by S 8 _ 2 , with d R b 2 = 7 from S 8 with d R a = 3 .
In Figure 1, the basic limit functions of S 8 (red), S 8 _ 1 (yellow), S 8 _ 2 (green) are illustrated.

4.3. New Subdivision Scheme Deduced from the C 2 -Scheme

In 2013, Cashman et al. [18] presented C-schemes, which generalized the Lane–Riesenfeld algorithm [2] by using the same operator of a local cubic interpolant to every four adjacent values to define each refining and smoothing stage. In this paper, for simplicity, we call the C-scheme with k smoothing steps the C k -scheme. The symbol of the C 2 -scheme is:
a ( z ) = ( 1 + z ) 6 32 z 3 · ( 1 + 4 z z 2 ) ( 1 + 10 z z 2 ) 2 128 z 3 = ( 1 + z 12 ) + 18 ( z + z 11 ) 54 ( z 2 + z 10 ) 310 ( z 3 + z 9 ) + 225 ( z 4 + z 8 ) + 2340 ( z 5 + z 7 ) + 3756 z 6 4096 z 6 ,
N = 12 ,
( a 0 , a 1 , a 2 , , a 12 )
= 1 4096 ( 1 , 18 , 54 , 310 , 225 , 2340 , 3756 , 2340 , 225 , 310 , 54 , 18 , 1 ) ,
d G a = 5, and d R a = 3 [18]. From Proposition 1, k = 1 , one scheme S b 1 with d R b 1 = 5 can be deduced from S a .
Using (9) and (8), we have
( a 0 , a 0 [ 1 ] , a 0 [ 2 ] , a 0 [ 3 ] , a 0 [ 4 ] ) = 1 4096 ( 1 , 6 , 72 , 224 , 1920 ) ,
and
a ( z ) = ( 1 z ) 12 + 6 ( 1 z ) 10 · z + 72 ( 1 z ) 8 · z 2 224 ( 1 z ) 6 · z 3 1920 ( 1 z ) 4 · z 4 4096 z 6 + 2 .
By (10) and (11), we obtain l 1 = a 0 [ 4 ] 2 = 1920 2 · 4096 = 15 64 , and ρ 1 ( z ) = 15 ( 1 z ) 4 64 z 2 , thus
b 1 ( z ) = a ( z ) · ( 1 + ρ 1 ( z ) ) = ( 1 + z ) 6 32 z 3 · ( 1 + 4 z z 2 ) ( 1 + 10 z z 2 ) 2 128 z 3 · 15 60 z + 154 z 2 60 z 3 + 15 z 4 64 z 2 ,
which is a new scheme S b 1 , denoted by C 2 _ 1 , with d R b 1 = 5 deduced from the C 2 -scheme with d R a = 3 . It should be noted that, although d G b 1 = d R b 1 = 5 , the new scheme C 2 _ 1 is not an interpolatory one. In Figure 2, the basic limit functions of C 2 (red), C 2 _ 1 (yellow) are illustrated.
Now, let us consider the examples for even symmetric subdivision schemes.

4.4. Subdivision Scheme Deduced from Uniform Quartic B-Spline Refinement

The symbol of uniform quartic B-spline refinement is
a ( z ) = ( 1 + z ) 5 16 z 3 ,
then N = 5 ,
m ( z ) = a ( z 2 ) · z = ( 1 + z 2 ) 5 16 z 5 = ( 1 + z 10 ) + 5 ( z 2 + z 8 ) + 10 ( z 4 + z 6 ) 16 z 5 ,
( m 0 , m 1 , m 2 , , m 10 ) = 1 16 ( 1 , 0 , 5 , 0 , 10 , 0 , 10 , 0 , 5 , 0 , 1 ) ,
d G a = 4, and d R a = 1. From Proposition 2, k = 1 , one scheme S b 1 with d R b 1 = 3 can be deduced from S a . Using (20) and (19), we have
( m 0 , m 0 [ 1 ] , m 0 [ 2 ] , m 0 [ 3 ] , m 0 [ 4 ] ) = 1 16 ( 1 , 10 , 40 , 80 , 80 ) ,
and
a ( z 2 ) · z = m ( z ) = ( 1 z ) 10 + 10 ( 1 z ) 8 · z + 40 ( 1 z ) 6 · z 2 + 80 ( 1 z ) 4 · z 3 + 80 ( 1 z ) 2 · z 4 16 z 5 + 2 .
Using (3), we have
b ( z 2 ) · z 2 = a ( z 2 ) · z 2 + a ( z 2 ) · z · ρ ( z 2 ) = ( 1 z ) 10 + 10 ( 1 z ) 8 · z + 40 ( 1 z ) 6 · z 2 + 80 ( 1 z ) 4 · z 3 + 80 ( 1 z ) 2 · z 4 16 z 5 + ( 1 z ) 10 + 10 ( 1 z ) 8 · z + 40 ( 1 z ) 6 · z 2 + 80 ( 1 z ) 4 · z 3 + 80 ( 1 z ) 2 · z 4 + 32 z 5 16 z 5 · ρ ( z 2 ) .
At this point, by (23) and (22), we obtain t 1 = m 0 [ 4 ] 8 = 5 8 , and ρ 1 ( z 2 ) = 5 8 · ( 1 z 2 ) 2 z 2 , thus
b 1 ( z ) = a ( z ) · ( 1 + ρ 1 ( z ) ) = ( 1 + z ) 5 16 z 3 · 1 5 ( 1 z ) 2 8 z = ( 1 + z ) 5 16 z 3 · 5 + 18 z 5 z 2 8 z .
That is, from the given scheme S a with d R a = 1 , we obtain S b 1 with d R b 1 = 3 , which is precisely the dual 4-point scheme proposed in [6].

4.5. New Subdivision Scheme S 7 _ 1 Deduced from Hormann–Sabin’s Family S 7

S 7 [8] has the symbol
a ( z ) = ( 1 + z ) 7 64 z 4 · 7 + 22 z 7 z 2 8 z = 7 ( 1 + z 9 ) 27 ( z + z 8 ) + 168 ( z 3 + z 6 ) + 378 ( z 4 + z 5 ) 512 z 5 .
Here N = 9 ,
( m 0 , m 1 , m 2 , , m 18 ) = 1 512 ( 7 , 0 , 27 , 0 , 0 , 0 , 168 , 0 , 378 , 0 , 378 , 0 , 168 , 0 , 0 , 0 , 27 , 0 , 7 ) ,
d G a = 6, and d R a = 3. From Proposition 2, k = 1 , one scheme S b 1 with d R b 1 = 5 can be deduced from S a . Using (20), we have
( m 0 , m 0 [ 1 ] , m 0 [ 2 ] , , m 0 [ 7 ] ) = 1 512 ( 7 , 126 , 972 , 4200 , 11088 , 18144 , 17472 , 8064 ) .
By (24) and (25), we obtain t 1 = m 0 [ 7 ] 32 = 63 128 , and ρ 1 ( z ) = 63 128 · ( 1 z ) 4 z 2 , thus
b 1 ( z ) = a ( z ) · ( 1 + ρ 1 ( z ) ) = ( 1 + z ) 7 64 z 4 · 7 + 22 z 7 z 2 8 z · 63 252 z + 506 z 2 252 z 3 + 63 z 4 128 z 2 .
That is, from the given scheme S a with d R a = 3 , S b 1 , denoted by S 7 _ 1 , we obtain d R b 1 = 5 .
In Figure 3, the basic limit functions of S 7 (red), S 7 _ 1 (yellow) are illustrated.

4.6. New Subdivision Scheme Deduced from the C 3 -Scheme

The symbol of the C 3 -scheme is
a ( z ) = ( 1 + z ) 7 · ( 1 + 4 z z 2 ) · ( 1 + 10 z z 2 ) 3 65536 z 8 .
Here, N = 15 ,
( m 0 , m 1 , m 2 , , m 30 ) = 1 65536 ( 1 , 0 , 27 , 0 , 207 , 0 , 13 , 0 , 3519 , 0 , 3051 , 0 , 19639 , 0 , 52299 , 0 , 52299 , 0 , 19639 , 0 , 3051 , 0 , 3519 , 0 , 13 , 0 , 207 , 0 , 27 , 0 , 1 ) ,
d G a = 6, and d R a = 3 [18]. From Proposition 2, k = 1 , one scheme S b 1 with d R b 1 = 5 can be deduced from S a .
Using (20), we have
( m 0 , m 0 [ 1 ] , , m 0 [ 13 ] ) = 1 65536 ( 1 , 30 , 378 , 2548 , 9384 , 14256 , 21920 , 128448 , 155520 , 159488 , 559104 , 276480 , 483328 , 540672 ) .
By (24) and (25), we obtain t 1 = m 0 [ 13 ] 32 = 540672 32 · 65536 = 33 128 , and ρ 1 ( z 2 ) = 33 128 · ( 1 z 2 ) 4 z 4 , thus
b 1 ( z ) = a ( z ) · ( 1 + ρ 1 ( z ) ) = ( 1 + z ) 7 · ( 1 + 4 z z 2 ) · ( 1 + 10 z z 2 ) 3 65536 z 8 · 33 132 z + 326 z 2 132 z 3 + 33 z 4 128 z 2 .
Hence, from the given scheme S a with d R a = 3 , S b 1 , denoted by C 3 _ 1 , we obtain d R b 1 = 5 . In Figure 4, the basic limit functions of C 3 (red), C 3 _ 1 (yellow) are illustrated.
Remark 5. 
If the original scheme S a is a uniform B-spline refinement, then, using the computational formulas provided in Table 1 and Table 3, it can be verified that the new schemes deduced are exactly the rest of the members of the family of pseudo-splines, which S a is included in. For example, if S a is the primal pseudo-spline S a 4 0 , with d G a = 7 and d R a = 1, then we can compute via Table 1 that S b k ,   k = 1 , 2 , 3 , are exactly pseudo-splines S a 4 l ,   l = 1 , 2 , 3 , respectively.
Remark 6. 
Note that S 8 in [8] is actually pseudo-spline S a 4 1 , but S 8 _ 1 is not pseudo-spline S a 4 2 . Although S 8 _ 1 and S a 4 2 both have degree-5 polynomial reproduction, the support of S 8 _ 1 is wider than S a 4 2 . This is because S 8 _ 1 is obtained from S 8 = S a 4 1 , while S a 4 2 is deduced directly from S a 4 0 , whose support is narrower.
In Table 5, we list comparisons of the reproduction properties of deduced subdivision schemes with the corresponding properties of original subdivision schemes included in this section. Additionally, we illustrate the polynomial reproduction capabilities (red) of S 8 _ 2 with 1-step subdivision for Legendre polynomials (blue) from degree 2 to degree 7 from initial points (yellow) in Figure 5.

5. Conclusions and Future Work

We have introduced a direct method for constructing subdivision schemes reproducing high-degree polynomials. The method is demonstrated with odd symmetric subdivision and even symmetric subdivision. The advantages of our method are mainly twofold, i.e., the constants of l k and t k can be computed rapidly and the form of the generating function of deduced subdivision is determinate. We will consider a similar method for polynomial reproduction of tensor-product and non-tensor-product schemes as our future work.

Author Contributions

Conceptualization, J.S. and J.T.; methodology, J.S. and J.T.; software, L.Z.; validation, J.T.; formal analysis, J.S. and J.T.; resources, J.S. and L.Z.; data curation, J.S. and L.Z.; writing—original draft preparation, J.S.; writing—review and editing, J.T.; visualization, J.S. and L.Z.; supervision, J.T.; project administration, J.T.; funding acquisition, J.S. and J.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, under Grant Nos. 12001151, 62172135.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dyn, N.; Levin, D.; Gregory, J.A. A 4-point interpolatory subdivision scheme for curve design. Comput. Aided Geom. Des. 1987, 4, 257–268. [Google Scholar] [CrossRef]
  2. Lane, J.M.; Riesenfeld, R.F. A theoretical development for the computer generation and display of piecewise polynomial surfaces. IEEE Trans. Pattern Anal. Mach. Intell. 1980, 2, 35–46. [Google Scholar] [CrossRef] [PubMed]
  3. Dyn, N.; Hormann, K.; Sabin, M.A.; Shen, Z. Polynomial reproduction by symmetric subdivision schemes. J. Approx. Theory 2008, 155, 28–42. [Google Scholar] [CrossRef]
  4. Levin, A. Polynomial generation and quasi-interpolation in stationary non-uniform subdivision. Comput. Aided Geom. Des. 2003, 20, 41–60. [Google Scholar] [CrossRef]
  5. Conti, C.; Hormann, K. Polynomial reproduction for univariate subdivision schemes of any arity. J. Approx. Theory 2011, 163, 413–437. [Google Scholar] [CrossRef]
  6. Dyn, N.; Floater, M.S.; Hormann, K. A C2 Four-Point Subdivision Scheme with Fourth Order Accuracy and Its Extensions. In Mathematical Methods for Curves and Surfaces: Tromsø; Nashboro Press: Brentwood, TN, USA, 2005; pp. 145–156. [Google Scholar]
  7. Choi, S.W.; Lee, B.G.; Lee, Y.J.; Yoon, J. Stationary subdivision schemes reproducing polynomials. Comput. Aided Geom. Des. 2006, 23, 351–360. [Google Scholar] [CrossRef]
  8. Hormann, K.; Sabin, M.A. A family of subdivision schemes with cubic precision. Comput. Aided Geom. Des. 2008, 25, 41–52. [Google Scholar] [CrossRef]
  9. Charina, M.; Conti, C. Polynomial reproduction of multivariate scalar subdivision schemes. J. Comput. Appl. Math. 2013, 240, 51–61. [Google Scholar] [CrossRef]
  10. Deng, C.; Ma, W. Efficient evaluation of subdivision schemes with polynomial reproduction property. J. Comput. Appl. Math. 2016, 294, 403–412. [Google Scholar] [CrossRef]
  11. Jeong, B.; Yoon, J. Construction of Hermite subdivision schemes reproducing polynomials. J. Math. Anal. Appl. 2017, 451, 565–582. [Google Scholar] [CrossRef]
  12. Conti, C.; Hüning, S. An algebraic approach to polynomial reproduction of Hermite subdivision schemes. J. Comput. Appl. Math. 2019, 349, 302–315. [Google Scholar] [CrossRef]
  13. Hüning, S. Polynomial reproduction of Hermite subdivision schemes of any order. Math. Comput. Simul. 2020, 176, 195–205. [Google Scholar] [CrossRef]
  14. Dong, B.; Shen, Z. Pseudo-splines, wavelets and framelets. Appl. Comput. Harmon. Anal. 2007, 22, 78–104. [Google Scholar] [CrossRef]
  15. Conti, C.; Deng, C.; Hormann, K. Symmetric four-directional bivariate pseudo-spline symbols. Comput. Aided Geom. Des. 2018, 60, 10–17. [Google Scholar] [CrossRef]
  16. Conti, C.; Romani, L. Algebraic conditions on non-stationary subdivision symbols for exponential polynomial reproduction. J. Comput. Appl. Math. 2011, 236, 543–556. [Google Scholar] [CrossRef]
  17. Romani, L. From approximating subdivision schemes for exponential splines to high-performance interpolating algorithms. J. Comput. Appl. Math. 2009, 224, 383–396. [Google Scholar] [CrossRef]
  18. Cashman, T.J.; Hormann, K.; Reif, U. Generalized Lane-Riesenfeld algorithms. Comput. Aided Geom. Des. 2013, 30, 398–409. [Google Scholar] [CrossRef]
Figure 1. Basic limit functions of S 8 , S 8 _ 1 , S 8 _ 2 .
Figure 1. Basic limit functions of S 8 , S 8 _ 1 , S 8 _ 2 .
Symmetry 15 02202 g001
Figure 2. Basic limit functions of C 2 , C 2 _ 1 .
Figure 2. Basic limit functions of C 2 , C 2 _ 1 .
Symmetry 15 02202 g002
Figure 3. Basic limit functions of S 7 , S 7 _ 1 .
Figure 3. Basic limit functions of S 7 , S 7 _ 1 .
Symmetry 15 02202 g003
Figure 4. Basic limit functions of C 3 , C 3 _ 1 .
Figure 4. Basic limit functions of C 3 , C 3 _ 1 .
Symmetry 15 02202 g004
Figure 5. Polynomial reproduction capabilities (red) of S 8 _ 2 for Legendre polynomials (blue) from degree 2 to degree 7 with 1-step subdivision from initial points (yellow).
Figure 5. Polynomial reproduction capabilities (red) of S 8 _ 2 for Legendre polynomials (blue) from degree 2 to degree 7 with 1-step subdivision from initial points (yellow).
Symmetry 15 02202 g005
Table 1. Specific computational formulas of l k and ρ k ( z ) and corresponding d R b k for the case d R a = 1 when S a is odd symmetric.
Table 1. Specific computational formulas of l k and ρ k ( z ) and corresponding d R b k for the case d R a = 1 when S a is odd symmetric.
k ρ k ( z ) l k d R b k
1 l 1 · ( 1 z ) 2 z a 0 [ N 2 1 ] 2 3
2 ρ 1 ( z ) + l 2 · ( 1 z ) 4 z 2 l 1 a 0 [ N 2 1 ] + a 0 [ N 2 2 ] 2 5
M = d G a 1 2 ρ M 1 ( z ) + l M · ( 1 z ) 2 M z M l M 1 a 0 [ N 2 1 ] + + l 1 a 0 [ N 2 M + 1 ] + a 0 [ N 2 M ] 2 d G a
Table 2. Specific computational formulas of l k and ρ k ( z ) and corresponding d R b k for the case d R a = 3 when S a is odd symmetric.
Table 2. Specific computational formulas of l k and ρ k ( z ) and corresponding d R b k for the case d R a = 3 when S a is odd symmetric.
k ρ k ( z ) l k d R b k
1 l 1 · ( 1 z ) 4 z 2 a 0 [ N 2 2 ] 2 5
2 ρ 1 ( z ) + l 2 · ( 1 z ) 6 z 3 a 0 [ N 2 3 ] 2 7
3 ρ 2 ( z ) + l 3 · ( 1 z ) 8 z 4 l 1 a 0 [ N 2 2 ] + a 0 [ N 2 4 ] 2 9
M = d G a 3 2 ρ M 1 ( z ) + l M · ( 1 z ) 2 M + 2 z M + 1 l M 2 a 0 [ N 2 2 ] + + l 2 a 0 [ N 2 M + 2 ] + l 1 a 0 [ N 2 M + 1 ] + a 0 [ N 2 M 1 ] 2 d G a
Table 3. Specific computational formulas of t k and ρ k ( z ) and corresponding d R b k for the case d R a = 1 when S a is even symmetric.
Table 3. Specific computational formulas of t k and ρ k ( z ) and corresponding d R b k for the case d R a = 1 when S a is even symmetric.
k ρ k ( z ) t k d R b k
1 t 1 · ( 1 z ) 2 z m 0 [ N 1 ] 8 3
2 ρ 1 ( z ) + t 2 · ( 1 z ) 4 z 2 2 t 1 + 4 t 1 m 0 [ N 1 ] + m 0 [ N 2 ] 32 5
3 ρ 2 ( z ) + t 3 · ( 1 z ) 6 z 3 16 t 2 + ( t 1 + 16 t 2 ) m 0 [ N 1 ] + 4 t 1 m 0 [ N 2 ] + m 0 [ N 3 ] 128 7
Table 4. Specific computational formulas of t k and ρ k ( z ) and corresponding d R b k for the case d R a = 3 when S a is even symmetric.
Table 4. Specific computational formulas of t k and ρ k ( z ) and corresponding d R b k for the case d R a = 3 when S a is even symmetric.
k ρ k ( z ) t k d R b k
1 t 1 · ( 1 z ) 4 z 2 m 0 [ N 2 ] 32 5
2 ρ 1 ( z ) + t 2 · ( 1 z ) 6 z 3 16 t 1 + m 0 [ N 3 ] 128 7
3 ρ 2 ( z ) + t 3 · ( 1 z ) 8 z 4 2 t 1 + 96 t 2 + 16 t 1 m 0 [ N 2 ] + m 0 [ N 4 ] 512 9
4 ρ 3 ( z ) + t 4 · ( 1 z ) 10 z 5 24 t 2 + 512 t 3 + ( 8 t 1 + 64 t 2 ) m 0 [ N 2 ] + 16 t 1 m 0 [ N 3 ] + m 0 [ N 5 ] 2048 11
Table 5. Comparison of the reproduction properties of deduced subdivision schemes with the corresponding properties of original subdivision schemes. (I—interpolatory, A—approximating).
Table 5. Comparison of the reproduction properties of deduced subdivision schemes with the corresponding properties of original subdivision schemes. (I—interpolatory, A—approximating).
Original SchemePrimal/Dual d G a d R a New Scheme Deduced ( d R b )I/A
Cubic B-splinePrimal314-point scheme (3)I
S 8 Primal73 S 8 _ 1 (5); S 8 _ 2 (7)A; A
C 2 Primal53 C 2 _ 1 (5)A
Quartic B-splineDual41Dual 4-point scheme (3)A
S 7 Dual63 S 7 _ 1 (5)A
C 3 Dual63 C 3 _ 1 (5)A
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shi, J.; Tan, J.; Zhang, L. A Simple Method for Constructing Symmetric Subdivision Schemes with High-Degree Polynomial Reproduction. Symmetry 2023, 15, 2202. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15122202

AMA Style

Shi J, Tan J, Zhang L. A Simple Method for Constructing Symmetric Subdivision Schemes with High-Degree Polynomial Reproduction. Symmetry. 2023; 15(12):2202. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15122202

Chicago/Turabian Style

Shi, Jun, Jieqing Tan, and Li Zhang. 2023. "A Simple Method for Constructing Symmetric Subdivision Schemes with High-Degree Polynomial Reproduction" Symmetry 15, no. 12: 2202. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15122202

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop