# A Family of Even-Point Ternary Approximating Schemes

Authors: Abdul Ghaffar; Ghulam Mustafa
Journal: ISRN Applied Mathematics
 (2012)

---

## Abstract

We presented a general formula to generate the family of even-point
ternary approximating subdivision schemes with a shape parameter for describing curves. Some sufficient conditions forC0 to C7 continuity and
approximation order for certain ranges of parameter are discussed. The proposed even-point ternary schemes compare remarkably with existing even-point ternary schemes because they are able to generate limit functions with
higher smoothness and approximation order. In addition, we measured curvature and torsion that assist the quality of subdivided curves.


---

## Body

## 1. Introduction and Preliminaries

There are numbers of binary subdivision schemes in the literature. The interest in investigating arities higher than two has been started by Hassan et al. [1, 2]. Nowadays, we have numbers of ternary schemes introduced by [3–7] and so forth, But the research communities are still gaining interest in introducing schemes higher than three arities (i.e., quaternary, quinary, senary,…, n-ary). Mustafa and Khan [8] introduced a new 4-point C3 quaternary approximating subdivision scheme. Lian [9, 10] introduced 3-, 4-, 5-, and 6-point a-ary schemes. Lian [11] also offered 2m-point and (2m+1)-point interpolating a-ary schemes for curve design.The 2-, 3-,…, 6-point binary and ternary schemes are very common in the literature. The schemes involving convex combination of more or less than six points at coarse refinement level to insert a new point at next refinement level is introduced by Ko et al. [7]. They introduced (2n+2)- and (2n+4)-point binary schemes. Zheng et al. [12] investigated ternary interpolatory schemes with an odd number of control points, namely, (2n-1)-point ternary interpolatory subdivision scheme. They also investigated ternary even symmetric 2n-point [13] and p-ary [14] approximating subdivision scheme and presented the general ternary even symmetric 2n-point approximating subdivision rule and design alternative smooth ternary subdivision scheme of higher order. Mustafa and Rehman [15] presented general formulae for the mask of (2b+4)-point n-ary approximating as well as interpolating subdivision schemes for any integers b≥0 and n≥2.This motivates us to present the family of even-point ternary schemes with high smoothness and more degree of freedom for curve design. Proposed schemes not only provide the mask of even-point schemes but also generalize and unify several well-known schemes. Moreover, we measured curvature and torsion that can be used to describe the quality of curve. Also we compared plot of curvature and torsion, obtained by proposed schemes with the other existing schemes.A general compact form of univariate ternary subdivision schemeS, which maps polygon fk={fik}i∈ℤ to a refined polygon fk+1={fik+1}i∈ℤ is defined by
(1.1)fik+1=∑j∈ℤ‍a3j-ifjk,i∈ℤ,
where the set a={ai:i∈ℤ} of coefficients is called the mask at k-th level of refinement. A necessary condition for the uniform convergence of subdivision scheme (1.1) is that
(1.2)∑j∈ℤ‍a3j=∑j∈ℤ‍a3j+1=∑j∈ℤ‍a3j+2=1.
A subdivision scheme is uniformly convergent if for any initial data f0={fi0:i∈ℤ}, there exists a continuous function f such that for any closed interval I⊂ℝ(1.3)limk→∞supi∈3kI|fik-f(3-ki)|=0.
Obviously, f=S∞f0, introducing a symbol called Laurent polynomial
(1.4)a(z)=∑i∈ℤ‍aizi,
of the mask a={ai:i∈ℤ} which plays an efficient role to analyze the convergence and smoothness of subdivision scheme. From (1.2) and (1.4) the Laurent polynomial of convergent subdivision scheme satisfies
(1.5)a(e2iπ/3)=a(e4iπ/3)=0,a(1)=3.
This condition guarantees existence of a related subdivision scheme for the divided difference of the original control points and the existence of associated Laurent polynomial a(1)(z)(1.6)a(1)(z)=3z21+z+z2a(z).
The subdivision scheme S1 with Laurent polynomial a(1)(z) is related to scheme S with Laurent polynomial a(z) by the following theorem.Theorem 1.1 (see [1]).
LetS denote a subdivision scheme with Laurent polynomial a(z) satisfying (1.5). Then there exists a subdivision scheme S1 with the property
(1.7)Δfk=S1Δfk-1,
where fk=Skf0 and Δfk={(Δfk)i=3k(fi+1k-fik);i∈ℤ}. Furthermore, S is a uniformly convergent if and only if (1/3)S1 converges uniformly to zero function for all initial data f0, in the sense that
(1.8)limk→0(13S1)kf0=0.A schemeS1 satisfying (1.8) for all initial data f0 is termed “contractive.” The above theorem indicates that for any given scheme S, with mask “a” satisfying (1.2), we can prove the uniform convergence of S by deriving the mask of (1/3)S1 and computing ∥((1/3)S1)i∥∞ for i=1,2,3...,L, where L is the first integer for which ∥((1/3)S1)L∥∞<1. If such an L exists, then S converges uniformly. Since there are three rules for computing the values at next refinement level, so we define the norm as follows:
(1.9)∥S∥∞=max{∑j∈ℤ‍|a3j|,∑j∈ℤ‍|a3j+1|,∑j∈ℤ‍|a3j+2|},(1.10)∥(13Sn)L∥∞=max{∑j∈ℤ|bi+3Lj[n,L]|;i=0,1,2,…,3L-1},
where
(1.11)b[n,L](z)=13L∏j=0L-1‍a(n)(z3j),(1.12)a(n)(z)=(3z21+z+z2)a(n-1)(z)=(3z21+z+z2)na(z),n⩾1.The paper is organized as follows: the family of even-point ternary approximating scheme and analysis of two even-point ternary schemes are presented in Section2. Basic properties of even-point ternary schemes are discussed in Section 3. Comparison with existing even-point ternary schemes is shown in Section 4. A few remarks and future work constitute Section 5.


## 2. The Even-Point Ternary Approximating Schemes

Here we offer a general formula for even-point ternary approximating subdivision schemes with one parameter in the form of Laurent polynomial(2.1)a[λ+1](z)=13λ+n-1(1+z+z2)λ+n×{(112+ω)+(512-ω)z+(512-ω)z2+(112+ω)z3},
where ω is shape parameter, λ=2n+1,n⩾0 and depends on n.Despite the fact that we can generate(λ+1)-point ternary approximating schemes from (2.1) for any λ=2n+1, n⩾0, however for straightforwardness, we generated and discussed the smoothness of 4- and 6-point ternary schemes. The general proof regarding the support of even-point ternary schemes has been presented in this paper.


### 2.1.4-Point Ternary Scheme

From (2.1) for n=1, we get the following: Laurent polynomial for a 4-point ternary scheme
(2.2)a[4](z)=127{(112+ω)z0+(34+3ω)z1+(3512+5ω)z2+(294+3ω)z3+(514-3ω)z4+(674-9ω)z5+(674-9ω)z6+(514-3ω)z7+(294+3ω)z8(3512+5ω)z9+(34+3ω)z10+(112+ω)z11}.
From the above polynomial, we suggest the following 4-point ternary approximating scheme:
(2.3)f3ik+1=127[(3512+5ω)fi-1k+(674-9ω)fik+(294+3ω)fi+1k+(112+ω)fi+2k],f3i+1k+1=127[(34+3ω)fi-1k+(514-3ω)fik+(514-3ω)fi+1k+(34+3ω)fi+2k],f3i+2k+1=127[(112+ω)fi-1k+(294+3ω)fik+(674-9ω)fi+1k+(3512+5ω)fi+2k].
From (1.12) and (2.2), we have
(2.4)a[4](1)(z)=z29{(112+ω)+(23+2ω)z+(136+2ω)z2+(5312-ω)z3+(376-4ω)z4+(376-4ω)z5+(5312-ω)z6+(136+2ω)z7+(23+2ω)z8+(112+ω)z9},a[4](2)(z)=z43{(112+ω)z0+(712+ω)z1+3z2+(73-2ω)z3+(73-2ω)z4+3z5+(712+ω)z6+(112+ω)z7},a[4](3)(z)=(112+ω)+12z+(1112-ω)z2+(1112-ω)z3+12z4+(112+ω)z5,a[4](4)(z)=3z2{(112+ω)+(512-ω)z+(512-ω)z2+(112+ω)z3}.
If Sm[4] is the scheme corresponding to a[4](m), then for Cm-1 continuity, we require that a[4](m) satisfies (1.2), which it does and ∥((1/3)Sm[4])L∥∞<1, for m=1,2,3,4.Since by (1.10), for -71/24<ω<91/24, -11/6<ω<8/3, -13/12<ω<23/12, and -5/12<ω<7/12, we have
(2.5)∥13S1[4]∥∞=max{227|112+ω|+227|5312-ω|,127|23+2ω|+127|376-4ω|+127|136+2ω|}<1,∥13S2[4]∥∞=max{19|112+ω|+19|74-2ω|+29|712+ω|,23}<1,∥13S3[4]∥∞=max{13|112+ω|+13|1112-ω|,13}<1,∥13S4[4]∥∞=max{2|112+ω|,|512-ω|,|512-ω|}<1.
In order to obtain a C4 scheme, we substitute ω=1/12 in a[4](4)(z), we have
(2.6)a[4](4)(z)=z22(1+2z+2z2+z3).
Using (1.12), we get
(2.7)a[4](5)(z)=(3z21+z+z2)a4(z)=32z4{1+z}.
If Sm[4] is the scheme corresponding to a[4](m), then for Cm-1 continuity, we require that a[4](m) satisfy (1.2), which it does and ∥((1/3)Sm[4])L∥∞<1, for m=1,2,3,4.Form=5 and L=1, we have
(2.8)∥13S5[4]∥∞=13max{∑j∈ℤ‍|a3j|,∑j∈ℤ‍|a3j+1|,∑j∈ℤ‍|a3j+2|},=13max{32,32,0}=12<1.
Therefore (1/3)S5[4] is contractive, that is, ∥(1/3)S5[4]∥∞<1 and the scheme is C4.


### 2.2. 6-Point Ternary Scheme

From (2.1) for n=2, we get the following 6-point ternary approximating scheme:
(2.9)f3ik+1=1729[(173+20ω)fi-2k+(3713+56ω)fi-1k+(437312-127ω)fik+(256912+7ω)fi+1k+(25312+43w)fi+2k+(112+ω)fi+3k],f3i+1k+1=1729[(6w+1)fi-2k+(2314+63ω)fi-1k+(12234-69ω)fik+(12234-69ω)fi+1k+(2314+63ω)fi+2k+(1+6ω)fi+3k],f3i+2k+1=1729[(112+ω)fi-2k+(25312+43ω)fi-1k+(256912+7ω)fik+(437312-127ω)fi+1k+(3713+56ω)fi+2k+(173+20ω)fi+3k].
The continuity of this scheme can be computed in a similar way as we did for the 4-point ternary scheme.From above discussion, we reach the conclusion shown in Table1.Table 1
The order of continuity of proposed 4-point and 6-point ternary schemes for certain ranges of parameter.
Scheme
Parameter
Continuity
Scheme
Parameter
Continuity
4-point
-
71
24
<
ω
<
91
24
C
0
-
329
60
<
ω
<
1489
204
C
0
-
11
6
<
ω
<
8
3
C
1
-
59
12
<
ω
<
509
12
C
1
-
13
12
<
ω
<
23
12
C
2
-
149
36
<
ω
<
175
36
C
2
-
5
12
<
ω
<
7
12
C
3
-
71
24
<
ω
<
91
24
C
3
ω
=
1
12
C
4
6-point
-
11
6
<
ω
<
8
3
C
4
-
13
12
<
ω
<
23
12
C
5
-
5
12
<
ω
<
7
12
C
6
ω
=
1
12
C
7Remark 2.1.
(i)
Forn=1, ω=-1/4+9u, it becomes 4-point ternary subdivision scheme [14].
(ii)
If we setn=1, ω=-1/4+9w, it becomes 4-point ternary subdivision scheme [13].
(iii)
In casen=2, ω=-11/12+(81/7)u, we get mask of 6-point ternary subdivision scheme [13].


## 2.1.4-Point Ternary Scheme

From (2.1) for n=1, we get the following: Laurent polynomial for a 4-point ternary scheme
(2.2)a[4](z)=127{(112+ω)z0+(34+3ω)z1+(3512+5ω)z2+(294+3ω)z3+(514-3ω)z4+(674-9ω)z5+(674-9ω)z6+(514-3ω)z7+(294+3ω)z8(3512+5ω)z9+(34+3ω)z10+(112+ω)z11}.
From the above polynomial, we suggest the following 4-point ternary approximating scheme:
(2.3)f3ik+1=127[(3512+5ω)fi-1k+(674-9ω)fik+(294+3ω)fi+1k+(112+ω)fi+2k],f3i+1k+1=127[(34+3ω)fi-1k+(514-3ω)fik+(514-3ω)fi+1k+(34+3ω)fi+2k],f3i+2k+1=127[(112+ω)fi-1k+(294+3ω)fik+(674-9ω)fi+1k+(3512+5ω)fi+2k].
From (1.12) and (2.2), we have
(2.4)a[4](1)(z)=z29{(112+ω)+(23+2ω)z+(136+2ω)z2+(5312-ω)z3+(376-4ω)z4+(376-4ω)z5+(5312-ω)z6+(136+2ω)z7+(23+2ω)z8+(112+ω)z9},a[4](2)(z)=z43{(112+ω)z0+(712+ω)z1+3z2+(73-2ω)z3+(73-2ω)z4+3z5+(712+ω)z6+(112+ω)z7},a[4](3)(z)=(112+ω)+12z+(1112-ω)z2+(1112-ω)z3+12z4+(112+ω)z5,a[4](4)(z)=3z2{(112+ω)+(512-ω)z+(512-ω)z2+(112+ω)z3}.
If Sm[4] is the scheme corresponding to a[4](m), then for Cm-1 continuity, we require that a[4](m) satisfies (1.2), which it does and ∥((1/3)Sm[4])L∥∞<1, for m=1,2,3,4.Since by (1.10), for -71/24<ω<91/24, -11/6<ω<8/3, -13/12<ω<23/12, and -5/12<ω<7/12, we have
(2.5)∥13S1[4]∥∞=max{227|112+ω|+227|5312-ω|,127|23+2ω|+127|376-4ω|+127|136+2ω|}<1,∥13S2[4]∥∞=max{19|112+ω|+19|74-2ω|+29|712+ω|,23}<1,∥13S3[4]∥∞=max{13|112+ω|+13|1112-ω|,13}<1,∥13S4[4]∥∞=max{2|112+ω|,|512-ω|,|512-ω|}<1.
In order to obtain a C4 scheme, we substitute ω=1/12 in a[4](4)(z), we have
(2.6)a[4](4)(z)=z22(1+2z+2z2+z3).
Using (1.12), we get
(2.7)a[4](5)(z)=(3z21+z+z2)a4(z)=32z4{1+z}.
If Sm[4] is the scheme corresponding to a[4](m), then for Cm-1 continuity, we require that a[4](m) satisfy (1.2), which it does and ∥((1/3)Sm[4])L∥∞<1, for m=1,2,3,4.Form=5 and L=1, we have
(2.8)∥13S5[4]∥∞=13max{∑j∈ℤ‍|a3j|,∑j∈ℤ‍|a3j+1|,∑j∈ℤ‍|a3j+2|},=13max{32,32,0}=12<1.
Therefore (1/3)S5[4] is contractive, that is, ∥(1/3)S5[4]∥∞<1 and the scheme is C4.


## 2.2. 6-Point Ternary Scheme

From (2.1) for n=2, we get the following 6-point ternary approximating scheme:
(2.9)f3ik+1=1729[(173+20ω)fi-2k+(3713+56ω)fi-1k+(437312-127ω)fik+(256912+7ω)fi+1k+(25312+43w)fi+2k+(112+ω)fi+3k],f3i+1k+1=1729[(6w+1)fi-2k+(2314+63ω)fi-1k+(12234-69ω)fik+(12234-69ω)fi+1k+(2314+63ω)fi+2k+(1+6ω)fi+3k],f3i+2k+1=1729[(112+ω)fi-2k+(25312+43ω)fi-1k+(256912+7ω)fik+(437312-127ω)fi+1k+(3713+56ω)fi+2k+(173+20ω)fi+3k].
The continuity of this scheme can be computed in a similar way as we did for the 4-point ternary scheme.From above discussion, we reach the conclusion shown in Table1.Table 1
The order of continuity of proposed 4-point and 6-point ternary schemes for certain ranges of parameter.
Scheme
Parameter
Continuity
Scheme
Parameter
Continuity
4-point
-
71
24
<
ω
<
91
24
C
0
-
329
60
<
ω
<
1489
204
C
0
-
11
6
<
ω
<
8
3
C
1
-
59
12
<
ω
<
509
12
C
1
-
13
12
<
ω
<
23
12
C
2
-
149
36
<
ω
<
175
36
C
2
-
5
12
<
ω
<
7
12
C
3
-
71
24
<
ω
<
91
24
C
3
ω
=
1
12
C
4
6-point
-
11
6
<
ω
<
8
3
C
4
-
13
12
<
ω
<
23
12
C
5
-
5
12
<
ω
<
7
12
C
6
ω
=
1
12
C
7Remark 2.1.
(i)
Forn=1, ω=-1/4+9u, it becomes 4-point ternary subdivision scheme [14].
(ii)
If we setn=1, ω=-1/4+9w, it becomes 4-point ternary subdivision scheme [13].
(iii)
In casen=2, ω=-11/12+(81/7)u, we get mask of 6-point ternary subdivision scheme [13].


## 3. Basic Properties of the Schemes

In this section, we discuss approximation order and support of basic limit function of even-point ternary approximating schemes.


### 3.1. Approximation Order

Here we only find the approximation order of proposed 4-point ternary scheme. The approximation order of other even-point ternary schemes can be computed in a similar fashion.Theorem 3.1.
A 4-point ternary approximating scheme has approximating order 4 forω∈(-5/12<ω<7/12) and 5 for ω=5/12.
A 6-point ternary approximating scheme has approximating order 7 forω∈(-5/12<ω<7/12) and 8 for ω=1/12.Proof.
We carry out this result by taking our origin the middle of an original span with ordinate…,(-5)n,(-3)n,(-1)n,1n,3n,5n,….
Ify=xn, then we have
(3.1)[y]=…,a1(-5)n+a2(-3)n+a3(-1)n+a4(1)n+a5(-5)n+a6(-3)n+a6(-1)n+a5(1)n+a4(-5)n+a3(-3)n+a2(-1)n+a1(1)n+a1(-3)n+a2(-1)n+a3(1)n+a4(3)n+a5(-3)n+a6(-1)n+a6(1)n+a5(3)n+⋯++a5(-1)n+a6(1)n+a6(3)n+a5(5)n+a4(-1)n+a3(1)n+a2(3)n+a1(5)n,…,
wherea1=(1/27)(35/12+5ω), a2=(1/27)(67/4-9ω), a3=(1/27)(29/4+3ω), a4=(1/27)(1/12+ω), a5=(1/27)(3/4+3ω), a6=(1/27)(51/3-3ω).
Ify=x1, then
(3.2)[y]=…,-83,-2,-43,-23,2,23,43,2,83,…[δy]=…,23,23,23,23,23,23,23,…,[δ2y]=0,
where δ represents the differences of the vertices.
Ify=x2, then
(3.3)[y]=…,779+169ω,499+169ω,299+169ω,179+169ω,139+169ω,179+169ω,299+169ω,499+169ω,779+169ω,….
Taking further differences, we get [δ3y]=0.
Ify=x3, then
(3.4)[y]=…,-82427-1289ω,-503-323ω,-22027-649ω,-8627-329ω,0,8627+329ω,22027+649ω,503+323ω,82427+12827ω,….
This implies that [δ4y]=0.
Ify=x4, then
(3.5)[y]=…,318727+8009ω,5059+5449ω,65927+32ω,899+1609ω,499+1609ω,899+1609ω,65927+32ω,5059+5449ω,318727+8009ω,….
By taking further differences, we have
(3.6)[δ5y]=…,16027-1289ω,16027-1289ω,0,16027-1289ω,16027-1289ω,…,[δ5y]=0,atω=512.
Since by [8], 4-point ternary scheme has cubic precision for ω∈(-5/12<ω<7/12) and quartic at ω=5/12 then 4-point ternary scheme has approximating order 4 for ω∈(-5/12<ω<7/12) and 5 for ω=5/12. The proof of other part is similar.


### 3.2. Support of Basic Limit Function

The basic limit functions,φ[λ+1]=(S[λ+1])∞fi0 of proposed (λ+1)-point ternary approximating schemes, for λ=2n+1, n=1,2 are presented in Figures 1(a) and 1(b). The following theorem is related to the support of basic limit functions of even point ternary schemes.The basic limit functions of proposed schemes atω=1/12.
(a)
4-point scheme
(b)
6-point schemeTheorem 3.2.
The basic limit functionsφ[λ+1] of proposed (λ+1)-point ternary approximating schemes have support width s=(2(λ+n)+3)/2, where λ=2n+1, n⩾0, which implies that it vanishes outside the interval [-(2(λ+n)+3)/4,(2(λ+n)+3)/4].Proof.
The support width “s” of the basic limit functions can be determined by computing how for the effect of the nonzero vertex f00 will propagate along by. As the mask of (λ+1) scheme is 3(λ+1)-long sequence by centering it on that vertex, the distances to the last of its left and right nonzero coefficients are equal to λ+n+2 and λ+n+1, respectively. At the first subdivision step we see that the vertices on the left and right sides of f01 at (λ+n+2)/3 and (λ+n+1)/3 are the furthest nonzero new vertices. At each refinement, the distances on both sides are reduced by the factor 1/3. At the next step of the scheme this will propagate along by ((λ+n+2)/3)×1/3 on the left and ((λ+n+1)/3)×1/3 on the right. Hence after k subdivision steps the furthest nonzero vertex on the left will be at (λ+n+2)(1/3+1/32+⋯+1/3k)=((λ+n+2)/3)(∑j=0k-1(1/3j)) and on the right will be at (λ+n+1)(1/3+1/32+⋯+1/3k)=((λ+n+1)/3)(∑j=0k-1(1/3j)). So the total support width is ((λ+n+2)/3)(∑j=0∞(1/3j))+((λ+n+1)/3)(∑j=0∞(1/3j))=(2(λ+n)+3)/2.


## 3.1. Approximation Order

Here we only find the approximation order of proposed 4-point ternary scheme. The approximation order of other even-point ternary schemes can be computed in a similar fashion.Theorem 3.1.
A 4-point ternary approximating scheme has approximating order 4 forω∈(-5/12<ω<7/12) and 5 for ω=5/12.
A 6-point ternary approximating scheme has approximating order 7 forω∈(-5/12<ω<7/12) and 8 for ω=1/12.Proof.
We carry out this result by taking our origin the middle of an original span with ordinate…,(-5)n,(-3)n,(-1)n,1n,3n,5n,….
Ify=xn, then we have
(3.1)[y]=…,a1(-5)n+a2(-3)n+a3(-1)n+a4(1)n+a5(-5)n+a6(-3)n+a6(-1)n+a5(1)n+a4(-5)n+a3(-3)n+a2(-1)n+a1(1)n+a1(-3)n+a2(-1)n+a3(1)n+a4(3)n+a5(-3)n+a6(-1)n+a6(1)n+a5(3)n+⋯++a5(-1)n+a6(1)n+a6(3)n+a5(5)n+a4(-1)n+a3(1)n+a2(3)n+a1(5)n,…,
wherea1=(1/27)(35/12+5ω), a2=(1/27)(67/4-9ω), a3=(1/27)(29/4+3ω), a4=(1/27)(1/12+ω), a5=(1/27)(3/4+3ω), a6=(1/27)(51/3-3ω).
Ify=x1, then
(3.2)[y]=…,-83,-2,-43,-23,2,23,43,2,83,…[δy]=…,23,23,23,23,23,23,23,…,[δ2y]=0,
where δ represents the differences of the vertices.
Ify=x2, then
(3.3)[y]=…,779+169ω,499+169ω,299+169ω,179+169ω,139+169ω,179+169ω,299+169ω,499+169ω,779+169ω,….
Taking further differences, we get [δ3y]=0.
Ify=x3, then
(3.4)[y]=…,-82427-1289ω,-503-323ω,-22027-649ω,-8627-329ω,0,8627+329ω,22027+649ω,503+323ω,82427+12827ω,….
This implies that [δ4y]=0.
Ify=x4, then
(3.5)[y]=…,318727+8009ω,5059+5449ω,65927+32ω,899+1609ω,499+1609ω,899+1609ω,65927+32ω,5059+5449ω,318727+8009ω,….
By taking further differences, we have
(3.6)[δ5y]=…,16027-1289ω,16027-1289ω,0,16027-1289ω,16027-1289ω,…,[δ5y]=0,atω=512.
Since by [8], 4-point ternary scheme has cubic precision for ω∈(-5/12<ω<7/12) and quartic at ω=5/12 then 4-point ternary scheme has approximating order 4 for ω∈(-5/12<ω<7/12) and 5 for ω=5/12. The proof of other part is similar.


## 3.2. Support of Basic Limit Function

The basic limit functions,φ[λ+1]=(S[λ+1])∞fi0 of proposed (λ+1)-point ternary approximating schemes, for λ=2n+1, n=1,2 are presented in Figures 1(a) and 1(b). The following theorem is related to the support of basic limit functions of even point ternary schemes.The basic limit functions of proposed schemes atω=1/12.
(a)
4-point scheme
(b)
6-point schemeTheorem 3.2.
The basic limit functionsφ[λ+1] of proposed (λ+1)-point ternary approximating schemes have support width s=(2(λ+n)+3)/2, where λ=2n+1, n⩾0, which implies that it vanishes outside the interval [-(2(λ+n)+3)/4,(2(λ+n)+3)/4].Proof.
The support width “s” of the basic limit functions can be determined by computing how for the effect of the nonzero vertex f00 will propagate along by. As the mask of (λ+1) scheme is 3(λ+1)-long sequence by centering it on that vertex, the distances to the last of its left and right nonzero coefficients are equal to λ+n+2 and λ+n+1, respectively. At the first subdivision step we see that the vertices on the left and right sides of f01 at (λ+n+2)/3 and (λ+n+1)/3 are the furthest nonzero new vertices. At each refinement, the distances on both sides are reduced by the factor 1/3. At the next step of the scheme this will propagate along by ((λ+n+2)/3)×1/3 on the left and ((λ+n+1)/3)×1/3 on the right. Hence after k subdivision steps the furthest nonzero vertex on the left will be at (λ+n+2)(1/3+1/32+⋯+1/3k)=((λ+n+2)/3)(∑j=0k-1(1/3j)) and on the right will be at (λ+n+1)(1/3+1/32+⋯+1/3k)=((λ+n+1)/3)(∑j=0k-1(1/3j)). So the total support width is ((λ+n+2)/3)(∑j=0∞(1/3j))+((λ+n+1)/3)(∑j=0∞(1/3j))=(2(λ+n)+3)/2.


## 4. Comparison and Application

In order to show the performance of the proposed schemes, we compare continuity, support, approximation order, and shape of limit curves. We also discuss curvature and torsion.Comparison of various 4- and 6-point ternary schemes are given in Tables2 and 3. Figure 2 shows visual comparison of 4- and 6-point ternary interpolating and approximating schemes with the proposed ternary approximating schemes.Table 2
Comparison of 4-point ternary schemes.
Scheme
Type
Support
Order
C
n
Ternary 4-point [3]
Interpolating
5
3
2
Ternary 4-point [4]
Interpolating
5
3
2
Ternary 4-point [4]
Interpolating
5
4
1
Ternary 4-point [5, 9]
Interpolating
5
3
1
Ternary 4-point [7]
Approximating
5.5
4
2
Ternary 4-point [14]
Approximating
5
4
3
Ternary 4-point [15]
Approximating
5.5
4
2
Ternary 4-point proposed
Approximating
5.5
5
4Table 3
Comparison of 6-point ternary schemes.
Scheme
Type
Support
Order
C
n
Ternary 6-point [5, 9]
Interpolating
8
5
2
Ternary 6-point [6]
Interpolating
8
4
2
Ternary 6-point proposed
Approximating
8.5
8
7Comparison: bold solid continuous curves are generated by proposed 4-point and 6-point ternary approximating schemes (a) Ko et al. [7] and proposed 4-point ternary schemes (b) [5, 9] and proposed 4-point ternary scheme (c) Hassan and Dodgson [1] and proposed 4-point ternary schemes (d) [5, 9] and proposed 6-point ternary scheme (e) Khan and mustafa [6] and proposed 6-point ternary schemes.
(a)
4-point
(b)
4-point
(c)
4-point
(d)
6-point
(e)
6-point


### 4.1. Curvature and Torsion

The quality of subdivided curves can be assessed quantitatively by measuring curvature and torsion, as functions of cumulative chord length. Curvature is the amount by which a geometric object deviates from being flat, or straight in the case of a line, while in the elementary geometry of space curves, torsion measures the rate at which a twisted curve tends to depart from its osculating plane. When the torsion is zero, osculating plane never changes, and we have a plane curve. We used the method described in [16] to determine the curvature and torsion. The appearance of curvature and torsion obtained by our schemes is measured and compared with the plot or appearance generated by other schemes. The examples in Figure 3 show the results of curvature and torsion, together with their control polygon. The appearance of curvature of the schemes [2, 3, 5, 9] is quite similar and plot of curvature of the schemes [4, 7, 14] is different. Figure 4 shows detailed view of space curve, curvature and torsion generated by existing schemes [2–5, 7, 9] and proposed schemes. Overall one can observe that the plot of space curve, curvature, and torsion of proposed schemes are more smoothing than other existing schemes.Comparison of existing and proposed 4-point schemes sampled from control polygon by using the 5th iteration. The results are shown on the left together with their corresponding curvature on the right.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)Left: control data sampled from a space curve and four iteration of 4-point subdivision scheme. The corresponding curvature and torsion are shown in the center and right column, respectively.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)
(n)
(o)


## 4.1. Curvature and Torsion

The quality of subdivided curves can be assessed quantitatively by measuring curvature and torsion, as functions of cumulative chord length. Curvature is the amount by which a geometric object deviates from being flat, or straight in the case of a line, while in the elementary geometry of space curves, torsion measures the rate at which a twisted curve tends to depart from its osculating plane. When the torsion is zero, osculating plane never changes, and we have a plane curve. We used the method described in [16] to determine the curvature and torsion. The appearance of curvature and torsion obtained by our schemes is measured and compared with the plot or appearance generated by other schemes. The examples in Figure 3 show the results of curvature and torsion, together with their control polygon. The appearance of curvature of the schemes [2, 3, 5, 9] is quite similar and plot of curvature of the schemes [4, 7, 14] is different. Figure 4 shows detailed view of space curve, curvature and torsion generated by existing schemes [2–5, 7, 9] and proposed schemes. Overall one can observe that the plot of space curve, curvature, and torsion of proposed schemes are more smoothing than other existing schemes.Comparison of existing and proposed 4-point schemes sampled from control polygon by using the 5th iteration. The results are shown on the left together with their corresponding curvature on the right.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)Left: control data sampled from a space curve and four iteration of 4-point subdivision scheme. The corresponding curvature and torsion are shown in the center and right column, respectively.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
(j)
(k)
(l)
(m)
(n)
(o)


## 5. Conclusion

The family of even-point approximating schemes for curve design has been established. The 4- and 6-point ternary schemes introduced by Zheng et al. [13, 14] are special cases of our proposed even-point approximating schemes. Smoothness and approximation order of proposed 4- and 6-point ternary schemes have been discussed. Support of family of even-point ternary schemes has been computed in general. It has been shown that proposed schemes are much better than existing ternary schemes in the sense of smoothness, approximation order, curvature, and torsion. The family of odd-point ternary approximating schemes will be studied in detail in the forthcoming paper.



---

*Source: 197383-2012-11-04.xml*