Next Article in Journal
Higher-Order Asymptotic Numerical Solutions for Singularly Perturbed Problems with Variable Coefficients
Next Article in Special Issue
HOME: 3D Human–Object Mesh Topology-Enhanced Interaction Recognition in Images
Previous Article in Journal
Fuzzy Quality Evaluation and Analysis Model for Improving the Quality of Unleaded Gasoline to Reduce Air Pollution
Previous Article in Special Issue
Some New n-Point Ternary Subdivision Schemes without the Gibbs Phenomenon
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On a Nonlinear Three-Point Subdivision Scheme Reproducing Piecewise Constant Functions

1
Ecole Supérieure des Sciences Appliquées-Alger, Algiers 16001, Algeria
2
Laboratory of EDPNL & HM of ENS Kouba, Kouba 16308, Algeria
3
Departamento de Matemática Aplicada y Estadística, Universidad Politécnica de Cartagena, 30202 Cartagena, Spain
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(15), 2790; https://doi.org/10.3390/math10152790
Submission received: 19 June 2022 / Revised: 26 July 2022 / Accepted: 29 July 2022 / Published: 5 August 2022
(This article belongs to the Special Issue Computer-Aided Geometric Design)

Abstract

:
In this article, a nonlinear binary three-point non-interpolatory subdivision scheme is presented. It is based on a nonlinear perturbation of the three-point subdivision scheme: A new three point approximating C 2 subdivision scheme, where the convergence and the stability of this linear subdivision scheme are analyzed. It is possible to prove that this scheme does not present Gibbs oscillations in the limit functions obtained. The numerical experiments show that the linear scheme is stable even in the presence of jump discontinuities. Even though, close to jump discontinuities, the accuracy is loosed. This order reduction is equivalent to the introduction of some diffusion. Diffusion is a good property for subdivision schemes when the discontinuities are numerical, i.e., they appear when discretizing a continuous function close to high gradients. On the other hand, if the initial control points come from the discretization of a piecewise continuous function, it can be interesting that the subdivision scheme produces a piecewise continuous limit function. For instance, in the approximation of conservation laws, real discontinuities appear as shocks in the solution. The nonlinear modification introduced in this work allows to attain this objective. As far as we know, this is the first subdivision scheme that appears in the literature with these properties.

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 C 2 regularity. They analyze the following binary subdivision scheme:
( S w f k ) 2 n = a f n 1 k + b f n k + c f n + 1 k , ( S w f k ) 2 n + 1 = c f n 1 k + b f n k + a f n + 1 k ,
where { f i k } is a set of initial control points with c = 1 2 w 2 , b = 1 2 ( 1 + 2 w 2 w 2 ) and a = 1 2 ( w 1 ) 2 and w [ 0 , 1 2 ] . 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 [ x i , x i + 1 ] . Due to the convex hull property of B-splines, this scheme cannot present Gibbs phenomenon [34]. The authors prove, in [31], that, for w = 1 4 , the scheme presents a limit function with C 2 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 R + : = ( 0 , ) , R + 2 : = R + × R + , and consider the function M : R + 2 R + . 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):
( S w f k ) 2 n = f n k a d f n k + c d f n + 1 k , ( S w f k ) 2 n + 1 = f n k c d f n k + a d f n + 1 k ,
where d f n k = f n k f n 1 k denotes the first order differences. The nonlinear scheme S N L that we propose is given by:
  • If | d f n k | | d f n + 1 k | ,
    ( S N L f k ) 2 n = f n k + ( c + a ) d f n + 1 k 2 a M ( d f n k , d f n + 1 k ) , ( S N L f k ) 2 n + 1 = f n k + ( c + a ) d f n + 1 k 2 c M ( d f n k , d f n + 1 k ) .
  • If | d f n k | < | d f n + 1 k | ,
    ( S N L f k ) 2 n = f n k ( c + a ) d f n k + 2 c M ( d f n k , d f n + 1 k ) , ( S N L f k ) 2 n + 1 = f n k ( c + a ) d f n k + 2 a M ( d f n k , d f n + 1 k ) .
For M a m e a n ( x , y ) = x + y 2 , 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 S N L , in the following definition we summarize the most important properties that the function M must verify (see [23] for more details).
Definition 1. 
( M 1 )
The function M is a mean on R + 2 , and thus:
( M 1 a )
M ( x , y ) max { x , y } .
( M 1 b )
M ( x , y ) min { x , y } .
( M 1 c )
M ( x , x ) = x .
( M 2 )
M ( x h , x + h ) = x + O ( h 2 ) as h 0 uniformly with respect to x
on compact subintervals on R + .
( M 3 )
M ( t x , t y ) = t M ( x , y ) ( x , y , t R + ) .
( M 4 )
M ( x + o ( h ) , y + o ( h ) ) = M ( x , y ) + o ( h ) .
( M 5 )
M ( x , y ) x + y 2 ( x , y R + ) .
( M 6 )
M ( x , y ) γ min ( x , y ) . ( x , y R + ) for some constant γ > 1 .
( M 7 )
If x = O ( 1 ) , y = O ( 1 ) , | y x | = O ( h ) and x y > 0 , then
x + y 2 M ( x , y ) = O ( h 2 ) .
( M 8 )
for all x , x , y , y R + there exist non-negative numbers α and β (depending on these points) such that
M ( x , y ) M ( x , y ) = α ( x x ) + β ( y y ) ,
and α + β γ with a constant γ > 1 not depending on x , x , y , y .
As a consequence of ( M 5 ) , the authors presented, in [23], the following proposition:
Proposition 1. 
Suppose that M is continuous on R + 2 . If ( M 5 ) holds, then the extension of M by
M ( x , y ) : = sign ( x ) M ( | x | , | y | ) if   x y > 0 , 0 if   x y 0 ,
is continuous on R 2 .
The proof can be found in Proposition 2 of [23]. We use this extension to R 2 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 a = { a j } j Z with finite support is defined as beginning with an initial sequence of finite data f 0 = { f i 0 } i J 0 . New values are obtained through refinement at level k + 1 , denoted by f k + 1 = { f i k + 1 } i J k + 1 . The subdivision scheme is the maximal set obtained by applying the rule:
( S a f k ) i = f i k + 1 : = j J k a i 2 j f j k , then , J k + 1 = 2 J k + s u p p ( a ) = { i Z : i = 2 j + γ , j J k , γ s u p p ( a ) } .
There are two rules to define the points on the level k + 1 :
( S a f k ) 2 i = f 2 i k + 1 = γ Z a 2 γ f i γ k , i Z , ( S a f k ) 2 i + 1 = f 2 i + 1 k + 1 = γ Z a 2 γ + 1 f i γ k , i Z .
We can represent these rules using an algebraic formalism in terms of z-transforms. The symbol of the mask a = { a j } j Z is defined as a ( z ) = j Z a j z j . If a [ k ] is denoted as the k iterated symbol (see [3]), then, it is known that a [ k ] ( z ) = l = 1 k a ( z 2 l 1 ) . Therefore, if 0 l < 2 k , then
f 2 k i + l k = ( S a k f 0 ) 2 k i + l = ( S a [ k ] f 0 ) 2 k i + l = γ Z a 2 k γ + l [ k ] f i γ 0 ,
with
a j [ k ] = i Z a i [ k 1 ] a j 2 i ,
being a j [ 1 ] = a j , j . As in [19], in this paper, we impose to the coefficients of the mask the following conditions:
  • The mask a = { a j } j Z is defined as:
    1 a j 1 , U j U + V , a U , a U + V 0 , a j = 0 , j > U + V or j < U ,
    with U < 0 and V fixed integers with V 2 and U + V > 0 .
  • The scheme S a is convergent. Therefore, one necessary condition is
    γ Z a 2 γ = γ Z a 2 γ + 1 = 1 .
  • For any f C n , n 2 and any h > 0 , with
    f i 0 = f ( i h ) , i Z , if the scheme has primal parametrization , f i 0 = f i 1 2 h , i Z , if the scheme has dual parametrization ,
    and then
    max x R | ( S a f 0 ) ( x ) f ( x ) | K h n ,
    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 0 ξ h , let f be any function defined by
x ξ , f ( x ) = f ( x ) , f C n ( ] , ξ ] ) , x ξ , f ( x ) = f + ( x ) , f + C n ( [ ξ , + [ ) ,
with n 2 and f ( ξ ) > f + ( ξ ) . Let S a be a univariate stationary subdivision scheme with:
λ l [ k ] ( i ) = τ i a 2 k τ + l [ k ] , i < 0 , 0 , i = 0 , τ i a 2 k τ + l [ k ] , i > 0 ,
0 l < 2 k , being a [ k ] defined in Equation (8). Then, if λ l [ k ] ( i ) 0 , i , k ; 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 | x | max { | U 1 2 | , | U + V 2 + 1 | } h , then
| f ( x ) ( S a f h ) ( x ) | = O ( h n ) ,
with n 2 .
2 . 
If | x | max { | U 1 2 | , | U + V 2 + 1 | } h , there exists α h = O ( h ) such that
f 1 , h α h f + ( h ) α h ( S a f h ) ( x ) f ( 0 ) + α h = f 0 , h + α h .
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 if
f l 1 ( Z ) , S f L 1 ( R ) s u c h   t h a t lim j + | | A S j f S f | | L 1 ( R ) = 0 ,
where A S j denotes the piecewise affine interpolation of S j f at the grid points.
Definition 3. 
A convergent subdivision scheme is stable if
C < + s u c h   t h a t f 0 , g 0 l 1 ( Z ) | | S f S g | | L 1 ( R ) C | | f 0 g 0 | | l 1 ( Z ) .
Let us first rewrite the scheme S N L (3) and (4) as a particular perturbation of the linear subdivision scheme S L (projection), which is defined as,
( S L f k ) 2 n = ( S L f k ) 2 n + 1 = f n k .
We can see that the nonlinear scheme in (3) and (4) can be written in terms of the linear scheme S L . For all f l 1 ( Z ) , we have
S N L f = S L f + F ( d f ) ,
where ( S L f ) 2 n = ( S L f ) 2 n + 1 = f n , and F is the nonlinear operator
F ( d f ) 2 n = ( c + a ) d f n + 1 2 a M ( d f n , d f n + 1 ) , if | d f n | | d f n + 1 | , ( c + a ) d f n + 2 c M ( d f n , d f n + 1 ) , if | d f n | < | d f n + 1 | .
F ( d f ) 2 n + 1 = ( c + a ) d f n + 1 2 c M ( d f n , d f n + 1 ) , if | d f n | | d f n + 1 | , ( c + a ) d f n + 2 a M ( d f n , d f n + 1 ) , if | d f n | < | d f n + 1 | .
We also need to recall the following result established in [14]. Let S N L be a subdivision scheme defined by
f l 1 ( Z ) , n Z ( S N L f ) n = ( S f ) n + F ( δ f ) n ,
where F is a nonlinear operator defined on l 1 ( Z ) , δ is a linear and continuous operator on l 1 ( Z ) and S is a linear subdivision scheme convergent in l 1 ( Z ) .
Theorem 2. 
If S is a linear convergent subdivision scheme and if S N L , F and δ verify
K > 0 s u c h   t h a t d l 1 ( Z ) | | F ( d ) | | l 1 ( Z ) K | | d | | l 1 ( Z ) ,
L > 0 , c < 1 s u c h   t h a t f l 1 ( Z ) | | δ S N L L ( f ) | | l 1 ( Z ) c | | δ f | | l 1 ( Z ) ,
then the subdivision scheme S N L defined by (16) is convergent.
Using Theorem 2, we can prove the following theorem:
Theorem 3. 
The nonlinear subdivision scheme S N L is l 1 ( Z ) convergent for all w [ 0 , 1 2 ] .
Proof. 
The perturbation F defined in (14), (15) can be bounded using that | M ( x , y ) | x + y 2 .
| F ( d f k ) 2 n | ( a + 2 c ) | d f n k | + ( 2 a + c ) | d f n + 1 k | , | F ( d f k ) 2 n + 1 | ( 2 a + c ) | d f n k | + ( a + 2 c ) | d f n + 1 k | .
Then
F ( d f k ) l 1 ( Z ) = h k n Z | F ( d f k ) n | 6 ( a + c ) d f k l 1 ( Z ) ,
that satisfies precisely hypothesis (17).
Now, we consider hypothesis (18), which is related to the contraction of first order differences ( d f k ) .
We have 16 different cases (see the Appendix A at the end of the paper).
In all of the cases, we obtain:
n Z , k N , C n k ] 0 , 2 [ : | d S f 2 n 2 k | + | d S f 2 n 1 k | + | d S f 2 n k | + | d S f 2 n + 1 k | + | d S f 2 n + 2 k | C n k | d f n k | + other terms .
Thus, for all f l ( Z ) :
d ( S N L ( f k ) l 1 = h k + 1 n Z | d ( S N L ( f k ) n | 1 + 2 a + 2 c 2 d f k l 1 , 2 w 2 2 w + 1 2 d f k l 1 .
If w [ 0 , 1 2 ] , then 2 w 2 2 w + 1 2 < 1 .
Thus, hypothesis (18) is satisfied, and S N L is l 1 ( Z ) convergent. □
Remark 1. 
Other terms in the last theorem means that they do not depend on d f n k .
In this paper, we prove the stability only at the regularity parts using the l ( Z ) norm; thus, we also have to prove the convergence at these parts using the l ( Z ) norm.
Theorem 4. 
( l ( Z ) convergence)
If F and δ verify:
M > 0 s u c h   t h a t d l ( Z ) F ( δ ) l M δ l ( Z ) .
f l ( Z ) δ S N L ( f ) l ( Z ) 1 2 δ f l ( Z ) + O ( h ) .
Then, the subdivision scheme S N L is uniformly convergent in l ( Z ) .
Proof. 
Let φ be the limit function of the linear subdivision S. j > 0 , f l ( Z ) , we consider the continuous function f j ( . ) as:
f j ( x ) = n f n j φ ( 2 j x n ) x R
where f j is the sequence S N L j f . We have
f j + 1 ( x ) f j ( x ) = n f n j + 1 φ ( 2 j + 1 x n ) n f n j φ ( 2 j x n ) ,
Using lemma 2.3 in [29], the definition of S N L and hypotheses (20) and (21), we find:
f j + 1 ( x ) f j ( x ) = n ( f n j + 1 ( S f j ) n ) φ ( 2 j + 1 x n ) . f j + 1 f j L M δ ( S N L f k 1 ) l , M [ 1 2 δ ( f k 1 ) l , + o ( h k 1 ) ] , M 2 k δ ( f 0 ) l + M i = 0 k 1 o ( h i ) 2 k 1 i M 2 k δ ( f 0 ) l + M 2 k 1 i = 0 k 1 2 i o ( h i ) M 2 k δ ( f 0 ) l + M k 2 k 1 o ( h 0 )
Then, ( f k ( . ) ) is a Cauchy Sequence in l ( Z ) and, hence, converges uniformly to a continuous function. Using Proposition II.1 of [40], we conclude that the subdivision scheme S N L is l ( Z ) convergent. □
Theorem 5. 
Let g C 1 ( [ a i , b i ] ) , i (this function can be discontinuous at some points),
f 0 = g ( n h ) n Z . Then, the nonlinear subdivision scheme S N L is l ( Z ) convergent in [ a i , b i ] , i .
Proof. 
The inequalities (A1) imply that F ( d f k ) l ( 3 a + 3 c ) d f k l .
In the regularity zones, we have that d f k h k = f + o ( h k ) and then d f k + 1 = d f k 2 + o ( h ) . Thus, hypothesis (21) is satisfied, and S N L is l convergent at the regularity zones. □
Theorem 6. 
Let g a continuous function and strictly monotonic or constant in [ a , b ] R , f 0 = g ( n h ) n Z . If the linear subdivision scheme S w is C 1 convergent in [ a , b ] , then:
k > k 0 , | ( f w k ) n ( f N L k ) n | = o ( h k 2 ) , in [ a , b ] .
Proof. 
Let k 0 the first integer when we define the nonlinear subdivision scheme as:
  • If | d ( f w ) n k 0 | | d ( f w ) n + 1 k 0 | ,
    ( f N L ) 2 n k 0 + 1 = ( S N L f k 0 ) 2 n = ( f w ) n k 0 + ( c + a ) d ( f w ) n + 1 k 0 2 a M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) , ( f N L ) 2 n + 1 k 0 + 1 = ( S N L f k 0 ) 2 n + 1 = ( f w ) n k 0 + ( c + a ) d ( f w ) n + 1 k 0 2 c M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) .
  • If | d ( f w ) n k 0 | < | d ( f w ) n + 1 k 0 | ,
    ( f N L ) 2 n k 0 + 1 = ( S N L f k 0 ) 2 n = ( f w ) n k 0 ( c + a ) d ( f w ) n + 1 k 0 + 2 c M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) , ( f N L ) 2 n + 1 k 0 + 1 = ( S N L f k 0 ) 2 n + 1 = ( f w ) n k 0 ( c + a ) d ( f w ) n + 1 k 0 + 2 a M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) .
Using the properties ( M 4 ) and ( M 7 ) , we find
  • If | d ( f w ) n k 0 | | d ( f w ) n + 1 k 0 | ,
    | ( f w k 0 + 1 ) 2 n ( f N L k 0 + 1 ) 2 n | = 2 a | d ( f w ) n k 0 + d ( f w ) n + 1 k 0 2 M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) | = o ( h 2 ) , | ( f w k 0 + 1 ) 2 n + 1 ( f N L k 0 + 1 ) 2 n + 1 | = 2 c | d ( f w ) n k 0 + d ( f w ) n + 1 k 0 2 M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) | = o ( h 2 ) .
  • If | d ( f w ) n k 0 | < | d ( f w ) n + 1 k 0 | ,
    | ( f w k 0 + 1 ) 2 n ( f N L k 0 + 1 ) 2 n | = 2 c | d ( f w ) n k 0 + d ( f w ) n + 1 k 0 2 M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) | = o ( h 2 ) , | ( f w k 0 + 1 ) 2 n + 1 ( f N L k 0 + 1 ) 2 n + 1 | = 2 a | d ( f w ) n k 0 + d ( f w ) n + 1 k 0 2 M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) | = o ( h 2 ) .
Using an induction strategy (see Appendix A at the end of the paper), we arrive at
k > k 0 , | ( f w k ) n ( f N L k ) n | = o ( h k 2 ) , in [ a , b ] .
 □
Corollary 1. 
Under the same hypothesis of the above theorem, the nonlinear subdivision scheme S N L is convergent in l ( [ a , b ] ) .

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, g l ( Z ) ,
S f S g L ( R ) C f g l ( Z ) .
In order to prove the stability at the smooth zones, we present, only for theoretical reasons, the following subdivision scheme, which we denote by S ¯ :
S ¯ f n k = f n k , | f n k f n 1 k | = O ( 1 ) or | f n k f n + 1 k | = O ( 1 ) , ( S L f k ) n , otherwise .
Due to the fact that the linear subdivision scheme S L is convergent with C 2 regularity, S ¯ is convergent; moreover, if the initial points come from the discretization of a piecewise continuous function, the limit function of the scheme S ¯ has C 2 regularity at the regularity zones.
Theorem 7. 
Let g , g ˜ C 1 ( [ a i , b i ] ) , i ( these functions can be discontinuous at some points), h 0 = x n + 1 0 x n 0 n . f 0 and f ˜ 0 are two sequences defined as f 0 = g ( n h ) n Z , f ˜ 0 = g ˜ ( n h ) n Z . If C such that h 0 C f 0 f ˜ 0 l ( Z ) , then S N L is stable in [ a i , b i ] , i .
Proof . 
We start by proving the following formula
f 0 l ( Z ) , h 0 > 0 , S N L f 0 S ¯ N L f 0 l ( Z ) = O ( h 0 )
Using the definition of S N L , S ¯ N L we find:
  • In the fist iteration f 1 f ¯ 1 l ( Z ) = S N L f 0 S ¯ N L f 0 l ( Z ) = O ( h 0 ) .
  • In the next iterations, we use the stability of the linear scheme:
    S N L f 1 S ¯ N L f ¯ 1 l ( Z ) S N L f 1 S ¯ N L f 1 l ( Z ) + S ¯ N L f 1 S ¯ N L f ¯ 1 l ( Z ) . O ( h 1 ) + O ( h 0 ) S N L f 0 S ¯ N L f 0 l ( Z ) C i h i = O ( h 0 )
  • Using the stability of the linear subdivision, we obtain the result:
    S N L ( f 0 ) S N L ( f ˜ 0 ) l ( Z ) S N L ( f 0 ) S ¯ N L ( f 0 ) + S ¯ N L ( f 0 ) S ¯ N L ( f ˜ 0 ) + S ¯ N L ( f ˜ 0 ) S N L ( f ˜ 0 ) O ( h 0 ) + C 1 f 0 f ˜ 0 l ( Z ) ( C + C 1 ) f 0 f ˜ 0 l ( Z ) .
 □

5. Numerical Experiments

In Section 3, we studied that the contraction of the differences in the l 1 ( Z ) norm implies convergence. If the contraction of the differences occurs in the l ( Z ) norm, the limit function is at least continuous. As it occurs in the l 1 ( Z ) norm, the limit function can be discontinuous. In fact, the linear version of the scheme has C 2 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 w = 1 / 4 .
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,
( x , y ) R 2 H ( x , y ) : = x y x + y ( sgn ( x y ) + 1 ) ,
with sgn ( x ) = 1 if x 0 and sgn ( x ) = 1 if x < 0 .
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 ( x , y ) R 2 , the harmonic mean H ( x , y ) satisfies
1 . 
H ( x , y ) = H ( y , x ) .
2 . 
H ( x , y ) = 0 i f x y 0 .
3 . 
H ( x , y ) = H ( x , y ) .
4 . 
H ( x , y ) = s i g n ( x ) + s i g n ( y ) 2 m i n ( | x | , | y | ) 1 + x y x + y .
5 . 
| H ( x , y ) | max ( | x | , | y | ) .
6 . 
| H ( x , y ) | 2 min ( | x | , | y | ) .
7 . 
For x , y > 0 , min ( x , y ) H ( x , y ) x + y 2 .
8 . 
If x = O ( 1 ) , y = O ( 1 ) , | y x | = O ( h ) and x y > 0 , then
x + y 2 H ( x , y ) = O ( h 2 ) .
9 . 
For all x , x , y , y R + , there exist non-negative numbers α and β (depending on these points) such that
H ( x , y ) H ( x , y ) = α ( x x ) + β ( y y )
and α + β γ with a constant γ = 2 not depending on x , x , y , y .
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 S N L defined in (34) for M equal to the arithmetic mean or the harmonic mean (22) and w = 1 / 4 . Then, the subdivision schemes for the differences of order k = 1 and k = 2 associated with S N L (which can be derived due to the specific definition of S N L ), allow to estimate the regularity constant β k for k = 1 , 2 using different number of subdivision scales j:
β k = log 2 2 k | | ( S k j + 1 f ) n + 1 ( S k j + 1 f ) n | | | | ( S k j f ) n + 1 ( S k j f ) n | | .
This expression provides an estimate for β 1 and β 2 such that the limit functions belong to C 1 + β 1 and C 2 + β 2 . In Table 1, we present some numerical estimations of the regularity constant for the nonlinear scheme defined in (34) with M equal to the harmonic mean and for the linear scheme defined in (2). The table was obtained using from j = 5 to j = 10 levels of subdivision, taking as initial data the function f ( x ) = e x in [ 1 , 1 ] 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., C 2 , independently of the mean that we choose. Thus, the nonlinear perturbation introduced by the substitution of the arithmetic mean by the harmonic mean in (34), 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 S N L 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),
g ( x ) = sin ( 5 x ) , 0 x 0.5 , 1 + sin 2 ( 3 x ) , 0.5 < x 1 .
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 O ( h ) close to the discontinuity.
Example 2. 
Let us continue with the function in (24),
g ( x ) = 3 x 6 x 4 2 x 2 3 x , 0 x 0.5 , 3 x 6 + x 4 + x 2 + 3 x , 0.5 < x 1 ,
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),
g ( x ) = 2 , 0 x 0.5 , 5 , 0.5 < x 1 .
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 w = 1 / 4 in (34) with M equal to the harmonic mean (22), the nonlinear scheme is convergent in the l 1 ( Z ) 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 O ( h ) 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 l ( Z ) 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 (34). Fixing the number of scales of subdivision j, we refine the grid spacing and check the l ( Z ) norm of the error. We define the order of accuracy of the reconstruction as,
o r d e r = l o g 2 e i e i + 1 ,
with e i the error obtained with a grid spacing h and e i + 1 the error obtained with a grid spacing h / 2 . We use the discrete l norm,
| | f j | | = max n Z { | f n j | } .
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 (34). 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,
f ( x ) = cos ( π x ) cos ( π y ) , if ( x + 1 2 ) 2 + ( y 1 2 ) 2 < 1 1 cos ( π x ) sin ( π y ) , if ( x + 1 2 ) 2 + ( y 1 2 ) 2 1 .
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.

Author Contributions

Conceptualization, S.A.; Formal analysis, S.Z.; Methodology, J.R.; Supervision, S.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by project 20928/PI/18 (Proyecto financiado por la Comunidad Autónoma de la Región de Murcia a través de la convocatoria de Ayudas a proyectos para el desarrollo de investigación científica y técnica por grupos competitivos, incluida en el Programa Regional de Fomento de la Investigación Científica y Técnica (Plan de Actuación 2018) de la Fundación Séneca-Agencia de Ciencia y Tecnología de la Región de Murcia) and by the national research project PID2019-108336GB-I00 (MINECO/FEDER).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof of Theorem 3. 
The perturbation F defined in (14), (15) can be bounded using that | M ( x , y ) | x + y 2 .
| F ( d f k ) 2 n | ( a + 2 c ) | d f n k | + ( 2 a + c ) | d f n + 1 k | , | F ( d f k ) 2 n + 1 | ( 2 a + c ) | d f n k | + ( a + 2 c ) | d f n + 1 k | .
Then,
F ( d f k ) l 1 ( Z ) = h k n Z | F ( d f k ) n | 6 ( a + c ) d f k l 1 ( Z ) ,
that satisfies precisely hypothesis (17).
Now, we consider hypothesis (18), which is related to the contraction of first order differences ( d f k ) .
We have 16 different cases:
Case A :   | df n k | > | df n + 1 k | , | df n 1 k | > | df n k | case A 1 :   d f n k > 0 , d f n + 1 k > 0 , d f n 1 k > 0 , case A 2 :   d f n k > 0 , d f n + 1 k > 0 , d f n 1 k < 0 , case A 3 :   d f n k > 0 , d f n + 1 k < 0 , d f n 1 k > 0 , case A 4 :   d f n k > 0 , d f n + 1 k < 0 , d f n 1 k < 0 .
Case B : | df n k | > | df n + 1 k | , | df n 1 k | < | df n k | case B 1 :   d f n k > 0 , d f n + 1 k > 0 , d f n 1 k > 0 , case B 2 :   d f n k > 0 , d f n + 1 k > 0 , d f n 1 k < 0 , case B 3 :   d f n k > 0 , d f n + 1 k < 0 , d f n 1 k > 0 , case B 4 :   d f n k > 0 , d f n + 1 k < 0 , d f n 1 k < 0 . Case C :   | df n k | < | df n + 1 k | , | df n 1 k | > | df n k | case C 1 :   d f n k > 0 , d f n + 1 k > 0 , d f n 1 k > 0 , case C 2 :   d f n k > 0 , d f n + 1 k > 0 , d f n 1 k < 0 , case C 3 :   d f n k > 0 , d f n + 1 k < 0 , d f n 1 k > 0 , case C 4 :   d f n k > 0 , d f n + 1 k < 0 , d f n 1 k < 0 . Case D :   | df n k | < | df n + 1 k | , | df n 1 k | < | df n k | case D 1 :   d f n k > 0 , d f n + 1 k > 0 , d f n 1 k > 0 , case D 2 :   d f n k > 0 , d f n + 1 k > 0 , d f n 1 k < 0 , case D 3 :   d f n k > 0 , d f n + 1 k < 0 , d f n 1 k > 0 , case D 4 :   d f n k > 0 , d f n + 1 k < 0 , d f n 1 k < 0 .
Case A1: considering the cases where each d f n k appears—that is,
| ( d S f k ) 2 n 2 | = | d f n 1 k + ( a + c ) ( d f n k d f n 1 k ) 2 a M ( d f n 1 k , d f n k ) + 2 c M ( d f n 2 k , d f n 1 k ) | d f n 1 k 2 a M ( d f n 1 k , d f n k ) + ( a + c ) d f n k + | 2 c M ( d f n 2 k , d f n 1 k ) ( a + c ) d f n 1 k | | d S f 2 n 1 k | = 2 ( a c ) M ( d f n 1 k , d f n k ) . | ( d S f ) 2 n | = d f n k 2 a M ( d f n k , d f n + 1 k ) + 2 c M ( d f n 1 k , d f n k ) ( a + c ) d f n k + ( a + c ) d f n + 1 k . | d S f 2 n + 1 k | = 2 ( a c ) M ( d f n k , d f n + 1 k ) . | d S f 2 n + 2 k | = d f n + 1 k + ( a + c ) ( d f n + 2 k d f n + 1 k ) 2 a M ( d f n + 1 k , d f n + 2 k ) + 2 c M ( d f n k , d f n + 1 k ) 2 c M ( d f n k , d f n + 1 k ) + | d f n + 1 k + ( a + c ) ( d f n + 2 k d f n + 1 k ) 2 a M ( d f n + 1 k , d f n + 2 k ) |
Then,
| d S f 2 n 2 k | + | d S f 2 n 1 k | + | d S f 2 n k | + | d S f 2 n + 1 k | + | d S f 2 n + 2 k | | d f n k | + other terms .
Using the same strategy, we have:
n Z , k N , C n k ] 0 , 2 [ : | d S f 2 n 2 k | + | d S f 2 n 1 k | + | d S f 2 n k | + | d S f 2 n + 1 k | + | d S f 2 n + 2 k | C n k | d f n k | + other terms .
where the different values of C n k are presented in Table A1.
Table A1. Values of C n k for the cases A , B , C and D considered in the proof.
Table A1. Values of C n k for the cases A , B , C and D considered in the proof.
Case C n k Case C n k Case C n k Case C n k
A 1 1 B 1 1 + 2 a + 2 c C 1 1 D 1 1 + 2 a c
A 2 1 B 2 1 + 2 a + 2 c C 2 1 a c D 2 1
A 3 1 B 3 1 + 2 a + 2 c C 3 1 a c D 3 1
A 4 1 B 4 1 + 2 a + 2 c C 4 1 2 ( a + c ) D 4 1
We can see that, for all f l ( Z ) ,
d ( S N L ( f k ) l 1 = h k + 1 n Z | d ( S N L ( f k ) n | 1 + 2 a + 2 c 2 d f k l 1 , 2 w 2 2 w + 1 2 d f k l 1 .
If w [ 0 , 1 2 ] , then 2 w 2 2 w + 1 2 < 1 .
Thus, hypothesis (18) is satisfied, and S N L is l 1 ( Z ) convergent. □
Proof of Theorem 6. 
Let k 0 the first integer when we define the nonlinear subdivision scheme as:
  • If | d ( f w ) n k 0 | | d ( f w ) n + 1 k 0 | ,
    ( f N L ) 2 n k 0 + 1 = ( S N L f k 0 ) 2 n = ( f w ) n k 0 + ( c + a ) d ( f w ) n + 1 k 0 2 a M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) , ( f N L ) 2 n + 1 k 0 + 1 = ( S N L f k 0 ) 2 n + 1 = ( f w ) n k 0 + ( c + a ) d ( f w ) n + 1 k 0 2 c M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) .
  • If | d ( f w ) n k 0 | < | d ( f w ) n + 1 k 0 | ,
    ( f N L ) 2 n k 0 + 1 = ( S N L f k 0 ) 2 n = ( f w ) n k 0 ( c + a ) d ( f w ) n + 1 k 0 + 2 c M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) , ( f N L ) 2 n + 1 k 0 + 1 = ( S N L f k 0 ) 2 n + 1 = ( f w ) n k 0 ( c + a ) d ( f w ) n + 1 k 0 + 2 a M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) .
Using the properties ( M 4 ) and ( M 7 ) , we find
  • If | d ( f w ) n k 0 | | d ( f w ) n + 1 k 0 | ,
    | ( f w k 0 + 1 ) 2 n ( f N L k 0 + 1 ) 2 n | = 2 a | d ( f w ) n k 0 + d ( f w ) n + 1 k 0 2 M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) | = o ( h 2 ) , | ( f w k 0 + 1 ) 2 n + 1 ( f N L k 0 + 1 ) 2 n + 1 | = 2 c | d ( f w ) n k 0 + d ( f w ) n + 1 k 0 2 M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) | = o ( h 2 ) .
  • If | d ( f w ) n k 0 | < | d ( f w ) n + 1 k 0 | ,
    | ( f w k 0 + 1 ) 2 n ( f N L k 0 + 1 ) 2 n | = 2 c | d ( f w ) n k 0 + d ( f w ) n + 1 k 0 2 M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) | = o ( h 2 ) , | ( f w k 0 + 1 ) 2 n + 1 ( f N L k 0 + 1 ) 2 n + 1 | = 2 a | d ( f w ) n k 0 + d ( f w ) n + 1 k 0 2 M ( d ( f w ) n k 0 , d ( f w ) n + 1 k 0 ) | = o ( h 2 ) .
Now, we suppose that, for some k, we have | ( f N L k ) n ( f w k ) n | = o ( h k 2 ) , and we prove that | ( f N L k + 1 ) n ( f w k + 1 ) n | = o ( h k 2 ) .
In the iteration k + 1 , the nonlinear subdivision scheme becomes
  • If | ( d f w ) n k | | ( d f w ) n + 1 k | ,
    ( f N L ) 2 n k + 1 = ( f N L ) n k + ( c + a ) d ( f N L ) n + 1 k 2 a M ( d ( f N L ) n k , d ( f N L ) n + 1 k ) , = ( f w ) n k + o ( h 2 ) + ( c + a ) d ( f w ) n + 1 k + o ( h 2 ) 2 a M ( d ( f w ) n k + + o ( h 2 ) , d ( f w ) n + 1 k + o ( h 2 ) ) , = ( f w ) n k + ( c + a ) d ( f w ) n + 1 k 2 a M ( d ( f w ) n k , d ( f w ) n + 1 k ) + o ( h 2 ) , ( f N L ) 2 n + 1 k + 1 = ( f N L ) n k + ( c + a ) d ( f N L ) n + 1 k 2 c M ( d ( f N L ) n k , d ( f N L ) n + 1 k ) , = ( f w ) n k + o ( h 2 ) + ( c + a ) d ( f w ) n + 1 k + o ( h 2 ) 2 c M ( d ( f w ) n k + o ( h 2 ) , d ( f w ) n + 1 k + o ( h 2 ) ) , = ( f w ) n k + ( c + a ) d ( f w ) n + 1 k 2 c M ( d ( f w ) n k , d ( f w ) n + 1 k ) + o ( h 2 ) .
    and
    | ( f w k + 1 ) 2 n ( f N L k + 1 ) 2 n | = 2 a | d ( f w ) n k + d ( f w ) n + 1 k 2 M ( d ( f w ) n k , d ( f w ) n + 1 k ) + o ( h 2 ) | = o ( h 2 ) , | ( f w k + 1 ) 2 n + 1 ( f N L k + 1 ) 2 n + 1 | = 2 c | d ( f w ) n k + d ( f w ) n + 1 k 2 M ( d ( f w ) n k , d ( f w ) n + 1 k ) + o ( h 2 ) | = o ( h 2 ) .
  • If | ( d f w ) n k | < | ( d f w ) n + 1 k | ,
    ( f N L ) 2 n k + 1 = ( f N L ) n k ( c + a ) d ( f N L ) n + 1 k + 2 c M ( d ( f N L ) n k , d ( f N L ) n + 1 k ) , = ( f w ) n k + o ( h 2 ) ( c + a ) d ( f w ) n + 1 k + o ( h 2 ) + 2 c M ( d ( f w ) n k + o ( h 2 ) , d ( f w ) n + 1 k + o ( h 2 ) ) , = ( f w ) n k ( c + a ) d ( f w ) n + 1 k + 2 c M ( d ( f w ) n k , d ( f w ) n + 1 k ) + o ( h 2 ) , ( f N L ) 2 n + 1 k + 1 = ( f N L ) n k ( c + a ) d ( f N L ) n + 1 k + 2 a M ( d ( f N L ) n k , d ( f N L ) n + 1 k ) . = ( f w ) n k + o ( h 2 ) ( c + a ) d ( f w ) n + 1 k + o ( h 2 ) + 2 a M ( d ( f w ) n k + o ( h 2 ) , d ( f w ) n + 1 k + o ( h 2 ) ) . = ( f w ) n k ( c + a ) d ( f w ) n + 1 k + 2 a M ( d ( f w ) n k , d ( f w ) n + 1 k ) + o ( h 2 ) .
    and
    | ( f w k + 1 ) 2 n ( f N L k + 1 ) 2 n | = 2 c | d ( f w ) n k + d ( f w ) n + 1 k 2 M ( d ( f w ) n k , d ( f w ) n + 1 k ) + o ( h 2 ) | = o ( h 2 ) , | ( f w k + 1 ) 2 n + 1 ( f N L k + 1 ) 2 n + 1 | = 2 a | d ( f w ) n k + d ( f w ) n + 1 k 2 M ( d ( f w ) n k , d ( f w ) n + 1 k ) + o ( h 2 ) | = o ( h 2 ) .
Thus,
k > k 0 , | ( f w k ) n ( f N L k ) n | = o ( h k 2 ) , in [ a , b ] .
 □

References

  1. Dyn, N.; Gregory, J.; Levin, D. A four-point interpolatory subdivision scheme for curve design. Comput. Aided Geom. Des. 1987, 4, 257–268. [Google Scholar] [CrossRef]
  2. Beccari, C.; Casciola, G.; Romani, L. An interpolating 4-point C2 ternary non-stationary subdivision scheme with tension control. Comput. Aided Geom. Des. 2007, 24, 210–219. [Google Scholar] [CrossRef]
  3. Dyn, N.; Levin, D. Subdivision Schemes in Geometric Modelling. Act. Num. 2002, 11, 73–144. [Google Scholar] [CrossRef]
  4. Hassan, M.F.; Ivrissimtzis, I.P.; Dodgson, N.A.; Sabin, M.A. An interpolating 4-point ternary stationary subdivision scheme. Comput. Aided Geom. Des. 2002, 19, 1–18. [Google Scholar] [CrossRef]
  5. Hassan, M.F.; Dodgson, N.A. Ternary and three-point univariate subdivision schemes. In Curve and Surface Fitting: Sant-Malo 2002; Cohen, A., Merrien, J.-L., Schumaker, L.L., Eds.; Nashboro Press: Brentwood, CA, USA, 2003; pp. 199–208. [Google Scholar]
  6. Jeon, M.; Han, D.; Park, K.; Choi, G. Ternary univariate curvature-preserving subdivision. J. Appl. Math. Comput. 2005, 18, 235–246. [Google Scholar] [CrossRef]
  7. Kwan, P.K.; Byung-Gook, L.; Gang, J.Y. A ternary 4-point approximating subdivision scheme. App. Math. Comp. 2007, 190, 1563–1573. [Google Scholar]
  8. Siddiqi, S.; Salam, W.; Rehan, K. Binary 3-point and 4-point non-stationary subdivision schemes using hyperbolic function. Appl. Math. Comput. 2015, 258, 120–129. [Google Scholar] [CrossRef]
  9. Siddiqi, S.; Younis, M. Construction of m-point binary approximating subdivision scheme. Appl. Math. Lett. 2013, 26, 337–343. [Google Scholar] [CrossRef]
  10. Siddiqi, S.; Rehan, K. Improved binary four point subdivision scheme and new corner cutting scheme. Comput. Math. Appl. 2010, 59, 2647–2657. [Google Scholar]
  11. Sharon, N.; Dyn, N. Bivariate interpolation based on univariate subdivision schemes. J. Approx. Theory 2012, 164, 709–730. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, H.; Qin, K. Improved ternary subdivision interpolation scheme. Tsinghua Sci. Technol. 2005, 10, 128–132. [Google Scholar] [CrossRef]
  13. Zheng, H.; Zhenglin, Y.; Zuoping, C.; Hongxing, Z. A controllable ternary interpolatory subdivision scheme. Int. J. CAD/CAM 2005, 5, 9. [Google Scholar]
  14. Amat, S.; Dadourian, K.; Liandrat, J. Analysis of a class of nonlinear subdivision schemes and associated multiresolution transforms. Adv. Comput. Math. 2011, 34, 253–277. [Google Scholar] [CrossRef]
  15. Amat, S.; Dadourian, K.; Liandrat, J. On a nonlinear subdivision scheme avoiding Gibbs oscillations and convergence towards Cs functions with s>1. Math. Comp. 2010, 80, 959–971. [Google Scholar] [CrossRef] [Green Version]
  16. Amat, S.; Dadourian, K.; Liandrat, J. On the convergence of various subdivision schemes using a perturbation theorem. In Curves and Surfaces Fitting: Avignon 2006; Nashboro Press: Brentwood, TN, USA, 2006; pp. 1–10. [Google Scholar]
  17. Amat, S.; Dadourian, K.; Liandrat, J. On a nonlinear 4-point ternary and interpolatory multiresolution scheme eliminating the Gibbs phenomenon. Int. J. Numer. Anal. Model. 2010, 7, 261–280. [Google Scholar]
  18. Amat, S.; Liandrat, J. On the stability of the PPH nonlinear multiresolution. Appl. Comp. Harm. Anal. 2005, 18, 198–206. [Google Scholar] [CrossRef] [Green Version]
  19. Amat, S.; Ruiz, J.; Trillo, J.C.; Yáñez, D.F. Analysis of the Gibbs phenomenon in stationary subdivision schemes. Appl. Math. Lett. 2018, 76, 157–163. [Google Scholar] [CrossRef]
  20. Cohen, A.; Dyn, N.; Matei, B. Quasilinear subdivision schemes with applications to ENO interpolation. Appl. Comp. Harm. Anal. 2003, 15, 89–116. [Google Scholar] [CrossRef] [Green Version]
  21. Dyn, N.; Hormann, K.; Mancinelli, C. Non-uniform interpolatory subdivision schemes with improved smoothness. Comput. Aided Geom. Des. 2022, 94, 102083. [Google Scholar] [CrossRef]
  22. Floater, M.S.; Michelli, C.A. Nonlinear stationary subdivision. In Approximation Theory: In Memory of A.K. Varna; Govil, N.K., Mohapatra, N., Nashed, Z., Sharma, A., Szabados, J., Eds.; CRC Press: Boca Raton, FL, USA, 1998; pp. 209–224. [Google Scholar]
  23. Guessab, A.; Moncayo, M.; Schmeisser, G. A class of nonlinear four-point subdivision schemes—Properties in terms of conditions. Adv. Comput. Math. 2012, 37, 151–190. [Google Scholar] [CrossRef]
  24. Harizanov, S.; Oswald, P. Stability of nonlinear subdivision and multiscale transforms. Constr. Approx. 2010, 31, 359–393. [Google Scholar] [CrossRef]
  25. Lian, J.A.; Wang, Y.; Yang, Y. Circular nonlinear subdivision schemes for curve design. Appl. Appl. Math. Int. J. 2009, 4, 1. [Google Scholar]
  26. Ortiz, P.; Trillo, J.C. Analysis of a New Nonlinear Interpolatory Subdivision Scheme on Quasi-Uniform Grids. Mathematics 2021, 9, 1320. [Google Scholar] [CrossRef]
  27. Oswald, P. Smoothness of Nonlinear Median-Interpolation Subdivision. Adv. Comput. Math. 2004, 20, 401–423. [Google Scholar]
  28. Deslauriers, G.; Dubuc, S. Symmetric iterative interpolation processes. Constr. Approx. 1989, 5, 49–68. [Google Scholar] [CrossRef]
  29. Dyn, N. Subdivision schemes in computer aided geometric design. In Advances in Numerical Analysis II., Subdivision algorithms and radial functions; Light, W.A., Ed.; Oxford University Press Prentice-Hall: Hoboken, NJ, USA, 1992; pp. 36–104. [Google Scholar]
  30. Rabia, H.; Ghulam, M. Construction and analysis of binary subdivision schemes for curves and surfaces originated from Chaikin points. Int. J. Anal. 2016, 2016, 1092476. [Google Scholar]
  31. Siddiqi, S.; Nadeem, A. A new three point approximating C2 subdivision scheme. Appl. Math. Lett. 2007, 20, 707–711. [Google Scholar] [CrossRef] [Green Version]
  32. Tan, J.; Sun, J.; Tong, G. A non-stationary binary three-point approximating subdivision scheme. Appl. Math. Comput. 2016, 276, 37–43. [Google Scholar] [CrossRef]
  33. Gottlieb, D.; Shu, C.-W. On the Gibbs phenomenon and its resolution. SIAM Rev. 1997, 39, 644–668. [Google Scholar] [CrossRef] [Green Version]
  34. Bartels, R.H.; Beatty, J.C.; Barsky, B.A. An Introduction to Splines for Use in Computer Graphics and Geometric Modeling. The Morgan Kaufmann: San Francisco, CA, USA, 1995. [Google Scholar]
  35. Amat, S.; Donat, R.; Liandrat, J.; Trillo, J.C. Analysis of a fully nonlinear multiresolution scheme for image processing. Found. Comput. Math. 2006, 6, 193–225. [Google Scholar]
  36. Amat, S.; Liandrat, J.; Ruiz, J.; Trillo, J.C. On a nonlinear mean and its application to image compression using multiresolution schemes. Numer. Algorithms 2016, 71, 729–752. [Google Scholar] [CrossRef]
  37. Amat, S.; Liandrat, J.; Ruiz, J.; Trillo, J.C. On a family of nonlinear cell-average multiresolution schemes for image processing: An experimental study. Math. Comput. Simul. 2015, 118, 30–48. [Google Scholar] [CrossRef]
  38. Amat, S.; Ruiz, J.; Trillo, J.C. Adaptive interpolation of images using a new nonlinear cell-average scheme. Math. Comput. Simul. 2012, 82, 1586–1596. [Google Scholar] [CrossRef]
  39. Amat, S.; Dadourian, K.; Liandrat, J.; Ruiz, J.; Trillo, J.C. On a class of L1-stable nonlinear cell-average multiresolution schemes. J. Comput. Appl. Math. 2010, 234, 1129–1139. [Google Scholar] [CrossRef] [Green Version]
  40. Dadourian, K. Schémas de Subdivision, Analyses Multirésolutions Non-Linéaires. Applications. Ph.D. Thesis, University of Provence—Aix-Marseille I, Aix-en-Provence, France, 2008. [Google Scholar]
  41. Amat, S.; Busquier, S.; Legaz, M.; Manzano, F.; Ruiz, J. Reciprocal polynomial extrapolation vs. richardson extrapolation for singular perturbed boundary problems. Numer. Algorithms 2012, 61, 631–647. [Google Scholar] [CrossRef]
  42. Kuijt, F. Convexity Preserving Interpolation: Nonlinear Subdivision and Splines. Ph.D. Thesis, University of Twente, Enschede, The Netherlands, 1998. [Google Scholar]
Figure 1. The results obtained for the function in (23) using 10 scales of subdivision. The top left presents the result obtained by the nonlinear version of the three points binary algorithm in (34). The plot at the top right shows a zoom around the discontinuity. At the bottom of the 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 reduce the diffusion introduced by the linear algorithm.
Figure 1. The results obtained for the function in (23) using 10 scales of subdivision. The top left presents the result obtained by the nonlinear version of the three points binary algorithm in (34). The plot at the top right shows a zoom around the discontinuity. At the bottom of the 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 reduce the diffusion introduced by the linear algorithm.
Mathematics 10 02790 g001
Figure 2. The results obtained for the function in (24) using 10 scales of subdivision. The top left presents the result obtained by the nonlinear version of the three points binary algorithm in (34). The plot at the top right shows a zoom around the discontinuity. At the bottom of the 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 reduce the diffusion that the linear algorithm introduces.
Figure 2. The results obtained for the function in (24) using 10 scales of subdivision. The top left presents the result obtained by the nonlinear version of the three points binary algorithm in (34). The plot at the top right shows a zoom around the discontinuity. At the bottom of the 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 reduce the diffusion that the linear algorithm introduces.
Mathematics 10 02790 g002
Figure 3. The results obtained for the function in (25) using 10 scales of subdivision. The top left presents the result obtained by the nonlinear version of the three points binary algorithm in (34). The plot at the top right shows a zoom around the discontinuity. At the bottom of the figure, we can observe the results obtained by the linear algorithm in (2). The nonlinear algorithm manages to reproduce piecewise constant functions, while the linear algorithm introduces diffusion.
Figure 3. The results obtained for the function in (25) using 10 scales of subdivision. The top left presents the result obtained by the nonlinear version of the three points binary algorithm in (34). The plot at the top right shows a zoom around the discontinuity. At the bottom of the figure, we can observe the results obtained by the linear algorithm in (2). The nonlinear algorithm manages to reproduce piecewise constant functions, while the linear algorithm introduces diffusion.
Mathematics 10 02790 g003
Figure 4. Plot of the low resolution data obtained from the function in (26).
Figure 4. Plot of the low resolution data obtained from the function in (26).
Mathematics 10 02790 g004
Figure 5. Top, subdivided data using w = 0.5 and one scale of subdivision. Bottom, subdivided data using w = 1 2 and one scale of subdivision. To the left, linear algorithm and to the right non linear algorithm.
Figure 5. Top, subdivided data using w = 0.5 and one scale of subdivision. Bottom, subdivided data using w = 1 2 and one scale of subdivision. To the left, linear algorithm and to the right non linear algorithm.
Mathematics 10 02790 g005
Table 1. Numerical estimation of the limit functions regularity C 1 + β 1 and C 2 + β 2 for S N L in (34) using the arithmetic mean and the harmonic mean. The number of subdivision scales is represented by j.
Table 1. Numerical estimation of the limit functions regularity C 1 + β 1 and C 2 + β 2 for S N L in (34) using the arithmetic mean and the harmonic mean. The number of subdivision scales is represented by j.
j5678910
β 1 Arithmetic9.997 × 10 01 1.000 × 10 + 00 1.000 × 10 + 00 1.000 × 10 e + 00 1.000 × 10 + 00 1.000 × 10 + 00
Harmonic7.648 × 10 01 9.450 × 10 01 9.773 × 10 01 9.987 × 10 01 9.999 × 10 01 1.000 × 10 + 00
β 2 Arithmetic4.150 × 10 01 4.150 × 10 01 4.150 × 10 01 4.150 × 10 01 4.150 × 10 01 4.150 × 10 01
Harmonic−2.430 × 10 01 −1.566 × 10 02 −3.817 × 10 03 −1.161 × 10 03 −6.646 × 10 05 −3.798 × 10 06
Table 2. Grid refinement analysis in the L norm for the numerical estimation of the order of approximation using j scales of subdivision. By i, we represent the number of initial control points that we use. We used the nonlinear scheme (with M equal to the harmonic mean in (34)) and the function in (24). For this piecewise continuous function, the greatest error is located at the discontinuity. Thus, the numerical errors presented in the table correspond to the errors at the discontinuity. We can see that they reduce as we refine the grid, reaching O ( h ) accuracy.
Table 2. Grid refinement analysis in the L norm for the numerical estimation of the order of approximation using j scales of subdivision. By i, we represent the number of initial control points that we use. We used the nonlinear scheme (with M equal to the harmonic mean in (34)) and the function in (24). For this piecewise continuous function, the greatest error is located at the discontinuity. Thus, the numerical errors presented in the table correspond to the errors at the discontinuity. We can see that they reduce as we refine the grid, reaching O ( h ) accuracy.
j123456
i e i log 2 e i e i + 1 e i log 2 e i e i + 1 e i log 2 e i e i + 1 e i log 2 e i e i + 1 e i log 2 e i e i + 1 e i log 2 e i e i + 1
640.21777-0.34569-0.4637-0.57448-0.6794-0.77926-
1280.10871.00250.174460.986560.236530.971190.296140.955960.353930.94080.410220.9257
2560.0542991.00140.0876190.993590.119410.986070.150280.978590.180540.971120.210340.96366
5120.0271361.00070.0439050.996870.059990.993150.0756940.989450.0911680.985740.106490.98204
10240.0135651.00040.0219760.998450.0300660.99660.0379850.994760.0458080.992920.0535750.99107
Table 3. Grid refinement analysis in the L norm for the numerical estimation of the order of approximation using j scales of subdivision. By i, we represent the number of initial control points that we use. We used the linear scheme (with M equal to the arithmetic mean in (34)) and the function in (24). For this piecewise continuous function, the greatest error is located at the discontinuity. Thus, the numerical errors presented in the table correspond to the errors at the discontinuity. We can see that they do not reduce as we refine the grid, reaching O ( 1 ) accuracy.
Table 3. Grid refinement analysis in the L norm for the numerical estimation of the order of approximation using j scales of subdivision. By i, we represent the number of initial control points that we use. We used the linear scheme (with M equal to the arithmetic mean in (34)) and the function in (24). For this piecewise continuous function, the greatest error is located at the discontinuity. Thus, the numerical errors presented in the table correspond to the errors at the discontinuity. We can see that they do not reduce as we refine the grid, reaching O ( 1 ) accuracy.
j123456
i e i log 2 e i e i + 1 e i log 2 e i e i + 1 e i log 2 e i e i + 1 e i log 2 e i e i + 1 e i log 2 e i e i + 1 e i log 2 e i e i + 1
641.0834-1.5284-1.7328-1.8316-1.8803-1.9045-
1281.07330.0134691.51370.0139871.71580.0142221.81350.0143291.86160.0143791.88560.014403
2561.06840.00667451.50640.00695061.70740.00707741.80450.00713581.85240.00716321.87630.0071764
5121.06590.00332211.50280.00346441.70320.00353011.80010.00356061.84780.00357491.87160.0035818
10241.06470.00165721.5010.00172951.70110.00176291.79780.00177841.84560.00178581.86930.0017893
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zouaoui, S.; Amat, S.; Busquier, S.; Ruiz, J. On a Nonlinear Three-Point Subdivision Scheme Reproducing Piecewise Constant Functions. Mathematics 2022, 10, 2790. https://doi.org/10.3390/math10152790

AMA Style

Zouaoui S, Amat S, Busquier S, Ruiz J. On a Nonlinear Three-Point Subdivision Scheme Reproducing Piecewise Constant Functions. Mathematics. 2022; 10(15):2790. https://doi.org/10.3390/math10152790

Chicago/Turabian Style

Zouaoui, Sofiane, Sergio Amat, Sonia Busquier, and Juan Ruiz. 2022. "On a Nonlinear Three-Point Subdivision Scheme Reproducing Piecewise Constant Functions" Mathematics 10, no. 15: 2790. https://doi.org/10.3390/math10152790

APA Style

Zouaoui, S., Amat, S., Busquier, S., & Ruiz, J. (2022). On a Nonlinear Three-Point Subdivision Scheme Reproducing Piecewise Constant Functions. Mathematics, 10(15), 2790. https://doi.org/10.3390/math10152790

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop