Next Article in Journal / Special Issue
Physics of the Chemical Asymmetry of the Cell Membrane: Implications in Gene Regulation and Pharmacology
Previous Article in Journal
New Nonlocal Symmetries of Diffusion-Convection Equations and Their Connection with Generalized Hodograph Transformation
Previous Article in Special Issue
Symmetry Breaking and Establishment of Dorsal/Ventral Polarity in the Early Sea Urchin Embryo
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Initial Symmetry on the Global Symmetry of One-Dimensional Legal Cellular Automata

Department of Earth and Planetary Sciences, Kobe University, Nada, 657-8501 Kobe, Japan
Symmetry 2015, 7(4), 1768-1779; https://doi.org/10.3390/sym7041768
Submission received: 10 July 2015 / Revised: 22 August 2015 / Accepted: 25 September 2015 / Published: 29 September 2015
(This article belongs to the Special Issue Symmetry and Asymmetry in Biology)

Abstract

:
To examine the development of pattern formation from the viewpoint of symmetry, we applied a two-dimensional discrete Walsh analysis to a one-dimensional cellular automata model under two types of regular initial conditions. The amount of symmetropy of cellular automata (CA) models under regular and random initial conditions corresponds to three Wolfram’s classes of CAs, identified as Classes II, III, and IV. Regular initial conditions occur in two groups. One group that makes a broken, regular pattern formation has four types of symmetry, whereas the other group that makes a higher hierarchy pattern formation has only two types. Additionally, both final pattern formations show an increased amount of symmetropy as time passes. Moreover, the final pattern formations are affected by iterations of base rules of CA models of chaos dynamical systems. The growth design formations limit possibilities: the ratio of developing final pattern formations under a regular initial condition decreases in the order of Classes III, II, and IV. This might be related to the difference in degree in reference to surrounding conditions. These findings suggest that calculations of symmetries of the structures of one-dimensional cellular automata models are useful for revealing rules of pattern generation for animal bodies.

1. Introduction

Animals have many kinds of pattern formation. Leopards living in closed habitats, such as forests, have a frequency of stripes that is much higher, more irregular, and more complex than those living in open habitats [1]. These body patterns reflect adaptations to the surrounding environment and, thus, pattern formation is substantially related to the conservation of species [1]. When we observe pattern formation, we read “texture”, which includes shape, among other factors. Texture is measured by terms including spectrum, entropy, and fractal dimension [2,3]. The method of measurement combining entropy and symmetry, the two-dimensional discrete Walsh analysis, can be applied to analyses of two-dimensional surfaces [3,4,5]. In previous studies, the two-dimensional discrete Walsh analysis has been used in theoretical studies with one-dimensional legal cellular automata. As a result, studies have revealed the theoretical process of pattern formation, which reduces entropy and increases specific symmetry (self-organization) [4], and the real biological process of pattern formation, which reduces entropy and increases specific symmetry. They have also revealed that the real biological process has phase transitions, as seen in calculations with the two-dimensional Ising model [5]. However, the calculations in previous studies [3,4,5] always start with the condition of a random state. Pattern formation of a biological or animal surface does not necessarily start from such a random state, but rather from a regular state. Thus, clarifying the theoretical process of regular pattern formation has the possibility of revealing general features and the real biological process of pattern formation.
The purpose of this study was to clarify the theoretical process of pattern formation, starting from a regular state. We applied a two-dimensional discrete Walsh analysis to one-dimensional legal cellular automata.
The structure of this paper is as follows. In Section 2, we explain the method of analysis and data on one-dimensional legal cellular automata (CAs). In Section 3, we explain the data of the regular initial conditions. In Section 4, we briefly review the discrete Walsh analysis for calculating the symmetry and symmetropy of CA patterns. In Section 5, we describe and discuss the results and consider the symmetry relationships between the initial conditions and final patterns. Section 6 provides our conclusions.

2. Methods

The Discrete Walsh Analysis

We briefly outline the discrete Walsh analysis here (Figure 1). For details of the mathematical procedures, see previous studies [3,4,6]. The Walsh function w a l ( r , x ) of order r and argument x can be represented over the interval 0 x < 1 as follows:
w a l ( r , x ) = i = 0 m 1 sgn [ ( cos 2 i π x ) r i ]
where r i = 0   or   1 , r = i = 0 m 1 r i 2 i , and m is the smallest positive integer such that 2 m > r . Dyadic addition of non-negative integers r and s is defined as r s = i = 0 f | r i s i | 2 i , where r = i = 0 f r i 2 i and s = i = 0 f s i 2 i . Consequently, the product of the two Walsh functions is given by w a l ( r , x ) w a l ( s , x ) = w a l ( r s , x ) . From this calculation, the Walsh function can be shown to form an orthonormal set. In this case, we can define a Walsh transform, similar to a Fourier transform.
Figure 1. (a) Examples of the two-dimensional discrete Walsh functions for m = n = 8 . Black represents +1, and white represents –1; and (b) four types of symmetry in terms of the discrete Walsh function.
Figure 1. (a) Examples of the two-dimensional discrete Walsh functions for m = n = 8 . Black represents +1, and white represents –1; and (b) four types of symmetry in terms of the discrete Walsh function.
Symmetry 07 01768 g001
The two-dimensional discrete Walsh function can be represented in matrix form as [ W m n ( i ,   j ) ] , where W m n ( i ,   j ) is the amount of the (m,n)th -order Walsh function in the ith row cell in the jth column.
This pattern can be written as [ x i j ] , where x i j is the amount of the gray level in the i t h row cell in the jth column and i , j = 0 , 1 , , N 1. If just two gray levels exist, for example, “black” and “white”, then x i j is typically represented by 1 and 0, respectively (Figure 1a). The two-dimensional discrete Walsh transform of the pattern [ x i j ] is given by:
a m n = 1 N 2 i = 0 N 1 j = 0 N 1 x i j W m n ( i ,   j )
where   m , n = 0 , 1 , , N 1. The functions a m n and ( a m n ) 2 are the two-dimensional Walsh spectrum and power spectrum, respectively.
The Walsh power spectrum can be normalized as follows:
P m n = ( a m n ) 2 K
where K = m = 0 N 1 n = 0 N 1 ( a m n ) 2 ( a 00 ) 2 . In this case, we obtain:
P m n = 1
where the sum is taken over all ordered pairs ( m , n ) for 0 m ,   n N 1 , except for ( m , n ) = ( 0 , 0 ) .
The spatial pattern is considered as an information source, consisting of dot patterns. The dot patterns emitted from the source are assumed to occur with the corresponding probabilities given by Equation (4).
Next, we consider the information entropy resulting from the symmetry of the pattern. Because the two-dimensional Walsh functions can be easily divided into four types of symmetry (Figure 1b), Equation (4) can be rewritten as:
i = 1 4 P i = 1
where vertical symmetric component:
P 1 = m = even , n = odd P m n
horizontal symmetric component:
P 2 = m = odd , n = even P m n
centro symmetric component:
P 3 = m = odd , n = odd P m n
and double symmetric component:
P 4 = m = even , n = even P m n
When the spatial pattern is regarded as an information source consisting of four types of symmetry, the corresponding probabilities are given by Equations (6)–(9). Applying the entropy function from information theory to these four symmetric components, we obtain:
S = i = 1 4 P i log 2 P i
As shown in Equations (5) and (10), S ranges from 0 to 2 bits. Since this entropy relates to the symmetry, it is called the “symmetropy” [6]. The symmetropy can be considered as a quantitative and objective measure of symmetry. If the amount of a certain component is larger than the amounts of the other three components, then the pattern is rich in symmetry related to the particular component. In this case, Equation (10) shows that S decreases. Conversely, if the amounts of the four components are almost equal to one another, then the pattern is poor in symmetry, and   S   increases. Previous studies have concentrated on symmetry [4]. Following this idea, in this study, we measured the symmetry of the initial condition and the final pattern to examine the effects on final patterns of the initial condition of symmetry in CA patterns.

3. Calculations

For simplicity, this paper uses only “classical” CAs according to the Wolfram classification [7]; that is, the one-dimensional legal CA that obeys the 32 possible legal totalistic rules. The totalistic rules involve nearest and next nearest neighbors, and each cell has two possible colors: white (0) or black (1). In addition, because the total number of neighborhoods of five cells is equal to 25, we can obtain the 32 possible legal totalistic rules. Wolfram [7,8] discovered that, although the patterns obtained with different rules differed in detail, they appeared to fall into four qualitative classes:
  • Class I: fixed point (rules 0, 4, 16, 32, 36, 48, 54, 60, and 62) (Figure 2a).
  • Class II: periodic (rules 8, 24, 40, 56, and 58) (Figure 2b).
  • Class III: chaotic (rules 2, 6, 10, 12, 14, 18, 22, 26, 28, 30, 34, 38, 42, 44, 46, and 50) (Figure 2c).
  • Class IV: complex (rules 20, 52) (Figure 2d).
Figure 2. Examples of final patterns in each class. (a) is an example of Class I, (b) of Class II, (c) Class III, and (d) is an example of Class IV.
Figure 2. Examples of final patterns in each class. (a) is an example of Class I, (b) of Class II, (c) Class III, and (d) is an example of Class IV.
Symmetry 07 01768 g002
In this paper, we do not use all codes in class I because the patterns are trivial. We measure the initial condition (26 × 1) and the final pattern. The number of cells is 26 × 26 because, to measure the number of symmetries, this number of cells is the minimum needed [4]. In the present study, we can use only two kinds of initial condition, horizontal symmetry and double symmetry, because we consider n = 2 in the discrete Walsh function Wmn (as in Figure 3).
Figure 3. An example of initial symmetry. (a) Double Symmetry is one example of D S 0 , and (b) Horizontal Symmetry is one example of H S 0 . In this study, the initial condition was limited to n = 2 because 21 = 2 = n. Thus, (a) shows the case of m = even, and (b) shows the case of m = odd.
Figure 3. An example of initial symmetry. (a) Double Symmetry is one example of D S 0 , and (b) Horizontal Symmetry is one example of H S 0 . In this study, the initial condition was limited to n = 2 because 21 = 2 = n. Thus, (a) shows the case of m = even, and (b) shows the case of m = odd.
Symmetry 07 01768 g003
First, we run two types of regular initial conditions; the horizontal symmetry with 32 conditions, abbreviated as H S 0 , and the double symmetry with 31 conditions, abbreviated as D S 0 . Then, under each initial condition, the amounts of symmetry were averaged for each class. We compared this result with the random initial condition in [4], abbreviated as R a 0 . Then, for each initial condition, we compared the relationship between the amount of symmetry in the initial condition and that in the final pattern. The method can be extended to 2D cellular automata, but requires a large number of analysis.

4. Results

4.1. Class

Figure 4 shows the amount of symmetropy for the three initial conditions ( R a 0 , D S 0 and H S 0 ) and the amount of symmetropy in each initial condition for the final patterns. The results show that under all initial conditions, the numbers of symmetries of the CA patterns correspond to three qualitative classes of CAs, identified as Classes II, III, and IV, according to Wolfram’s original classification [7,8]. Additionally, the result of the measurement of symmetropy under regular and random initial conditions shows a relational expression: Class II < Class IV < Class III. In all classes, the amount of symmetropy under R a 0 is larger than that under D S 0 and H S 0 . In R a 0 , the amount of symmetropy of all classes decreased as time passed. In contrast, in D S 0 and H S 0 , the amount of symmetropy of all classes increased as time passed. Thus, the final patterns of the regular initial condition ( D S 0 and H S 0 ) are more regular than those of the random initial conditions ( R a 0 ).

4.2. Four Types of Symmetry

Figure 4 shows increases or decreases in the amount of symmetry for four kinds of symmetries in initial conditions and in final patterns. In H S 0 , Classes II, III, and IV have all of the four kinds of symmetries: vertical symmetry, horizontal symmetry, centro symmetry, and double symmetry. However, in D S 0 , Classes II, III, and IV have only two kinds of symmetries: vertical symmetry and double symmetry (Figure 4).
Figure 4. Pie charts of the three types of initial condition ( H S 0 , D S 0 and R a 0 ). For H S 0 and R a 0 , all classes have all four types of symmetry: vertical symmetry, horizontal symmetry, centro symmetry, and double symmetry. In contrast, D S 0 with all classes had two types of symmetry: vertical symmetry and double symmetry. Regarding vertical symmetry and double symmetry, the relational expression was H S 0 < R a 0 < D S 0 , and regarding horizontal symmetry and centro symmetry, it was D S 0 < R a 0 < H S 0 .
Figure 4. Pie charts of the three types of initial condition ( H S 0 , D S 0 and R a 0 ). For H S 0 and R a 0 , all classes have all four types of symmetry: vertical symmetry, horizontal symmetry, centro symmetry, and double symmetry. In contrast, D S 0 with all classes had two types of symmetry: vertical symmetry and double symmetry. Regarding vertical symmetry and double symmetry, the relational expression was H S 0 < R a 0 < D S 0 , and regarding horizontal symmetry and centro symmetry, it was D S 0 < R a 0 < H S 0 .
Symmetry 07 01768 g004
In D S 0 , two kinds of symmetries, double symmetry and vertical symmetry, developed globally in the final formations. The design developed by D S 0 has a complicated pattern that lined up with a regular pattern (Figure 5a,b). Development was accompanied by a large increase in the amount of vertical symmetry. Thus, vertical symmetry contributed to the development of the spatial direction of the final patterns, whereas the double symmetry of D S 0 was preserved. In contrast, the design developed by H S 0 has all four kinds of symmetry (Figure 4), and the horizontal symmetry of H S 0 was not preserved.
Figure 5. Examples of final patterns with regular initial conditions. (a,b) belong to D S 0 , and (c,d) belong to H S 0 .
Figure 5. Examples of final patterns with regular initial conditions. (a,b) belong to D S 0 , and (c,d) belong to H S 0 .
Symmetry 07 01768 g005

4.3. The Ratio of Final Formation

The ratios of developed final patterns in each Class for DS 0 , Ra 0 , and HS 0 are shown in Table 1. The ratio of developed of final patterns is calculated by dividing the number of final pattern formation for each class (periodic or chaotic or complex) by a product of the number of initial condition (if H S 0 : 32, if D S 0 : 31) and the total number of classes (the sum rules of periodic, chaotic and complex: 23). For example, in the case of H S 0 ’s periodic in Table 1, the ratio of periodic is 12 is calculated from 89, 32, and 23 for the numbers of final pattern formation, initial condition, and the sum rules of periodic, chaotic and complex, respectively. The ratio of 12 is obtained as: 89(final patterns)/32rules/23(number of classes) × 100. According to Yamazaki et al. [4], R a 0   made final patterns using every rule and class; the ratio of formation was 100%. In contrast, for D S 0  and  H S 0 , the ratio of formation is less than 50% for all classes, and the developing final pattern ratios are 12% for Class II, 43%–44% for Class III, and 6% for Class IV (Table 1).
Table 1. The rate of final patterns of D S 0 and H S 0 .
Table 1. The rate of final patterns of D S 0 and H S 0 .
ClassHorizontal SymmetryDouble SymmetryRandom
Periodic1212100×
Chaotic4344100×
Complex66100×
In all classes, the rate of the forming pattern of R a 0 was larger than that of D S 0 and H S 0 . Comparing D S 0 with H S 0 , the rate of forming patterns of class III was larger than that of class II and class IV; also, in each class, the two types of D S 0 and H S 0 had the same rates. Data indicated by asterisks are from [4].

5. Discussion

5.1. Class

For all classes, the amount of symmetry for     H S 0 is larger than that for D S 0 , and the final pattern formations of D S 0 have a higher hierarchy than H S 0 ; in contrast, the final pattern formations of H S 0   had more broken regularity than H S 0   (Figure 5c,d). This result is interpreted as the double symmetry having a higher amount of symmetry than the horizontal symmetry, because D S 0 has an axis of symmetry with space and time, whereas H S 0 has an axis of symmetry with time alone.
From the measurement of symmetropy, developed pattern formations of H S 0 and   D S 0   have many types of symmetries that predominate as time passes. This result corresponds to the viewpoint of the Curie Symmetry Principle [4]. In contrast, the developed pattern formations of R a 0 show a situation in which one type of symmetry was predominant as time passed [4]; this result does not correspond to the viewpoint of the Curie Symmetry Principle [4]. Here, it is shown that the theoretical process of pattern formation reduced entropy and increased specific symmetry [4,5]; however, this process has been examined only from R a 0   [4,5]. In this study, the final pattern formations from regular initial conditions did not reduce entropy and increased specific symmetry. Thus, the Curie Symmetry Principle cannot be used to evaluate the proxy of pattern formation.
The final pattern formation of D S 0   is composed of the same scale forms, to stand in a line (Figure 5a,b). The pattern formations are supposed to be affected by “iteration” of the base rules of CAs, and this result is correlated with the previous studies, such as [9] that shows that the Walsh transform allows one to compute the correlation between the inputs and the outputs of the iterations of a cellular automaton. Generally, as systems are put under greater stress, the observed states lose symmetry and gain in complexity (both spatial and temporal). However, as the systems are further stressed, states with more symmetries and more structures appear. Moreover, symmetries and chaos—patterns and disorder—can coexist naturally within the same simple mathematical framework [10]. In this study, the natural coexistence of symmetry and chaos (iteration of CAs) is found in the final pattern formation of D S 0 ; this means a simple mathematical framework of a one-dimensional legal cellular automata. Complicated structure pattern formations are developed by a simple model; thus, the quantity of information needed to describe an object may be more than that needed to prescribe it. Thus, under an initial condition of strong regularity, theoretical phenomena revealed occurrences of iterations with regularity, and the concepts of the theoretical phenomena could be helpful to interpret the process and the circumstance of biological patterns with higher-complexity structures.

5.2. Four Types of Symmetry

For all classes, the amount of symmetropy for H S 0 is larger than that for D S 0 (Figure 4), because H S 0   has four types of symmetry, while D S 0   has only two types (Figure 4). This result coincides with that the final pattern of H S 0 , which shows steady state decay, and that of D S 0 , which shows switching of a steady state to a hierarchy.
H S 0   has four types of symmetry: vertical, double, horizontal, and centro symmetries. D S 0 has only two types: vertical and double symmetries. These results show that double and vertical symmetries are the basic ones for forming patterns, and whether they form random or regular patterns is affected by the amount of horizontal and centro symmetries. The direction of each symmetry shows that (1) vertical symmetry has an axis of symmetry with space, (2) horizontal symmetry has an axis of symmetry with time, (3) double symmetry has an axis of symmetry with space and time, and (4) centro symmetry has an axis of symmetry with space and time. Thus, combinations of double and vertical symmetries, which are the basic symmetries of pattern forming, have axes of symmetries with space and time. In contrast, the combination of horizontal and centro symmetries is a supplemental one, for pattern formation. Thus, the basic and supplemental symmetries have balance with space and time. Additionally, the initial condition constrains the combination of symmetries (the analysis decided the symmetry direction): vertical and double symmetries versus horizontal and centro symmetries.
We compared the amounts of the four types of symmetries among D S 0 ,   R a 0 , and  H S 0 . The amount of vertical symmetry and double symmetry increased in the order H S 0 < R a 0 < D S 0 , except for Class II’s vertical symmetry. Additionally, the amount of horizontal symmetry and centro symmetry increased in the order D S 0 < R a 0 < H S 0 .
In regular and random initial conditions, Class III has equal amount of symmetries of the different types of symmetries, but Classes II and IV have different amounts. Here, we explain biological patterns from these results. The organisms’ developed pattern formations can select the amounts of two types of symmetry (horizontal and centro) according to the circumstances, because both horizontal and centro symmetries have symmetry axes related to “time”. Here, time is related to the future. Thus, the combination of symmetries with wide varieties, as in this study, may provide the opportunity for understanding how to develop pattern formation (heredity or circumstances), and may be able to manage pattern formation, and to predict the future.

5.3. The Ratio of Final Formation

The growth design formations are limited in H S 0 and D S 0 , but there is no limit in R a 0 (Table 1). This means the variability of symmetries of R a 0 is larger than that of H S 0 and D S 0 . The numbers of developed final pattern rules for H S 0 and D S 0 are almost equal in every class, because the numbers of initial conditions are almost the same; 32 for H S 0 and 31 for D S 0 , and most of the H S 0 rules correspond to D S 0 ones.
The rules that did not develop final patterns under regular initial conditions need larger amounts and more symmetries. The ratio of rules developing final patterns to those with no developing final patterns was calculated for each Class (Table 1). The ratio of rules that developed final patterns under the regular initial condition decreased in the order of Class III, Class II, and Class IV. Thus, Class III does not need more of the amounts and types of symmetries, whereas Class IV requires much more in the way of amounts and many types of symmetries, and Class II requires moderate amounts and types of symmetry. We explain biological patterns from these theories. Class III is not influenced by circumstances, and allows easy development of final patterns, which means heredity. On the other hand, Class IV is influenced by circumstances, and Class II by circumstances and heredity. Animals having Class III biological patterns may select the circumstances by themselves, because their biological patterns are interpreted as genetic characteristics. The animals having Class IV biological patterns may adapt themselves to circumstances, because their biological patterns are interpreted as acquired characteristics. The animals having Class II biological patterns may be able to select a circumstance by themselves and may be able to adapt themselves to a circumstance, because their biological patterns are interpreted as random genetic characteristics and random circumstances. They may take an unplanned survival strategy.
Here, we discuss living animals. Although clams appear in Class III [11], their shell patterns are in heredity [12]. Additionally, the habitat of clams is decided by the height of ground level and drying time within 1 day, which has no relationship with pattern formation [13]. This might be supported by the hypothesis “The animals having Class III biological pattern formation may select circumstance by themselves, because their biological patterns are interpreted as genetic characters”. However, Neritina, a typical shellfish, appears in Classes II and IV [11]; their shell design pattern formation has been found along the structure of the growth lines [14]. Thus, Neritina may be affected by the circumstances around them. In addition, their shell patterns have been found to be of color bands [14]; thus, Neritina may have been affected by circumstance and heredity.

6. Conclusions

We examined the structure of one-dimensional legal cellular automata from the viewpoint of symmetry to reveal the rules of animals’ pattern generation. The following conclusions are derived. (1) Under all regular and random initial conditions, a relational expression, Class II <Class IV <Class III, was seen. In all classes, the pattern formation with regular initial conditions was more regular than that with random initial conditions. (2) We compared the pattern formation made with two types of regular initial conditions. The final pattern formations of H S 0 , which have more broken regularity, have four types of symmetries, whereas that of D S 0 has only two types of symmetries, which are of a higher hierarchy. Both regular initial conditions show situations in which many types of symmetry predominated as time passed. We showed that when pattern formation develops, amounts of symmetries and types of symmetries do not always decrease. (3) We found that growth design formations have limited possibilities with regular initial conditions; the results of calculating the ratio of developing final patterns decreases in the order of Class III, Class II, and Class IV.

Acknowledgements

We would like to give thanks to Masayuki Hyodo, Kazuhito Yamasaki, Takao Ubukata, Sunao Mochizuki and Yusuke Seto for valuable discussions and advices on the manuscript. This work was financially supported by a grant from the Japan Society for the Promotion of Science awarded to IT (15J00351).

Conflicts of Interest

The authors declares no conflict of interest.

References

  1. Allen, W.L.; Cuthill, I.C.; Scott-Samuel, N.E.; Baddeley, R. Why the leopard got its spots: Relating pattern development to ecology in felids. Proc. R Soc. B. 2010, 1373–1380. [Google Scholar] [CrossRef] [PubMed]
  2. Ubukata, T. Theoretical morphology of composite prismatic, fibrous prismatic and foliated microstructures in bivalves. Venus 2000, 59, 297–305. [Google Scholar]
  3. Nishiyama, Y.; Nanjo, K.Z.; Yamasaki, K. Geometrical minimum units of fracture patterns in two-dimensional space: Lattice and discrete Walsh functions. Phys. A 2008, 387, 6252–6262. [Google Scholar]
  4. Yamasaki, K.; Nanjo, K.Z.; Chiba, S. Symmetry and entropy of one-dimensional legal cellular automata. Complex Syst. 2012, 20, 352–361. [Google Scholar]
  5. Yamasaki, K.; Nanjo, K.Z.; Chiba, S. Symmetry and entropy of biological patterns: Discrete Walsh functions for 2D image analysis. BioSystems 2010, 103, 105–112. [Google Scholar] [CrossRef] [PubMed]
  6. Yodogawa, E. Symmetropy, an entropy-like measure of visual symmetry. Percept. Psychophys. 1982, 32, 230–240. [Google Scholar] [CrossRef] [PubMed]
  7. Wolfram, S. Universality and complexity, in cellular automata. Phys. D 1984, 10, 1–35. [Google Scholar] [CrossRef]
  8. Wolfram, S. A New Kind of Science; Wolfram Media, Inc.: Champaign, IL, USA, 2002. [Google Scholar]
  9. Martin, B. A Walsh exploration of Wolfram CA rules. In Proceedings of the 12th International Workshop on Cellular Automata, Hiroshima University, Higashi-Hiroshima, Japan, 12–15 September 2006; pp. 25–30.
  10. Stewart, I.; Golubitsky, M. Fearful Symmetry: Is God a Geometer? Blackwell Pub: Oxford, UK, 1992; p. 346. [Google Scholar]
  11. Schiff, J.L. Cellular Automata: Discrete View of the World; John Wiley & Sons Inc.: Hoboken, NJ, USA, 2008; p. 250. [Google Scholar]
  12. Peignon, J.M.; Gérard, A.; Naciri, Y.; Ledu, C.; Phélipot, P. Analyse du déterminisme de la coloration et de l'ornementation chez la palourde japonaise Ruditapes Philippinarum. Aqua. Liv. Res. 1995, 8, 181–185. (In French) [Google Scholar] [CrossRef]
  13. Akiyama, B.Y.; Saito, H.; Nanbu, R.; Tanaka, Y.; Kuwahara, H. The spatial distribution of Manila clam Ruditapes philippinarum associated with habitat environment in a sandy tidal flat on the coast of Matsunase, Mie Prefecture, Japan. Ecol. Civil. Eng. 2011, 14, 21–34. (In Japanese) [Google Scholar] [CrossRef]
  14. Habe, T. Gakken Picture Book of Shellfish; Gakken Holdings: Tokyo, Japan, 1975; p. 294. (In Japanese) [Google Scholar]

Share and Cite

MDPI and ACS Style

Tanaka, I. Effects of Initial Symmetry on the Global Symmetry of One-Dimensional Legal Cellular Automata. Symmetry 2015, 7, 1768-1779. https://doi.org/10.3390/sym7041768

AMA Style

Tanaka I. Effects of Initial Symmetry on the Global Symmetry of One-Dimensional Legal Cellular Automata. Symmetry. 2015; 7(4):1768-1779. https://doi.org/10.3390/sym7041768

Chicago/Turabian Style

Tanaka, Ikuko. 2015. "Effects of Initial Symmetry on the Global Symmetry of One-Dimensional Legal Cellular Automata" Symmetry 7, no. 4: 1768-1779. https://doi.org/10.3390/sym7041768

APA Style

Tanaka, I. (2015). Effects of Initial Symmetry on the Global Symmetry of One-Dimensional Legal Cellular Automata. Symmetry, 7(4), 1768-1779. https://doi.org/10.3390/sym7041768

Article Metrics

Back to TopTop