# Positivity Preserving Interpolation Using Rational Bicubic Spline

Authors: Samsul Ariffin Abdul Karim; Kong Voon Pang; Azizan Saaban
Journal: Journal of Applied Mathematics
 (2015)

---

## Abstract

This paper discusses the positivity preserving interpolation for positive surfaces data by extending theC
1 rational cubic spline interpolant of Karim and Kong to the bivariate cases. The partially blended rational bicubic spline has 12 parameters in the descriptions where 8 of them are free parameters. The sufficient conditions for the positivity are derived on every four boundary curves network on the rectangular patch. Numerical comparison with existing schemes also has been done in detail. Based on Root Mean Square Error (RMSE), our partially blended rational bicubic spline is on a par with the established methods.


---

## Body

## 1. Introduction

Shape preserving interpolation is a process that involves some mathematical derivation to visualize the two-dimensional (2D) or three-dimensional (3D) data sets in the form of piecewise curves or surfaces that satisfies some degree of smoothness. Notably the shape preserving method will be able to maintain the original shape of the data sets. For example, if the given data is positive, then the final resulting curves or surface will be positive everywhere and it should be visual pleasing for the computer graphics displays. Usually the collected data from some laboratory experiment may contain little unwanted noise. Common cubic spline interpolation is not able to produce the shape preserving interpolation for positive, monotone and convexity of the data sets. Even though the interpolating curves are very smooth, that is, withC
2 continuity, there may exist some uncharacteristic behavior such as wiggle that will result in the negative values along the whole intervals. This unwanted flaw must be eliminated before the data are processed by the user. Therefore for the shape preserving interpolation and practical designing purposes, cubic spline and cubic Hermite spline are not recommended. Thus many researchers have proposed the idea of how to preserve the positive data by using some rational spline function with C
1 or C
2 continuity.Some literature studies for positivity preserving are as follows. Abbas et al. [1, 2] have discussed the positivity preserving by using partially blended rational bicubic spline function and rational bicubic spline function, respectively. Asim and Brodlie [3], Butt and Brodlie [4], and Brodlie et al. [5] discussed the use of bicubic spline for positivity preserving by inserting one or two extra knots in the region or interval where the negativity of the interpolant is found. Asim et al. [6] have studied the use of modified quadratic Shepard method for positivity preserving of large data sets. Duan et al. [7] have proposed new bivariate rational interpolation based on function values. They utilized the rational cubic spline of the forms of cubic numerator and linear denominator. Meanwhile, Duan et al. [8, 9] investigated the bounded property, point control, and shape control of the bivariate rational cubic spline with linear denominator. Goodman et al. [10] discussed the constrained interpolation by using rational cubic with cubic denominator. Goodman [11] summarized all shape preserving criteria for the univariate cases. Hussain and Sarfraz [12] and M. Z. Hussain and M. Hussain [13] also studied the positivity preserving by utilizing the rational cubic spline of the forms cubic/cubic and cubic/quadratic, respectively. In the study by M. Z. Hussain and M. Hussain [13], there are no free parameters for the refining processes. M. Z. Hussain and M. Hussain [14] have developed the surface construction scheme for positivity preserving of positive scattered data arranged over triangular grids. They utilized the rational cubic spline of Tian et al. [15] without any free parameters. Hussain et al. [16] have developed shape preserving surfaces for positive and convex data using rational biquadratic splines. The degree of smoothness attained is C
0 as can clearly be seen from their graphical results. Ibraheem et al. [17] develop a positivity preserving scheme by using trigonometric rational functions which will increase the computation times. Karim and Kong [18] have proposed new rational cubic spline of the form cubic/quadratic with three parameters for local control of interpolating function. Positivity preserving by using this rational cubic spline has been addressed in the study by Karim and Kong [19]. Liu et al. [20] have studied the positivity and monotonicity preserving by using new rational biquartic spline by utilizing the extension cubic Ball basis function with two parameters. Saaban et al. [21] studied the positivity preserving by using quintic triangular Bezier patches. The surfaces are constructed by using convex combination of quintic triangular Bezier patches. They use the real data collected from the rainfall at various stations in West Peninsular of Malaysia. Sarfraz et al. [22] also have developed rational cubic spline scheme for positivity preserving of 3D data. Their scheme has 8 parameters without any free parameter. Furthermore Abbas et al. [1] indicated that Sarfraz et al. [22] works do not produce C
1 positive surfaces. The surfaces may only have C
0 continuity. Bastian-Walther and Schmidt [23] discussed the range restricted data interpolation for scattered data by using rational cubic spline of Gregory [24]. Wu et al. [25] studied the positivity preserving data approximation by using compactly supported radial basis functions (CSRBFs). Their scheme requires the solving of optimization problems in order to calculate the parameters that will guarantee that the positivity of the data is preserved. Peng et al. [26] have used the rational cubic spline of Hussain and Sarfraz [12] to construct the positive rational cubic Bezier interpolant for scattered positive data. Their methods require the modification of the first derivatives if the positivity of the rational interpolant is not met. The bicubic Hermite spline [27] interpolation cannot preserve the positivity of the surface data sets, since there maybe exist some negative values on the constructed surfaces. This paper is a continuation of our previous works. The univariate rational cubic spline interpolation in the studies by Karim and Kong [18, 19] is extended to the bivariate cases. Positivity preserving interpolation is our main object in this paper. We compare the performance between our rational bicubic spline and that of the work of Abbas et al. [1] by using Root Mean Square Error (RMSE).We identify several advantages of our proposed partially blended rational bicubic spline interpolation for positive data as follows.(i)
Our schemes have 12 parameters and 8 of them are free parameters; meanwhile, in the study by M. Z. Hussain and M. Hussain [13] the bivariate spline has 8 parameters and none of them are free parameters.
(ii)
Our schemes are applicable for both equally or unequally spaced data; meanwhile, the rational cubic spline of Duan et al. [7–9] requires that the data be equally spaced.
(iii)
The scheme can also be used whether the first derivative is given or not. Meanwhile, the works of Duan et al. [7, 9] require the true function values without the first derivative values.
(iv)
Our rational bicubic spline is on a par with the works of Abbas et al. [1, 2] and Hussain and Sarfraz [12]. Furthermore the presented scheme is easy to use and not involving any trigonometric functions compared with the work of Ibraheem et al. [17].
(v)
Whenγ
i
,
j
=
γ
^
i
,
j
=
0, our partially blended rational bicubic spline is reduced to the work of M. Z. Hussain and M. Hussain [13].
(vi)
The partially blended rational bicubic spline is unique and the free parameters in its description provide the flexibility in controlling the final shape of the positive surface data sets.
(vii)
Furthermore based on Root Mean Square Error (RMSE) our partially blended rational bicubic spline gives numerical results as good as the numerical results from the study by Abbas et al. [1].The remainder of the paper is organized as follows. Section 2 is devoted to the derivative estimation for 2D and 3D data by using Arithmetic Mean Method (AMM). Meanwhile, Section 3 reviews the univariate of rational cubic spline initiated by Karim and Kong [18, 19]. Positivity preserving for 2D data also is discussed in this section. The partially bicubic rational cubic spline interpolation is discussed in Section 4. Section 5 discusses the positivity preserving for 3D positive data sets by using the constructed partially bicubic rational cubic spline interpolation with numerical demonstrations with comparison to the existing schemes. Finally a summary and conclusions are given in Section 6.


## 2. Derivative Estimation

The first derivative that depends on the original data can be chosen by using Arithmetic Mean Method (AMM). The details are as follows:For 2D Data
(1)
d
0
=
Δ
0
+
Δ
0
-
Δ
1
h
0
h
0
+
h
1
,
d
n
=
Δ
n
-
1
+
Δ
n
-
1
-
Δ
n
-
2
h
n
-
1
h
n
-
1
+
h
n
-
2
.At the interior points,x
i, i
=
1,2
,
…
,
n
-
1, the values of d
i are given as (2)
d
i
=
h
i
-
1
Δ
i
+
h
i
Δ
i
-
1
h
i
-
1
+
h
i
.For 3D data
(3)
F
0
,
j
x
=
Δ
0
,
j
+
Δ
0
,
j
-
Δ
1
,
j
h
1
h
1
+
h
2
,
F
n
,
j
x
=
Δ
n
-
1
,
j
+
Δ
n
-
1
,
j
-
Δ
n
-
2
,
j
h
n
-
1
h
n
-
1
+
h
n
-
2
,
F
i
,
j
x
=
Δ
i
,
j
+
Δ
i
-
1
,
j
2
,
i
=
1,2
,
…
,
n
-
1
,
j
=
0,1
,
…
,
m
,
F
i
,
0
y
=
Δ
^
i
,
0
+
Δ
^
i
,
0
-
Δ
^
i
,
1
h
^
1
h
^
1
+
h
^
2
,
F
i
,
m
y
=
Δ
^
i
,
m
-
1
+
Δ
^
i
,
m
-
1
-
Δ
^
i
,
m
-
2
h
^
m
-
1
h
^
m
-
1
+
h
^
m
-
2
,
F
i
,
j
y
=
Δ
^
i
,
j
+
Δ
^
i
,
j
-
1
2
,
i
=
0,1
,
2
,
…
,
n
,
j
=
1,2
,
…
,
m
-
1
,
F
i
,
j
x
y
=
F
i
+
1
,
j
y
-
F
i
-
1
,
j
y
h
i
-
1
+
h
i
+
F
i
,
j
+
1
x
-
F
i
,
j
-
1
x
h
^
j
-
1
+
h
^
j
,
i
=
1,2
,
…
,
n
-
1
,
j
=
1,2
,
…
,
m
-
1
,where(4)
h
i
=
x
i
+
1
-
x
i
,
h
^
j
=
y
j
+
1
-
y
j
,
Δ
i
,
j
=
f
i
+
1
-
f
i
h
i
,
Δ
^
i
,
j
=
f
i
,
j
+
1
-
f
i
,
j
h
^
j
.


## 3. Review of Rational Cubic Spline Interpolant

In this section the rational cubic spline of Karim and Kong [18, 19] is discussed briefly before it extends to the bivariate cases.Given the set of data points{
(
x
i
,
f
i
)
,
i
=
0,1
,
…
,
n
} such that x
0
<
x
1
<
⋯
<
x
n, let h
i
=
x
i
+
1
-
x
i, Δ
i
=
(
f
i
+
1
-
f
i
)
/
h
i, and θ
=
(
x
-
x
i
)
/
h
i, where 0
≤
θ
≤
1. For x
∈
[
x
i
,
x
i
+
1
], i
=
0,1
,
2
,
…
,
n
-
1,(5)
s
x
≡
s
i
x
=
P
i
θ
Q
i
θ
,where(6)
P
i
θ
=
A
0
1
-
θ
3
+
A
1
θ
1
-
θ
2
+
A
2
θ
2
1
-
θ
+
A
3
θ
3
,
Q
i
θ
=
1
-
θ
2
α
i
+
θ
1
-
θ
2
α
i
β
i
+
γ
i
+
θ
2
β
i
.The rational function in (5) satisfies C
1 continuity as follows:(7)
s
x
i
=
f
i
,
s
x
i
+
1
=
f
i
+
1
,
s
i
1
x
i
=
d
i
,
s
i
1
x
i
+
1
=
d
i
+
1
.It can be shown that the unknowns A
i, i
=
0,1
,
2,3, are given as follows: (8)
A
0
=
α
i
f
i
,
A
1
=
2
α
i
β
i
+
α
i
+
γ
i
f
i
+
α
i
h
i
d
i
,
A
2
=
2
α
i
β
i
+
β
i
+
γ
i
f
i
+
1
-
β
i
h
i
d
i
+
1
,
A
3
=
β
i
f
i
+
1
,where s
i
(
1
)
(
x
) denotes derivative with respect to x and d
i denotes the derivative value which is given at the knot x
i, i
=
0,1
,
2
,
…
,
n. The parameters are chosen such that α
i
,
β
i
>
0, γ
i
≥
0. When α
i
=
β
i
=
1, γ
i
=
0, the rational cubic interpolant in (5) is reduced to cubic Hermite spline [19].Positivity Preserving. For the strictly positive set of data (
x
i
,
f
i
), i
=
0,1
,
…
,
n, x
0
<
x
1
<
⋯
<
x
n,(9)
f
i
>
0
,
i
=
0,1
,
…
,
n
.Karim and Kong [19] have developed the following results for strictly positive data sets given in (9).Theorem 1.
For strictly positive data defined in (9), the C
1 rational cubic interpolant defined over the interval [
x
0
,
x
n
] is positive if in each subinterval [
x
i
,
x
i
+
1
], i
=
0,1
,
…
,
n
-
1, the involving parameters α
i, β
i, and γ
i satisfy the following sufficient conditions:(10)
α
i
,
β
i
>
0
,
γ
i
=
λ
i
+
Max
⁡
0
,
-
α
i
h
i
d
i
+
2
β
i
+
1
f
i
f
i
,
h
h
h
h
h
h
h
l
h
β
i
h
i
d
i
+
1
-
2
α
i
+
1
f
i
+
1
f
i
+
1
,
λ
i
>
0
.


### 3.1. Numerical Demonstrations

In this section the positivity shape preserving usingC
1 rational cubic spline of Karim and Kong [19] is discussed. The positive data is taken from the study by Hussain et al. [28].Example 2.
See  Tables1 and 2 and Figure 1.Table 1
Positive data from the study by Hussain et al. [28].
i
0
1
2
3
x
i
0.0
1.0
1.70
1.80
f
i
0.25
1.0
11.10
25Table 2
Numerical results.
i
0
1
2
3
d
i
−7.296
8.796
123.429
154.570
Δ
i
0.75
14.429
139
α
i
0.5
0.5
0.5
β
i
0.5
0.5
0.5
γ
i
13.84
3.14
0.25Figure 1
Comparison of (a) cubic Hermite spline curve and (b) our rational cubic spline withα
i
=
β
i
=
0.5.
(a)
(b)


## 3.1. Numerical Demonstrations

In this section the positivity shape preserving usingC
1 rational cubic spline of Karim and Kong [19] is discussed. The positive data is taken from the study by Hussain et al. [28].Example 2.
See  Tables1 and 2 and Figure 1.Table 1
Positive data from the study by Hussain et al. [28].
i
0
1
2
3
x
i
0.0
1.0
1.70
1.80
f
i
0.25
1.0
11.10
25Table 2
Numerical results.
i
0
1
2
3
d
i
−7.296
8.796
123.429
154.570
Δ
i
0.75
14.429
139
α
i
0.5
0.5
0.5
β
i
0.5
0.5
0.5
γ
i
13.84
3.14
0.25Figure 1
Comparison of (a) cubic Hermite spline curve and (b) our rational cubic spline withα
i
=
β
i
=
0.5.
(a)
(b)


## 4. Rational Bicubic Spline Interpolant

The piecewise rational cubic spline in (5) is extended to bicubic partially blended rational function S
(
x
,
y
) over rectangular domain Ω
=
[
a
,
b
]
×
[
c
,
d
]. The partition of the intervals [
a
,
b
] and [
c
,
d
] is defined as π
:
a
=
x
0
<
x
1
<
⋯
<
x
n
=
b and π
¯
:
c
=
y
0
<
y
1
<
⋯
<
y
m
=
d. The partially blended rational bicubic function over each rectangular patch x
i
,
x
i
+
1
×
[
y
j
,
y
j
+
1
], i
=
0,1
,
…
,
n
-
1 and j
=
0,1
,
…
,
m
-
1, is defined as follows: (11)
S
x
,
y
=
-
A
F
B
T
,where (12)
F
=
0
S
x
,
y
j
S
x
,
y
j
+
1
S
x
i
,
y
S
x
i
,
y
j
S
x
i
,
y
j
+
1
S
x
i
+
1
,
y
S
x
i
+
1
,
y
j
S
x
i
,
y
j
+
1
,
A
=
-
1
a
0
θ
a
1
θ
,
B
=
-
1
b
0
ϕ
b
1
ϕ
,where(13)
a
0
θ
=
1
-
θ
2
1
+
2
θ
,
a
1
θ
=
θ
2
3
-
2
θ
,
b
0
ϕ
=
1
-
ϕ
2
1
+
2
ϕ
,
b
1
ϕ
=
ϕ
2
3
-
2
ϕ
,
θ
=
x
-
x
i
h
i
,
ϕ
=
y
-
y
j
h
^
j
,
h
i
=
x
i
+
1
-
x
i
,
h
^
j
=
y
j
+
1
-
y
j
.
S
(
x
,
y
j
),  S
(
x
,
y
j
+
1
), S
(
x
i
,
y
), and S
(
x
i
+
1
,
y
) are rational cubic functions (5) defined on the boundary of rectangular patch x
i
,
x
i
+
1
×
[
y
j
,
y
j
+
1
], i
=
0,1
,
…
,
n
-
1, j
=
0,1
,
…
,
m
-
1, and are given as follows: (14)
S
x
,
y
j
=
∑
i
=
0
3
1
-
θ
3
-
i
θ
i
A
i
q
1
θ
,With A
0
=
α
i
,
j
F
i
,
j, A
1
=
(
2
α
i
,
j
β
i
,
j
+
α
i
,
j
+
γ
i
,
j
)
F
i
,
j
+
α
i
,
j
h
i
F
i
,
j
x, A
2
=
(
2
α
i
,
j
β
i
,
j
+
β
i
,
j
+
γ
i
,
j
)
F
i
+
1
,
j
-
β
i
,
j
h
i
F
i
+
1
,
j
x, A
3
=
β
i
,
j
F
i
+
1
,
j, and q
1
(
θ
)
=
(
1
-
θ
)
2
α
i
,
j
+
(
2
α
i
,
j
β
i
,
j
+
γ
i
,
j
)
θ
(
1
-
θ
)
+
θ
2
β
i
,
j;(15)
S
x
,
y
j
+
1
=
∑
i
=
0
3
1
-
θ
3
-
i
θ
i
B
i
q
2
θ
,with B
0
=
α
i
,
j
+
1
F
i
,
j
+
1, B
1
=
(
2
α
i
,
j
+
1
β
i
,
j
+
1
+
α
i
,
j
+
1
+
γ
i
,
j
+
1
)
F
i
,
j
+
1
+
α
i
,
j
+
1
h
i
F
i
,
j
+
1
x, B
2
=
(
2
α
i
,
j
+
1
β
i
,
j
+
1
+
β
i
,
j
+
1
+
γ
i
,
j
+
1
)
F
i
+
1
,
j
+
1
-
β
i
,
j
+
1
h
i
F
i
+
1
,
j
+
1
x, B
3
=
β
i
,
j
F
i
+
1
,
j
+
1, and q
2
(
θ
)
=
(
1
-
θ
)
2
α
i
,
j
+
1
+
(
2
α
i
,
j
+
1
β
i
,
j
+
1
+
γ
i
,
j
+
1
)
θ
(
1
-
θ
)
+
θ
2
β
i
,
j
+
1;(16)
S
x
i
,
y
=
∑
i
=
0
3
1
-
ϕ
3
-
i
ϕ
i
C
i
q
3
ϕ
,with C
0
=
α
^
i
,
j
F
i
,
j, C
1
=
(
2
α
^
i
,
j
β
^
i
,
j
+
α
^
i
,
j
+
γ
^
i
,
j
)
F
i
,
j
+
α
^
i
,
j
h
^
j
F
i
,
j
y, C
2
=
(
2
α
^
i
,
j
β
^
i
,
j
+
β
^
i
,
j
+
γ
^
i
,
j
)
F
i
,
j
+
1
-
β
^
i
,
j
h
^
j
F
i
,
j
+
1
y, C
3
=
β
^
i
,
j
F
i
,
j
+
1, and q
3
(
ϕ
)
=
1
-
ϕ
2
α
^
i
,
j
+
(
2
α
^
i
,
j
β
^
i
,
j
+
γ
^
i
,
j
)
ϕ
(
1
-
ϕ
)
+
ϕ
2
β
^
i
,
j;(17)
S
x
i
+
1
,
y
=
∑
i
=
0
3
1
-
ϕ
3
-
i
ϕ
i
D
i
q
4
ϕ
,with D
0
=
α
^
i
+
1
,
j
F
i
+
1
,
j, D
1
=
(
2
α
^
i
+
1
,
j
β
^
i
+
1
,
j
+
α
^
i
+
1
,
j
+
γ
^
i
+
1
,
j
)
F
i
+
1
,
j
+
α
^
i
+
1
,
j
h
^
j
F
i
+
1
,
j
y, D
2
=
(
2
α
^
i
+
1
,
j
β
^
i
+
1
,
j
+
β
^
i
+
1
,
j
+
γ
^
i
+
1
,
j
)
F
i
+
1
,
j
+
1
-
β
^
i
+
1
,
j
h
^
j
F
i
+
1
,
j
+
1
y, D
3
=
β
^
i
,
j
F
i
+
1
,
j
+
1, and q
4
(
ϕ
)
=
1
-
ϕ
2
α
^
i
+
1
,
j
+
(
2
α
^
i
+
1
,
j
β
^
i
+
1
,
j
+
γ
^
i
+
1
,
j
)
ϕ
(
1
-
ϕ
)
+
ϕ
2
β
^
i
+
1
,
j.Theorem 3.
The partially blended rational bicubic splineS
(
x
,
y
) defined by (11) is C
1-continuous degree-seven piecewise rational surface.Proof.
The results trivially follow from Proposition 3 in the study by Casciola and Romani [29].


## 5. Positivity Preserving of 3D Positive Data

Let(
x
i
,
y
i
,
F
i
,
j
) be positive data defined over rectangular grid x
i
,
x
i
+
1
×
[
y
j
,
y
j
+
1
], i
=
0,1
,
…
,
n
-
1, j
=
0,1
,
…
,
m
-
1, such that (18)
F
i
,
j
>
0
,
∀
i
,
j
.Casciola and Romani [29] have developed the following results: bicubic partially blended surface patch inherits all the properties of network of boundary curves.Thus the partially bicubic surface patch in (21) is positive if boundary curves defined in (22)–(25) are positive; that is, S
(
x
,
y
j
)
>
0, S
(
x
,
y
j
+
1
)
>
0, S
(
x
i
,
y
)
>
0, and S
(
x
i
+
1
,
y
)
>
0, respectively. By simple algebraic manipulation we have the following theorem.Theorem 4.
The piecewise rational partially bicubic functionS
(
x
,
y
) in (21) defined over the rectangular mesh x
i
,
x
i
+
1
×
[
y
j
,
y
j
+
1
], i
=
0,1
,
…
,
n
-
1, j
=
0,1
,
…
,
m
-
1, preserves the positive surface data if the parameters satisfy the following sufficient conditions: (19)
α
i
,
j
>
0
,
α
i
,
j
+
1
>
0
,
β
i
,
j
>
0
,
β
i
,
j
+
1
>
0
,
α
^
i
,
j
>
0
,
α
^
i
+
1
,
j
>
0
,
β
^
i
,
j
>
0
,
β
^
i
+
1
,
j
>
0
,
γ
i
,
j
>
Max
⁡
0
,
-
α
i
,
j
h
i
F
i
,
j
x
+
2
β
i
,
j
+
1
F
i
,
j
F
i
,
j
,
h
h
h
h
h
h
β
i
,
j
h
i
F
i
+
1
,
j
x
-
2
α
i
,
j
+
1
F
i
+
1
,
j
F
i
+
1
,
j
,
γ
i
,
j
+
1
>
Max
⁡
0
,
-
α
i
,
j
+
1
h
i
F
i
,
j
+
1
x
+
2
β
i
,
j
+
1
+
1
F
i
,
j
+
1
F
i
,
j
+
1
,
h
h
h
h
h
h
β
i
,
j
+
1
h
i
F
i
+
1
,
j
+
1
x
-
2
α
i
,
j
+
1
+
1
F
i
+
1
,
j
+
1
F
i
+
1
,
j
+
1
,
γ
^
i
,
j
>
Max
⁡
0
,
-
α
^
i
,
j
h
^
j
F
i
,
j
y
+
2
β
^
i
,
j
+
1
F
i
,
j
F
i
,
j
,
h
h
h
h
h
h
β
^
i
,
j
h
^
j
F
i
,
j
+
1
y
-
2
α
^
i
,
j
+
1
F
i
,
j
+
1
F
i
,
j
+
1
,
γ
^
i
+
1
,
j
>
Max
⁡
0
,
-
α
^
i
+
1
,
j
h
^
j
F
i
+
1
,
j
y
+
2
β
^
i
+
1
,
j
+
1
F
i
+
1
,
j
F
i
+
1
,
j
,
h
h
h
h
h
h
β
^
i
+
1
,
j
h
^
j
F
i
+
1
,
j
+
1
y
-
2
α
^
i
+
1
,
j
+
1
F
i
+
1
,
j
+
1
F
i
+
1
,
j
+
1
.Proof.
From (14)–(17), S
(
x
,
y
j
)
>
0, if A
i
>
0, i
=
0,1
,
2,3. Now for α
i
,
j
>
0, β
i
,
j
>
0, and γ
i
,
j
>
0, S
(
x
,
y
j
)
>
0, if (20)
γ
i
,
j
>
-
α
i
,
j
h
i
F
i
,
j
x
+
2
β
i
,
j
+
1
F
i
,
j
F
i
,
j
,
γ
i
,
j
>
β
i
,
j
h
i
F
i
+
1
,
j
x
-
2
α
i
,
j
+
1
F
i
+
1
,
j
F
i
+
1
,
j
.Combining the conditions in (20) results in the following sufficient conditions for the positivity of (
x
,
y
j
):(21)
γ
i
,
j
>
Max
⁡
0
,
-
α
i
,
j
h
i
F
i
,
j
x
+
2
β
i
,
j
+
1
F
i
,
j
F
i
,
j
,
h
h
h
h
h
h
β
i
,
j
h
i
F
i
+
1
,
j
x
-
2
α
i
,
j
+
1
F
i
+
1
,
j
F
i
+
1
,
j
.Likewise the remaining three boundary curves S
(
x
,
y
j
+
1
), S
(
x
i
,
y
), and S
(
x
i
+
1
,
y
) are positive if the following sufficient conditions are satisfied:(22)
γ
i
,
j
+
1
>
-
α
i
,
j
+
1
h
i
F
i
,
j
+
1
x
+
2
β
i
,
j
+
1
+
1
F
i
,
j
+
1
F
i
,
j
+
1
,
(23)
γ
i
,
j
+
1
>
β
i
,
j
+
1
h
i
F
i
+
1
,
j
+
1
x
-
2
α
i
,
j
+
1
+
1
F
i
+
1
,
j
+
1
F
i
+
1
,
j
+
1
,
(24)
γ
^
i
,
j
>
-
α
^
i
,
j
h
^
j
F
i
,
j
y
+
2
β
^
i
,
j
+
1
F
i
,
j
F
i
,
j
,
(25)
γ
^
i
,
j
>
β
^
i
,
j
h
^
j
F
i
,
j
+
1
y
-
2
α
^
i
,
j
+
1
F
i
,
j
+
1
F
i
,
j
+
1
,
(26)
γ
^
i
+
1
,
j
>
-
α
^
i
+
1
,
j
h
^
j
F
i
+
1
,
j
y
+
2
β
^
i
+
1
,
j
+
1
F
i
+
1
,
j
F
i
+
1
,
j
,
(27)
γ
^
i
+
1
,
j
>
β
^
i
+
1
,
j
h
^
j
F
i
+
1
,
j
+
1
y
-
2
α
^
i
+
1
,
j
+
1
F
i
+
1
,
j
+
1
F
i
+
1
,
j
+
1
.Combining the conditions in (22) and (23), (24) and (25), and (26) and (27) leads us to the following:(28)
γ
i
,
j
+
1
>
Max
⁡
0
,
-
α
i
,
j
+
1
h
i
F
i
,
j
+
1
x
+
2
β
i
,
j
+
1
+
1
F
i
,
j
+
1
F
i
,
j
+
1
,
h
h
h
h
h
h
β
i
,
j
+
1
h
i
F
i
+
1
,
j
+
1
x
-
2
α
i
,
j
+
1
+
1
F
i
+
1
,
j
+
1
F
i
+
1
,
j
+
1
,
γ
^
i
,
j
>
Max
⁡
0
,
-
α
^
i
,
j
h
^
j
F
i
,
j
y
+
2
β
^
i
,
j
+
1
F
i
,
j
F
i
,
j
,
h
h
h
h
h
h
β
^
i
,
j
h
^
j
F
i
,
j
+
1
y
-
2
α
^
i
,
j
+
1
F
i
,
j
+
1
F
i
,
j
+
1
,
γ
^
i
+
1
,
j
>
Max
⁡
0
,
-
α
^
i
+
1
,
j
h
^
j
F
i
+
1
,
j
y
+
2
β
^
i
+
1
,
j
+
1
F
i
+
1
,
j
F
i
+
1
,
j
,
h
h
h
h
h
h
β
^
i
+
1
,
j
h
^
j
F
i
+
1
,
j
+
1
y
-
2
α
^
i
+
1
,
j
+
1
F
i
+
1
,
j
+
1
F
i
+
1
,
j
+
1
.Then from (21) and (28), sufficient conditions for the positivity of network of all four curves S
(
x
,
y
j
), S
(
x
,
y
j
+
1
), S
(
x
i
,
y
), and S
(
x
i
+
1
,
y
) are obtained. Thus the proof is complete.The sufficient conditions in (19) can be rewritten as(29)
α
i
,
j
>
0
,
α
i
,
j
+
1
>
0
,
β
i
,
j
>
0
,
β
i
,
j
+
1
>
0
,
α
^
i
,
j
>
0
,
α
^
i
+
1
,
j
>
0
,
β
^
i
,
j
>
0
,
β
^
i
+
1
,
j
>
0
,
γ
i
,
j
=
k
i
,
j
+
Max
⁡
0
,
-
α
i
,
j
h
i
F
i
,
j
x
+
2
β
i
,
j
+
1
F
i
,
j
F
i
,
j
,
h
k
i
,
j
+
h
h
h
h
h
β
i
,
j
h
i
F
i
+
1
,
j
x
-
2
α
i
,
j
+
1
F
i
+
1
,
j
F
i
+
1
,
j
,
γ
i
,
j
+
1
=
l
i
,
j
+
Max
⁡
0
,
-
α
i
,
j
+
1
h
i
F
i
,
j
+
1
x
+
2
β
i
,
j
+
1
+
1
F
i
,
j
+
1
F
i
,
j
+
1
,
h
h
h
h
h
h
h
β
i
,
j
+
1
h
i
F
i
+
1
,
j
+
1
x

-
2
α
i
,
j
+
1
+
1

F
i
+
1
,
j
+
1
F
i
+
1
,
j
+
1

,
γ
^
i
,
j
=
m
i
,
j
+
Max
⁡
0
,
-
α
^
i
,
j
h
^
j
F
i
,
j
y
+
2
β
^
i
,
j
+
1
F
i
,
j
F
i
,
j
,
h
k
i
,
j
+
h
h
l
h
h
h
β
^
i
,
j
h
^
j
F
i
,
j
+
1
y
-
2
α
^
i
,
j
+
1
F
i
,
j
+
1
F
i
,
j
+
1
,
γ
^
i
+
1
,
j
=
n
i
,
j
+
Max
⁡
0
,
-
α
^
i
+
1
,
j
h
^
j
F
i
+
1
,
j
y
+
2
β
^
i
+
1
,
j
+
1
F
i
+
1
,
j
F
i
+
1
,
j
,
h
h
h
h
l
h
h
h
β
^
i
+
1
,
j
h
^
j
F
i
+
1
,
j
+
1
y

-
2
α
^
i
+
1
,
j
+
1

F
i
+
1
,
j
+
1
F
i
+
1
,
j
+
1

,where k
i
,
j
>
0, l
i
,
j
>
0, m
i
,
j
>
0, and n
i
,
j
>
0.Theorem 5.
The partially blended rational bicubic splineS
(
x
,
y
) constructed by using shape parameters calculated from sufficient condition for positivity in (29) is C
1-continuous degree-seven piecewise rational surface with positivity preserving properties to the positive data sets.Proof.
The proof also trivially follows from the work of Casciola and Romani [29].


### 5.1. Algorithm for Computer Implementations

The algorithm is for the purpose of computer implementation of positivity preserving by using our partially rational bicubic spline presented in Theorem4.Input. The input is the positive data points (
x
i
,
y
j
,
F
i
,
j
), i
=
0,1
,
…
,
n and j
=
0,1
,
…
,
m.Output. The output is the positive surfaces.Step 1.
Calculate the first derivative valuesF
i
,
j
x, F
i
,
j
y, F
i
,
j
x
y by using AMM.Step 2.
Choose suitable values for free parametersα
i
,
j, β
i
,
j, α
i
,
j
+
1, β
i
,
j
+
1, α
^
i
,
j, α
^
i
+
1
,
j, β
^
i
,
j, β
^
i
+
1
,
j and k
i
,
j
,
l
i
,
j
,
m
i
,
j
,
n
i
,
j, and compute the values of parameters γ
i
,
j, γ
i
,
j
+
1, γ
^
i
,
j, γ
^
i
+
1
,
j by using Theorem 4.Step 3.
Construct the positive surfaces by substituting all the required parameters into the partially blended rational bicubic spline function (11).Repeat Steps2 and 3 to generate different positive surfaces through different positive data sets.


### 5.2. Numerical Demonstrations

In this section, numerical examples for positivity preserving by using partially blended rational bicubic spline interpolation of positive surface data will be discussed in detail. Four positive data sets are taken from the studies by Abbas et al. [1, 2] and Hussain and Sarfraz [12], respectively.Example 6.
Positive data from the following function is truncated to five decimal places [1]: (30)
F
1
x
,
y
=
e
-
(
x
2
+
y
2
)
/
15
sin
x
+
cos
⁡
y
+
0.33
,
0
≤
x
,
y
≤
6
.Figure 2 shows the default bicubic Hermite spline for the positive data given in Table 3. Figures 3(a) and 3(b) show the x
z-view and y
z-view for Figure 2, respectively. Figure 4 shows the positivity preserving by using the proposed rational bicubic spline with α
i
,
j
=
β
i
,
j
=
0.5, α
^
i
,
j
=
β
^
i
,
j
=
0.5. Meanwhile, Figures 5(a) and 5(b) show the x
z-view and y
z-view for Figure 4, respectively.Table 3
Positive surface data from functionF
1
(
x
,
y
).
y
/
x
0
2
4
6
0
1.33000
0.011261
0.10505
0.41710
2
1.79240
0.619300
0.39739
0.45990
4
0.41370
0.020814
0.16294
0.33635
6
0.39537
0.281670
0.30087
0.33560Figure 2
Bicubic Hermite surface.Figure 3
Bicubic Hermite surface: (a)x
z-view and (b) y
z-view for Figure 2, respectively.
(a)
(b)Figure 4
Positivity preserving by using our partially rational bicubic spline.Figure 5
Positivity preserving surface: (a)x
z-view and (b) y
z-view for Figure 4, respectively.
(a)
(b)Example 7.
Positive data from the following function is truncated to five decimal places [1]: (31)
F
2
x
,
y
=
sin
y
e
-
x
+
1
,
-
3
≤
x
,
y
≤
3
,
x
,
y
≠
0
.Figure 6 shows the default bicubic Hermite spline for the positive data given in Table 4. Figures 7(a) and 7(b) show the x
z-view and y
z-view for Figure 6, respectively. Figure 8 shows the positivity preserving by using the proposed rational bicubic spline with α
i
,
j
=
β
i
,
j
=
0.5, α
^
i
,
j
=
β
^
i
,
j
=
0.5. Meanwhile, Figures 9(a) and 9(b) show the x
z-view and y
z-view for Figure 8, respectively.Table 4
Positive surface data from functionF
2
(
x
,
y
).
y
/
x
−3
−2
−1
1
2
3
−3
1.536600
0.37930
0.055529
1.9445
1.62070
0.46344
−2
1.175100
0.19847
0.106150
1.8939
1.80150
0.82489
−1
0.044919
1.74900
0.589220
1.4108
0.25095
1.95510
1
0.107150
0.32885
0.640360
1.3596
1.67110
1.89290
2
0.605060
0.73262
0.865080
1.1349
1.26740
1.39490
3
0.851190
0.90059
0.950230
1.0498
1.09940
1.14880Figure 6
Bicubic Hermite surface.Figure 7
Bicubic Hermite surface: (a)x
z-view and (b) y
z-view for Figure 6, respectively.
(a)
(b)Figure 8
Positivity preserving by using our partially rational bicubic spline.Figure 9
Positivity preserving surface: (a)x
z-view and (b) y
z-view for Figure 8, respectively.
(a)
(b)Example 8.
Positive data from the following function is truncated to four decimal places [2]: (32)
F
3
x
,
y
=
e
-
x
2
+
e
-
2
y
2
+
0.04
,
-
3
≤
x
,
y
≤
3
.Figure 10 shows the default bicubic Hermite spline for the positive data given in Table 5. Figures 11(a) and 11(b) show the x
z-view and y
z-view for Figure 10, respectively. Figure 12 shows the positivity preserving by using the proposed rational bicubic spline with α
i
,
j
=
β
i
,
j
=
1.5, α
^
i
,
j
=
β
^
i
,
j
=
1.5. Meanwhile, Figures 13(a) and 13(b) show the x
z-view and y
z-view for Figure 12, respectively.Table 5
Positive surface data from functionF
3
(
x
,
y
).
y
/
x
−3
−2
−1
0
1
2
3
−3
0.0401
0.0404
0.1755
1.0401
0.1755
0.0404
0.0401
−2
0.0583
0.0586
0.1936
1.0583
0.1936
0.0586
0.0583
−1
0.4078
0.4082
0.5432
1.4079
0.5432
0.4082
0.4078
0
1.0400
1.0403
1.1753
2.0400
1.1753
1.0403
1.0400
1
0.4078
0.4082
0.5432
1.4079
0.5432
0.4082
0.4078
2
0.0583
0.0586
0.1936
1.0583
0.1936
0.0586
0.0583
3
0.0401
0.0404
0.1755
1.0401
0.1755
0.0404
0.0401Figure 10
Bicubic Hermite surface.Figure 11
Bicubic Hermite surface: (a)x
z-view and (b) y
z-view for Figure 10, respectively.
(a)
(b)Figure 12
Positivity preserving by using our partially rational bicubic spline.Figure 13
Positivity preserving surface: (a)x
z-view and (b) y
z-view for Figure 12, respectively.
(a)
(b)Example 9.
Positive data from the following function is truncated to four decimal places [12]: (33)
F
4
x
,
y
=
4
x
2
+
y
2
2
-
1
,
-
3
≤
x
,
y
≤
3
,
x
,
y
≠
0
.Figure 14 shows the default bicubic Hermite spline for the positive data given in Table 6. Figures 15(a) and 15(b) show the x
z-view and y
z-view for Figure 14, respectively. Figure 16 shows the positivity preserving by using the proposed rational bicubic spline with α
i
,
j
=
β
i
,
j
=
1.5, α
^
i
,
j
=
β
^
i
,
j
=
1.5. Meanwhile, Figures 17(a) and 17(b) show the x
z-view and y
z-view for Figure 16, respectively.Table 6
Positive surface data from functionF
4
(
x
,
y
).
y
/
x
−3
−2
−1
1
2
3
−3
0.0124
0.0238
0.0404
0.0404
0.0238
0.0124
−2
0.0238
0.0635
0.1667
0.1667
0.0635
0.0238
−1
0.0404
0.1667
1.3333
1.3333
0.1667
0.0404
1
0.0404
0.1667
1.3333
1.3333
0.1667
0.0404
2
0.0238
0.0635
0.1667
0.1667
0.0635
0.0238
3
0.0124
0.0238
0.0404
0.0404
0.0238
0.0124Figure 14
Bicubic Hermite surface.Figure 15
Bicubic Hermite surface: (a)x
z-view and (b) y
z-view for Figure 14, respectively.
(a)
(b)Figure 16
Positivity preserving by using our partially rational bicubic spline.Figure 17
Positivity preserving surface: (a)x
z-view and (b) y
z-view for Figure 16, respectively.
(a)
(b)Meanwhile, Figures18(a) and 18(b) show the positivity preserving by using the study by Abbas et al. [1] for positive data in Tables 3 and 4, respectively. Figure 19 shows the positivity preserving by using the study by Abbas et al. [2] for positive data in Table 5. Finally Figures 20 and 21 show the positive surfaces by using Hussain and Sarfraz [12] rational bicubic spline and Liu et al. [20] rational biquartic spline surfaces.Figure 18
Positivity preserving by using the study of Abbas et al. [1] for (a) positive data in Table 3 and (b) positive data in Table 4.
(a)
(b)Figure 19
Positivity preserving by using the study of Abbas et al. [2] for positive data in Table 5.Figure 20
Positivity preserving by using the study of Hussain and Sarfraz [12] for positive data in Table 6.Figure 21
Positivity preserving by using the study of Liu et al. [20] for positive data in Table 6.We compare the performance between our proposed partially blended rational bicubic spline interpolation and that of the work of Abbas et al. [1] by calculating Root Mean Square Error (RMSE). Table 7 summarized the results of RMSE estimation. It can be seen clearly that our rational bicubic spline interpolation is on a par with the work of Abbas et al. [1]. The value for RMSE for all four data sets is almost the same. Thus we conclude that the proposed partially blended rational bicubic spline interpolation produces the positive surfaces through positive data sets with smaller RMSE. Furthermore we believe that the computation is less by using our rational bicubic spline interpolation compared with the other rational bicubic spline interpolation of Hussain and Sarfraz [12] and Sarfraz et al. [22].Table 7
RMSE estimation.
Method
Data sets
Example6
Example7
Example8
Example9
Abbas et al. [1]
0.3900
0.0252
0.0011
0.0103
Ours
0.3947
0.0233
0.0012
0.0106


## 5.1. Algorithm for Computer Implementations

The algorithm is for the purpose of computer implementation of positivity preserving by using our partially rational bicubic spline presented in Theorem4.Input. The input is the positive data points (
x
i
,
y
j
,
F
i
,
j
), i
=
0,1
,
…
,
n and j
=
0,1
,
…
,
m.Output. The output is the positive surfaces.Step 1.
Calculate the first derivative valuesF
i
,
j
x, F
i
,
j
y, F
i
,
j
x
y by using AMM.Step 2.
Choose suitable values for free parametersα
i
,
j, β
i
,
j, α
i
,
j
+
1, β
i
,
j
+
1, α
^
i
,
j, α
^
i
+
1
,
j, β
^
i
,
j, β
^
i
+
1
,
j and k
i
,
j
,
l
i
,
j
,
m
i
,
j
,
n
i
,
j, and compute the values of parameters γ
i
,
j, γ
i
,
j
+
1, γ
^
i
,
j, γ
^
i
+
1
,
j by using Theorem 4.Step 3.
Construct the positive surfaces by substituting all the required parameters into the partially blended rational bicubic spline function (11).Repeat Steps2 and 3 to generate different positive surfaces through different positive data sets.


## 5.2. Numerical Demonstrations

In this section, numerical examples for positivity preserving by using partially blended rational bicubic spline interpolation of positive surface data will be discussed in detail. Four positive data sets are taken from the studies by Abbas et al. [1, 2] and Hussain and Sarfraz [12], respectively.Example 6.
Positive data from the following function is truncated to five decimal places [1]: (30)
F
1
x
,
y
=
e
-
(
x
2
+
y
2
)
/
15
sin
x
+
cos
⁡
y
+
0.33
,
0
≤
x
,
y
≤
6
.Figure 2 shows the default bicubic Hermite spline for the positive data given in Table 3. Figures 3(a) and 3(b) show the x
z-view and y
z-view for Figure 2, respectively. Figure 4 shows the positivity preserving by using the proposed rational bicubic spline with α
i
,
j
=
β
i
,
j
=
0.5, α
^
i
,
j
=
β
^
i
,
j
=
0.5. Meanwhile, Figures 5(a) and 5(b) show the x
z-view and y
z-view for Figure 4, respectively.Table 3
Positive surface data from functionF
1
(
x
,
y
).
y
/
x
0
2
4
6
0
1.33000
0.011261
0.10505
0.41710
2
1.79240
0.619300
0.39739
0.45990
4
0.41370
0.020814
0.16294
0.33635
6
0.39537
0.281670
0.30087
0.33560Figure 2
Bicubic Hermite surface.Figure 3
Bicubic Hermite surface: (a)x
z-view and (b) y
z-view for Figure 2, respectively.
(a)
(b)Figure 4
Positivity preserving by using our partially rational bicubic spline.Figure 5
Positivity preserving surface: (a)x
z-view and (b) y
z-view for Figure 4, respectively.
(a)
(b)Example 7.
Positive data from the following function is truncated to five decimal places [1]: (31)
F
2
x
,
y
=
sin
y
e
-
x
+
1
,
-
3
≤
x
,
y
≤
3
,
x
,
y
≠
0
.Figure 6 shows the default bicubic Hermite spline for the positive data given in Table 4. Figures 7(a) and 7(b) show the x
z-view and y
z-view for Figure 6, respectively. Figure 8 shows the positivity preserving by using the proposed rational bicubic spline with α
i
,
j
=
β
i
,
j
=
0.5, α
^
i
,
j
=
β
^
i
,
j
=
0.5. Meanwhile, Figures 9(a) and 9(b) show the x
z-view and y
z-view for Figure 8, respectively.Table 4
Positive surface data from functionF
2
(
x
,
y
).
y
/
x
−3
−2
−1
1
2
3
−3
1.536600
0.37930
0.055529
1.9445
1.62070
0.46344
−2
1.175100
0.19847
0.106150
1.8939
1.80150
0.82489
−1
0.044919
1.74900
0.589220
1.4108
0.25095
1.95510
1
0.107150
0.32885
0.640360
1.3596
1.67110
1.89290
2
0.605060
0.73262
0.865080
1.1349
1.26740
1.39490
3
0.851190
0.90059
0.950230
1.0498
1.09940
1.14880Figure 6
Bicubic Hermite surface.Figure 7
Bicubic Hermite surface: (a)x
z-view and (b) y
z-view for Figure 6, respectively.
(a)
(b)Figure 8
Positivity preserving by using our partially rational bicubic spline.Figure 9
Positivity preserving surface: (a)x
z-view and (b) y
z-view for Figure 8, respectively.
(a)
(b)Example 8.
Positive data from the following function is truncated to four decimal places [2]: (32)
F
3
x
,
y
=
e
-
x
2
+
e
-
2
y
2
+
0.04
,
-
3
≤
x
,
y
≤
3
.Figure 10 shows the default bicubic Hermite spline for the positive data given in Table 5. Figures 11(a) and 11(b) show the x
z-view and y
z-view for Figure 10, respectively. Figure 12 shows the positivity preserving by using the proposed rational bicubic spline with α
i
,
j
=
β
i
,
j
=
1.5, α
^
i
,
j
=
β
^
i
,
j
=
1.5. Meanwhile, Figures 13(a) and 13(b) show the x
z-view and y
z-view for Figure 12, respectively.Table 5
Positive surface data from functionF
3
(
x
,
y
).
y
/
x
−3
−2
−1
0
1
2
3
−3
0.0401
0.0404
0.1755
1.0401
0.1755
0.0404
0.0401
−2
0.0583
0.0586
0.1936
1.0583
0.1936
0.0586
0.0583
−1
0.4078
0.4082
0.5432
1.4079
0.5432
0.4082
0.4078
0
1.0400
1.0403
1.1753
2.0400
1.1753
1.0403
1.0400
1
0.4078
0.4082
0.5432
1.4079
0.5432
0.4082
0.4078
2
0.0583
0.0586
0.1936
1.0583
0.1936
0.0586
0.0583
3
0.0401
0.0404
0.1755
1.0401
0.1755
0.0404
0.0401Figure 10
Bicubic Hermite surface.Figure 11
Bicubic Hermite surface: (a)x
z-view and (b) y
z-view for Figure 10, respectively.
(a)
(b)Figure 12
Positivity preserving by using our partially rational bicubic spline.Figure 13
Positivity preserving surface: (a)x
z-view and (b) y
z-view for Figure 12, respectively.
(a)
(b)Example 9.
Positive data from the following function is truncated to four decimal places [12]: (33)
F
4
x
,
y
=
4
x
2
+
y
2
2
-
1
,
-
3
≤
x
,
y
≤
3
,
x
,
y
≠
0
.Figure 14 shows the default bicubic Hermite spline for the positive data given in Table 6. Figures 15(a) and 15(b) show the x
z-view and y
z-view for Figure 14, respectively. Figure 16 shows the positivity preserving by using the proposed rational bicubic spline with α
i
,
j
=
β
i
,
j
=
1.5, α
^
i
,
j
=
β
^
i
,
j
=
1.5. Meanwhile, Figures 17(a) and 17(b) show the x
z-view and y
z-view for Figure 16, respectively.Table 6
Positive surface data from functionF
4
(
x
,
y
).
y
/
x
−3
−2
−1
1
2
3
−3
0.0124
0.0238
0.0404
0.0404
0.0238
0.0124
−2
0.0238
0.0635
0.1667
0.1667
0.0635
0.0238
−1
0.0404
0.1667
1.3333
1.3333
0.1667
0.0404
1
0.0404
0.1667
1.3333
1.3333
0.1667
0.0404
2
0.0238
0.0635
0.1667
0.1667
0.0635
0.0238
3
0.0124
0.0238
0.0404
0.0404
0.0238
0.0124Figure 14
Bicubic Hermite surface.Figure 15
Bicubic Hermite surface: (a)x
z-view and (b) y
z-view for Figure 14, respectively.
(a)
(b)Figure 16
Positivity preserving by using our partially rational bicubic spline.Figure 17
Positivity preserving surface: (a)x
z-view and (b) y
z-view for Figure 16, respectively.
(a)
(b)Meanwhile, Figures18(a) and 18(b) show the positivity preserving by using the study by Abbas et al. [1] for positive data in Tables 3 and 4, respectively. Figure 19 shows the positivity preserving by using the study by Abbas et al. [2] for positive data in Table 5. Finally Figures 20 and 21 show the positive surfaces by using Hussain and Sarfraz [12] rational bicubic spline and Liu et al. [20] rational biquartic spline surfaces.Figure 18
Positivity preserving by using the study of Abbas et al. [1] for (a) positive data in Table 3 and (b) positive data in Table 4.
(a)
(b)Figure 19
Positivity preserving by using the study of Abbas et al. [2] for positive data in Table 5.Figure 20
Positivity preserving by using the study of Hussain and Sarfraz [12] for positive data in Table 6.Figure 21
Positivity preserving by using the study of Liu et al. [20] for positive data in Table 6.We compare the performance between our proposed partially blended rational bicubic spline interpolation and that of the work of Abbas et al. [1] by calculating Root Mean Square Error (RMSE). Table 7 summarized the results of RMSE estimation. It can be seen clearly that our rational bicubic spline interpolation is on a par with the work of Abbas et al. [1]. The value for RMSE for all four data sets is almost the same. Thus we conclude that the proposed partially blended rational bicubic spline interpolation produces the positive surfaces through positive data sets with smaller RMSE. Furthermore we believe that the computation is less by using our rational bicubic spline interpolation compared with the other rational bicubic spline interpolation of Hussain and Sarfraz [12] and Sarfraz et al. [22].Table 7
RMSE estimation.
Method
Data sets
Example6
Example7
Example8
Example9
Abbas et al. [1]
0.3900
0.0252
0.0011
0.0103
Ours
0.3947
0.0233
0.0012
0.0106


## 6. Discussions

The proposed partially blended rational bicubic spline with 12 parameters has been used for positivity preserving of 3D positive data. The scheme has been tested on 4 different data sets. The bicubic Hermite spline does not have the capability to preserve the positivity of the data without changing the original data points. Figures2, 6, 10, and 14 show that there exist some negative values when the data sets are interpolated by using bicubic Hermite spline. From Figures 4, 8, 12, and 16 it can be seen clearly that the proposed scheme gives comparable results with existing schemes such as those of Abbas et al. [1, 2], Hussain and Sarfraz [12], and Liu et al. [20]. Similar to the work of Abbas et al. [1], our scheme also is partially blended bicubic rational spline. Meanwhile, Abbas et al. [2] and Hussain and Sarfraz [12] have used rational bicubic spline interpolation. Table 7 summarized the RMSE values comparison between the proposed partially blended rational bicubic spline interpolation and that of the work of Abbas et al. [1].


## 7. Summary and Conclusions

In this paper the rational cubic spline of Karim and Kong [18, 19] has been extended to the bivariate cases. The partially blended bicubic rational spline has 12 parameters in the descriptions and 8 of them are free parameters. These free parameters can be used to change the shape of the final surfaces of the positive data. Graphical comparison between the proposed schemes and existing schemes also has been done. From all numerical results, we concluded that the proposed partially blended bicubic rational spline works well for all tested data sets. Future works will be developing the rational bicubic spline for scattered data interpolation without the need to modify the first derivative. This will improve the work of Peng et al. [26] that requires the modification of the first derivative in which positivity preserving is not met. Work on monotonicity and convexity preservation by using the proposed partially blended rational bicubic spline interpolation is underway.



---

*Source: 572768-2015-03-22.xml*