Next Article in Journal
An Efficient Class of Traub-Steffensen-Like Seventh Order Multiple-Root Solvers with Applications
Previous Article in Journal
Leaf Counting with Multi-Scale Convolutional Neural Network Features and Fisher Vector Coding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Modeling of a Class of Symmetrical Islamic Design

College of Sciences, Department of Mathematics, Research Group MASEP, University of Sharjah, P.O. Box 27272, Sharjah 61174, UAE
Symmetry 2019, 11(4), 517; https://doi.org/10.3390/sym11040517
Submission received: 19 March 2019 / Revised: 3 April 2019 / Accepted: 8 April 2019 / Published: 10 April 2019

Abstract

:
In this paper, we present a new model for simulating an interesting class of Islamic design. Based on periodic sequences on the one-dimensional manifolds, and from emerging numbers, we construct closed graphs with edges on the unit circle. These graphs build very nice shapes and lead to a symmetrical class of geometric patterns of so-called Islamic design. Moreover, we mathematically characterize and analyze some convergence properties of the used up-down sequences. Finally, four planar type of patterns are simulated.

1. Motivation

The art in the Islamic culture is very rich in spiritual and meditative messages. It uses complex geometric patterns and shapes; see Figure 1. For instance, in wood for ornamental patterns, in bricks of buildings, in brass for decorations, paper, tiling, plaster, glass, etc. Islamic design can be found mainly in three distinct geometrical patterns, namely rectangular Kufic fonts, arabesque shapes, and polygonal designs [1]. The largest class of Islamic design employs, in general, complex and symmetrical polygons, where, in each polygonal graph, the starting design is a cell in the form of a symmetrical shape. The creation of the symmetrical rectangular, triangular, and circular patterns leads to beautiful ornamental patterns [2,3,4,5]. The mathematical modeling of these cells is based on iterative techniques for generating beautiful and colored patterns. Different methods were used for seeking explicit relationships between a simple single cell and the corresponding complex multiplicity, either for constructing symmetrical groups for infinite periodic patterns form the 17 classes of the two-dimensional crystallographic groups [6] or by using planar geometrical transformations, translations, rotations, reflections, and glide reflections for rearranging and space filling. The modeling of symmetries in nature and sciences brought together the interest of many scientists [7,8,9,10] for developing several techniques. In [11,12], the authors used the logistic mapping for generating and describing pictures created by merging symmetry and chaos.
In [13], we studied the existence and the symmetry of the multiple barycenters, especially in the unit circle considered as one-dimensional manifold. The generated vertices are, for special values of n, the so-called multiple barycenters on S 1 . Moreover, since mapping circular patterns on curved surfaces in a three-dimensional space is in general a nontrivial task—see [14,15,16]—we propose a new algorithm for generating planar and isometric patterns. We construct and characterize vertices on the unit circle similar to the multiple barycenters using repeated iterations of emerged numbers. Thereafter, connecting these barycenters as vertices with finite number of edges describes and classifies special periodic sequences. These sequences are the main source for constructing symmetrical designs, especially Islamic design.
This paper is organized as follows: In Section 2, we introduce the periodic sequences on the unit circle, and we present necessary analytical properties. Based on the results presented in the previous section, we define and construct, in Section 3, the Islamic design. In Section 4, we present some numerical simulations of different types of designs and patterns. Conclusions and future works are listed in Section 5.

2. Periodic Sequences on the Unit Circle

Let a be a real number in [ 0 , 1 ] and S 1 be the unit circle considered as one-dimensional manifold. For a given starting value x a ( 0 ) = α , we define the real iterative sequence x a : I N S 1 with values on S 1 modulo 2 π by:
x a ( 0 ) = α , x a ( n ) : = x a ( n - 1 ) + a π if n odd , x a ( n - 1 ) - ( 1 - a ) π if n even .
The up-down behavior of the sequence x a ( n ) [ mod 2 π ] will be observed as emerging numbers on the unit circle, considered as a one-dimensional manifold (see Appendix A). In this example, we compute some special cases of the sequence (1) for the initial guess x a ( 0 ) = 0 . The following table presents the first eight values of the up-down sequence:
From our experiments, we observed that the sequence x a ( n ) is periodic for special values of the parameter a only. Therefore, we assign to the clockwise rotation a negative direction (negative integer k Z - ) and we assign the positive direction to any rotation in the counterclockwise direction (i.e.,  k Z + ). Furthermore, we are interested to characterizing the minimal number of rotations k Z till the first-come-back x a ( n ) = x a ( 0 ) . In the following, let us denote by X a ( I N ) : = { x a ( n ) | n I N } , the range of the up-down sequence x a for a given parameter a [ 0 , 1 ] .
Definition 1.
A sequence ( x a ) n 0 is called Alternating sequence if the corresponding range set X a ( I N ) contains exactly two elements. Otherwise, we call it a Loop-sequence.
Lemma 1.
The sequence x a is an Alternating sequence for a = 0 , 0 . 5 and 1 and a Loop-sequence for a [ 0 , 1 ] { 0 , 0 . 5 , 1 } .
Proof. 
From Formula (1), we have
x a ( 2 n + 1 ) - x a ( 2 n ) = a π ,
x a ( 2 n ) - x a ( 2 n - 1 ) = ( a - 1 ) π ,
x a ( 2 n + 1 ) - x a ( 2 n - 1 ) = ( 2 a - 1 ) π .
  • If a = 0 , then from the identities (2)–(4) we get
    x 0 ( 2 n + 1 ) = x 0 ( 2 n ) ,
    x 0 ( 2 n ) x 0 ( 2 n - 1 ) .
    Thus, the range of x 0 contains exactly two elements: X 0 ( I N ) = { x 0 ( 2 n - 1 ) , x 0 ( 2 n ) } , for n I N . Thus, it is an Alternating sequence.
  • If a = 1 , then from the identities (2)–(4) we get
    x 1 ( 2 n ) = x 1 ( 2 n - 1 ) ,
    x 1 ( 2 n ) x 1 ( 2 n + 1 ) .
    Hence, the range of x 1 contains exactly two elements: X 1 ( I N ) = { x 1 ( 2 n ) , x 1 ( 2 n + 1 ) } , for n I N . It follows that x 1 is an Alternating sequence.
  • If a = 0 . 5 , then from the identities (2)–(4) we get
    x 0 . 5 ( 2 n + 1 ) = x 0 . 5 ( 2 n - 1 ) ,
    x 0 . 5 ( 2 n ) x 0 . 5 ( 2 n - 1 ) .
    Thus, the range of x 0 . 5 contains exactly two elements: X 0 . 5 ( I N ) = { x 0 . 5 ( 2 n - 1 ) , x 0 . 5 ( 2 n ) } , for  n I N . Thus, it is an Alternating sequence.
  • If a ( 0 , 0 . 5 ) or a ( 0 . 5 , 1 ) , from the identities (2)–(4), the range of x a contains at least three distinct elements X a ( I N ) = { x a ( 2 n - 1 ) , x a ( 2 n ) , x a ( 2 n + 1 ) , } , for n I N . Thus, it is a Loop-sequence. □
In the bullet diagram below, we illustrate the results of the lemma above:
  • For a = 0 . 00 , we needed one rotation in the counterclockwise direction and four iteration to reach the first-come-back x 0 ( 4 ) = 0 for n = 4 and k = - 1 (see Table 1 and Figure 2a).
  • For a = 0 . 25 , we needed one rotation in the counterclockwise direction and eight iteration to reach the first-come-back x 0 . 25 ( 8 ) = 0 for n = 8 and k = - 1 (see Table 1 and Figure 2b).
  • For a = 0 . 50 , we consider this case as a critical case, while there is a come-back to the origin x 0 . 5 ( 0 ) = 0 without rotating a complete round. it is in fact an oscillatory behavior (see Table 1 and Figure 2c).
  • For a = 0 . 75 we needed one rotation in the positive clockwise direction and eight iteration to reach the first-come-back X 0 . 75 ( 8 ) = 0 for n = 8 and k = 1 (see Table 1 and Figure 2d).
  • For a = 1 . 00 we needed one rotation in the clockwise direction and four iteration to reach the first-come-back x 1 ( 4 ) = 0 for k = 1 (see Table 1 and Figure 2e).
  • For a = 0 . 95 we needed one rotation in the counterclockwise direction and four iteration to reach the first-come-back x 0 . 95 ( 40 ) = 0 for n = 40 and k = 9 (see Table 1 and Figure 2f).
It should be stressed that for a > 0 . 5 , the up-down sequence rotate in the counterclockwise direction (i.e., positive direction k > 0 ) and if a < 0 . 5 it rotates in the clockwise direction (negative direction k < 0 ). In the following, we are looking for a relationship between a and the minimal index such that x a ( n ) = 0 .
It is important to note that the graphs in the Figure 2b,d,f are closed and lead to the construction symmetric areas. While the path trajectory in the Figure 2a,c,e represents a segment. Moreover, in the Figure 2c, the up-down sequence oscillates between the origin and π / 2 .
According to the results in the example above, it should be stressed that for a certain value of a one can find the pair ( k , n ) satisfying the first-come-back property. In the same context, the following theorem characterize the parameter a I R using, in somehow, Archimedean principle applied to a circle. In other words, we determine a necessary and sufficient condition for constructing the minimal number of loops k Z till the first-come-back n N :
Theorem 1.
For any rational number a [ 0 , 1 ] Q and any initial value x a ( 0 ) , there exists n such that x a ( n ) x a ( 0 ) mod 2 π and if a [ 0 , 1 ] Q , we have x a ( n ) x a ( 0 ) for all n N .
Proof. 
Let us assume that x a ( 0 ) = 0 . If x a ( 0 ) 0 , then x ˜ 0 = x a - α = 0 . Therefore, we restrict the proof for to the case x a ( 0 ) = 0 and we distinguish between the following two cases:
Case 1: 
n even In this case, we have
x a ( n ) = x a ( 0 ) + n 2 a π + n 2 ( a - 1 ) π = x a ( 0 ) + n 2 ( 2 a - 1 ) π
Let us consider the number the integer k as the number of complete rotations. From (11) it follows that
2 k = n 2 ( 2 a - 1 ) ,
thus, the parameter a, could be written as
a = 1 2 4 k n + 1 .
Case 2: 
n odd means that n - 1 is even
x a ( n ) = x a ( n - 1 ) + a π = x a ( 0 ) + n - 1 2 a π + n - 1 2 ( a - 1 ) π + a π = x a ( 0 ) + n - 1 2 ( 2 a - 1 ) + a π = x a ( 0 ) + n 2 ( 2 a - 1 ) - 1 2 π
Similar to the previous case, from (11) we have
2 k = n 2 ( 2 a - 1 ) - 1 2 .
Thus, the parameter a, could be written as follows:
a = 1 2 4 k - 1 n + 1 .
Hence, if a [ 0 , 1 ] Q , we conclude that there exists k Z and n I N such that either (12) or (14) admits a rational solution. It is clear if a [ 0 , 1 ] Q , neither (12) nor (14) are solvable for such k Z and n N . □
In general, we distinguish between the odd and the even come-back events, therefore we define the following mappings.
Definition 2.
According to Theorem 1, we define the come-back mappings (for construction of the parameter a) as
e ( k , n ) : Z × 2 N [ 0 , 1 ] ; e ( k , n ) : = 1 2 4 k n + 1 .
and
o ( k , n ) : Z × ( 2 N + 1 ) [ 0 , 1 ] ; o ( k , n ) : = 1 2 4 k - 1 n + 1 .
Lemma 2.
Let k Z and n be an even natural number, then the following identities hold
1.
e ( - k , n ) + e ( k , n ) = 1 .
2.
o ( - k , n ) + o ( k , n ) = 1 - 1 2 n .
Proof. 
For 1 using the definition of the mapping e,
e ( - k , n ) = 1 2 4 ( - k ) n + 1 = - 1 2 4 k n - 1 = - 1 2 4 k n + 1 - 2 = 1 - e ( k , n ) .
For 2: similar to 1, we can easily prove that the following identity is true:
o ( - k , n ) = 1 - 1 2 n - o ( k , n )
The geometrical interpretation of the results (17) and (19) leads to an interesting observation about the up-down behavior of the sequence x a , namely if the first-come-back n is even, making k loops in the positive direction or in the negative one have the same path (i.e., the same range). While if the first-come-back n is odd, the direction affects the range of the up-down sequence. However, we could admit similar result in the even case for large n only.
Lemma 3.
Let us consider k Z and n in N , then the following identities hold
1.
e ( k , n ) = e ( p k , p n ) for all p I N .
2.
o ( k , n ) = o p k - p - 1 4 , p n for all odd p N such that ( p - 1 ) / 4 N .
Proof. 
For 1 using the definition of the mapping e,
e ( p k , p n ) = 1 2 4 ( p k ) p n + 1 = e ( k , n ) .
For 2 we have
o ( k , n ) = 1 2 4 k - 1 n + 1 = 1 2 4 p k - p p n + 1 = 1 2 4 p k - 1 - ( p - 1 ) p n + 1 = 1 2 4 ( p k - ( p - 1 ) / 4 ) - 1 p n + 1 = o p k - ( p - 1 ) / 4 , p n .
The results in Lemma 3 mean that the range of x o ( k , n ) and x e ( p k , p n ) are identical for all p N but x o ( k , n ) and x o ( p k , p n ) are identical only for all odd p N such that ( p - 1 ) / 4 N .
Theorem 2.
If j e = min j | j ( 2 a - 1 ) / 2 I N and j o = min j | ( j ( 2 a - 1 ) + a ) / 2 I N , then the first-come-back n f of x a is given by
n f = min 2 j e , 2 j o + 1 .
Proof. 
Suppose that x a ( 1 ) = 0 , then the 1st and the 2nd values of the sequence x a are given by
x a ( 1 ) = x a ( 0 ) + a π = a π , x a ( 2 ) = a π + ( a - 1 ) π = ( 2 a - 1 ) π .
Per recurrence, we deduce the 2 j t h and the ( 2 j + 1 ) t h values of the sequence x a :
x a ( 2 j ) = j ( 2 a - 1 ) π ,
x a ( 2 j + 1 ) = ( j ( 2 a - 1 ) + a ) π .
Thus, we can define the first come-back, if at least one of the following conditions is satisfied
j e ( 2 a - 1 ) = 2 k e ; for k e I N ( j o ( 2 a - 1 ) + a ) = 2 k o ; for k o I N
Since per definition the natural numbers k e and k o represent the number of loops needed till the first-come-back, numbers
k e = j e ( 2 a - 1 ) 2 and k o = j o ( 2 a - 1 ) + a 2 .
Thus, according to (22) and (23), we can determine the first come-back as a function of the minimal number of loops needed as follows:
n f = min 2 j e , 2 j o + 1 .
According to the result of the theorem above, we remarked that the range of the sequence x a contains exactly n f elements if n f = 2 j e . However, if n f = 2 j o + 1 , then the range of the sequence x a contains at least n f element.
(i)
If n f = 2 k e , then n f is an even number, it follows
x a ( n f + 1 ) = x a ( n f ) + a π = 0 + a π = x a ( 1 ) .
By this way, we construct all previous elements and x a ( 2 ) , x a ( 3 ) , , x a ( n f ) .
(ii)
If n f = 2 k o + 1 , then n f is an odd number, it follows
x a ( n f + 1 ) = x a ( n f ) + ( a - 1 ) π = 0 + ( a - 1 ) π x a ( 1 ) .
Therefore, we can construct other elements x a ( n f + 1 ) , which are not necessary in the set x a ( 0 ) , , x a ( n f + 1 ) . To make the following observations clear, we consider the following examples.
Definition 3.
We define the following indices:
(i) 
The first-come-back index n f is defined as n f = min { n | x a ( n ) = x a ( 0 ) } .
(ii) 
If Range ( x a ) : = { x a ( 0 ) , , x a ( m ) } , then we denote by n r = m , the maximal index in the range of the sequence R ( x a ) .
(iii) 
The cardinal number of the range R ( x a ) of the sequence x a is defined by n c : = Range ( x a ) .
Example 1.
For a = 4 / 5 , we compute the values of the first 15 elements of the sequence x 4 / 5 :We deduce the indices n f , n c and n r (see below Table 2). Similarly, we compute the same indices for different parameters a. We also must mention that an interesting open question could be deduced from the results in Table 2 and Table 3. Namely, the relationship between the three indices n f , n c and n r .

3. Construction of Geometric Patterns in Islamic Design

In general, the Islamic design employs complex polygons with vertices on a circle. In our case, these polygons represent connected graphs G ( a , α ) with vertices V ( a , α ) and edges E ( a , α ) , where α represents the initial value of the sequence x a α . For a given parameter a and an initial value α , the vertices set is defined as spatial points in I R 2 by:
V ( a , α ) : = p a α ( 0 ) , p a α ( 1 ) , , p a α ( n r ) ,
where p a α ( t ) = p ( x a α , t ) = cos ( x a α ( t ) ) , sin ( x a α ( t ) ) I R 2 for t = 0 , 1 , , n r .
The set of edges E ( a , α ) of the closed graph is
E ( a , α ) : = p a α ( 0 ) , p a α ( 1 ) , p a α ( 1 ) , p a α ( 2 ) , , p a α ( n r - 1 ) , p a α ( n r ) .
where n r denotes the range number and a is rational number in [ 0 , 1 ] .
Definition 4 (single cell).
A cell C ( q , r , a , α ) defining a design is a graph G ( a , α ) with colored sub-areas generated by a sequence x a α on a circle with radius r > 0 and a center q = ( q x , q y ) I R 2 . We denote by D ( I R 2 ) , the set of all possible geometric patterns of the Islamic designs in I R 2 . (see Figure 2b,d,f).
Throughout this paper, starting from a single cell, we construct special class of Islamic designs and patterns. Any geometrical transformation of a single cell leads to the formation of symmetrical patterns. Therefore, geometric patterns in the Islamic design, at least in the context of this work, is a repeated and/or modified planar transformation of a single cell.
Definition5 (Planar transformations).
For a given v I R , we define the vertical translation mapping f v as
f v : D ( I R 2 ) D ( I R 2 ) f ( 0 , v ) ( C ( ( q x , q y ) , r , a , α ) ) = C ( ( q x , q y + v ) , r , a , α ) .
The horizontal translation mapping f h is defined by
f h : D ( I R 2 ) D ( I R 2 ) f ( h , 0 ) ( C ( ( q x , q y ) , r , a , α ) ) = C ( ( q x + h , q y ) , r , a , α ) .
The diagonal translation mapping f ( v , h ) is defined by
f ( v , h ) : D ( I R 2 ) D ( I R 2 ) f ( v , h ) ( C ( ( q x , q y ) , r , a , α ) ) = C ( ( q x + h , q y + v ) , r , a , α ) ,
where h 0 and v 0 .
The Design (space filling) mapping F ( m v , m h ) is defined by
F ( m v , m h ) : D ( I R 2 ) D ( I R 2 ) F ( v , h ) ( C ( q , r , a , α ) ) = h = 1 m h v = 1 m v f ( K v v , K h h ) ( C ( ( q x , q y ) , r , a , α ) ) ,
where m v , m h I N , K v , K h are positive constants.
The radial translation mapping f s is defined by
f s : D ( I R 2 ) D ( I R 2 ) f s ( C ( q , r , a , α ) ) = C ( q , s r , a , α ) .
To construct designs according to all constructions presented in this paper, we suggest the following mathematical definition of geometric patterns in Islamic designs:
Definition6 (Geometric Patterns in Islamic Design).
A Circular Islamic Design D ( m v , m h , C ( q , r , a , α ) ) with m v × m h I N 2 elements is a symmetrical repeated placement of circular cells of the form C ( q , r , a , α ) .
Please note that any single cell from D ( I R 2 ) represents the class of geometric patterns in Islamic design. Moreover, any composite cell could be considered as a single cell, thus it is also an element of the Circular Islamic design. The mathematical construction of the circle-based class of geometric patterns in Islamic design is then achieved by using the geometrical translational and transformation mappings given by (29)–(33). More details are presented in the following section.

4. Numerical Simulation Islamic Geometric Patterns

In this work, we discuss four classes of circular Islamic designs, namely:
  • Modeling of Circular Islamic Design using single cells (see Appendix B).
  • Modeling of Circular Islamic Design using composite cells.
  • Modeling of Circular Islamic Design using circular cells.
  • Modeling of Circular Islamic Design using overlapping cells.

4.1. Modeling with Single Cells

The first class of Islamic designs presented in this work is achieved by arranging single symmetrical cells in vertical and horizontal forms. For any even number n, representing the first-come-back and any integer k representing the number of loops needed, we construct a sequence x a α using the Formula (12). Thus, we can generate Islamic designs using single cells using vertical and horizontal translation as follows:
D ( 1 , 1 , C ( q , r , a , α ) ) = f ( K v , K h ) C ( q , r , a , α ) ,
see Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 (left and middle figures). To generate m v cells in the vertical direction and m h cells in the horizontal direction, we use the following algorithmic translations:
D ( m v , m h , C ( q , r , a , α ) ) = h = 1 m h v = 1 m v f ( K v , K h ) C ( q , r , a , α ) .
where the distance needed for constructing non-overlapping shapes are given by K v for the vertical direction and by K h for the horizontal direction.
For even number n and for different values of k, according to theorem 1, we construct different symmetrical single cells. We show the resulting shapes in the first column of Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7. In the second column, we plot a 2 × 2 design and in the third column we present the result of 4 × 4 designs. Moreover, between the 2 × 2 design, we observe the presence of the contour of a new design; see the second column.
It is important to note that changing the initial value α leads to other designs, especially with totally different contours between the single cells. Moreover, it should be stressed that all our simulations are presented as the output of mathematical function without any modifications. For more interesting definition and historical evolution of the Islamic arts we refer to [1]. In Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7, we show five different designs using even numbers n and different number of loops: (a) e ( 1 , 8 ) , (b) e ( 1 , 12 ) , (c)  e ( 3 , 34 ) , (d) e ( 5 , 44 ) and (e) e ( 12 , 78 ) .

4.2. Modeling with Composite Cells

The second class of Islamic designs is achieved by arranging four single non-symmetrical cells in vertical and horizontal forms. In this case, the number n is odd. Because of the non-symmetry of the cells, we can construct a composite cell using four identical cells, using four different initial values of the same sequence x a α z given the Formula (12), where α z is given by x a α z ( 0 ) = α z . The governing function generating symmetrical shapes using composite cells is defined by
D ( 1 , 1 , C ( q z , r , a , ( α z ) 1 z 4 ) ) = z = 1 4 C ( q z , r , a , α z ) ,
where x a α z ( 0 ) = α z , r a constant radius of the circle and n = 0 , , n r .
After the construction of the composite single cells (36) for a satisfying the theorem 1, we use different initial values α z for z = 1 , 2 , 3 , 4 , for instance α 1 = π / 4 , α 1 = 3 π / 4 , α 1 = 5 π / 4 and α 1 = 7 π / 4 and similar to the first case, we use vertical and horizontal translation to simulate the patterns showed in Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12, second and third column. For m v cells in the vertical direction and m h cells in the horizontal direction, we use the following translations; see the third column in Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12. To generate m v cells in the vertical direction and m h cells in the horizontal direction, we use the following algorithmic translations:
D m v , m h , D ( 1 , 1 , C ( q z , r , a , ( α z ) ) ) = h = 1 m h v = 1 m v f ( K v , K h ) D ( 1 , 1 , C ( q z , r , a , ( α z ) ) ) ,
where the distance needed for constructing non-overlapping four cells shapes are given by K v for the vertical directions and by K h for the horizontal direction. For an odd number n and for different values of k, according to theorem 1, we construct different non-symmetrical single cells using odd numbers n and different number of loops: (a) o ( 1 , 9 ) , (b) o ( 1 , 15 ) , (c) o ( 2 , 17 ) , (d) o ( 3 , 27 ) and (e) o ( 4 , 35 ) . In the first column of the Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12, we show the resulting shapes of sequences with odd numbers and different initial values. In the second column, we plot the corresponding m v × m h = 2 × 2 composite cell. In the third column, we present the result of 4 × 4 designs. Moreover, between the 2 × 2 designs, we remarked the presence new shapes; see the second column in Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12.

4.3. Modeling with Circular Composite Cells

The third class of Islamic designs is the composition of symmetrical cells in circular forms and a centered design. For even numbers n z , z = 1 , 2 , , H , we construct composite design using single cells only. Because of the non-symmetry of the cells, we can construct a composite cell using four identical cells, using four different initial values of the same sequence x a α z given by Formula (12), where α z is given by x a α z ( 0 ) = α z . This could be clearly seen in the first row of Figure 13. The governing function generating symmetrical composite cell is defined by:
D ( 1 , 1 , C ( q z , ( r z ) , a , ( α z ) ) ) = C ( q 0 , r 0 , a 0 , α z ) + z = 1 H C ( q z , r z , a z , α z ) ,
where H I N refers to the number of surrounded circles, q 0 = ( 0 , 0 ) is the center of the centered design, the coordinates q z = K cos 2 z π H , K sin 2 z π H represent the centers of surrounded circles and r z is the corresponding radius for z = 1 , 2 , , H . Similar to the first and second cases, we use vertical and horizontal translation to simulate the patterns showed in Figure 13.
D ( m v , m h , C ( q z , ( r z ) , ( a z ) , ( α z ) ) ) = h = 1 m h v = 1 m v f ( K v , K h ) D ( 1 , 1 , C ( q z , ( r z ) , ( a z ) , ( α z ) ) ) ,
where m v represents the number of composite cells in the vertical direction and m h the number of composite cells in the horizontal direction (see the third column in Figure 13).
In the first row of the Figure 13, we plot the design using H = 5 , m v = m h = 2 . For the case (a), we have a 0 = e ( 15 , 34 ) and a 1 = e ( 5 , 78 ) . In the second row, we simulate a different design using H = 6 , m v = m h = 2 . For the case (b), a 0 = e ( 3 , 20 ) and a 1 = e ( 1 , 26 ) . In the third row, case (c), we use H = 10 , m v = m h = 2 , a 0 = e ( 5 , 44 ) and a 1 = e ( 3 , 34 ) .
For these three designs, we remarked the presence of three different contours between the 2 × 2 circles (see the second row of Figure 13).

4.4. Modeling with Overlapping Cells

The fourth class of Islamic designs is constructed by using overlapping single cells with the same center and increasing radius. Because of the symmetry of the cells, the number n should be an even number. Thus, we can construct overlapping cells with non-necessary identical sequences x a z α z with different initial values α z and parameters a z . The governing function generating symmetrical overlapping cells is defined by:
D ( 1 , 1 , C ( q 0 , r z , a z , α z ) ) = s = 1 H f s ( C ( q 0 , r z , a z , α ) ) = C ( q 0 , s r z , a z , α z ) ,
where q 0 = ( 0 , 0 ) is the center of the second centered design, the coordinates, s r z > 0 . To construct more overlapping designs, we apply the radial transformation mapping several times. For z = 1 , 2 , , H , the governing function generating this type of designs is given as follows:
D ( m v , m h , C ( q z , ( r z ) , ( a z ) , ( α z ) ) ) = h = 1 m h v = 1 m v f ( s ) D ( 1 , 1 , C ( q z , s r z , a z , α z ) ) ,
where for z = 1 , 2 , , H , the positive decreasing sequence ( r z ) represents the radius of the overlapping circles. In Figure 14, we see six different overlapping designs for H = 4 , α z = 0 . In Figure 15, we present repeated overlapping designs using different colors with different parameters, (see the titles of the Figure 15). For Figure 14, the following entries are implemented:
a.
e 1 ( 1 , 20 ) = e 2 ( 1 , 20 ) = e 3 ( 1 , 20 ) = e 4 ( 1 , 20 ) .
b.
e 1 ( 5 , 42 ) , e 2 ( 5 , 46 ) , e 3 ( 5 , 74 ) and e 4 ( 5 , 74 ) .
c.
e 1 ( 3 , 42 ) , e 2 ( 1 , 32 ) , e 3 ( 4 , 42 ) and e 4 ( 3 , 64 ) .
d.
e 1 ( 5 , 66 ) , e 2 ( 1 , 44 ) , e 3 ( 3 , 64 ) and e 4 ( 4 , 64 ) .
e.
e 1 ( 1 , 14 ) , e 2 ( 5 , 42 ) , e 3 ( 4 , 62 and e 4 ( 3 , 40 ) .
f.
e 1 ( 3 , 42 ) , e 2 ( 4 , 36 ) , e 3 ( 5 , 68 ) and e 4 ( 1 , 44 ) .

5. Concluding Remarks and Future Outlook

Using the geometrical properties of connected and closed graphs with vertices on the unit circle, we presented a new mathematical approach for modeling an interesting class of geometric patterns in Islamic design. The symmetries of vertices on the one-dimensional manifold leads to the construction of two-dimensional planar patterns, which are in the literature well known as a part of Islamic design. The resulting designs define beautiful fine arts and are, in general, expressed using basic mathematical transformations. Moreover, we believe that the class of geometric patterns presented in this work could be seen as a starting point for advanced studies in modeling of other classes of Islamic design using advanced techniques in Graph Theory, in modeling of more complex designs, or in characterizing the symmetry phenomenon in nature. In light of the simple techniques used here, the present geometrical framework can also serve as a valuable exercise in mathematics education. Furthermore, a generalization of the proposed work to curved surface inputs is an interesting avenue for future research works, in particular for developing similar and easy models to contribute to making mathematics more attractive for students as well as for non-mathematician researchers.

Funding

This work is supported by the University of Sharjah (UOS), Research group MASEP, RGS-2018-2019, No. 1802144062-P and Research project No. 1802144069.

Acknowledgments

I would like to thank Prof. Madjid Merabti and Prof. Abdelaziz Soufyane for their fruitful and valuable comments. In addition, I would like to thank the anonymous reviewers for their comments and suggestions.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. The Sequence x a

%% The sequence function  sequenceID.m
function f = sequenceID(n,x0,a)
ff = zeros(n,1);
ff(1) = x0;
for k = 1:n
     p=mod(k,2);
if(p==1)
     ff(k+1) = ff(k)+a*pi;
else
     ff(k+1) = ff(k)+(a-1)*pi ;
end
end
f=ff;
%---------------------------------------------

Appendix B. The Main Code for Generating Islamic Geometric Patterns

function ID=IslamicDesignE(k,N,x0,mv,mh)
     % N should even, for the odd case use the thm formula
     % construction the parameter of the sequence
     a=((4*k)/N+1)/2;
     % construction of the sequence
     X=sequenceID(N,x0,a);
     % adapting the distance between the cells
     tx=2; ty=1.85;
   % generating designs
   for r=1:mv
   for t=1:mh
A=tx*t+cos(X);
B=ty*r+sin(X);
   hold on; fill(A,B,’k’);
   hold on; plot(A,B,’r’);
     axis off; axis equal;
     end
   end
ID=a;
%---------------------------------------------
Compile: IslamicDesignE(1,20,0,5,5)

References

  1. Abas, S.J.; Salman, A.S. Symmetries of Islamic Geometrical Patterns; World Scientific Publishing Company: Singapore, 1994. [Google Scholar]
  2. Denny, W.B. The Ceramics of the Mosque of Rüstem Pasha and the Environment of Change; Garland Publishers, Inc.: New York, NY, USA, 1977. [Google Scholar]
  3. Grube, E.J. Islamic Pottery; Faber & Faber, Ltd.: London, UK, 1976. [Google Scholar]
  4. Humbert, C. Islamic Ornamental Design; Faber & Faber, Ltd.: London, UK, 1980. [Google Scholar]
  5. van der Hoop, A.N.J.Th.àTh. Indonesian Ornamental Design; A.C. Nix & Co.: Bandoeng, Indonesia, 1949. [Google Scholar]
  6. Tennant, R.F. Subgroup Lattices for Crystallografic Groups. Vis. Math. Art Sci. Electron. J. 2002, 4, 4. [Google Scholar]
  7. Arnone, S.; Falcolini, C.; Moauro, F.; Siccardi, M. On Numbers in Different Bases: Symmetries and a Conjecture. Exp. Math. 2017, 26, 197–209. [Google Scholar] [CrossRef]
  8. Chossat, P.; Golubitsky, M. Symmetry increasing bifurcations of chaotic attractors. Physica D 1988, 32, 423–436. [Google Scholar] [CrossRef]
  9. Hodgson, C.D.; Weeks, J.R. Symmetries, Isometries and Length Spectra of Closed Hyperbolic Three-Manifolds. Exp. Math. 1994, 3, 261–274. [Google Scholar] [CrossRef]
  10. Weyl, H. Symmetry; Princeton University Press: Princeton, NJ, USA, 1952. [Google Scholar]
  11. Field, M.; Golubitsky, M. Symmetry in Chaos: A Search for Pattern in Mathematics, Art, and Nature; SIAM: Philadelphia, PA, USA, 2009. [Google Scholar]
  12. Metzler, W.; Beau, W.; Frees, W.; Ueberla, A. Symmetry and Self-Similarity with Coupled Logistic Maps. Zeitschrift für Naturforschung 1987, 42, 310–318. [Google Scholar] [CrossRef]
  13. Zahri, M.; Khallaf, N. Some Results on Riemannean Multiple Barycenters. J. Appl. Math. Phys. 2014, 2, 503–509. [Google Scholar] [CrossRef]
  14. Barton, M.; Shi, L.; Kilian, M.; Wallner, J.; Pottmann, H. Circular arc snakes and kinematic surface generation. In Computer Graphics Forum; Blackwell Publishing Ltd.: Oxford, UK, 2013; Volume 32, pp. 1–10. [Google Scholar]
  15. Bo, P.; Pottmann, H.; Kilian, M.; Wang, W.; Wallner, J. Circular Arc Structures; Association for Computing Machinery: New York, NY, USA, 2011; Volume 30, p. 101. [Google Scholar]
  16. Song, X.; Aigner, M.; Chen, F.; Jüttler, B. Circular spline fitting using an evolution process. J. Comput. Appl. Math. 2009, 231, 423–433. [Google Scholar] [CrossRef]
Figure 1. Examples of geometric patterns in Islamic design.
Figure 1. Examples of geometric patterns in Islamic design.
Symmetry 11 00517 g001
Figure 2. The path of the periodic sequence x a ( n ) on S 1 . (a) a = 0 ; k = - 1 ; n = 4 ; (b) a = 0 . 25 ; k = - 1 ; n = 8 ; (c) a = 0 . 5 ; k = 0 ; n = 2 ; (d) a = 0 . 75 ; k = 1 ; n = 8 ; (e) a = 1 ; k = 1 ; n = 3 ; (f) a = 0 . 7 ; k = 1 ; n = 10 .
Figure 2. The path of the periodic sequence x a ( n ) on S 1 . (a) a = 0 ; k = - 1 ; n = 4 ; (b) a = 0 . 25 ; k = - 1 ; n = 8 ; (c) a = 0 . 5 ; k = 0 ; n = 2 ; (d) a = 0 . 75 ; k = 1 ; n = 8 ; (e) a = 1 ; k = 1 ; n = 3 ; (f) a = 0 . 7 ; k = 1 ; n = 10 .
Symmetry 11 00517 g002
Figure 3. Single cells Patterns for k = 1 and n = 8 .
Figure 3. Single cells Patterns for k = 1 and n = 8 .
Symmetry 11 00517 g003
Figure 4. Single cells Patterns for k = 1 and n = 12 .
Figure 4. Single cells Patterns for k = 1 and n = 12 .
Symmetry 11 00517 g004
Figure 5. Single cells Patterns for k = 3 and n = 34 .
Figure 5. Single cells Patterns for k = 3 and n = 34 .
Symmetry 11 00517 g005
Figure 6. Single cells Patterns for k = 5 and n = 44 .
Figure 6. Single cells Patterns for k = 5 and n = 44 .
Symmetry 11 00517 g006
Figure 7. Single cells Patterns for k = 12 and n = 78 .
Figure 7. Single cells Patterns for k = 12 and n = 78 .
Symmetry 11 00517 g007
Figure 8. Composite cells Patterns for k = 1 and n = 9 .
Figure 8. Composite cells Patterns for k = 1 and n = 9 .
Symmetry 11 00517 g008
Figure 9. Composite cells Patterns for k = 1 and n = 15 .
Figure 9. Composite cells Patterns for k = 1 and n = 15 .
Symmetry 11 00517 g009
Figure 10. Composite cells Patterns for k = 2 and n = 17 .
Figure 10. Composite cells Patterns for k = 2 and n = 17 .
Symmetry 11 00517 g010
Figure 11. Composite cells Patterns for k = 3 and n = 27 .
Figure 11. Composite cells Patterns for k = 3 and n = 27 .
Symmetry 11 00517 g011
Figure 12. Composite cells Patterns for k = 4 and n = 35 .
Figure 12. Composite cells Patterns for k = 4 and n = 35 .
Symmetry 11 00517 g012
Figure 13. Circular composite Patterns.
Figure 13. Circular composite Patterns.
Symmetry 11 00517 g013aSymmetry 11 00517 g013b
Figure 14. Overlapping Patterns for ( k 0 , k 1 , k 2 , k 3 ) loops and ( n 0 , n 1 , n 2 , n 3 ) come-backs.
Figure 14. Overlapping Patterns for ( k 0 , k 1 , k 2 , k 3 ) loops and ( n 0 , n 1 , n 2 , n 3 ) come-backs.
Symmetry 11 00517 g014
Figure 15. Repeated overlapping patterns for ( k 0 , k 1 , k 2 , k 3 ) loops and ( n 0 , n 1 , n 2 , n 3 ) come-backs.
Figure 15. Repeated overlapping patterns for ( k 0 , k 1 , k 2 , k 3 ) loops and ( n 0 , n 1 , n 2 , n 3 ) come-backs.
Symmetry 11 00517 g015
Table 1. Value of the up-down sequence using different values of the parameter a.
Table 1. Value of the up-down sequence using different values of the parameter a.
a x a ( 0 ) x a ( 1 ) x a ( 2 ) x a ( 3 ) x a ( 4 ) x a ( 5 ) x a ( 6 ) x a ( 7 ) x a ( 8 ) x a ( 9 )
000 - π - π - 2 π - 2 π - 3 π - 3 π - 4 π - 4 π
1 4 0 1 4 π - 2 4 π - 1 4 π - 4 4 π - 3 4 π - 6 4 π - 5 4 π - 8 4 π - 7 4 π
1 2 0 1 2 π 0 1 2 π 0 1 2 π 0 1 2 π 0 1 2 π
3 4 0 3 4 π 2 4 π 5 4 π 4 4 π 7 4 π 6 4 π 9 4 π 8 4 π 11 4 π
10 π π 2 π 2 π 3 π 3 π 4 π 4 π 5 π
Table 2. Value of the up-down sequence using different values of the parameter a.
Table 2. Value of the up-down sequence using different values of the parameter a.
n x 4 5 (n)Remarks
00Initial Element in the range: 1
1 4 5 π New Element: 2
2 3 5 π New Element: 3
3 7 5 π New Element: 4
4 6 5 π New Element: 5
5 10 5 π = 2 π 0 First-come-back and Existent Element n f = 5
6 9 5 π New Element: 6
7 3 5 π Existent Element x a ( 7 ) = x a ( 2 )
8 2 5 π New Element: 7
9 6 5 π Existent Element x a ( 9 ) = x a ( 4 )
10 π New Element: 8
11 9 5 π Existent Element x a ( 11 ) = x a ( 6 )
12 8 5 π New Element: 9
13 2 5 π Existent Element x a ( 13 ) = x a ( 8 )
14 1 5 π New Element: n c = 10 but n r = 14
Table 3. Examples of values of n f , n c and n r using different values of the parameters a.
Table 3. Examples of values of n f , n c and n r using different values of the parameters a.
Index/a01 1 2 3 4 1 4 7 10 4 5 3 4 2 3 5 7
n f 4388105897
n c 2228810108614
n r 4228810148920

Share and Cite

MDPI and ACS Style

Zahri, M. Mathematical Modeling of a Class of Symmetrical Islamic Design. Symmetry 2019, 11, 517. https://doi.org/10.3390/sym11040517

AMA Style

Zahri M. Mathematical Modeling of a Class of Symmetrical Islamic Design. Symmetry. 2019; 11(4):517. https://doi.org/10.3390/sym11040517

Chicago/Turabian Style

Zahri, Mostafa. 2019. "Mathematical Modeling of a Class of Symmetrical Islamic Design" Symmetry 11, no. 4: 517. https://doi.org/10.3390/sym11040517

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