![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
The subdivision scheme is a valuable tool for designing shapes and representing geometry in computer-aided geometric design. It has excellent geometric properties, such as fractals and adjustable shape. In this research paper, we explore the generation of fractal curves using a novel binary 8-point interpolatory subdivision scheme with two parameters. We analyse different properties of the proposed scheme, including convergence, special cases, and fractals. Additionally, we demonstrate through various examples the relationship between the shape parameters and the fractal behaviour of the resulting curve. Our research also identifies a specific range of shape parameters that can effectively produce fractal curves. The findings of this study provide a fast and efficient method for generating fractals, as demonstrated by numerous examples. Modelling examples show that the 8-point interpolatory scheme can enhance the efficiency of computer design for complex models.
1. Introduction
Subdivision is a well-organized technique for producing curves and surfaces in Computer Aided Geometric Design (CAGD). The subdivision technique is mathematical by its nature thus very simple to implement on computer systems and consequently has a wide range of applications in various industries. Some examples are animation i.e film making, aircrafts i.e ship making, utensils making, furniture industry, machinery production, the vehicles industry and jewelry making. In general, subdivision schemes can be divided into two categories: Interpolatory schemes and Approximating schemes. Interpolatory schemes get better shape control while approximating schemes have better smoothness. The interpolating subdivisions schemes are more gorgeous then approximating schemes because of their interpolating property.
The history of subdivision schemes started when the first subdivision scheme was presented by Rham [Citation1]. This scheme was a recursive corner-cutting piecewise linear approximation to generate - continuous smooth curves, which opened new dimensions in the field of geometric modelling. Ko et al. [Citation2] reconstructed the mask of interpolating symmetric subdivision systems,
-point binary interpolating systems and the 4-point ternary interpolating scheme by using the essential condition for smoothness and symmetry. Rehan and Sabri [Citation3] introduced a combined ternary four-point subdivision scheme. Laurent polynomial method is used to evaluate the continuity and observed that this scheme generates interpolating
and approximating
,
,
limiting curves. Ghaffar et al [Citation4] introduced a new non-stationary
-point subdivision scheme that satisfies some of the mentioned properties are curvature, torsion monotonicity shape preservation and continuity. Ghaffar et al [Citation5] introduced the non-stationary binary
-point subdivision schemes for integer
. The Lagrange polynomial method is used to establish the schemes that reproduce polynomial degree
. The importance of these schemes is proved with the help of different examples.
In last few decades, computational geometry by using SSs have attained much attention. Mustafa et al. [Citation6–8] devised parametric schemes using a new type of Lagrange polynomial that employs a single shape parameter. Ghaffar et al. [Citation9–12] created continuity conditions to represent complex, free-form unified Shape Splines (SSs) with one shape control parameter. These methods allow for the creation of intricate shapes through the use of shape control unified SSs. These approaches inherit the strengths of non-parametric SSs and resolve the issue of shape control in SSs by employing tension control shape parameters. The proposed methods use varied techniques to construct composite curves with higher smoothness conditions, using simple and straightforward computations with linear polynomials instead of trigonometric functions. In 2020, Ashraf et al. [Citation13–15] introduced a new method for generalizing SSs with a single shape parameter. Shahzad et al. [Citation16] developed a numerical algorithm to calculate the subdivision depth of binary SSs with shape parameters.
Hussain et al. [Citation17] proposed a subdivision algorithm with varying arity that explores properties such as continuity, regularity, error bound, and limit curve shape. Ashraf et al. [Citation18, Citation19] presented a generalized algorithm to develop a class of approximating binary subdivision schemes. Mustafa et al. [Citation20] introduced a numerical approach for solving second-order singularly perturbed boundary value problems. Bari et al. [Citation21] unified the interpolating nonstationary subdivision scheme with a tension control parameter to construct a link between nonstationary SSs and applied geometry.
Fractals are the patterns that are infinitely complex patterns. These are made up by repeatedly using a basic operation in an ongoing feedback iterative process. Nature is filled with hundreds of fractals varieties. Trees and their leaves, mountains ranges, systems of rivers, snowflakes, crystals, sea-shells, sea-stars, clouds, lighting, fruits like pineapple, vegetables like cauliflower or broccoli, blood vessels, and many other examples like these exhibit formation of fractals. The self-similarity of fractals makes it an attractive factor for creative work. From paintings, sculptures, textiles to graphic designing, animations and in every domain of artistic work is affected by fractals.
Doo and Sabin [Citation22] discussed the functioning of the limits surface using eigenvalues of a set of matrices, defined by a recursive division construction. Carpenter [Citation23], generated fractal curves and surfaces by using a simple set of methods and presented the idea of the statistical subdivision. Mandelbrot [Citation24] shed light on fractals that exists in nature's beauty and devised fractal geometry that studied self-similar and irregular geometric shapes and has an in-depth presentation at arbitrarily small scales. Soon after that, several procedures had started to formulate to create fractals. Bouboulis and Dalla [Citation25] derived the fractal interpolation surfaces from its respective functions and proved some interesting properties of them. Zheng et al. [Citation26] used a new mechanism to evaluates the fractal characteristics of binary four point and ternary three point interpolating subdivision schemes. Later on a relationship is found between shape parameter and fractal behaviour of the boundary curves of the both schemes respectively. The production of fractal curves and surfaces are also discussed. The three point ternary interpolating scheme with two shape parameters that was devised by Hassan is examined by Zheng et al. [Citation27].
Many examples are elaborated by author that demonstrates the true direction for creation of fractal curves and surfaces. By using this method fractal range of this scheme is also obtained. Mustafa et al. [Citation28] discussed the engineering images by the fractal proprieties of binary 6-point interpolating scheme. A relationship between the fractal dimension of the fractal limit curve and shape parameter is also found and show that this scheme is best for the modelling of fractal antennas, garari's, bearing and rock etc. Peng et al. [Citation29] have conducted a study on the fractal properties of a ternary 4-point rational interpolation SS. Meanwhile, Bari et al. [Citation30] have proposed a method for constructing fractal curves using a 7-point binary approximating SS with one shape parameter. On the other hand, Yao et al. [Citation31] have discussed the fractal and convexity analysis of the quaternary 4-point SS.
On the basis of the above mention literature of subdivision schemes, in this paper we are going to construct ‘a new 8- point interpolating subdivision scheme’ which is the generalized form of schemes presented by Weiseman [Citation32] and Deslauriers and Dubuc [Citation33] and the fractal behaviour of this proposed subdivision scheme is investigated and analysed.
2. Preliminaries
Here, we present some definitions and theorems which will be used in coming section.
Definition 2.1
Binary Subdivision Scheme
A general form of univariate binary subdivision scheme S which maps polygon to a refined polygon
is defined by
(1)
(1) Where the set
of coefficients is called mask/stencil of subdivision scheme.
Definition 2.2
Laurent Polynomial
The Laurent polynomial of the mask of the scheme (Equation1
(1)
(1) ) is defined as
Theorem 2.1
Suppose S be a convergent Subdivision Scheme with mask then, the mask should satisfy:
(2)
(2)
Theorem 2.2
Suppose S indicate a binary subdivision scheme with mask and the i-th order divided difference
with mask
satisfies:
(3)
(3) If there exists a smallest positive integer L satisfying
, then the binary subdivision scheme is
continuous. In particular, when:
3. Construction
The introducing rule of 8-point interpolating subdivision scheme is denoted as
(4)
(4) Laurent Polynomial of the above 8- point interpolating technique is
by using symmetric condition, we have
now the above Laurent Polynomial of this technique can be written as
by using necessary condition
, we have
then
(5)
(5) Laurent polynomial of
is
It is easy to note that
is always true.Laurent polynomial of
is
by using the necessary condition for
, that is
(6)
(6) Now putting
and
into equations (Equation6
(6)
(6) ), then we get
Now putting the values of
and
in equation (Equation5
(5)
(5) ), then
so the required new 8-pint interpolating subdivision scheme is
(7)
(7)
4. Convergence
This section presents the condition of convergence and smoothness of limit curves generated by scheme (Equation7(7)
(7) ) and discussed the special cases of this scheme. The mask of suggested 8-point interpolating subdivision scheme that generates polynomial
can be written as:
To prove
, it is required that
must satisfied the (Equation2
(2)
(2) ), which it does and
, we have
Now
This 8-point subdivion scheme is uniformly convergent and generates
limit curves. To prove
, we require that
satisfy the (Equation2
(2)
(2) ), which it does and
, we have
Now
This 8-point subdivion scheme generates
limit curves. To prove
, we require that
satisfy the (Equation2
(2)
(2) ), which it does and
, we have
Now
This 8-point subdivion scheme generates
limit curves.
4.1. Special cases
It can be seen that there are certain traditional subdivision techniques that are special cases of said scheme here.
by letting
and
, we get 4-point interpolating subdivision scheme.
this is called Deslauries and Dubuc [Citation33] 4-point interpolating scheme and this scheme is
continuous.
by letting
and
, we get 6-point interpolation subdivision scheme.
this is Deslauries and Dubuc [Citation33] 6-point interpolating scheme and this scheme produces
limit curves.
by letting
, we get 6-point interpolation subdivision scheme.
this is called Weissman [Citation32] 6-point interpolating sundivision scheme, that generates
and
continues curves at
and
respectively.
by letting
and
, we get 8-point interpolation subdivision scheme.
this is Deslauriers and Dubuc [Citation33] 8-point interpolating scheme is
continues.
5. Generation of fractal curves
The 8-point subdivion scheme (Equation7(7)
(7) ) is uniformly convergent and generates
limit curves. let two choosen fixed control points
and
after k subdivision iterations, where ∀
and
. For next subdivision iteration the effect of the parameters must be examined on the sum of all the smaller edges between two points. The effect between two initial control points
and
is investigated for simplicity's sake. As per the interpolating property
,
. we have
and
.since
so we have
so we can write as
(8)
(8) for the solution of the edge vectors
and
, we use
,
,we can write it as
further
The vector
can be written as
(9)
(9) Putting i = 0 in equation (Equation7
(7)
(7) )
(10)
(10) Putting i = −1 in equation (Equation7
(7)
(7) )
(11)
(11) from Equations (Equation10
(10)
(10) ) and (Equation11
(11)
(11) )
using the dyadic parametrization, we can write the terms as
and
now using the even rule we get
and
.Therefore
here we have, linear differential equation of second order
(12)
(12) The characteristic equation is
(13)
(13) by solving equation (Equation13
(13)
(13) ), we get
and
,
when
and
The solution of equation (Equation12
(12)
(12) ) is
(14)
(14) where
and
When
and
the solution of the equation (Equation14
(14)
(14) ) will be
and
When
and
then solution of equation (Equation14
(14)
(14) ) is
where
and
since
Since by dyadric parametrization
,
and
then
Then by taking
,
and
, we get
we know
⇒
, then
⇒
(15)
(15) Now we will find the solution of equation (Equation15
(15)
(15) ). The characteristic equation is
The solution of the above equation is
and
and
,
,
so equation (Equation15
(15)
(15) ) can be write as
(16)
(16) The solution of the equation (Equation16
(16)
(16) ) is
(17)
(17) where
and
are given as
From equation (Equation8
(8)
(8) ) and equation (Equation17
(17)
(17) ), we have
(18)
(18)
Theorem 5.1
For and
, the limit curve of the 8-point interpolating scheme is a fractal curve.
Proof.
From Equations (Equation17(17)
(17) ) and (Equation18
(18)
(18) ), it can be concluded by mathematical induction that
small edge vectors between
and
, after k subdivision level can be shown as:
(19)
(19) for
, where
, i = 1, 2, 3, 4. the equation (Equation19
(19)
(19) ) is the linear combination of eqs (Equation17
(17)
(17) ) and (Equation18
(18)
(18) ). By taking
and
. Let
denote the length of a vector v and
, then we obtained
tends to ∞ as k tends to ∞ Hence, when k approaches to infinity, the sum of the length of all small edges between initial control points
and
has no bound, after the k subdivision levels, so the limit curve of 8-point interpolating scheme is a fractal curve when
and
.
6. Fractal curves
Introducing two shape parameters in the generation of fractal curves using a new binary 8-point interpolatory subdivision scheme can offer several advantages. These parameters provide additional degrees of freedom, allowing for more flexibility and control over the characteristics of the generated fractals. Here are some potential advantages:
Two shape parameters provide greater control over the interpolation process. Each parameter can be independently adjusted, allowing for more fine-tuned control over specific aspects of the fractal curve, such as detail, smoothness, or overall shape.
With two shape parameters, it becomes possible to control both local and global features of the fractal curves independently. One parameter may influence local details, while the other could govern broader patterns, giving the user a more nuanced control over the aesthetics.
The dual parameter system allows for the customization of the fractal curves to achieve specific visual effects. For instance, one parameter could control the density of fine details, while the other might influence the overall structure or coarseness of the curve.
Different applications or artistic requirements may demand various characteristics in fractal curves. The two shape parameters make the subdivision scheme adaptable to a broader range of scenarios, providing a versatile tool for generating diverse fractal patterns.
The fractal dimension of a curve can be influenced by the shape parameters. Having two parameters allows for a more precise adjustment of the fractal dimension, enabling the generation of curves with specific dimensionality characteristics.
A dual parameter system enhances aesthetic control, enabling users to achieve visually appealing fractal curves with a higher level of precision. This is particularly valuable in artistic and design applications where the appearance of the fractals is a key consideration.
The use of two shape parameters supports dynamic and interactive design processes. Users can iteratively adjust the parameters and observe real-time changes in the fractal curves, facilitating an explorative and intuitive design experience.
The presence of two shape parameters encourages creativity and exploration in fractal curve design. Users can experiment with various combinations of parameter values, leading to the discovery of novel and interesting fractal patterns.
Complex fractal patterns often require nuanced control over multiple features. Two shape parameters provide a mechanism to navigate the complexity and intricacy of the fractal curve space, facilitating the generation of more sophisticated patterns.
Users can balance trade-offs between conflicting design goals by adjusting the two parameters. For example, one parameter might control detail while the other manages smoothness, allowing for a balance that suits specific preferences or requirements.
Here some examples of fractal curves are given by using the new binary 8-point interpolating subdivision scheme.
Figure shows three fractal curves are produced by applying the 8-point binary interpolating subdivision scheme. Figure (i) represents the initial control polygon whereas Figure (ii-iv) shows the results obtained by setting the values of and
after first, fourth and tenth iteration respectively.
Figure 1. Fractals:(i) Elaborate the initial polygon and (ii-iv) elaborates the fractal curves at first, forth and tenth subdivision level of the scheme (Equation7(7)
(7) ) at
respectively.
![Figure 1. Fractals:(i) Elaborate the initial polygon and (ii-iv) elaborates the fractal curves at first, forth and tenth subdivision level of the scheme (Equation7(7) f2ik+1=fik,f2i+1k+1=−ηfi−3k+(μ+5η)fi−2k−(3μ+9η+116)fi−1k+(2μ+5η+916)fik+(2μ+5η+916)fi+1k−(3μ+9η+116)fi+2k+(μ+5η)fi+3k−ηfi+4k.(7) ) at η=−0.0328,μ=−0.0328 respectively.](/cms/asset/ce07eb40-02de-4bfd-93bc-20a9b5cdbd05/gipe_a_2367718_f0001_oc.jpg)
Figure explains that the two types fractal curves generated by implementing the 8-point interpolating binary subdivision scheme after first, third, fifth, seventh, 9th, 11th and 13th subdivision level of the scheme along with initial polygon by using (Equation7(7)
(7) ) at
and
respectively. It is observed that the nature of fractal curves are totally different on the positive and negative range of parameters.
Figure 2. Fractals:(i) shows the initial polygon and (ii-xvi) shows the fractal curves at first, third, fifth, seventh, 9th, 11th and 13th subdivision level of the scheme (Equation7(7)
(7) ) at
respectively.
![Figure 2. Fractals:(i) shows the initial polygon and (ii-xvi) shows the fractal curves at first, third, fifth, seventh, 9th, 11th and 13th subdivision level of the scheme (Equation7(7) f2ik+1=fik,f2i+1k+1=−ηfi−3k+(μ+5η)fi−2k−(3μ+9η+116)fi−1k+(2μ+5η+916)fik+(2μ+5η+916)fi+1k−(3μ+9η+116)fi+2k+(μ+5η)fi+3k−ηfi+4k.(7) ) at η=−0.0576,μ=0.0576 respectively.](/cms/asset/9e94768e-5f01-45f5-9966-28637f586834/gipe_a_2367718_f0002_oc.jpg)
Figure shows three fractal curves generated by using 8-point interpolating scheme after first and tenth level respectively. Figure (i) depicts their initial polygon, Figure (ii) shows the result obtained with and
after first iteration and Figure (iii) shows the result obtained with
and
after tenth iteration.
Figure 3. (i) shows the initial polygon and Figure (ii-iii) shows the fractal curves at first and tenth subdivision level of the scheme (Equation7(7)
(7) ) at
respectively.
![Figure 3. (i) shows the initial polygon and Figure 4(ii-iii) shows the fractal curves at first and tenth subdivision level of the scheme (Equation7(7) f2ik+1=fik,f2i+1k+1=−ηfi−3k+(μ+5η)fi−2k−(3μ+9η+116)fi−1k+(2μ+5η+916)fik+(2μ+5η+916)fi+1k−(3μ+9η+116)fi+2k+(μ+5η)fi+3k−ηfi+4k.(7) ) at η=0.029,μ=0.029 respectively.](/cms/asset/30e49c67-aac5-4bc8-8453-01e737c2e5f9/gipe_a_2367718_f0003_oc.jpg)
Figure 4. Fractals:(i) shows the initial polygon and (ii-iv) shows the fractal curves at first, forth and tenth subdivision level of the scheme (Equation7(7)
(7) ) at
respectively.
![Figure 4. Fractals:(i) shows the initial polygon and (ii-iv) shows the fractal curves at first, forth and tenth subdivision level of the scheme (Equation7(7) f2ik+1=fik,f2i+1k+1=−ηfi−3k+(μ+5η)fi−2k−(3μ+9η+116)fi−1k+(2μ+5η+916)fik+(2μ+5η+916)fi+1k−(3μ+9η+116)fi+2k+(μ+5η)fi+3k−ηfi+4k.(7) ) at η=−0.029,μ=−0.029 respectively.](/cms/asset/3abfc7cb-1acf-4bd1-a4d2-059675020313/gipe_a_2367718_f0004_oc.jpg)
Figure (i) shows the initial polygon for a bearing and Figure (ii-iv) shows the actual bearing generated by setting the values of after first, fifth and tenth level respectively. The fractal curves in Figure (ii-iii) are just like the fractal antenna of Weissman [Citation32] binary 6-point interpolating subdivision scheme that is generated with one parameter at
.
Figure 5. Terrain like fractal curves by applying the 8-point binary interpolating subdivision scheme. (i) shows the initial polygon and (ii-iii) shows the fractal curves at fifth and ninth subdivision level of the scheme (Equation7(7)
(7) ) at
respectively.
![Figure 5. Terrain like fractal curves by applying the 8-point binary interpolating subdivision scheme. (i) shows the initial polygon and (ii-iii) shows the fractal curves at fifth and ninth subdivision level of the scheme (Equation7(7) f2ik+1=fik,f2i+1k+1=−ηfi−3k+(μ+5η)fi−2k−(3μ+9η+116)fi−1k+(2μ+5η+916)fik+(2μ+5η+916)fi+1k−(3μ+9η+116)fi+2k+(μ+5η)fi+3k−ηfi+4k.(7) ) at η=−0.03,μ=−0.03 respectively.](/cms/asset/197c55e9-8f36-4405-ac3f-3d118d8a6a93/gipe_a_2367718_f0005_ob.jpg)
Figure shows the terrain like fractal curves produced by applying the proposed 8-point interpolating binary scheme after fifth and ninth iteration. Figure 6(i) depicts the initial control polygon. Figure 6(ii-iii) which shows the results obtained by letting the values of respectively. From the above examples we can show that given the control polygon we can generate fractal curves by choosing the appropriate values of parameters according to theorem (5.1).
7. Conclusion
In this paper, A new 8-point binary interpolating subdivision scheme is constructed by using two parameters η and μ. Laurent polynomial is used to analysed the continuity of scheme. Fractal properties of this scheme with two parameters are discussed and relationship between the parameters and the limit fractal curve of 8-point binary interpolating subdivision scheme is also presented. The advantages of incorporating two shape parameters in a binary 8-point interpolatory subdivision scheme are rooted in the enhanced control, adaptability, and creative possibilities they offer for generating diverse and aesthetically pleasing fractal curves.
Disclosure statement
No potential conflict of interest was reported by the author(s).
References
- Rham Gde. Unpeu de mathématiques à propos d'unecourbe plane. Elemente Der Math. 1947;2:73–76.
- Ko KP, Lee BG, Yoon GJ. A study on the mask of interpolatory symmetric subdivision schemes. Appl Math Comput. 2007;187(2):609–621. doi: 10.1016/j.amc.2006.08.089
- Rehan K, Sabri MA. A combined ternary 4-point subdivision scheme. Appl Math Comput. 2016;276:278–283. doi: 10.1016/j.amc.2015.12.016
- Ghaffar A, Ullah Z, Bari M, et al. A new class of 2 m-point binary non-stationary subdivision schemes. Adv Differ Equ. 2019;2019(1):1–19. doi: 10.1186/s13662-019-2264-4
- Ghaffar A, Ullah Z, Bari M, et al. Family of odd point non-stationary subdivision schemes and their applications. Adv Differ Equ. 2019;2019(1):1–20. doi: 10.1186/s13662-019-2105-5
- Mustafa G, Khan F, Ghaffar A. The m-point approximating subdivision scheme. Lobachevskii J Math. 2009;30(2):138–145. doi: 10.1134/S1995080209020061
- Mustafa G, Ghaffar A, Khan F. The odd-point ternary approximating schemes. Am J Comput Math. 2011;01(02):111.118 doi: 10.4236/ajcm.2011.12011
- Mustafa G, Ghaffar A, Aslam M. A subdivision-regularization framework for preventing over fitting of data by a model. Appl Appl Math Int J. 2013;8(1):178–190.
- Ghaffar A, Mustafa G, Qin K. Unification and application of 3-point approximating subdivision schemes of varying arity. Open J Appl Sci. 2012;02(04):48–52. doi: 10.4236/ojapps.2012.24B012
- Ghaffar A, Mustafa G, Qin K. The 4-point a-ary approximating subdivision scheme. Open J Appl Sci. 2013;03(01):106–111. doi: 10.4236/ojapps.2013.31B1022
- Ghaffar A, Bari M, Ullah Z, et al. A new class of 2 q-point nonstationary subdivision schemes and their applications. Mathematics. 2019;7(7):639. doi: 10.3390/math7070639
- Ghaffar A, Iqbal M, Bari M, et al. Construction and application of nine-tic B-spline tensor product SS. Mathematics. 2019;7(8):675. doi: 10.3390/math7080675
- Ashraf P, Nawaz B, Baleanu D, et al. Analysis of geometric properties of ternary four-point rational interpolating subdivision scheme. Mathematics. 2020;8(3):338. doi: 10.3390/math8030338
- Ashraf P, Ghaffar A, Baleanu D, et al. Shape-Preserving properties of a relaxed four-Point interpolating subdivision scheme. Mathematics. 2020;8(5):806. doi: 10.3390/math8050806
- Ashraf P, Sabir M, Ghaffar A, et al. Shape-preservation of the four-point ternary interpolating non-stationary subdivision scheme. Front Phys. 2020;7:241. doi: 10.3389/fphy.2019.00241
- Shahzad A, Khan F, Ghaffar A, et al. A novel numerical algorithm to estimate the subdivision depth of binary subdivision schemes. Symmetry. 2020;12(1):66. doi: 10.3390/sym12010066
- Hussain SM, Rehman AU, Baleanu D, et al. Generalized 5-point approximating subdivision scheme of varying arity. Mathematics. 2020;8(4):474. doi: 10.3390/math8040474
- Ashraf P, Mustafa G, Ghaffar A, et al. Unified framework of approximating and interpolatory subdivision schemes for construction of class of binary subdivision schemes. J Funct Spaces. 2020;2020. 1 19. doi: 10.1155/2020/6677778
- Ashraf P, Mustafa G, Khan HA, et al. A shape-preserving variant of Lane-Riesenfeld algorithm. AIMS Math. 2021;6(3):2152–2170. doi: 10.3934/math.2021131
- Mustafa G, Ejaz ST, Baleanu D, et al. A subdivision-based approach for singularly perturbed boundary value problem. Adv Differ Equ. 2020;2020(1):1–20. doi: 10.1186/s13662-020-02732-8
- Bari M, Mustafa G, Ghaffar A, et al. Construction and analysis of unified 4-point interpolating nonstationary subdivision surfaces. Adv Differ Equ. 2021;2021(1):1–17. doi: 10.1186/s13662-021-03234-x
- Doo D, Sabin M. Behaviour of recursive division surfaces near extraordinary points. Comput-Aided Des. 1978;10(6):356–360. doi: 10.1016/0010-4485(78)90111-2
- Carpenter LC. Computer rendering of fractal curves and surfaces. ACM SIGGRAPH Comput Graphics. 1980;14(SI):9–9. doi: 10.1145/988447.988449
- Mandelbrot BB, Mandelbrot BB. The fractal geometry of nature. (Vol. New York: WH freeman; 1982. doi: 10.1002/esp.3290080415
- Bouboulis P, Dalla L. Fractal interpolation surfaces derived from fractal interpolation functions. J Math Anal Appl. 2007;336(2):919–936. doi: 10.1016/j.jmaa.2007.01.112
- Zheng H, Ye Z, Lei Y, et al. Fractal properties of interpolatory subdivision schemes and their application in fractal generation. Chaos Solit Fractals. 2007;32(1):113–123. doi: 10.1016/j.chaos.2005.10.075
- Zheng H, Ye Z, Chen Z, et al. Fractal range of a 3-point ternary interpolatory subdivision scheme with two parameters. Chaos Solit Fractals. 2007;32(5):1838–1845. doi: 10.1016/j.chaos.2005.12.017
- Mustafa G, Bari M, Jamil S. Engineering images designed by fractal subdivision scheme. SpringerPlus. 2016;5(1):1–11. doi: 10.1186/s40064-015-1659-2
- Peng K, Tan J, Li Z, et al. Fractal behavior of a ternary 4point rational interpolation subdivision scheme. Math Comput Appl. 2018;23(4):65. doi: 10.3390/mca23040065
- Bari M, Mustafa G, Hussain A. Fractal generation by 7-Point binary approximating subdivision scheme. Int J Appl Phys Math. 2018;8(4):105–111. doi: 10.17706/ijapm.2018.8.4.105-111
- Yao SW, Ashraf P, Ghaffar A, et al. Fractal and convexity analysis of the quaternary 4-point scheme and its applications. Fractals. 2023;31(10):1–14. doi: 10.1142/S0218348X23400881
- Weissman A. A 6-point interpolatory subdivision scheme for curve design. University of Tel-Aviv. 1989.
- Deslauriers G, Dubuc S. Symmetric iterative interpolation processes. In Constructive approximation (pp. 49–68). Springer, Boston, MA. 1989.