1. Introduction
Subdivision schemes have attracted a lot of attention for applications in CAGD since Dyn et al. introduced, in [
1], a four point interpolatory subdivision scheme in 1987. These schemes can be considered an efficient way to generate curves and surfaces from an initial set of control points. In past years, a great number of works have been published on linear subdivision schemes. See, for example, Refs. [
2,
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13] and the references therein for an incomplete list of relevant publications. Nonlinear subdivision schemes have also attracted great attention. See, for example, Refs. [
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27] and the references therein for an incomplete list of publications about the subject.
Subdivision schemes can be classified into stationary or nonstationary in determining their dependence on the level
k of application of the scheme. If the scheme is stationary, the masks have an invariant expression independently of the level of subdivision
k. Otherwise, the scheme is said to be nonstationary. It is possible to find many stationary subdivision schemes in the state of the art literature. Among them, some binary schemes have been also proposed in the past years [
1,
28,
29,
30,
31,
32]. All of them suffer from diffusion close to jump discontinuities and some of them present Gibbs phenomenon [
33].
As mentioned in the abstract, if the control points are obtained from the sampling of a continuous function with high gradients, numerical discontinuities can appear due to low sampling frequencies that do not manage to capture the regularity of the function at those high gradients. If this is the case, the diffusion introduced in the limit function is not a problem and can even be considered a good quality of the subdivision scheme. If the control points are obtained from the sampling of a piecewise continuous function, the discontinuities should also appear in the limit function.
Thus, the main objective of this paper is to introduce a new nonlinear, noninterpolatory, stationary and binary family of subdivision schemes converging to a piecewise continuous function and that do not oscillate close to jump discontinuities. In fact, close to the discontinuities, the new family of schemes do not introduce diffusion. The properties of the family of subdivision schemes will be analyzed. In particular, convergence and stability of the family of subdivision schemes are also established. As far as we know, this is the first time that a binary, stable and convergent subdivision scheme without Gibbs oscillations close to jump discontinuities and reproducing piecewise constant functions has appeared in the literature.
The paper is organized as follows. In
Section 2, we introduce a new nonlinear, noninterpolatory family of subdivision schemes based on the one studied in [
31]. Writing the scheme as a perturbation of a linear scheme and establishing a contractivity property of this perturbation, we deduce the convergence (studied in
Section 3) and stability (studied in
Section 4) of the subdivision scheme, that, due to the nonlinear nature of the scheme, is not a consequence of the convergence.
Section 5 is devoted to numerical examples.
Section 5.1 analyzes the numerical regularity of some members of the family of schemes presented.
Section 5.2 presents some experiments of subdivision for piecewise smooth functions.
Section 5.3 is dedicated to the analysis of the order of accuracy of the schemes close to the discontinuities. Finally,
Section 6 presents our conclusions.
2. A Nonlinear Family of Binary Subdivision Schemes
This section is devoted to the introduction of a nonlinear binary subdivision scheme as a perturbation of the linear scheme introduced in [
31]. In this work, S. Shahid and A. Nadeem prove that the scheme that they proposed converges to a limit function with
regularity. They analyze the following binary subdivision scheme:
where
is a set of initial control points with
,
and
and
. The theory presented in [
19] assures the absence of Gibbs oscillations close to jump discontinuities.
Thus, the only drawback that the previous linear noninterpolatory scheme can present is that diffusion appears in the limit function when the scheme starts from data that might arise from the sampling of a discontinuous function. The mentioned linear subdivision scheme is based on a quadratic B-spline parametrized in the interval
. Due to the convex hull property of B-splines, this scheme cannot present Gibbs phenomenon [
34]. The authors prove, in [
31], that, for
, the scheme presents a limit function with
regularity. This result can also be expected because the scheme is based on a binary B-spline.
In order to obtain a nonlinear family of schemes from (
1), let us define
and consider the function
. The nonlinear scheme we analyze in this paper follows [
18,
35], and it is based on the substitution of the arithmetic mean with a
M function. We first propose a new formulation of the scheme (
1):
where
denotes the first order differences. The nonlinear scheme
that we propose is given by:
If
,
If
,
For
, we recover the original linear scheme in (
1). In the numerical experiments, we use the harmonic mean, properly introduced in
Section 5.
In [
36,
37,
38,
39], we use a similar technique to adapt subdivision schemes to the presence of discontinuities. Before analyzing in detail the properties of the new scheme
, in the following definition we summarize the most important properties that the function
M must verify (see [
23] for more details).
Definition 1. The function M is a mean on , and thus:
as uniformly with respect to x
on compact subintervals on
.
. for some constant
If , , and , then for all there exist non-negative numbers α and β (depending on these points) such that and with a constant not depending on
As a consequence of
, the authors presented, in [
23], the following proposition:
Proposition 1. Suppose that M is continuous on . If holds, then the extension of M byis continuous on . The proof can be found in Proposition 2 of [
23]. We use this extension to
in the rest of the paper.
For the linear family of schemes presented in this section, the elimination of the Gibbs phenomenon in the presence of discontinuities is not needed to be studied, as the family of schemes present positive masks that assure the absence of oscillations close to discontinuities [
19]. Indeed, in [
19], we prove a theorem that assures the absence of Gibbs oscillations for stationary subdivision schemes with positive masks. They follow the notation used by Dyn and Levin in [
3], where a univariate stationary subdivision scheme with a mask
with finite support is defined as beginning with an initial sequence of finite data
. New values are obtained through refinement at level
, denoted by
. The subdivision scheme is the maximal set obtained by applying the rule:
There are two rules to define the points on the level
:
We can represent these rules using an algebraic formalism in terms of z-transforms. The symbol of the mask
is defined as
. If
is denoted as the
k iterated symbol (see [
3]), then, it is known that
. Therefore, if
, then
with
being
,
. As in [
19], in this paper, we impose to the coefficients of the mask the following conditions:
The mask
is defined as:
with
and
V fixed integers with
and
.
The scheme
is convergent. Therefore, one necessary condition is
For any
,
and any
, with
and then
being
K a constant, which depends on
f, and it does not depend on
h.
With this notation, the theorem that assures the absence of Gibbs phenomenon for subdivision schemes with positive masks is the following.
Theorem 1. Given , let f be any function defined bywith and . Let be a univariate stationary subdivision scheme with:, being defined in Equation (8). Then, if , ; and if h is sufficiently small (meaning an h for which the differences between adjacent data at smooth zones are significantly smaller than the differences close to discontinuities), we have: - 1 .
If , then with .
- 2 .
If , there exists such that
The interested reader can refer to [
19] for the proof of this theorem.
3. Convergence of the New Family of Subdivision Schemes
Let us recall the following definitions that will be useful for our purposes,
Definition 2. A binary subdivision scheme S is said to be convergent ifwhere denotes the piecewise affine interpolation of at the grid points. Definition 3. A convergent subdivision scheme is stable if Let us first rewrite the scheme
(
3) and (
4) as a particular perturbation of the linear subdivision scheme
(projection), which is defined as,
We can see that the nonlinear scheme in (
3) and (
4) can be written in terms of the linear scheme
. For all
, we have
where
, and
F is the nonlinear operator
We also need to recall the following result established in [
14]. Let
be a subdivision scheme defined by
where
F is a nonlinear operator defined on
,
is a linear and continuous operator on
and
S is a linear subdivision scheme convergent in
.
Theorem 2. If S is a linear convergent subdivision scheme and if , F and δ verifythen the subdivision scheme defined by (16) is convergent. Using Theorem 2, we can prove the following theorem:
Theorem 3. The nonlinear subdivision scheme is convergent for all .
Proof. The perturbation
F defined in (
14), (
15) can be bounded using that
.
Then
that satisfies precisely hypothesis (
17).
Now, we consider hypothesis (
18), which is related to the contraction of first order differences
.
We have 16 different cases (see the
Appendix A at the end of the paper).
In all of the cases, we obtain:
Thus, for all
:
If , then .
Thus, hypothesis (
18) is satisfied, and
is
convergent. □
Remark 1. Other terms in the last theorem means that they do not depend on .
In this paper, we prove the stability only at the regularity parts using the norm; thus, we also have to prove the convergence at these parts using the norm.
Theorem 4. ( convergence)
If F and δ verify:Then, the subdivision scheme is uniformly convergent in . Proof. Let
be the limit function of the linear subdivision
S.
, we consider the continuous function
as:
where
is the sequence
. We have
Using lemma 2.3 in [
29], the definition of
and hypotheses (
20) and (
21), we find:
Then,
is a Cauchy Sequence in
and, hence, converges uniformly to a continuous function. Using Proposition II.1 of [
40], we conclude that the subdivision scheme
is
convergent. □
Theorem 5. Let , (this function can be discontinuous at some points),
. Then, the nonlinear subdivision scheme is convergent in , .
Proof. The inequalities (A1) imply that .
In the regularity zones, we have that
and then
. Thus, hypothesis (
21) is satisfied, and
is
convergent at the regularity zones. □
Theorem 6. Let g a continuous function and strictly monotonic or constant in , . If the linear subdivision scheme is convergent in , then: Proof. Let the first integer when we define the nonlinear subdivision scheme as:
If
,
If
,
Using the properties and , we find
If
,
If
,
Using an induction strategy (see
Appendix A at the end of the paper), we arrive at
□
Corollary 1. Under the same hypothesis of the above theorem, the nonlinear subdivision scheme is convergent in .
4. Stability of the New Family of Subdivision Schemes
Let us start by recalling the following definition,
Definition 4. A convergent nonlinear subdivision operator is called stable if there exists a constant C such that, for every pair of initial data f, , In order to prove the stability at the smooth zones, we present, only for theoretical reasons, the following subdivision scheme, which we denote by
:
Due to the fact that the linear subdivision scheme is convergent with regularity, is convergent; moreover, if the initial points come from the discretization of a piecewise continuous function, the limit function of the scheme has regularity at the regularity zones.
Theorem 7. Let , ( these functions can be discontinuous at some points), . and are two sequences defined as . If such that , then is stable in , .
Proof . We start by proving the following formula
Using the definition of
,
we find:
□
5. Numerical Experiments
In
Section 3, we studied that the contraction of the differences in the
norm implies convergence. If the contraction of the differences occurs in the
norm, the limit function is at least continuous. As it occurs in the
norm, the limit function can be discontinuous. In fact, the linear version of the scheme has
limit functions. The nonlinear version of the scheme should keep this property at smooth zones. In this section, we investigate the validity of these assertions through numerical experiments. In all the experiments that we perform, we set
.
In order to obtain a nonlinear version of the scheme, we particularize the mean
M used in the previous sections to the harmonic mean, defined as,
with
if
and
if
.
The most important properties of this mean are included in the following proposition (see [
18] for more details) and are related to the ones presented in Definition 1.
Proposition 2. For all , the harmonic mean satisfies
- 1 .
.
- 2 .
.
- 3 .
.
- 4 .
- 5 .
.
- 6 .
.
- 7 .
For , .
- 8 .
If , , and , then - 9 .
For all , there exist non-negative numbers α and β (depending on these points) such that and with a constant not depending on
Remark 2. Property 2 of Proposition 2 implies that the use of the harmonic mean (22) produces order reduction when the arguments have different sign. This problem can be solved translating both arguments (such that both are positive or negative), obtaining the mean and then translating back the result by the same amount. A typical value for the translation is two times the minimum of the absolute value of the arguments, as proposed in [41]. This technique was used in all the experiments in order to avoid order reduction close to the discontinuities. 5.1. Numerical Regularity
Following [
42], the regularity of a limit function can be evaluated numerically. We use the scheme
defined in (
3–
4) for
M equal to the arithmetic mean or the harmonic mean (
22) and
. Then, the subdivision schemes for the differences of order
and
associated with
(which can be derived due to the specific definition of
), allow to estimate the regularity constant
for
using different number of subdivision scales
j:
This expression provides an estimate for
and
such that the limit functions belong to
and
. In
Table 1, we present some numerical estimations of the regularity constant for the nonlinear scheme defined in (
3–
4) with
M equal to the harmonic mean and for the linear scheme defined in (
2). The table was obtained using from
to
levels of subdivision, taking as initial data the function
in
discretized with 16 equidistant points.
From this table, we can see that the numerical estimation of the regularity constant for the nonlinear scheme is close to the one obtained for the linear scheme [
31], i.e.,
, independently of the mean that we choose. Thus, the nonlinear perturbation introduced by the substitution of the arithmetic mean by the harmonic mean in (
3–
4), has a weak influence on the regularity at smooth zones.
5.2. Subdivision of Univariate Functions
In this subsection, we analyze the performance of the scheme when subdividing univariate functions using the arithmetic mean and the harmonic mean. In all cases, we start from an initial data vector of length equal to 128, and we apply 10 scales of subdivision.
Example 1. Let us start by the function in (23),This is a piecewise continuous function with a jump discontinuity. Figure 1 top left presents the limit function obtained by the nonlinear version of the three points binary algorithm in (3-4). The hollow green circles represent the original sampled function, and the blue points are the function obtained using subdivision. Figure 1 at the top right presents a zoom around the discontinuity. At the bottom of the same figure, we can observe the results obtained by the linear algorithm in (2). None of the results show Gibbs oscillations; however, the nonlinear algorithm manages to remove the diffusion close to the discontinuity. The numerical effect that we can appreciate is not due to diffusion but to the reduction of the order of accuracy to close to the discontinuity. Example 2. Let us continue with the function in (24),which is a piecewise polynomial with a jump discontinuity. Figure 2 shows the results obtained for the nonlinear algorithm (top) and for the linear one (bottom). The conclusions are similar: none of the algorithms introduce Gibbs oscillations, and the nonlinear algorithm removes the diffusion close to the discontinuity. Example 3. Let us finish with the piecewise constant function shown in (25),In this experiment, our aim is to show that the new non linear scheme is capable of reproducing piecewise constant functions. In Figure 3, we show the result obtained for the nonlinear scheme (top) and for the linear one (bottom). It is clear that the new new nonlinear algorithm reproduces piecewise constant functions while the linear algorithm introduces diffusion. 5.3. Analysis of the Order of Approximation Close to the Discontinuities
To say that a scheme is
convergent means that the numerical approximation of the data becomes closer to the true data as the grid is refined. In
Section 3, we analyzed that choosing
in (
3–
4) with
M equal to the harmonic mean (
22), the nonlinear scheme is convergent in the
norm. As mentioned before, this means that the limit function can be discontinuous. In this case, the mentioned convergence implies that the
diffusion in the numerical results obtained by the nonlinear algorithm in the previous subsection results from approximation errors and is due to the fact that the nonlinear scheme loses its order of accuracy to
close to the discontinuities.
In this subsection, we investigate whether this hypothesis is true through a grid refinement analysis. If the errors that we observe close to the discontinuity are due to the diffusion introduced by the scheme, they should not disappear as the grid is refined. If we start from data obtained from the discretization of a piecewise continuous function, we know that the greatest errors will be located close to the discontinuity, and thus we can choose the norm of the error to perform a grid refinement analysis.
In this way, we are sure that the conclusions that we reach about the order of accuracy of the scheme will be limited by the numerical artifacts that we observe close to the discontinuity.
Table 2 shows a grid refinement analysis performed for the piecewise continuous function in (
24) using
M equal to the harmonic mean in (
3–
4). Fixing the number of scales of subdivision
j, we refine the grid spacing and check the
norm of the error. We define the order of accuracy of the reconstruction as,
with
the error obtained with a grid spacing
h and
the error obtained with a grid spacing
. We use the discrete
norm,
In
Table 2 and
Table 3, we represent, by
i, the number of initial control points that we use and, by
j, the number of subdivision scales.
The results presented in
Table 2 show that the nonlinear scheme works like a zoom of first order near the discontinuities and preserving jumps. The table shows that the artifacts that we can see close to the discontinuity in the results obtained by the nonlinear scheme are not due to diffusion. Instead, they are approximation errors of first order.
Table 3 shows the results obtained by the linear algorithm, i.e., with
M equal to the arithmetic mean in (
3–
4). We can see how the errors do not reduce as the grid is refined. In this case, the linear algorithm is affected by diffusion. At the moment, we are working in high order schemes in order to increase the accuracy of the nonlinear schemes presented in this work.
5.4. 2D Experiments: Subdivision of Bivariate Functions
Let us consider the next bivariate function,
This is represented in
Figure 4. The subdivision of this data is represented in
Figure 5.
6. Conclusions
In this paper, we analyzed the behavior of a new family of stationary nonlinear, noninterpolatory subdivision schemes in the presence of strongly varying data and, in particular, the possible apparition of diffusion. A way to construct subdivision schemes with improved diffusion close to jump discontinuities was introduced. The convergence and stability of the new family was also analyzed. The components of the family do not present Gibbs oscillations as the masks are considered positive.
The numerical experiments presented agree with the theoretical results obtained. As far as we know, this is the first time that a binary, stable and convergent subdivision scheme without Gibbs oscillations close to jump discontinuities and reproducing piecewise constant functions has appeared in the literature.