A Variant Cubic Exponential B-Spline Scheme with Shape Control

: This paper presents a variant scheme of the cubic exponential B-spline scheme, which, with two parameters, can generate curves with different shapes. This variant scheme is obtained based on the iteration from the generation of exponentials and a suitably chosen function. For such a scheme, we show its C 2 -convergence and analyze the effect of the parameters on the shape of the generated curves and also discuss its convexity preservation. In addition, a non-uniform version of this variant scheme is derived in order to locally control the shape of the generated curves. Numerical examples are given to illustrate the performance of the new schemes in this paper.


Introduction
Subdivision can be used as an effetive tool to generate smooth curves/surfaces starting from a given set of initial control points. Due to their advantages such as numerical stability and refinability, subdivision schemes have been widely used in fields such as geometric modeling [1], wavelets and framelets [2], isogeometric analysis [3], and multigrid [4]. Generally speaking, subdivision can be classified into stationary and non-stationary schemes. Compared with stationary schemes, non-stationary subdivision depend on the recursion level and thus have more flexibility to generate different curves. Recall that the stationary subdivision are the schemes whose subdivision rules are independent of the recursion level while the non-stationary subdivision are the ones whose subdivision rules depend on the recursion level. This means that, in the subdivision process, as the recursion level changes, the subdivision rules of the stationary schemes keep the same while the subdivision rules of the non-stationary schemes change.
Therefore, there have been amounts of works on the construction of non-stationary subdivision. Beccari. et al. [5] derived the non-stationary 4-pt interpolatory scheme reproducing conics. Dyn et al. [6] investigated the non-stationary interpolatory schemes reproducing exponential polynomials. Fang [7] gave a family of generalized curve subdivision generating different kinds of curves. All of these schemes can be seen as obtained from the view of generating exponentials and thus can produce curves such as conics and helix. For other references on this kind of non-stationary schemes, please refer to [8][9][10][11][12][13] and the references therein. For such schemes, the interpolatory ones are usually more powerful than the approximating ones in controlling the shape of the limit curves while the approximating ones can have higher smoothness order than the interpolatory ones.
In fact, besides such non-stationary subdivision, there is another kind of non-stationary schemes, which can be seen as obtained based on some fixed point iteration. For this kind of schemes, Beccari [14] constructed a ternary interpolatory non-stationary scheme generating various kinds of curves. Tan et al. [15] presented a binary approximating nonstationary scheme with a similar property. Compared with the ones obtained in the view of generating exponentials, the approximating ones based on iterations can be as powerful as the interpolatory ones in controlling the shape of the limit curves without reducing the smoothness order. In fact, similar to [14,15], Zhang et al. [9] also proposed a modified cubic exponential B-spline scheme and a non-stationary Catmull-Clark scheme controlling the shape of the curves and surfaces freely. Yet, the analysis of the parameters on the shape of the generated curves/surfaces and the convexity preservation is absent. Note that the analysis of these properties is also important in applications, since it is better that we know how to adjust the parameters to generate the curves with the needed shape. In fact, for the convexity preservation of non-stationary schemes, Akram et al. [16] analyzed the shape preservation of the non-stationary 4-pt interpolatory scheme [5].
Therefore, in this paper, continuing our work in [9], we base on the iteration from the generation of exponentials and a suitable function of this iteration, and give a variant scheme of the cubic exponential B-spline scheme, which can generate different kinds of curves, including the conics. For this scheme, we analyze its C l -convergence and show that it has the same smoothness order as the cubic exponential B-spline scheme. Then, we give the limit position of the initial control points and analyze the convexity preservation. This analysis can show us how to control the shape of the limit curves and preserve the convexity of the initial control polygon by adjusting the parameters. Furthermore, a nonuniform version of this variant scheme is also given to locally control the shape of the limit curves. Numerous examples are given to illustrate the performance of the new schemes.
The rest of this paper is organized as follows. Section 2 gives a brief review of nonstationary subdivision. In Section 3, we present the variant cubic exponential B-spline scheme and show that it can indeed generate curves with different shapes. Section 4 is devoted to the analysis of this variant cubic exponential B-spline scheme, including the smoothness analysis, the limit position of the initial control points and the convexity preservation. In Section 5, we move a further step and give a discussion on the local control of the shape of the limit curves to generate more kinds of curves. Section 6 concludes this paper.

Preliminaries
In this section, we present some basic knowledge about the non-stationary subdivision schemes, which is needed in the rest of this paper.
Let l(Z) denote the linear space of real sequences indexed by Z and l 0 (Z) denote the linear space of real sequences with finite support. Given an initial data sequence q 0 = {q 0 i , i ∈ Z}, we consider the following binary non-stationary subdivision where S a k is the k-level subdivision operator, the sequence a k = {a k i , i ∈ Z} is the k-level mask with finite support, and the k-level symbol is the Laurent polynomial a k (z) = ∑ i∈Z a k i z i . We denote this non-stationary subdivision by {S a k } k≥0 .
By attaching q k i , i ∈ Z to the parameter values i 2 k , k ≥ 0, we say the scheme {S a k } k≥0 is C 0 -convergent if, for the initial data sequence q 0 ∈ l 0 (Z), there exists a function f q 0 ∈ C 0 (R) satisfying with K being any compact set in R, and f q 0 = 0 for at least some initial data sequence q 0 . If For the convergence and smoothness of the non-stationary subdivision in this paper, we recall the following definitions.

Definition 1 ([17]).
A non-stationary scheme {S u k } k≥0 with the k-level mask u k is said to be asymptotically similar to the stationary scheme S u with the mask u, if the k-level mask u k and the mask u have the same support U and satisfy

Definition 2 ([18]).
A non-stationary scheme {S u k } k≥0 with the k-level symbol u k (z) is said to satisfy the approximate sum rules of order r + 1 if Now we cite the result on the smoothness of non-stationary subdivision for the C lconvergence of the variant cubic exponential B-spline scheme in this paper as follows.
Theorem 1 ([17]). Assume that the non-stationary scheme {S u k } k≥0 satisfies the approximate sum rules of order r + 1, and is asymptotically similar to a convergent stationary scheme S u , which is C r -convergent. Then the non-stationary scheme {S u k } k≥0 is C r -convergent.

The Variant Cubic Exponential B-Spline Scheme
The main purpose of this section is to derive the variant cubic exponential B-spline scheme in this paper. To this purpose, we start from reviewing the cubic exponential B-spline [19] and the non-stationary 4-pt interpolatory scheme [5].
As is known, the cubic exponential B-spline scheme [19] and the non-stationary 4-point interpolatory scheme [5] can be characterized by their k-level masks 2 .
It can be seen that, v k and v k+1 satisfies the iteration , then, the k-level masks in (1) and (2) can be rewritten as respectively. Both schemes generate conics but the cubic exponential B-spline scheme is C 2 convergent while the non-stationary 4-pt interpolatory scheme is C 1 convergent. In addition, as seen in [9], for h 0 (x) and h 1 (x), combining with the iteration in (3), we have and It follows that, with the functions h 0 (x) and h 1 (x) behaving differently as shown in (4) and (5), the non-stationary 4-pt interpolatory scheme can generate curves with more kinds of shapes than the exponential cubic B-spline scheme [9]. In this way, similar to [9], based on the iteration in (3), we choose a function such that it satisfies (4) and (5) with i = 1. Here, we can choose h(x) as the following one Now with the iteration in (3) and the function h(x), we can derive the following k-level mask We denote the corresponding scheme by {S a k } k≥0 , which is just the desired variant cubic exponential B-spline scheme. It can be seen that, when β = 0, h(x) reduces to h 0 (x) and the scheme {S a k } k≥0 reduces to the cubic exponential B-spline and, when β > 0, h(x) behaves like h 1 (x) as in (4) and (5), and thus the scheme {S a k } k≥0 can generate curves with different shapes. Figures 1 and 2 show the curves generated by the scheme {S a k } k≥0 with different choices of v 0 and β. It can be seen from these two figures that, with different choices of the parameters v 0 and β, the scheme {S a k } k≥0 can indeed generate curves with different shapes. In particular, when β = 0, the scheme {S a k } k≥0 reduces to the cubic exponential B-spline scheme and Figures 1 and 2 show the generation of circles with v 0 = cos(π/3) and β = 0.

Analysis of the Variant Cubic Exponential B-Spline Scheme
This section is devoted to the analysis of the properties of the scheme {S a k } k≥0 , including the C l -convergence, the limit position of the initial control points and the convexity preservation.

Smoothness Analysis
We now analyze the C l -convergence of the variant scheme {S a k } k≥0 . In fact, this can be done similar to [9]. Yet, for completeness, we prove the following result.

Proof.
The k-level symbol of the scheme {S a k } k≥0 can be written as Recall that for the sequence {v k } k≥0 , we have lim k→∞ v k = 1, for v 0 ∈ (−1, ∞). Therefore, lim k→∞ h(v k+1 ) = 1 8 and the scheme {S a k } k≥0 is, by Definition 1, asymptotically similar to the cubic B-spline scheme, which is C 2 -convergent. If we can show that the scheme {S a k } k≥0 satisfies the approximate sum rules of order 3, then by Theorem 1, the scheme {S a k } k≥0 is C 2 -convergent.
In fact, from (7), it can be computed that a k (1) = 2 and a k (−1) = D 1 a k (−1) = 0, D 2 a k (−1) = 8h(v k+1 ) − 1. Thus, for µ k and δ k in Definition 2, it follows that, for η ∈ N 0 , In this way, following [9], by Definition 2 and the function h(x) in (6), there exists a constant c 1 independent of k, such that Therefore, according to Definition 2, we only have to show Let v k+1 = ϕ(v k ) = v k 2 + 1 2 with ϕ(x) = x 2 + 1 2 , then ϕ(1) = 1 and thus 1 is the fixed point of the iteration v k+1 = ϕ(v k ). Then, as shown in [9], for v 0 ∈ (−1, ∞), with 0 < L < 1 and c 2 is a constant independent of k. Thus, where c 3 is a constant independent of k. Therefore, we can conclude that the scheme {S a k } k≥0 satisfies the approximate sum rules of order 3 and, thus, is C 2 -convergent.

The Limit Position of the Initial Control Points
Now we give the limit position of the initial control points, by which we can discuss the effect of the parameters on the shape of the limit curves. To this purpose, we give the following result.

Theorem 3.
In the subdivision process of the scheme {S a k } k≥0 , the limit position of the initial control point P 0 i is Proof. For the point P 0 i , according to the subdivision rules of the scheme {S a k } k≥0 , after one step of subdivision, it is moved to Similarly, after another step of subdivision, the point P 1 2i is moved to This can be seen in Figure 3. Generally, we can show that after k + 1 steps of subdivision using the scheme {S a k } k≥0 , we actually have For the displacement ∆P k 2 k i , note that combining with the subdivision rules of the scheme {S a k } k≥0 , it follows that Generally, we have, for k ∈ N, Therefore, after k + 1 steps of subdivision, the initial control point P 0 i is moved to Recall that as k tends to ∞, v k tends to 1. Thus h(v k ) tends to 1 8 , 1 2 − 2h(v k ) tends to 1 4 . Together with h(x) in (6), it can be shown that the term Then it can seen that ( 1 2 − 2h(v s+1 )) converges as k → ∞. Therefore, the limit position of P 0 i , which can be denoted by P ∞ i , is which is just (8). Now, denote the displacement between the initial control point P 0 i and its limit position by Then Theorem 3 actually implies the effect of the parameters v 0 and β on the shape of the limit curve. In fact, from Theorem 3, it can be seen that, for given β, with a larger v 0 , each term of ∆ gets smaller and thus the limit point P ∞ i gets closer to P 0 i , which makes the limit curve shrinks to the initial control polygon. Similarly, for given v 0 > 1, with a larger β, the limit point P ∞ i also gets closer to P 0 i and, thus, the limit curve shrinks to the initial control polygon. For v 0 = 1, (8) actually becomes which actually implies the limit position of the initial control points using the cubic B-spline scheme. In addition, for given −1 < v 0 < 1, when β increases, the limit point P ∞ i gets far away from P 0 i , which implies the big change of the shape of the limit curve. The above analysis can be seen in Figure 4, which shows the curves generated by the scheme {S a k } k≥0 with different choices of the parameters v 0 and β.

Convexity Preservation
Now we analyze the convexity preservation of the scheme {S a k } k≥0 and show how to choose the parameters to preserve the convexity of the initial control polygon.
Given a set of control points {p k i } i∈Z , i.e., these points are equally spaced, and suppose they are strictly convex, which means that the second divided difference d 0 Then for the k level control points, they are also equally spaced, meaning that ) be the second order divided differences and let r k Then, we have the following result.

Theorem 4. Suppose the initial control points {p
Proof. We first show d k i > 0. In fact, according to the subdivision rules of the scheme {S a k } k≥0 , it follows that Similarly, we have Therefore, since d 0 i > 0 and h(v k ) > 0, by induction, it follows that d k i > 0. Now we show R k < λ by showing r k i , 1 In fact, from (9) and (10), we have Therefore, it follows that and r k+1 from (11) and (12), we have r k+1 Note that r k i , 1 r k i < λ also implies 1 λ < r k i . Then, similar to (11) and (12), if (13) is valid, we have   Figure 5, it can be seen that, although the scheme {S a k } k≥0 can generate curves with different shapes, with certain choice of v 0 and β (see Figure 5 with v 0 = −0.35 and β = 2 ), the curve generated by the scheme {S a k } k≥0 may not be convex. However, Theorem 4 implies that we can rechoose β to guarantee that the limit curve is convex. This can be seen in Figure 5 with v 0 = −0.35 and β = 1.2, 0.8, 0.1, 0.

Discussion on the Local Control
Now we move a further step and give a non-uniform version of the scheme {S a k } k≥0 to locally control the shape of the limit curves.
In fact, there have been works [5,21,22] on the local control, which assigned a parameter to each edge of the initial control polygon. Different from these schemes, the work in [9] assigned a parameter to each initial control point, which can also locally control the limit curves.
Here, following [9], we also assign each initial control point a parameter β i and then we can derive a non-uniform version of the scheme {S a k } k≥0 . Following the work in [9], we can also show the C 2 -convergence of this non-uniform scheme. Figure 6 shows the curves generated by the non-uniform version of the variant scheme {S a k } k≥0 with v 0 = −0.8 and different choices of β i assigned to the i-th point. From Figure 6, it can be seen that by suitably adjusting the parameters assigned to the initial control points, we can indeed locally control the shape of the limit curve. We point out that, similar to the above non-uniform scheme, we can also assign a parameter v 0 i to each initial control point to derive another non-uniform scheme in order to locally control the shape of the limit curves.

Conclusions
This paper presented a variant cubic exponential B-spline scheme, which can generate curves with different shapes, including the conics. Such a scheme is obtained based on the iteration from the generation of exponential polynomials and a suitable function. For such a scheme, we verified its C 2 -convergence and gave the limit position of the initial control points to show the effect of the parameters on the shape of the limit curves. Apart from these properties, we also discussed the convexity preservation of this variant scheme and showed how to preserve the convexity of the initial control polygon by adjusting the parameters. In addition, a non-uniform generalization of this variant scheme is also given, which can locally control the shape of the limit curves.