Next Article in Journal
Existence of Weak Solutions for a New Class of Fractional p-Laplacian Boundary Value Systems
Next Article in Special Issue
Shape-Preserving Properties of a Relaxed Four-Point Interpolating Subdivision Scheme
Previous Article in Journal
Magnifiers in Some Generalization of the Full Transformation Semigroups
Previous Article in Special Issue
Analysis of Geometric Properties of Ternary Four-Point Rational Interpolating Subdivision Scheme
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalized 5-Point Approximating Subdivision Scheme of Varying Arity

by
Sardar Muhammad Hussain
1,
Aziz Ur Rehman
1,
Dumitru Baleanu
2,3,4,
Kottakkaran Sooppy Nisar
5,*,
Abdul Ghaffar
6,7,* and
Samsul Ariffin Abdul Karim
8
1
Department of Mathematical Sciences, BUITEMS, Quetta 87300, Pakistan
2
Department of Mathematics, Cankaya University, 06790 Ankara, Turkey
3
Institute of Space Sciences, 077125 Magurele-Bucharest, Romania
4
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40447, Taiwan
5
Department of Mathematics, College of Arts and Sciences, Prince Sattam bin Abdulaziz University, Wadi Aldawaser 11991, Saudi Arabia
6
Informetrics Research Group, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
7
Faculty of Mathematics & Statistics, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
8
Fundamental and Applied Sciences Department and Centre for Smart Grid Energy Research (CSMER), Institute of Autonomous System, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, Seri Iskandar 32610, Perak DR, Malaysia
*
Authors to whom correspondence should be addressed.
Mathematics 2020, 8(4), 474; https://doi.org/10.3390/math8040474
Submission received: 22 February 2020 / Revised: 25 March 2020 / Accepted: 27 March 2020 / Published: 31 March 2020
(This article belongs to the Special Issue Discrete and Computational Geometry)

Abstract

:
The Subdivision Schemes (SSs) have been the heart of Computer Aided Geometric Design (CAGD) almost from its origin, and various analyses of SSs have been conducted. SSs are commonly used in CAGD and several methods have been invented to design curves/surfaces produced by SSs to applied geometry. In this article, we consider an algorithm that generates the 5-point approximating subdivision scheme with varying arity. By applying the algorithm, we further discuss several properties: continuity, Hölder regularity, limit stencils, error bound, and shape of limit curves. The efficiency of the scheme is also depicted with assuming different values of shape parameter along with its application.

1. Introduction

Computer Aided Geometric Design (CAGD) deals with studies of curves and surfaces used in computer graphics, data structure, and computational algebra. In CAGD, geometric shapes are related to the mathematical representations that satisfy approximation and interpolation properties of curves and surfaces. Surface modeling is one of the important studies in the fields of CAGD and computer graphics. It links mathematical sciences with computer science and engineering such as the animation industry, automotive and industrial design, aerospace, mechanical engineering, and numerical computing. Subdivision is an interesting subject and one of the common tools in CAGD, which provides an elegant way for the description of curves and surfaces modeling. Initially, Rham [1] worked on subdivision schemes and made a scheme which generates a function with the first derivative. Similarly, Chaikin started work and used subdivision scheme to design a curve [2]. Subdivision schemes gained importance when scientists generalized the tensor product in an arbitrary topology. Doo and Catmull used the subdivision schemes to establish surface design and control meshes in an arbitrary topology [3,4]. Deslauriers and Dubuc formed a 4-point scheme [5]. Later, Dyn et al. [6] generalized the scheme of Dubuc and Deslauriers, known as the butterfly scheme, which is based on approximated schemes. Cai used a 4-point scheme with non-uniform control points to calculate convergence and error estimation. He illustrated that the curves and surfaces generated from 4-point schemes gave better results [7]. Hassan et al. [8,9] worked on arity and number of control points, whereas Mustafa and Xuefeng [10] worked on the scheme of Bajaj with new a parameter which controls the shape of models and gave more flexibility to design a model over the soft and rough mesh network. Similarly, Siddiqui and Ahmad [11] presented a 6-point subdivision scheme that gives better smoothness. Moreover, Hormann and Sabin [12] produced a family of subdivision schemes to calculate support size, Hölder regularity, precision set, and degree of polynomial curve. Khan and Mustafa [13] calculated an interpolating 6-point subdivision scheme for complex eigenvalues as well as worked on an approximating 4-point subdivision scheme. They showed that their scheme has higher smoothness and small support size as compared to other 4-point schemes [14]. Mustafa et al. [15] worked on the m-point approximating subdivision schemes and illustrated that their schemes have higher smoothness as compared to other subdivision schemes. Siddiqi and Rehan [16] worked on a 4-point binary scheme to generate the family of curves. They introduced a scheme for C 1 continuity to generate a curve called corner cutting. Mustafa et al. [17] further worked on odd-point ternary approximating subdivision schemes and developed a formula to generalize them. Later, Ghaffar et al. [18] considered 3-point approximating subdivision schemes and observed that the given approach is more universal and is applied to schemes of arbitrary arity. Ghaffar et al. [19] introduced a general formula for 4-point a-ary approximating subdivision scheme for curve designing for any arity a 2 .
In addition, Mustafa et al. [20] worked over odd point ternary families of approximating subdivision schemes and showed that their schemes have high smoothness. They also worked on subdivision regularization, in which they showed that unified frame work can work well for both curve fitting and noise removal. They generalized unified families of interpolating subdivision schemes of 2 n -point and ( 2 n 1 ) -point p-ary which generate Lagrange’s polynomial for n 2 and p 3 , presented in [21]. In 2013, Younus and Siddiqi [22] established an algorithm based on Quaternary-point for ( m > 1 ) approximating subdivision scheme which has high smoothness and small support. Rehan et al. [23] discussed the continuity of a new class of 3-point ternary schemes and generated limiting curves using the proposed schemes. They also proposed a 4-point ternary scheme which creates C 0 interpolating and C 1 , C 2 , C 3 approximating limiting curves, described in [24]. For other recent work on this topic, we may refer to [25,26,27,28,29] and references therein.
The above-mentioned literature shows limited knowledge about the arity of the SSs. This motivated us to construct a unified 5-point approximating SS of varying arity with the shape controlling parameter. To show the performance of the schemes, we analyze the geometric properties such as continuity, Hölder regularity, and Limit stencils. Moreover, the limit curves with the specific value of shape control parameter w are depicted by the significant application of derived conditions on the initial data. The rest of the paper is organized as follows. The preliminaries regarding SSs are presented in Section 2. In Section 3, we analyze the geometric properties of the proposed schemes. The results and discussion are presented in Section 4. Some example are considered in this section to show the efficiency of the schemes. Finally, the concluding remarks are given in final section.

2. Preliminaries

In this section, we recall some well known concepts and basic results.
Definition 1.
A curve which is generated by applying a subdivision operator repeatedly to a given polygon is known as subdivision curve.
Definition 2.
If the mask of the scheme is similar for all points of the control polygon, the scheme is known as stationary.
Definition 3.
If the mask of the scheme is not similar for all points of the control polygon, the scheme is known as non-stationary.
Definition 4.
In an approximating scheme, every data point that belongs to a function generated at stage k does not belong at stage k + 1 .
Definition 5.
In an interpolating scheme, every data point belongs to a function at both stages/levels k and k + 1 .
Theorem 1.
[30] S a scheme converges iff the scheme S c is contractive. S c is contractive if c l < 1 for some l > 0 with c l = max j | c k 2 j l l | : 0 k < 2 l , where c j l are the coefficients of the scheme S c l with Laurent polynomial c l ( x ) = c ( x ) c ( x 2 ) c ( x 2 l 1 ) .
Theorem 2.
[30] If S b converges, then the limit curves of the scheme S a with Laurent polynomial a ( x ) = ( 1 + x 2 ) m b ( x ) are C m continuous, where S b is the scheme for the mth divided differences.
Theorem 3.
[30] The scheme S a with Laurent polynomial α ( x ) = 1 + x 2 m b ( x ) generates limit curves with Hölder regularity r m log 2 b l l for any l.
An a-ary scheme is said to be linear if it generates level k + 1 from level k with linear combination of control points, that is for all ‘k’ and ‘j’, there exists sets of real numbers known as masks a k = { a k } such that
λ j k + 1 = i Z a j n i k λ i k .
If the mask of the scheme is independent of k, then the scheme has finite support. Similarly, if the mask is independent of ‘j’, that is each refinement rule operates in the same way at all locations, then the scheme is known as uniform.
A general formula for the mask of the proposed scheme is defined as
α 5 a ( x ) = 1 2 a 4 1 x a 1 x 4 i = 0 4 4 i ω i
and
i = 0 4 a 8 4 i ω i = a ; where   ω j = ω 3 j   and   j = 0 , 1 , 2 ,
where w is called the shape controlling parameter and it is used to control the shape of the control of the polygon.

3. The 5-Point Approximating Schemes

This section consists of different 5-point approximating schemes together with the properties: convergence criteria, continuity, Hölder regularity, and limit stencils.

3.1. 5-Point Binary Approximating Scheme

By substituting a = 2 into Equations (1) and (2), we can get the scheme in the form
λ 2 i j + 1 = 1 256 ( 4 + 5 ω ) λ i 2 j + 84 λ i 1 j + ( 140 14 ω ) λ i j + ( 28 + 8 ω ) λ i + 1 j + ω λ i + 2 j , λ 2 i + 1 j + 1 = 1 256 ω λ i 2 j + ( 28 + 8 ω ) λ i 1 j + ( 140 14 ω ) λ i j + 84 λ i + 1 j + ( 4 + 5 ω ) λ i + 2 j ,
called 5-point binary approximating scheme.

3.1.1. Convergence Criteria

The mask of the binary scheme using Equation (3) may be written as
a = 1 256 ω , ( 4 + 5 ω ) , ( 28 + 8 ω ) , 84 , ( 140 14 ω ) , ( 140 14 ω ) , 84 , ( 28 + 8 ω ) , ( 4 + 5 ω ) , ω .
The even and odd stencil of the above scheme may be written as ( ω , ( 28 + 8 ω ) , ( 140 14 ω ) , 84 , ( 4 + 5 ω ) ) 256 and the sum of the coefficients may be written as
{ ω + ( 28 + 8 ω ) + ( 140 14 ω ) + 84 + ( 4 + 5 ω ) } 256 = 1 ,
which shows the convergence condition of the scheme.
The Laurent polynomial a ( x ) of Equation (3) takes the form
a ( x ) = 1 256 ω + ( 4 + 5 ω ) x + ( 28 + 8 ω ) x 2 + 84 x 3 + ( 140 14 ω ) x 4 + ( 140 14 ω ) x 5     + 84 x 6 + ( 28 + 8 ω ) x 7 + ( 4 + 5 ω ) x 8 + ω x 9 .
After simplification, Equation (5) becomes
a ( x ) = 1 + x 2 7 b ( x ) ,
where
b ( x ) = ω x 2 + 4 2 ω x + ω 2 .

3.1.2. Continuity

To find continuity of the scheme, Equation (6) gives the Laurent polynomial of the form
a ( x ) = 1 + x 2 b 1 ( x ) ,
where
b 1 ( x ) = ω x 8 128 + ω 32 + 1 32 x 7 + ω 32 + 3 16 x 6 + 15 32 ω 32 x 5     + 5 8 5 ω 64 x 4 + 15 32 ω 32 x 3 + ω 32 + 3 16 x 2 + ω 32 + 1 32 x + ω 128 .
C 0 continuity of the scheme S a analogous to a ( x ) , b 1 ( x ) should be convergent, where b 1 ( x ) should satisfy Theorem 1 with the given condition 1 2 S b 1 < 1 . From Theorem 1, for 9 < ω < 72 / 5 , we extract
1 2 S b 1 = 1 2 max ω 32 + 1 32 + 2 ω 32 15 32   ω 128 + ω 64 + 3 32 + 5 ω 128 5 16 < 1 .
The given condition satisfies Theorem 1, thus it must satisfy Theorem 2. This shows that the 5-point scheme is C 0 continuous.
For C 1 continuity, Equation (6) takes the form
b 2 ( x ) = 2 1 + x b 1 ( x ) ,
or
b 2 ( x ) = ω x 7 64 + 3 ω 64 + 1 16 x 6 + ω 64 + 5 16 x 5 + 5 ω 64 + 5 8 x 4     + 5 ω 64 + 5 8 x 3 + ω 64 + 5 16 x 2 + 3 ω 64 + 1 16 x + ω 64 ,
which satisfies Theorem 1 with the given condition 1 2 S b 2 < 1 . From Theorem 1, for 9 < ω < 72 / 5 , we extract
1 2 S b 2 = 1 2 max ω 64 + ω 64 + 5 16 + 5 ω 64 5 8 + 3 ω 64 + 1 16 < 1 ,
which shows that the scheme is C 1 continuous.
For C 2 continuity, Equation (6) may be written as
b 3 ( x ) = 2 1 + x b 2 ( x ) ,
or
b 3 ( x ) = 1 32 ω x 6 + ω 16 + 1 8 x 5 + ω 32 + 1 2 x 4 + ω 8 + 3 4 x 3     + ω 32 + 1 2 x 2 + ω 16 + 1 8 x + ω 32 ,
which satisfies Theorem 1 with the given condition 1 2 S b 3 < 1 . From Theorem 1, for 6 < ω < 10 , we can get
1 2 S b 3 = 1 2 max 1 16 ω + 2 ω 32 1 2 + 2 ω 16 + 1 8 + ω 8 3 4 < 1 .
Hence, the scheme is C 2 continuous.
For C 3 continuity, Equation (6) takes the form
b 4 ( x ) = 2 1 + x b 3 ( x ) ,
or
b 4 ( x ) = 1 16 ω x 5 + ω 16 + 1 4 x 4 + ω 8 + 3 4 x 3 + ω 8 + 3 4 x 2 + ω 16 + 1 4 x + ω 16 ,
which satisfies Theorem 1 with the given condition 1 2 S b 4 < 1 . From Theorem 1, for 6 < ω < 10 , we extract
1 2 S b 4 = 1 2 max 1 16 ω + ω 8 3 4 + ω 16 + 1 4 < 1 ,
which shows that the scheme is C 3 continuous.
For C 4 continuity, Equation (6) may be written as
b 5 ( x ) = 2 1 + x b 4 ( x ) ,
or
b 5 ( x ) = 1 8 ω x 4 + 1 2 x 3 + ω 4 + 1 x 2 + x 2 + ω 8 ,
with the given condition 1 2 S b 5 < 1 . Thus, from Theorem 1, for 2 < ω < 6 , we can get
1 2 S b 5 = 1 2 max 1 4 ω + ω 4 1 < 1 .
This satisfies Theorem 1, thus it must satisfy Theorem 2. Thus, the scheme is C 4 continuous.
For C 5 continuity, Equation (6) may be written as
b 6 ( x ) = 2 1 + x b 5 ( x ) ,
or
b 6 ( x ) = 1 4 ω x 3 + ω 4 + 1 x 2 + ω 4 + 1 x + ω 4 ,
with the given condition 1 2 S b 6 < 1 . From Theorem 1, for 2 < ω < 6 , we extract
1 2 S b 6 = 1 2 max 1 4 ω + ω 4 1 < 1 ,
which shows that the scheme is C 5 continuous.
For C 6 continuity, Equation (6) may be written as
b 7 ( x ) = 2 1 + x b 6 ( x ) ,
or
b 7 ( x ) = 1 2 ω x 2 + ω + 2 x + ω 2 ,
with the given condition 1 2 S b 7 < 1 . If we extract for 0 < ω < 2
1 2 S b 7 = 1 2 max ω < 1 ,
The given condition satisfies Theorem 1, thus it must satisfy Theorem 2, which shows C 6 continuity of the scheme.
Similarly, for C 7 continuity we substitute ω = 1 in Equation (6) and get
b 7 ( x ) = 1 2 ( x + 1 ) ( x + 1 ) ,
this implies that
b 8 ( x ) = 2 1 + x b 7 ( x ) .
or
b 8 ( x ) = ( x + 1 ) .
From Theorem 2, we have
1 2 S b 8 = 1 2 max 1 , 1 < 1 .
Therefore, the scheme is C 7 continuous.

3.1.3. Hölder Regularity

We use Theorem 3 to find Hölder regularity of the scheme. The Laurent polynomial of the binary scheme using Equation (3) may be written as
a ( x ) = 1 256 ω + ( 4 + 5 ω ) x + ( 28 + 8 ω ) x 2 + 84 x 3 + ( 140 14 ω ) x 4 + ( 140 14 ω ) x 5     + 84 x 6 + ( 28 + 8 ω ) x 7 + ( 4 + 5 ω ) x 8 + ω x 9 .
or
a ( x ) = 1 + x 2 7 b ( x ) ,
where
b ( x ) = ω x 2 + 4 2 ω x + ω 2 .
If ω = 1 , we can get b = 2 . Using Theorem 3 with m = 7 and b = 2 , we have r 7 log 2 ( 2 ) = 6.

3.1.4. Limit Stencils

The matrix form of the scheme using Equation (3) for w = 1 has the form
λ 3 j λ 2 j λ 1 j λ 0 j λ 1 j λ 2 j λ 3 j λ 4 j = 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 λ 3 j 1 λ 2 j 1 λ 1 j 1 λ 0 j 1 λ 1 j 1 λ 2 j 1 λ 3 j 1 λ 4 j 1
and the local subdivision matrix
X = 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 ,
which shows that the size of the invariant neighborhood is 8. After simplification, matrix X has eigenvalues λ = 1 , 1 2 , 1 4 , 1 8 , 1 16 , 1 32 , 1 64 , 1 128 and eigenvectors
1 1 2 1 4 1 8 1 16 1 32 1 64 1 128 1 1 1 1 1 1 1 1 127 967 11 23 1 7 29 363 2 7 5 469 5 7 1 33 967 3 23 1 21 47 3267 0 9 469 3 7 1 7 967 1 23 1 35 1 363 1 35 1 67 1 7 1 7 967 1 23 1 35 1 363 1 35 1 67 1 7 1 33 967 3 23 1 21 47 3267 0 9 469 3 7 1 127 967 11 23 1 7 29 363 2 7 5 469 5 7 1 1 1 1 1 1 1 1 1 .
Let
Y = 1 1 1 1 1 1 1 1 127 967 11 23 1 7 29 363 2 7 5 469 5 7 1 33 967 3 23 1 21 47 3267 0 9 469 3 7 1 7 967 1 23 1 35 1 363 1 35 1 67 1 7 1 7 967 1 23 1 35 1 363 1 35 1 67 1 7 1 33 967 3 23 1 21 47 3267 0 9 469 3 7 1 127 967 11 23 1 7 29 363 2 7 5 469 5 7 1 1 1 1 1 1 1 1 1 .
Thus, the decomposition of the local subdivision matrix X has the form X = Y Y 1 .
9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 = 1 1 1 1 1 1 1 1 127 967 11 23 1 7 29 363 2 7 5 469 5 7 1 33 967 3 23 1 21 47 3267 0 9 469 3 7 1 7 967 1 23 1 35 1 363 1 35 1 67 1 7 1 7 967 1 23 1 35 1 363 1 35 1 67 1 7 1 33 967 3 23 1 21 47 3267 0 9 469 3 7 1 127 967 11 23 1 7 29 363 2 7 5 469 5 7 1 1 1 1 1 1 1 1 1
1 0 0 0 0 0 0 0 0 1 2 0 0 0 0 0 0 0 0 1 4 0 0 0 0 0 0 0 0 1 8 0 0 0 0 0 0 0 0 1 16 0 0 0 0 0 0 0 0 1 32 0 0 0 0 0 0 0 0 1 64 0 0 0 0 0 0 0 0 1 128 967 5760 6769 5760 2901 640 18 , 373 5760 18 , 373 5760 2901 640 6769 5760 967 5760 23 2880 253 576 483 320 1127 576 1127 576 483 320 253 576 23 2880 1 8 7 8 21 8 35 8 35 8 21 8 7 8 1 8 363 1120 363 224 3267 1120 363 224 363 224 3267 1120 363 224 363 1120 7 144 161 144 7 16 665 144 665 144 7 16 161 144 7 144 469 1440 469 1440 469 160 2345 288 2345 288 469 160 469 1440 469 1440 1 1440 119 1440 119 160 245 288 245 288 119 160 119 1440 1 1440 1 40 , 320 247 40 , 320 477 4480 15 , 619 40 , 320 15 , 619 40 , 320 477 4480 247 40 , 320 1 40 , 320 .
Using diagonalization = X 1 A X , where as ∧ represents diagonal matrix and A = X X 1 . In addition, X = Y Y 1 implies that X j = Y j Y 1 and = X 1 A X implies that A j = X j X 1 .
Since the ∧ is diagonal matrix and also the power of a diagonal matrix is equal to the power of each diagonal element. Therefore,
j = ( 1 ) j 0 0 0 0 0 0 0 0 1 2 j 0 0 0 0 0 0 0 0 1 4 j 0 0 0 0 0 0 0 0 1 8 j 0 0 0 0 0 0 0 0 1 16 j 0 0 0 0 0 0 0 0 1 32 j 0 0 0 0 0 0 0 0 1 64 j 0 0 0 0 0 0 0 0 1 128 j
Since X j = Y j Y 1 , by substituting the eigen decomposition of X, we get λ j = X λ j 1 = X ( X λ j 2 ) = X 2 λ j 2 … up to λ j = X j λ 0 , or λ j = ( Y j Y 1 ) λ 0 . Taking lim j λ = lim j X j λ 0 = X λ 0 or λ = Y Y 1 λ 0 .
λ 3 j 1 λ 2 j λ 1 j λ 0 j λ 1 j λ 2 j λ 3 j λ 4 j 1 = 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 0 0 0 0 9 256 21 64 63 128 9 64 1 256 0 0 0 1 256 9 64 63 128 21 64 9 256 ( 1 ) j 0 0 0 0 0 0 0 0 1 2 j 0 0 0 0 0 0 0 0 1 4 j 0 0 0 0 0 0 0 0 1 8 j 0 0 0 0 0 0 0 0 1 16 j 0 0 0 0 0 0 0 0 1 32 j 0 0 0 0 0 0 0 0 1 64 j 0 0 0 0 0 0 0 0 1 128 j
967 5760 6769 5760 2901 640 18 , 373 5760 18 , 373 5760 2901 640 6769 5760 967 5760 23 2880 253 576 483 320 1127 576 1127 576 483 320 253 576 23 2880 1 8 7 8 21 8 35 8 35 8 21 8 7 8 1 8 363 1120 363 224 3267 1120 363 224 363 224 3267 1120 363 224 363 1120 7 144 161 144 7 16 665 144 665 144 7 16 161 144 7 144 469 1440 469 1440 469 160 2345 288 2345 288 469 160 469 1440 469 1440 1 1440 119 1440 119 160 245 288 245 288 119 160 119 1440 1 1440 1 40 , 320 247 40 , 320 477 4480 15 , 619 40 , 320 15 , 619 40 , 320 477 4480 247 40 , 320 1 40 , 320 λ 3 j 1 λ 2 j 1 λ 1 j 1 λ 0 j 1 λ 1 j 1 λ 2 j 1 λ 3 j 1 λ 4 j 1 .
After simplification, we can get
λ 3 j 1 λ 2 j λ 1 j λ 0 j λ 1 j λ 2 j λ 3 j λ 4 j 1 = 1 127 967 33 967 7 967 7 967 33 967 127 967 1 1 127 967 33 967 7 967 7 967 33 967 127 967 1 1 127 967 33 967 7 967 7 967 33 967 127 967 1 1 127 967 33 967 7 967 7 967 33 967 127 967 1 1 127 967 33 967 7 967 7 967 33 967 127 967 1 1 127 967 33 967 7 967 7 967 33 967 127 967 1 1 127 967 33 967 7 967 7 967 33 967 127 967 1 1 127 967 33 967 7 967 7 967 33 967 127 967 1 λ 3 j 1 λ 2 j 1 λ 1 j 1 λ 0 j 1 λ 1 j 1 λ 2 j 1 λ 3 j 1 λ 4 j 1 ,
which shows that the limit stencils are stable/constant.

3.2. The 5-Point Ternary Approximating Scheme

Substituting a = 3 into Equations (1) and (2), the 5-point ternary approximating scheme may be written as
λ 3 i j + 1 = 1 1296 ( 28 + 13 ω ) λ i 2 j + ( 480 4 ω ) λ i 1 j + ( 684 30 ω ) λ i j + ( 104 + 20 ω ) λ i + 1 j + ω λ i + 2 j , λ 3 i + 1 j + 1 = 1 1296 ( 4 + 5 ω ) λ i 2 j + ( 260 + 16 ω ) λ i 1 j + ( 768 42 ω ) λ i j + ( 260 + 16 ω ) λ i + 1 j + ( 4 + 5 ω ) λ i + 2 j , λ 3 i + 2 j + 1 = 1 1296 ω λ i 2 j + ( 104 + 20 ω ) λ i 1 j + ( 684 30 ω ) λ i j + ( 480 4 ω ) λ i + 1 j + ( 28 + 13 ω ) λ i + 2 j ,
β = 1 1296 ω , ( 4 + 5 ω ) , ( 28 + 13 ω ) , ( 104 + 20 ω ) , ( 260 + 16 ω ) , ( 480 4 ω ) , ( 684 30 ω ) ,     ( 768 42 ω ) , ( 684 30 ω ) , ( 480 4 ω ) , ( 260 + 16 ω ) , ( 104 + 20 ω ) , ( 28 + 13 ω ) , ( 4 + 5 ω ) , ω
and the Laurent polynomial β ( x ) of Equation (12) is
β ( x ) = 1 1296 ω + 4 + 5 ω x + 28 + 13 ω x 2 + 104 + 20 ω x 3 + 260 + 16 ω x 4 + 480 4 ω x 5     + 684 30 ω x 6 + 768 42 ω x 7 + 684 30 ω x 8 + 480 4 ω x 9 + 260 + 16 ω x 10     + 104 + 20 ω x 11 + 28 + 13 ω x 12 + 4 + 5 ω x 13 + x 14 ω
or
β ( x ) = 1 1296 1 + x + x 2 5 ( 1 + x ) 2 ω x 2 2 ω x + ω + 4 x .

3.2.1. Continuity

To find continuity of the scheme, further simplification of Equation (15) gives the Laurent polynomial in the form
β ( x ) = 1 + x + x 2 3 c 1 ( x ) ,
where
c 1 ( x ) = ω x 12 432 + ω 108 + 1 108 x 11 + ω 54 + 1 18 x 10 + ω 54 + 19 108 x 9     + 10 x 8 27 + ω 36 + 61 108 x 7 + ω 24 + 35 54 x 6 + ω 36 + 61 108 x 5     + 10 x 4 27 + ω 54 + 19 108 x 3 + ω 54 + 1 18 x 2 + ω 108 + 1 108 x + ω 432 .
To check C 0 continuity of the scheme S β analogous to β ( x ) , c 1 ( x ) should be convergent, where c 1 ( x ) should satisfy Theorem 1 with the given condition 1 3 S c 1 < 1 . From Theorem 1, for 292 9 < ω < 356 9 , we extract
1 3 S c 1 = 1 3 max ω 216 + 2 ω 54 + 19 108 + ω 24 35 54 ,   ω 108 + 1 108 + 10 27 + ω 54 + 1 18     + ω 36 61 108   ω 108 + 1 108 + 10 27 + ω 54 + 1 18 + ω 36 61 108 < 1 ,
which satisfies Theorem 1. Since the given condition satisfies Theorem 1, it must satisfy Theorem 2. This shows that the 5-point scheme is C 0 continuous.
Now, for C 1 continuity, Equation (15) may be written as
c 2 ( x ) = 3 1 + x + x 2 c 1 ( x ) ,
or
c 2 ( x ) = ω x 10 144 + ω 48 + 1 36 x 9 + ω 36 + 5 36 x 8 + ω 144 + 13 36 x 7 + 5 ω 144 + 11 18 x 6     + ω 18 + 13 18 x 5 + 5 ω 144 + 11 18 x 4 + ω 144 + 13 36 x 3 + ω 36 + 5 36 x 2 + ω 48 + 1 36 x     + ω 144 .
To check C 1 continuity, the scheme S β analogous to β ( x ) , c 2 ( x ) should be convergent, where c 2 ( x ) should satisfy Theorem 1 with the given condition 1 3 S c 2 < 1 . From Theorem 1, for 23 < ω < 31 , we extract
1 3 S c 2 = 1 3 max ω 144 + ω 144 + 13 36 + 5 ω 144 11 18 + ω 48 + 1 36 , 2 ω 36 + 5 36 + ω 18 13 18 < 1 ,
which satisfies Theorem 1. This shows that the 5-point scheme is C 1 continuous.
For C 2 continuity, Equation (15) may be written as
c 3 ( x ) = 3 1 + x + x 2 c 2 ( x ) ,
or
c 3 ( x ) = 1 48 ω x 8 + ω 24 + 1 12 x 7 + ω 48 + 1 3 x 6 + ω 24 + 2 3 x 5     + ω 12 + 5 6 x 4 + ω 24 + 2 3 x 3 + ω 48 + 1 3 x 2 + ω 24 + 1 12 x + ω 48 .
To check C 2 continuity, the scheme S β analogous to β ( x ) , c 3 ( x ) should be convergent, where c 3 ( x ) should satisfy Theorem 1 with the given condition 1 3 S c 3 < 1 . From Theorem 1, for 14 < ω < 22 , we extract
1 3 S c 3 = 1 3 max 1 48 ω + ω 24 2 3 + ω 48 + 1 3 , 2 ω 24 + 1 12 + ω 12 5 6 < 1 .
Since the given condition satisfies Theorem 1, it must satisfy Theorem 2. This means that the scheme is C 2 continuous.
For C 3 continuity, Equation (15) may be written as
c 4 ( x ) = 3 1 + x + x 2 c 3 ( x ) ,
or
c 4 ( x ) = 1 16 ω x 6 + ω 16 + 1 4 x 5 + ω 16 + 3 4 x 4 + ω 8 + 1 x 3 + ω 16 + 3 4 x 2     + ω 16 + 1 4 x + ω 16 .
To check C 3 continuity, the scheme S β analogous to β ( x ) , c 4 ( x ) should be convergent, where c 4 ( x ) should satisfy Theorem 1 with the given condition 1 3 S c 4 < 1 . From Theorem 1, for 8 < ω < 16 , we extract
1 3 S c 4 = 1 3 max 1 8 ω + ω 8 1 , ω 16 + 1 4 + ω 16 3 4 < 1 .
This shows that the scheme is C 3 continuous.
For C 4 continuity, Equation (15) may be written as
c 5 ( x ) = 3 1 + x + x 2 c 4 ( x ) ,
or
c 5 ( x ) = 3 16 ω x 4 + 3 4 x 3 + 3 8 ω + 3 2 x 2 + 3 4 x + 3 16 ω .
To check C 4 continuity, the scheme S β analogous to β ( x ) , c 5 ( x ) should be convergent, where c 5 ( x ) should satisfy Theorem 1 with the given condition 1 3 S c 5 < 1 . From Theorem 1, for 4 < ω < 12 , we extract
1 3 S c 5 = 1 3 max 1 16 ω < 3 4 , 3 8 + ω 3 2 , 1 3 < 1 ,
which satisfies Theorem 1. Thus, the scheme is C 4 continuous.
To check C 5 continuity we substitute ω = 4 3 in Equation (15) and get
c 5 ( x ) = 1 4 x 4 + 3 4 x 3 + x 2 + 3 4 x + 1 4 ,
which implies that
c 6 ( x ) = 3 1 + x + x 2 c 5 ( x ) .
or
c 6 ( x ) = 3 4 ( 1 + 2 x + x 2 ) .
Since from Theorem 2, we have
1 3 S b 8 = 1 3 max 1 4 , 2 4 < 1 .
Therefore, the scheme is C 5 continuous.

3.2.2. Hölder’s Regularity

To find Hölder regularity, we use Theorem 3. The Laurent polynomial of the ternary scheme using Equation (12) may be written as
β ( x ) = 1 1296 ω + 4 + 5 ω x + 28 + 13 ω x 2 + 104 + 20 ω x 3 + 260 + 16 ω x 4 + 480 4 ω x 5     + 684 30 ω x 6 + 768 42 ω x 7 + 684 30 ω x 8 + 480 4 ω x 9 + 260 + 16 ω x 10     + 104 + 20 ω x 11 + 28 + 13 ω x 12 + 4 + 5 ω x 13 + x 14 ω ,
β ( x ) = 1 1296 1 + x + x 2 5 c ( x ) ,
where
c ( x ) = ( 1 + x ) 2 ω x 2 2 ω x + ω + 4 x .
If ω = 4 3 , we can get b = 8 . Using Theorem 3 with m = 5 and b = 8 , we have r 5 log 2 ( 8 ) = 4.09691.

3.2.3. Limit Stencils

The matrix form of the scheme using Equation (12) has the form
λ 2 j λ 1 j λ 0 j λ 1 j λ 2 j λ 3 j = 17 486 89 243 161 324 49 486 1 972 0 2 243 211 972 89 162 211 972 2 243 0 1 972 49 486 161 324 89 243 17 486 0 0 17 486 89 243 161 324 49 486 1 972 0 2 243 211 972 89 162 211 972 2 243 0 1 972 49 486 161 324 89 243 17 486 λ 2 j 1 λ 1 j 1 λ 0 j 1 λ 1 j 1 λ 2 j 1 λ 3 j 1
and the local subdivision matrix
X = 17 486 89 243 161 324 49 486 1 972 0 2 243 211 972 89 162 211 972 2 243 0 1 972 49 486 161 324 89 243 17 486 0 0 17 486 89 243 161 324 49 486 1 972 0 2 243 211 972 89 162 211 972 2 243 0 1 972 49 486 161 324 89 243 17 486 ,
which shows the size of the invariant neighborhood is 6. After simplification, matrix X has eigenvalues λ = 1 , 1 3 , 1 9 , 1 27 , 1 81 , 1 243 , and eigenvectors.
1 1 3 1 9 1 27 1 81 1 243 1 1 1 1 1 1 27 91 27 365 661 3805 1957 24 , 283 1 3 5 5 91 23 365 313 3805 283 24 , 283 1 1 5 5 91 23 365 313 3805 283 24 , 283 1 1 5 27 91 27 365 661 3805 1957 24 , 283 1 3 5 1 1 1 1 1 1 .
After simplification on the same manner as presented in Section 3.1.4, we can get
λ 2 j λ 1 j λ 0 j λ 1 j λ 2 j λ 3 j = 1 27 91 5 91 5 91 27 91 1 1 27 91 5 91 5 91 27 91 1 1 27 91 5 91 5 91 27 91 1 1 27 91 5 91 5 91 27 91 1 1 27 91 5 91 5 91 27 91 1 1 27 91 5 91 5 91 27 91 1 λ 2 j 1 λ 1 j 1 λ 0 j 1 λ 1 j 1 λ 2 j 1 λ 3 j 1 ,
which shows that the limit stencils are stable/constant.

3.3. The 5-Point Quaternary Approximating Scheme

By substituting a = 4 into Equations (1) and (2), we can get the scheme in the form
λ 3 i j + 1 = 1 4096 ( 104 + 25 ω ) λ i 2 j + ( 1608 12 ω ) λ i 1 j + ( 2104 50 ω ) λ i j + ( 280 + 36 ω ) λ i + 1 j     + ω λ i + 2 j , λ 3 i + 1 j + 1 = 1 4096 ( 28 + 13 ω ) λ i 2 j + ( 1064 + 20 ω ) λ i 1 j + ( 2400 74 ω ) λ i j + ( 600 + 36 ω ) λ i + 1 j     + ( 4 + 5 ω ) λ i + 2 j , λ 3 i + 2 j + 1 = 1 4096 ( 4 + 5 ω ) λ i 2 j + ( 600 + 36 ω ) λ i 1 j + ( 2400 74 ω ) λ i j + ( 1064 + 20 ω ) λ i + 1 j     + ( 28 + 13 ω ) λ i + 2 j , λ 3 i + 3 j + 1 = 1 4096 ω λ i 2 j + ( 280 + 36 ω ) λ i 1 j + ( 2104 50 ω ) λ i j + ( 1608 12 ω ) λ i + 1 j     + ( 104 + 25 ω ) λ i + 2 j .
The mask of the quaternary scheme using Equation (18) takes the form
γ = 1 1296 ω , ( 4 + 5 ω ) , ( 28 + 13 ω ) , ( 104 + 13 ω ) , ( 104 + 25 ω ) , ( 280 + 36 ω ) , ( 600 + 36 ω ) ,     ( 1064 + 20 ω ) , ( 1608 12 ω ) , ( 2104 50 ω ) , ( 2400 74 ω ) , ( 2400 74 ω ) , ( 2104 50 ω ) ,     ( 1608 12 ω ) , ( 1064 + 20 ω ) , ( 600 + 36 ω ) , ( 280 + 36 ω ) , ( 104 + 25 ω ) , ( 28 + 13 ω ) ,     ( 4 + 5 ω ) , ω .
γ ( x ) = 1 4096 ω + 4 + 5 ω x + 28 + 13 ω x 2 + 104 + 25 ω x 3 + 280 + 36 ω x 4     + 600 + 36 ω x 5 + 1064 + 20 ω x 6 + 1608 12 ω x 7 + 2104 50 ω x 8     + 2400 74 ω x 9 + 2400 74 ω x 10 + 2104 50 ω x 11 + 1608 12 ω x 12     + 1064 + 20 ω x 13 + 600 + 36 ω x 14 + 280 + 36 ω x 15 + 104 + 25 ω x 16     + 28 + 13 ω x 17 + 4 + 5 ω x 18 + ω x 19 ,
or
γ ( x ) = 1 4096 x 2 + 1 5 ( 1 + x ) 7 ω x 2 2 ω x + ω + 4 x .
Hence,
γ ( x ) = 1 4096 1 + x + x 2 + x 3 5 ( 1 + x ) 2 ω x 2 2 ω x + ω + 4 x .

3.3.1. Continuity

To find continuity of the scheme, further simplification of Equation (21) gives the Laurent polynomial in the form
γ ( x ) = 1 + x + x 2 + x 3 4 e 1 ( x ) ,
where
e 1 ( x ) = ω x 16 1024 + 1 256 + ω 256 x 15 + ω 128 + 3 128 x 14 + 19 256 + 3 ω 256 x 13 + 11 64 + 3 ω 256 x 12     + 81 256 + ω 256 x 11 + 61 128 ω 128 x 10 + 155 256 5 ω 256 x 9 + 21 32 13 ω 512 x 8     + 155 256 5 ω 256 x 7 + 61 128 ω 128 x 6 + 81 256 + ω 256 x 5 + 11 64 + 3 ω 256 x 4     + 19 256 + 3 ω 256 x 3 + ω 128 + 3 128 x 2 + 1 256 + ω 256 x + ω 1024 .
To check C 0 continuity, the scheme S γ analogous to γ ( x ) , e 1 ( x ) should be convergent, where e 1 ( x ) should satisfy Theorem 1 with the given condition 1 4 S e 1 < 1 . From Theorem 1, for 3775 51 < ω < 4415 51 , we extract
1 4 S e 1 = 1 4 max ω 1024 + 2 3 ω 256 + 11 64 + 13 ω 512 21 32 , 1 256 + ω 256 + ω 256 + 81 256     + 5 ω 256 155 256 + 19 256 + 3 ω 256 , 2 3 128 + ω 128 + 2 ω 128 61 128 < 1 ,
which satisfies Theorem 1. Since the given condition satisfies Theorem 1, it must satisfy Theorem 2. This shows that 5-point scheme is C 0 continuous.
Now, for C 1 continuity, Equation (21) may be written as
e 2 ( x ) ) = 1 + x + x 2 + x 3 4 e 1 ( x ) ,
or
e 2 ( x ) = ω x 13 256 + 1 64 + 3 ω 256 x 12 + 5 64 + ω 64 x 11 + ω 64 + 13 64 x 10 + 25 64 + ω 256 x 9     + 19 32 5 ω 256 x 8 + 23 32 ω 32 x 7 + 23 32 ω 32 x 6 + 19 32 5 ω 256 x 5     + 25 64 + ω 256 x 4 + ω 64 + 13 64 x 3 + 5 64 + ω 64 x 2 + 1 64 + 3 ω 256 x + ω 256 .
To check C 1 continuity, the scheme S γ analogous to γ ( x ) , e 2 ( x ) should be convergent, where e 2 ( x ) should satisfy Theorem 1 with the given condition 1 4 S e 2 < 1 . From Theorem 1, for 57 < ω < 71 , we extract
1 4 S e 1 = 1 4 max 1 64 + 3 ω 256 + 19 32 + 5 ω 256 + 25 64 + ω 256 + ω 256 , 5 64 + ω 64     + 23 32 + ω 32 + ω 64 + 13 64 < 1 .
Since the given condition satisfies Theorem 1, it must satisfy Theorem 2. This means that the scheme is C 1 continuous.
For C 2 continuity, Equation (21) may be written as
e 3 ( x ) = 1 + x + x 2 + x 3 4 e 2 ( x ) ,
or
e 3 ( x ) = ω x 10 64 + ω 32 + 1 16 x 9 + ω 64 + 1 4 x 8 + 1 2 x 7 + ω 32 + 3 4 x 6 + ω 16 + 7 8 x 5     + ω 32 + 3 4 x 4 + 1 2 x 3 + ω 64 + 1 4 x 2 + ω 32 + 1 16 x + ω 64 .
To check C 2 continuity, the scheme S γ analogous to γ ( x ) , e 3 ( x ) should be convergent, where e 3 ( x ) should satisfies Theorem 1 with the given condition 1 4 S e 3 < 1 . From Theorem 1, for 26 < ω < 38 , we extract
1 4 S e 3 = 1 4 max ω 64 + ω 32 3 4 + ω 64 + 1 4 , 2 ω 32 + 1 16 + ω 16 7 8 < 1 ,
which satisfies Theorem 1. This shows that the scheme is C 2 continuous.
For C 3 continuity, Equation (21) may be written as
e 4 ( x ) = 1 + x + x 2 + x 3 4 e 3 ( x ) ,
or
e 4 ( x ) = 1 16 ω x 7 + ω 16 + 1 4 x 6 + ω 16 + 3 4 x 5 + ω 16 + 1 x 4 + ω 16 + 1 x 3     + ω 16 + 3 4 x 2 + ω 16 + 1 4 x + ω 16 .
To check C 3 continuity, the scheme S γ analogous to γ ( x ) , e 4 ( x ) should be convergent, where e 4 ( x ) should satisfy Theorem 1 with the given condition 1 4 S e 4 < 1 . From Theorem 1, for 24 < ω < 36 , we extract
1 4 S e 4 = 1 4 max 1 16 ω + ω 16 1 , ω 16 + 1 4 + ω 16 3 4 < 1 .
Since the given condition satisfies Theorem 1, it must satisfy Theorem 2. Thus, the scheme is C 3 continuous.
For C 4 continuity, Equation (21) may be written as
e 5 ( x ) = 1 + x + x 2 + x 3 4 e 4 ( x ) ,
or
e 5 ( x ) = 1 4 ω x 4 + x 3 + ω 2 + 2 x 2 + x + ω 4 .
To check C 4 continuity, the scheme S γ analogous to γ ( x ) , e 5 ( x ) should be convergent, where e 5 ( x ) should satisfy Theorem 1 with the given condition 1 4 S e 5 < 1 . From Theorem 1, for 4 < ω < 8 , we extract
1 4 S e 5 = 1 4 max 1 2 ω , ω 2 2 < 1 ,
which satisfies Theorem 1. Since the given condition satisfies Theorem 1, it must satisfy Theorem 2. Hence, the given scheme is C 4 continuous.
Now, for C 5 continuity, Equation (21) may be written as
e 6 ( x ) = 1 + x + x 2 + x 3 4 e 5 ( x ) .
Hence,
e 6 ( x ) = 1 2 ( 1 + x ) .
Since, from Theorem 2, we have
1 4 S e 6 = 1 4 max 1 2 , 1 2 < 1 ,
the scheme is C 5 continuous.

3.3.2. Hölder’s Regularity

To find Hölder regularity, we use Theorem 3. The Laurent polynomial of the quinary scheme using Equation (18) may be written as
γ ( x ) = 1 4096 ω + 4 + 5 ω x + 28 + 13 ω x 2 + 104 + 25 ω x 3 + 280 + 36 ω x 4     + 600 + 36 ω x 5 + 1064 + 20 ω x 6 + 1608 12 ω x 7 + 2104 50 ω x 8     + 2400 74 ω x 9 + 2400 74 ω x 10 + 2104 50 ω x 11 + 1608 12 ω x 12     + 1064 + 20 ω x 13 + 600 + 36 ω x 14 + 280 + 36 ω x 15 + 104 + 25 ω x 16     + 28 + 13 ω x 17 + 4 + 5 ω x 18 + ω x 19 ,
γ ( x ) = 1 4096 1 + x + x 2 + x 3 5 e ( x ) ,
where
e ( x ) = ( 1 + x ) 2 ω x 2 2 ω x + ω + 4 x .
If ω = 2 , we can get b = 8 . Using Theorem 3 with m = 5 and b = 8 , we have r 5 log 2 ( 8 ) = 4.09691.

3.3.3. Limit Stencil

The matrix form of the scheme using Equation (18) has the form
λ 2 j λ 1 j λ 0 j λ 1 j λ 2 j λ 3 j = 77 2048 99 256 501 1024 11 128 1 2048 0 27 2048 69 256 563 1024 21 128 7 2048 0 7 2048 21 128 563 1024 69 256 27 2048 0 1 2048 11 128 501 1024 99 256 77 2048 0 0 77 2048 99 256 501 1024 11 128 1 2048 0 27 2048 69 256 563 1024 21 128 7 2048 λ 2 j 1 λ 1 j 1 λ 0 j 1 λ 1 j 1 λ 2 j 1 λ 3 j 1
and the local subdivision matrix
X = 77 2048 99 256 501 1024 11 128 1 2048 0 27 2048 69 256 563 1024 21 128 7 2048 0 7 2048 21 128 563 1024 69 256 27 2048 0 1 2048 11 128 501 1024 99 256 77 2048 0 0 77 2048 99 256 501 1024 11 128 1 2048 0 27 2048 69 256 563 1024 21 128 7 2048 ,
which shows the size of the invariant neighborhood is 6. After simplification, matrix X has eigenvalues λ = 1 , 1 4 , 1 16 , 1 64 , 1 256 , 1 1024 and eigenvectors
1 1 4 1 16 1 64 1 256 1 1024 1 45 , 643 796 , 637 13 17 3679 9911 286 , 979 1 , 391 , 129 94 169 1 11 , 123 796 , 637 7 17 7 583 52 , 171 1 , 391 , 129 19 169 1 7081 796 , 637 1 17 137 9911 26 , 369 1 , 391 , 129 11 169 1 9905 796 , 637 5 17 35 901 28 , 201 1 , 391 , 129 4 169 1 34 , 309 796 , 637 11 17 163 901 59 , 519 1 , 391 , 129 64 169 1 1 1 1 1 1 .
After simplification on the same manner as presented in Section 3.1.4, we can get
λ 2 j λ 1 j λ 0 j λ 1 j λ 2 j λ 3 j = 1 1 1 1 1 1 45 , 643 796 , 637 45 , 643 796 , 637 45 , 643 796 , 637 45 , 643 796 , 637 45 , 643 796 , 637 45 , 643 796 , 637 13 17 13 17 13 17 13 17 13 17 13 17 3679 9911 3679 9911 3679 9911 3679 9911 3679 9911 3679 9911 286 , 979 1 , 391 , 129 286 , 979 1 , 391 , 129 286 , 979 1 , 391 , 129 286 , 979 1 , 391 , 129 286 , 979 1 , 391 , 129 286 , 979 1 , 391 , 129 94 169 94 169 94 169 94 169 94 169 94 169 λ 2 j 1 λ 1 j 1 λ 0 j 1 λ 1 j 1 λ 2 j 1 λ 3 j 1 ,
which shows that the limit stencils are stable/constant.

4. Results and Discussion

This section consists of three major parametric effects of the schemes presented by Equations (3), (12), and (18).

4.1. Error Bound

This section presents the error between control polygon and limit curve after kth subdivision level of 5-point binary, ternary, and quaternary subdivision schemes using different values mentioned in Table 1, Table 2 and Table 3 by applying the approach of Hashmi [31]. The error is minimum over the interval ω [ 0 , 8 ] , ω [ 0 , 15 ] , and ω [ 0 , 25 ] for binary, ternary, and quaternary, respectively, and increases on both sides of the intervals. In Table 1, Table 2 and Table 3, it is observed that increases in the arity of the schemes decrease the error of the proposed schemes. Figure 1, Figure 2 and Figure 3 illustrate graphical representation of error. Moreover, the proposed computational cost decreases by increasing the arity of subdivision schemes. Therefore, our experiments show that higher arity scheme are better than the lower arity schemes in the sense of computational cost and error bounds.

4.2. Continuity

This section describes the effects of parameters for the schemes in Equations (3), (12), and (18). The order of continuity and effects of parameters ω of the schemes are shown in Table 4, Table 5 and Table 6, respectively. This can easily be found over the parametric intervals using the approach of Hassan [8].

4.3. Shapes of Limit Curves

The parametric effect and continuity of the limit curve of the schemes are shown in Figure 4, Figure 5 and Figure 6, respectively. These figures illustrate the role of free parameters when 5-point binary, ternary, and quaternary approximating schemes are applied on discrete data point. One can see the looseness/tightness of the limit curves in Figure 4, Figure 5 and Figure 6 when the parameter values change.

5. Conclusions

In this work, we introduce the family of 5-point schemes which depict the representation of a wide variety of shapes with high smoothness (continuity) and less computational cost (processing time). These properties are useful in computer aided geometric design and geometric modeling. We apply Laurent polynomial to analyze our schemes. The shape parameter ω makes it able to provide different results along with its applications.

Author Contributions

Conceptualization, K.S.N.; Formal analysis, D.B.; Methodology, D.B.; Software, S.M.H., A.U.R., K.S.N., and A.G.; Supervision, K.S.N.; Writing—original draft, S.M.H., A.U.R., and A.G.; and Writing—review and editing, K.S.N. and S.A.A.K.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Rham, G. Un peu de mathématiques à propos d’une courbe plane. Elem. Math. 1947, 2, 73–76. [Google Scholar]
  2. Chaikin, G.M. An algorithm for high-speed curve generation. Comput. Gr. Imag. Process. 1974, 3, 346–349. [Google Scholar] [CrossRef]
  3. Doo, D.; Sabin, M. Behaviour of recursive division surfaces near extraordinary points. Comput. Aided Des. 1978, 10, 356–360. [Google Scholar] [CrossRef]
  4. Catmull, E.; Clark, J. Recursively generated B-spline surfaces on arbitrary topological meshes. Comput. Aided Des. 1978, 10, 350–355. [Google Scholar] [CrossRef]
  5. Deslauriers, G.; Dubuc, S. Symmetric iterative interpolation processes. In Constructive Approximation; DeVore, R.A., Saff, E.B., Eds.; Springer: Boston, MA, USA, 1989; pp. 49–68. [Google Scholar]
  6. Dyn, N.; Levine, D.; Gregory, J.A. A butterfly subdivision scheme for surface interpolation with tension control. ACM Trans. Gr. 1990, 9, 160–169. [Google Scholar] [CrossRef]
  7. Cai, Z. Convergence, error estimation and some properties of four-point interpolation subdivision scheme. Comput. Aided Geom. Des. 1995, 12, 459–468. [Google Scholar]
  8. Hassan, M.F.; Ivrissimitzis, I.P.; Dodgson, N.A.; Sabin, M.A. An interpolating 4-point C2 ternary stationary subdivision scheme. Comput. Aided Geom. Des. 2002, 19, 1–18. [Google Scholar] [CrossRef]
  9. Hassan, M.; Dodgson, N.A. Ternary and Three-point Univariate Subdivision Schemes; No. UCAM-CL-TR-520; Cambridge Computer Laboratory: Cambridge, UK, September 2001. [Google Scholar]
  10. Mustafa, G.; Liu, X. A subdivision scheme for volumetric models. Appl. Math. J. Chin. Univ. 2005, 20, 213–224. [Google Scholar] [CrossRef]
  11. Siddiqi, S.S.; Ahmad, N. A C6 approximating subdivision scheme. Appl. Math. Lett. 2008, 21, 722–728. [Google Scholar] [CrossRef] [Green Version]
  12. Hormann, K.; Sabin, M.A. A family of subdivision schemes with cubic precision. Comput. Aided Geom. Des. 2008, 25, 41–52. [Google Scholar] [CrossRef]
  13. Faheem, K.; Mustafa, G. Ternary six-point interpolating subdivision scheme. Lobachevskii J. Math. 2008, 29, 153–163. [Google Scholar] [CrossRef]
  14. Mustafa, G.; Khan, F. A new 4-point quaternary approximating subdivision scheme. Abstr. Appl. Anal. 2009, 2009. [Google Scholar] [CrossRef]
  15. Mustafa, G.; Khan, F.; Ghaffar, A. The m-point approximating subdivision scheme. Lobachevskii J. Math. 2009, 30, 138–145. [Google Scholar] [CrossRef]
  16. Siddiqi, S.S.; Rehan, K. Modified form of binary and ternary 3-point subdivision schemes. Appl. Math. Comput. 2010, 216, 970–982. [Google Scholar] [CrossRef]
  17. Mustafa, G.; Ghaffar, A.; Khan, F. The odd-point ternary approximating schemes. Am. J. Comput. Math. 2011, 1, 111–118. [Google Scholar] [CrossRef] [Green Version]
  18. Ghaffar, A.; Mustafa, G.; Qin, K. Unification and application of 3-point approximating subdivision schemes of varying arity. Open J. Appl. Sci. 2012, 2, 48–52. [Google Scholar] [CrossRef]
  19. Ghaffar, A.; Mustafa, G.; Qin, K. The 4-point 3-ary approximating subdivision scheme. Open J. Appl. Sci. 2013, 3, 106–111. [Google Scholar] [CrossRef]
  20. Mustafa, G.; Ghaffar, A.; Aslam, M. A subdivision-regularization framework for preventing over fitting of data by a model. AAM 2013, 8, 178–190. [Google Scholar]
  21. Mustafa, G.; Ashraf, P.; Deng, J. Generalized and unified families of interpolating subdivision schemes. Numer. Math. Theory Method. Appl. 2014, 7, 193–213. [Google Scholar] [CrossRef]
  22. Siddiqi, S.S.; Younis, M. The Quaternary Interpolating Scheme for Geometric Design. Int. Sch. Res. Not. 2013, 2013. [Google Scholar] [CrossRef]
  23. Rehan, K.; Siddiqi, S.S. A Family of Ternary Subdivision Schemes for Curves. Appl. Math. Comput. 2015, 270, 114–123. [Google Scholar] [CrossRef]
  24. Rehan, K.; Sabri, M.A. A combined ternary 4-point subdivision scheme. Appl. Math. Comput. 2016, 276, 278–283. [Google Scholar] [CrossRef]
  25. Ashraf, P.; Sabir, M.; Ghaffar, A.; Nisar, K.S.; Khan, I. Shape-Preservation of Ternary Four-point Interpolating Non-stationary Subdivision Scheme. Front. Phys. 2020, 7. [Google Scholar] [CrossRef]
  26. Ghaffar, A.; Ullah, Z.; Bari, M.; Nisar, K.S.; Al-Qurashi, M.M.; Baleanu, D. A new class of 2m-point binary non-stationary subdivision schemes. Adv. Differ. Equ. 2019, 2019, 325. [Google Scholar] [CrossRef]
  27. Ghaffar, A.; Ullah, Z.; Bari, M.; Nisar, K.S.; Baleanu, D. Family of odd point non-stationary subdivision schemes and their applications. Adv. Differ. Equ. 2019, 2019, 1–20. [Google Scholar] [CrossRef]
  28. Ghaffar, A.; Bari, M.; Ullah, Z.; Iqbal, M.; Nisar, K.S.; Baleanu, D. A New Class of 2q-Point Nonstationary Subdivision Schemes and Their Applications. Mathematics 2019, 7, 639. [Google Scholar] [CrossRef] [Green Version]
  29. Ghaffar, A.; Iqbal, M.; Bari, M.; Muhammad Hussain, S.; Manzoor, R.; Sooppy Nisar, K.; Baleanu, D. Construction and Application of Nine-Tic B-Spline Tensor Product SS. Mathematics 2019, 7, 675. [Google Scholar] [CrossRef] [Green Version]
  30. Dyn, N.; Iske, A.; Quak, E.; Floater, M.S. Tutorials on Multiresolution in Geometric Modelling, Summer School Lecture Notes Series: Mathematics and Visualization; Springer Science & Business Media: Berlin, Germany, 2002. [Google Scholar]
  31. Mustafa, G.; Hashmi, M.S. Subdivision depth computation for n-ary subdivision curves/surfaces. Vis. Comput. 2010, 26, 841–851. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Error bounds of binary scheme (Equation (3)).
Figure 1. Error bounds of binary scheme (Equation (3)).
Mathematics 08 00474 g001
Figure 2. Error bounds of ternary scheme (Equation (12)).
Figure 2. Error bounds of ternary scheme (Equation (12)).
Mathematics 08 00474 g002
Figure 3. Error bounds of Quaternary scheme (Equation (18)).
Figure 3. Error bounds of Quaternary scheme (Equation (18)).
Mathematics 08 00474 g003
Figure 4. Parameter effects on the limit curves with initial polygon of binary scheme.
Figure 4. Parameter effects on the limit curves with initial polygon of binary scheme.
Mathematics 08 00474 g004
Figure 5. Parameter effects on the limit curves with initial polygon of ternary scheme.
Figure 5. Parameter effects on the limit curves with initial polygon of ternary scheme.
Mathematics 08 00474 g005
Figure 6. Parameter effect on the limit curves with initial polygon of quaternary scheme.
Figure 6. Parameter effect on the limit curves with initial polygon of quaternary scheme.
Mathematics 08 00474 g006
Table 1. Binary scheme error bounds.
Table 1. Binary scheme error bounds.
k123456
ω [ 0 , 8 ] 0.1750000.0875000.0437500.0218750.0109380.005469
ω = 9 0.2398150.1386430.0801530.0463380.0267890.015488
ω = 1 0.1862900.0960560.0495290.0255380.0131680.006790
ω = 2 0.2250000.11250000.0562500.02815000.01407500.00703850
Table 2. Ternary scheme error bounds.
Table 2. Ternary scheme error bounds.
k123456
ω [ 0 , 15 ] 0.0833330.0277780.0092590.0030860.0010290.000343
ω = 16 0.0880500.0304370.0105210.0036370.0012570.000435
ω = 1 0.0844960.0284260.0095630.0032170.0010820.000364
ω = 2 0.0856700.0290850.0098750.0033530.0011380.000386
Table 3. Quaternary scheme error bounds.
Table 3. Quaternary scheme error bounds.
k123456
ω [ 0 , 25 ] 0.0541670.0135420.0033850.0008460.0002120.000053
ω = 8 0.0880500.0304370.0105210.0036370.0012570.000435
ω = 1 0.0544490.0136650.0034300.0008610.0002160.000054
ω = 2 0.0547320.0137900.0034740.0008750.0002210.000056
Table 4. Continuity order of binary scheme (Equation (3)).
Table 4. Continuity order of binary scheme (Equation (3)).
SchemeParameterContinuitySchemeParameterContinuity
Binary 9 < ω < 72 5 C 0 Binary 2 < ω < 6 C 4
Binary 9 < ω < 72 5 C 1 Binary 2 < ω < 6 C 5
Binary 6 < ω < 10 C 2 Binary 0 < ω < 2 C 6
Binary 6 < ω < 10 C 3 Binary ω = 1 C 7
Table 5. Continuity order of ternary scheme (Equation (12)).
Table 5. Continuity order of ternary scheme (Equation (12)).
SchemeParameterContinuity
Ternary 2929 9 < ω < 356 9 C 0
Ternary 23 < ω < 31 C 1
Ternary 14 < ω < 22 C 2
Ternary 8 < ω < 16 C 3
Ternary 4 < ω < 12 C 4
Ternary ω = 4 3 C 5
Table 6. Continuity order of Quaternary scheme (Equation (18)).
Table 6. Continuity order of Quaternary scheme (Equation (18)).
SchemeParameterContinuity
Quaternary 3775 51 < ω < 4415 51 C 0
Quaternary 57 < ω < 71 C 1
Quaternary 26 < ω < 38 C 2
Quaternary 24 < ω < 36 C 3
Quaternary 4 < ω < 8 C 4
Quaternary ω = 2 C 5

Share and Cite

MDPI and ACS Style

Hussain, S.M.; Rehman, A.U.; Baleanu, D.; Nisar, K.S.; Ghaffar, A.; Abdul Karim, S.A. Generalized 5-Point Approximating Subdivision Scheme of Varying Arity. Mathematics 2020, 8, 474. https://doi.org/10.3390/math8040474

AMA Style

Hussain SM, Rehman AU, Baleanu D, Nisar KS, Ghaffar A, Abdul Karim SA. Generalized 5-Point Approximating Subdivision Scheme of Varying Arity. Mathematics. 2020; 8(4):474. https://doi.org/10.3390/math8040474

Chicago/Turabian Style

Hussain, Sardar Muhammad, Aziz Ur Rehman, Dumitru Baleanu, Kottakkaran Sooppy Nisar, Abdul Ghaffar, and Samsul Ariffin Abdul Karim. 2020. "Generalized 5-Point Approximating Subdivision Scheme of Varying Arity" Mathematics 8, no. 4: 474. https://doi.org/10.3390/math8040474

APA Style

Hussain, S. M., Rehman, A. U., Baleanu, D., Nisar, K. S., Ghaffar, A., & Abdul Karim, S. A. (2020). Generalized 5-Point Approximating Subdivision Scheme of Varying Arity. Mathematics, 8(4), 474. https://doi.org/10.3390/math8040474

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