Further Properties of Quantum Spline Spaces

We construct q-B-splines using a new form of truncated power functions. We give basic properties to show that q-B-splines form a basis for quantum spline spaces. On the other hand, we derive algorithmic formulas for 1/q-integration and 1/q-differentiation for q-spline functions. Moreover, we show a way to find the polynomial pieces on each interval of a q-spline function.


Introduction
B-splines were constructed by Lobachevsky as convolutions of certain probability distributions in the early 19th century. Spline functions are used in various numerical analysis areas like interpolation, approximation, computer aided geometric design, numerical solutions of differential equations, etc. The modern theory of spline approximation was started by Schoenberg in 1946, when he used B-splines for statistical data smoothing, see [1]. De Boor gave the recurrence relation for B-splines in [2]. Gordon and Reisenfield formally introduced B-splines into computer aided design in [3].
Mangasarian and Schumaker introduced discrete splines, h-splines, to solve discrete analogoues of minimization problems in a Banach space (see [4]). These discrete splines of degree n are defined on a subset of real line of the form [a, b] h = {a, a + h, . . . , a + Nh}, b = a + Nh, whose knot sequence is in [a, b] and their polynomial pieces agree at the knots up to the order n − 1 of forward differences with step size h instead of derivatives. In [5], q-splines which allow us to model tolerances, jumps and quantum leaps in the derivatives at the joins, were defined recursively based on a q-analogue of the de Boor algorithm. After giving certain properties, they defined blossoms for q-B-splines. In [6], fundamental formulas of classical B-splines were extended to q-B-splines. The q-splines are piecewise polynomials whose q-derivatives up to some order agree at the joins. A recent study relates q-B-splines with the q-Peano kernels of divided differences and solves a best approximation problem in the space of quantum integrable functions, see [7].
Let t 0 < t 1 < · · · < t m be the knots, and let n be a nonnegative integer and q = 0 be a fixed real. A q-spline function of degree n having knots t 0 < t 1 < · · · < t m is a function S such that (i) S is a polynomial of degree up to n on each interval [t i−1 , t i ) , i = 1, 2, . . . , m.
(ii) S is quantum continuous of order n − 1 at the knots.
Then, S is a continuous piecewise polynomial of degree at most n and quantum continuous of order n − 1 at the knots. Here, "quantum continuous of order n − 1 at the knots" means the quantum derivatives D k 1/q , k = 0, . . . , n − 1 of adjacent polynomial segments agree at the knots. The rest of this paper is organized as follows. In Section 2, we begin with definitions and theorems in q-calculus concerning this work. In Section 3, we give some properties of q-B-splines and give proofs which are not stated explicitly in [5]. In Section 4, we find a new basis for the q-spline spaces with a given knot sequence and degree. We define q-B-splines in a different manner from [5,6] but similar to that of Curry and Schoenberg [8] in which the truncated power function played a significant role. Our approach is based on a certain q-truncated power function rather than the recursive definition of q-B-splines. We show how we can obtain quantum derivatives and quantum integrals of a given q-spline function. Furthermore, we find the polynomial pieces of a q-spline by using quantum derivatives.

Preliminaries
For a fixed parameter q = 1, the q-derivatives are defined by Indeed q-derivatives are approximations to classical derivatives and, if f is a differentiable function, then lim For polynomials, it follows from the definition of the q-derivative that where the q-integers [n] q are defined by Furthermore, the q-factorial is defined by The next definition states the q-analogues of a classical definite integral; for details, one may see [9].

Properties of q-B-Splines
The q-B-splines are introduced first by Simeonov and Goldman in [5]. We give from [5,6] the important properties of q-B-splines that will be used in this work. Throughout the paper N k,n (t; q) refers to q-B-spline basis functions.
where (x − t) n,q Property 4. For any integer n 2, we have Property 5. Partition of unity:

Remark 1.
The classical B-spline N k,n (t) is nonnegative for each k, n and all t ∈ [t k , t k+n+1 ]. However, Figure 1 shows that q-B-splines N k,n (t; q) may be nonnegative depending on the value of parameter q. The graphs of N 0,2 (t; q) with the knot sequence (1, 2, 3, 4) and q = 2, q = 1.5, q = 1.1 respectively, whose first order quantum derivatives agree at knots.

Remark 2.
If q = 1, it is obvious that q-B-splines become classical B-splines and they mimic the classical counterpart when q is near 1. The following Lemmas 1, 2, and Proposition 1 are noted in [5]. Here, we prove them by using induction. We use the notation N q k,n in the place of N k,n (t; q) to emphasize its dependence on the initial parameter value q. Lemma 1. For n 1, the q-B-splines N q k,n belong to continuity class q − C n−1 (R) that is quantum derivatives of N q k,n agree up to the order n − 1 at the joins.
The following result shows that the sequence of q-B-splines of the same degree is linearly independent in a single interval of its support.
Proof. When n = 0, it is obvious that the set is linearly independent. Let n 1 and assume that the Lemma 2 holds for n − 1. Let S = n ∑ k=0 c i+k N q i+k,n and suppose that the q-spline S restricted to the interval (t i+n , t i+n+1 ), S|(t i+n , t i+n+1 ) = 0. By Equation (3), we have Since N i+n+1,n−1 (t; q) = 0 and N i,n−1 (t; q) = 0 on (t i+n , t i+n+1 ), by induction hypothesis {N q i+1,n−1 , . . . , N q i+n,n−1 } is linearly independent on the interval (t i+n , t i+n+1 ). Therefore, we must have that all the coefficients in Equation (4) are zero, that is, c i = c i+1 = · · · = c i+n , say all of them equal to the value c. Thus, by the partition unity property of q-B-splines, we have S(t; q) = c on (t i+n , t i+n+1 ). Hence, by the assumption S|(t i+n , t i+n+1 ) = 0, we have c = 0 which shows that the set of q-B-splines {N q i,n , N q i+1,n , . . . , N q i+n,n } is linearly independent on an interval (t i+n , t i+n+1 ).
From the previous lemma, the set {N q i−n,n , N q i−n+1,n , . . . , N q i,n } is linearly independent on (t i , t i+1 ). Hence, the coefficients c k = 0 for i − n k i in Equation (5). If all the c k 's are zero, then there is nothing to prove. Suppose, on the contrary, that not all the c k 's are zero. Let j be the index such that c j = 0. Assume 0 j m − 1 and (t j , t j+1 ) ⊂ (t 0 , t m ). For any t ∈ (t j , t j+1 ), we obtain which contradicts S|(t j , t j+1 ) = 0. Therefore, all the c k 's are zero and this implies that the set {N q −n,n , N q −n+1,n , . . . , N q m−1,n } is linearly independent on (t 0 , t m ).

The Quantum Spline Space S n,q m
Let S n,q m denote the space of the q-spline functions which are quantum continuous up to the order n − 1, and n denotes the degree of polynomial pieces, m + 1 is the number of knots in the knot sequence, and q is a nonzero initial parameter.
The next theorem shows that q-B-splines form a basis for the quantum spline space S n,q m of degree n with the knot sequence t 0 < t 1 < · · · < t m .
We note that, although the work [5] investigates broadly q-B-splines via blossoming, the following theorem and its proof were not mentioned.

Theorem 2. A basis for the space
Consequently, the dimension of S Since in the interval [t 0 , t 1 ] all the truncated powers vanish, S n m (t; q) is a polynomial of degree n, say β 0 . Thus, we have p 0 (t; q) = n ∑ i=0 a i t i which determines all the coefficients a i . In the interval [t 1 , t 2 ], S n m (t; q) is another polynomial, say β 1 . According to the quantum continuity at the knots, we have at t 1 , When we repeat the same argument at the other internal knots t 2 , t 3 , . . . , t m−1 , we obtain the q-spline in Equation (6). Thus, any spline of degree n on the interval [t 0 , t m ] with intermediate knots t 1 , t 2 , . . . , t m−1 may be written as a sum of multiples of n + m functions 1, t, t 2 , . . . , t n , (t − t 1 ; q) n + , (t − t 2 ; q) n + , . . . , (t − t m−1 ; q) n + .
Since these functions are linearly independent, they form a basis for the this spline space and hence the dimension of the space S n,q m is n + m. Furthermore, it follows from Proposition 1 that q-B-splines N q i,n |[t 0 , t m ] with −n i m − 1 form a basis for the space S n,q m .
In [5], q-B-splines are constructed by using the de Boor algorithm. In the following theorem, we construct q-B-splines using properties of truncated power functions depending on q. Namely, we give explicit expression of q-B-spline basis functions in terms of linear combinations of q-truncated power functions.
with respect to the normalization ∑ k N k,n (t; q) = 1.
Proof. Let {Ψ j : j = 1, 2, . . . , n + m} be a basis in the expression such that each function {Ψ j (t; q) : j = 1, 2, . . . , n + m} is identically zero over a large part of the range a t b. Consider an element of the space S(n, t 0 , . . . , t m ) that is zero on the intervals [t 0 , t k ] and [t r , t m ], but that is nonzero on (t k , t r ), where 0 < k < r < m. If Φ is such a function, it can be expressed in the form where the parameters µ j have to satisfy the condition r ∑ j=k µ j (q n−1 t − t j ) · · · (qt − t j )(t − t j ) = 0, t r t b since (t − t j ) + = 0 for j = k, k + 1, . . . , r and a t t k .
By rearranging the terms and using the properties of Lagrange polynomial interpolation, we have r ∑ j=k µ j t i j = 0, i = 0, 1, . . . , n.
If r k + n + 1, then equations in (7) have a nonzero solution. If r = k + n + 1, then the coefficients are where c is a nonzero constant. Thus, we can conclude that By applying the normalization constraint ∑ k N k,n (t; q) = 1, we obtain the q-B-splines A consequence of the last expression is the following: Proof. This follows from the following property of divided differences by replacing f by the truncated power function . We have shown that q-B-splines form a basis for q-splines. Now, let us investigate how we can find the quantum derivatives D j 1/q of the q-spline functions. The following algorithms allow us to store, evaluate, and manipulate q-splines on a computer easily, see [10]. In addition, we show that 1/q-derivatives of q-splines are again q-spline functions.
where a (j) Proof. For j = 1, we have This shows that Equation (11) is true for j = 1. Repeating the same argument shows that the formula (11) is true for each j = 1, 2, . . . , n − 1. Now, we derive a formula for computing the 1/q-integral of a given q-spline function.
Theorem 5. Let S be a given q-spline function such that S(t; q) = Proof. Define a functionS byS This is a quantum spline of degree n + 1 and soS can be written as Using (10) and (11) for all t ∈ [a, b] gives whereã −m−1 = 0. It follows by comparing the appropriate coefficients of the basis that Therefore,ã i = i ∑ j=−n a j t j+m+1 − t j [n + 1] q for i = −n, −n + 1, . . . , k and this completes the proof.
In the next proposition, we demonstrate a way to find the polynomials on each interval of a q-spline function. In certain circumstances, we may have to evaluate the q-spline function S at a large number of points, then it is advantageous to determine the polynomial pieces at each interval P j (t; q) = S| [t j ,t j+1 ) for each j. Proposition 2. Let S(t; q) be a q-spline function such that and P j be the polynomial restricted to S| [t j ,t j+1 ) for j = 0, . . . , k. Then, P j (t; q) = n ∑ r=0 1 [r] q ! D r 1/q S(t j ; q) (q r−1 t − t j ) · · · (qt − t j )(t − t j ), t ∈ [t j , t j+1 ) Proof. Since S(t; q) is a polynomial of degree n on the interval [t j , t j+1 ) for each j = 0, 1, . . . , k, we can write where λ rj for r = 0, 1, . . . , n and j = 0, 1, . . . , k are the unknowns and (t − t j ; q) r = (q r−1 t − t j ) · · · (qt − t j )(t − t j ).
Then, it is easily seen that, by taking 1/q-derivative of S(t; q) r times, substituting t = t j in Equation (10) gives λ rj = 1 [r] q ! D r 1/q S(t j ; q).

Conclusions
In this work, we derived a new way to compute q-B-spline N k,n (t; q) and algorithms for q-derivatives and q-integrals of a q-spline function, and found the polynomial pieces on a specified single interval. These formulas will be useful in numerical methods incorporating a finite difference scheme with q-B-splines. For CAGD purposes, we take the parameter values q near 1 in q-B-spline curves so as to mimic the behavior of their classical counterpart. Besides solving a best approximation problem, not only q-B-splines but also h-splines share certain fundamental properties of the classical B-splines. In a future work, we will study how we can solve partial q-differential equations by using q-splines. We will introduce box q-splines as well as their rational counterpart and then investigate them in solving certain PDE problems towards isogeometric analysis.