Lenia

Last updated January 24, 2024

An animation showing the movement of a glider in Lenia. Peek 2021-10-12 22-29.gif — An animation showing the movement of a glider in Lenia.

Lenia is a family of cellular automata created by Bert Wang-Chak Chan.^[1]^[2]^[3] It is intended to be a continuous generalization of Conway's Game of Life, with continuous states, space and time. As a consequence of its continuous, high-resolution domain, the complex autonomous patterns ("lifeforms" or "spaceships") generated in Lenia are described as differing from those appearing in other cellular automata, being "geometric, metameric, fuzzy, resilient, adaptive, and rule-generic".^[1]

Rules

Iterative updates

Let ${\mathcal {L}}$ be the lattice or grid containing a set of states $S^{\mathcal {L}}$ . Like many cellular automata, Lenia is updated iteratively; each output state is a pure function of the previous state, such that

\Phi (A^{0})=A^{\Delta t},\Phi (A^{\Delta t})=A^{2\Delta t},\ldots ,\Phi (A^{t})=A^{t+\Delta t},\ldots

where $A^{0}$ is the initial state and $\Phi :S^{\mathcal {L}}\rightarrow S^{\mathcal {L}}$ is the global rule, representing the application of the local rule over every site $\mathbf {x} \in {\cal {L}}$ . Thus $\Phi ^{N}(A^{t})=A^{t+N\Delta t}$ .

If the simulation is advanced by $\Delta t$ at each timestep, then the time resolution $T={\frac {1}{\Delta t}}$ .

State sets

Let $S=\{0,1,\ldots ,P-1,P\}$ with maximum $P\in \mathbb {Z}$ . This is the state set of the automaton and characterizes the possible states that may be found at each site. Larger $P$ correspond to higher state resolutions in the simulation. Many cellular automata use the lowest possible state resolution, i.e. $P=1$ . Lenia allows for much higher resolutions. Note that the actual value at each site is not in $[0,P]$ but rather an integer multiple of $\Delta p={\frac {1}{P}}$ ; therefore we have $A^{t}(\mathbf {x} )\in [0,1]$ for all $\mathbf {x} \in {\mathcal {L}}$ . For example, given $P=4$ , $\mathbf {A} ^{t}(\mathbf {x} )\in \{0,0.25,0.75,1\}$ .

Neighborhoods

Mathematically, neighborhoods like those in Game of Life may be represented using a set of position vectors in $\mathbb {R} ^{2}$ . For the classic Moore neighborhood used by Game of Life, for instance, ${\mathcal {N}}=\{-1,0,1\}^{2}$ ; i.e. a square of size 3 centered on every site.

In Lenia's case, the neighborhood is instead a ball of radius $R$ centered on a site, ${\mathcal {N}}=\{\mathbf {x} \in {\mathcal {L}}:\lVert \mathbf {x} \rVert _{2}\leq R\}$ , which may include the original site itself.

Note that the neighborhood vectors are not the absolute position of the elements, but rather a set of relative positions (deltas) with respect to any given site.

Local rule

There are discrete and continuous variants of Lenia. Let $\mathbf {x}$ be a vector in $\mathbb {R} ^{2}$ within ${\mathcal {L}}$ representing the position of a given site, and ${\mathcal {N}}$ be the set of sites neighboring $\mathbf {x}$ . Both variations comprise two stages:

Using a convolution kernel $\mathbf {K$ to compute the potential distribution $\mathbf {U} ^{t}(\mathbf {x} )=\mathbf {K} *\mathbf {A} ^{t}(\mathbf {x} )$ .
Using a growth mapping $G:[0,1]\rightarrow [-1,1]$ to compute the final growth distribution $\mathbf {G} ^{t}(\mathbf {x} )=G(\mathbf {U} ^{t}(\mathbf {x} ))$ .

Once $\mathbf {G} ^{t}$ is computed, it is scaled by the chosen time resolution $\Delta t$ and added to the original state value:

\mathbf {A} ^{t+\Delta t}(\mathbf {x} )={\text{clip}}(\mathbf {A} ^{t}+\Delta t\;\mathbf {G} ^{t}(\mathbf {x} ),\;0,\;1)

Here, the clip function is defined by $\operatorname {clip} (u,a,b):=\min(\max(u,a),b)$ .

The local rules are defined as follows for discrete and continuous Lenia:

{\begin{aligned}\mathbf {U} ^{t}(\mathbf {x} )&={\begin{cases}\sum _{\mathbf {n} \in {\mathcal {N}}}\mathbf {K(n)} \mathbf {A} ^{t}(\mathbf {x} +\mathbf {n} )\Delta x^{2},&{\text{discrete Lenia}}\\\int _{\mathbf {n} \in {\mathcal {N}}}\mathbf {K(n)} \mathbf {A} ^{t}(\mathbf {x} +\mathbf {n} )dx^{2},&{\text{continuous Lenia}}\end{cases}}\\\mathbf {G} ^{t}(\mathbf {x} )&=G(\mathbf {U} ^{t}(\mathbf {x} ))\\\mathbf {A} ^{t+\Delta t}(\mathbf {x} )&={\text{clip}}(\mathbf {A} ^{t}(\mathbf {x} )+\Delta t\;\mathbf {G} ^{t}(\mathbf {x} ),\;0,\;1)\end{aligned}}

Kernel generation

There are many ways to generate the convolution kernel $\mathbf {K}$ . The final kernel is the composition of a kernel shell $K_{C}$ and a kernel skeleton $K_{S}$ .

For the kernel shell $K_{C}$ , Chan gives several functions that are defined radially. Kernel shell functions are unimodal and subject to the constraint $K_{C}(0)=K_{C}(1)=0$ (and typically $K_{C}\left({\frac {1}{2}}\right)=1$ as well). Example kernel functions include:

K_{C}(r)={\begin{cases}\exp \left(\alpha -{\frac {\alpha }{4r(1-r)}}\right),&{\text{exponential}},\alpha =4\\(4r(1-r))^{\alpha },&{\text{polynomial}},\alpha =4\\\mathbf {1} _{\left[{\frac {1}{4}},{\frac {3}{4}}\right]}(r),&{\text{rectangular}}\\\ldots ,&{\text{etc.}}\end{cases}}

Here, $\mathbf {1} _{A}(r)$ is the indicator function.

Once the kernel shell has been defined, the kernel skeleton $K_{S}$ is used to expand it and compute the actual values of the kernel by transforming the shell into a series of concentric rings. The height of each ring is controlled by a kernel peak vector $\beta =(\beta _{1},\beta _{2},\ldots ,\beta _{B})\in [0,1]^{B}$ , where $B$ is the rank of the parameter vector. Then the kernel skeleton $K_{S}$ is defined as

K_{S}(r;\beta )=\beta _{\lfloor Br\rfloor }K_{C}(Br{\text{ mod }}1)

The final kernel $\mathbf {K} (\mathbf {n} )$ is therefore

\mathbf {K} (\mathbf {n} )={\frac {K_{S}(\lVert \mathbf {n} \rVert _{2})}{|K_{S}|}}

such that $\mathbf {K}$ is normalized to have an element sum of $1$ and $\mathbf {K} *\mathbf {A} \in [0,1]$ (for conservation of mass). $|K_{S}|=\textstyle \sum _{\mathcal {N}}\displaystyle K_{S}\,\Delta x^{2}$ in the discrete case, and $\int _{N}K_{S}\,dx^{2}$ in the continuous case.

Growth mappings

The growth mapping $G:[0,1]\rightarrow [-1,1]$ , which is analogous to an activation function, may be any function that is unimodal, nonmonotonic, and accepts parameters $\mu ,\sigma \in \mathbb {R}$ . Examples include

G(u;\mu ,\sigma )={\begin{cases}2\exp \left(-{\frac {(u-\mu )^{2}}{2\sigma ^{2}}}\right)-1,&{\text{exponential}}\\2\cdot \mathbf {1} _{[\mu \pm 3\sigma ]}(u)\left(1-{\frac {(u-\mu )^{2}}{9\sigma ^{2}}}\right)^{\alpha }-1,&{\text{polynomial}},\alpha =4\\2\cdot \mathbf {1} _{[\mu \pm \sigma ]}(u)-1,&{\text{rectangular}}\\\ldots ,&{\text{etc.}}\end{cases}}

where $u$ is a potential value drawn from $\mathbf {U} ^{t}$ .

Game of Life

The Game of Life may be regarded as a special case of discrete Lenia with $R=T=P=1$ . In this case, the kernel would be rectangular, with the function

K_{C}(r)=\mathbf {1} _{\left[{\frac {1}{4}},{\frac {3}{4}}\right]}(r)+{\frac {1}{2}}\mathbf {1} _{\left[0,{\frac {1}{4}}\right)}(r)

and the growth rule also rectangular, with $\mu =0.35,\sigma =0.07$ .

Patterns

By varying the convolutional kernel, the growth mapping and the initial condition, over 400 "species" of "life" have been discovered in Lenia, displaying "self-organization, self-repair, bilateral and radial symmetries, locomotive dynamics, and sometimes chaotic nature".^[7] Chan has created a taxonomy for these patterns.^[1]

Related work

Cellular automata as a convolutional neural network.

Other works have noted the strong similarity between cellular automata update rules and convolutions. Indeed, these works have focused on reproducing cellular automata using simplified convolutional neural networks. Mordvintsev et al. investigated the emergence of self-repairing pattern generation.^[9] Gilpin found that any cellular automaton could be represented as a convolutional neural network, and trained neural networks to reproduce existing cellular automata^[8]

In this light, cellular automata may be seen as a special case of recurrent convolutional neural networks. Lenia's update rule may also be seen as a single-layer convolution (the "potential field" $\mathbf {K}$ ) with an activation function (the "growth mapping" $G$ ). However, Lenia uses far larger, fixed, kernels and is not trained via gradient descent.

External links

Related Research Articles

In mathematics, convolution is a mathematical operation on two functions that produces a third function. The term convolution refers to both the result function and to the process of computing it. It is defined as the integral of the product of the two functions after one is reflected about the y-axis and shifted. The integral is evaluated for all values of shift, producing the convolution function. The choice of which function is reflected and shifted before the integral does not change the integral result. Graphically, it expresses how the 'shape' of one function is modified by the other.

In mathematical physics and mathematics, the Pauli matrices are a set of three $2 \times 2$ complex matrices that are Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.

In mathematical analysis, the Dirac delta function, also known as the unit impulse, is a generalized function on the real numbers, whose value is zero everywhere except at zero, and whose integral over the entire real line is equal to one. Since there is no function having this property, to model the delta "function" rigorously involves the use of limits or, as is common in mathematics, measure theory and the theory of distributions.

In probability theory, the central limit theorem (CLT) states that, under appropriate conditions, the distribution of a normalized version of the sample mean converges to a standard normal distribution. This holds even if the original variables themselves are not normally distributed. There are several versions of the CLT, each applying in the context of different conditions.

The stress–energy tensor, sometimes called the stress–energy–momentum tensor or the energy–momentum tensor, is a tensor physical quantity that describes the density and flux of energy and momentum in spacetime, generalizing the stress tensor of Newtonian physics. It is an attribute of matter, radiation, and non-gravitational force fields. This density and flux of energy and momentum are the sources of the gravitational field in the Einstein field equations of general relativity, just as mass density is the source of such a field in Newtonian gravity.

Noether's theorem states that every continuous symmetry of the action of a physical system with conservative forces has a corresponding conservation law. This is the first of two theorems proven by mathematician Emmy Noether in 1915 and published in 1918. The action of a physical system is the integral over time of a Lagrangian function, from which the system's behavior can be determined by the principle of least action. This theorem only applies to continuous and smooth symmetries of physical space.

In physics, a Langevin equation is a stochastic differential equation describing how a system evolves when subjected to a combination of deterministic and fluctuating ("random") forces. The dependent variables in a Langevin equation typically are collective (macroscopic) variables changing only slowly in comparison to the other (microscopic) variables of the system. The fast (microscopic) variables are responsible for the stochastic nature of the Langevin equation. One application is to Brownian motion, which models the fluctuating motion of a small particle in a fluid.

<span class="mw-page-title-main">Fokker–Planck equation</span> Partial differential equation

In statistical mechanics and information theory, the Fokker–Planck equation is a partial differential equation that describes the time evolution of the probability density function of the velocity of a particle under the influence of drag forces and random forces, as in Brownian motion. The equation can be generalized to other observables as well. The Fokker-Planck equation has multiple applications in information theory, graph theory, data science, finance, economics etc.

In mathematics, a Gaussian function, often simply referred to as a Gaussian, is a function of the base form

Linear elasticity is a mathematical model of how solid objects deform and become internally stressed due to prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

In mathematics, the Radon transform is the integral transform which takes a function f defined on the plane to a function Rf defined on the (two-dimensional) space of lines in the plane, whose value at a particular line is equal to the line integral of the function over that line. The transform was introduced in 1917 by Johann Radon, who also provided a formula for the inverse transform. Radon further included formulas for the transform in three dimensions, in which the integral is taken over planes. It was later generalized to higher-dimensional Euclidean spaces and more broadly in the context of integral geometry. The complex analogue of the Radon transform is known as the Penrose transform. The Radon transform is widely applicable to tomography, the creation of an image from the projection data associated with cross-sectional scans of an object.

In signal processing, cross-correlation is a measure of similarity of two series as a function of the displacement of one relative to the other. This is also known as a sliding dot product or sliding inner-product. It is commonly used for searching a long signal for a shorter, known feature. It has applications in pattern recognition, single particle analysis, electron tomography, averaging, cryptanalysis, and neurophysiology. The cross-correlation is similar in nature to the convolution of two functions. In an autocorrelation, which is the cross-correlation of a signal with itself, there will always be a peak at a lag of zero, and its size will be the signal energy.

In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.

The covariant formulation of classical electromagnetism refers to ways of writing the laws of classical electromagnetism in a form that is manifestly invariant under Lorentz transformations, in the formalism of special relativity using rectilinear inertial coordinate systems. These expressions both make it simple to prove that the laws of classical electromagnetism take the same form in any inertial coordinate system, and also provide a way to translate the fields and forces from one frame to another. However, this is not as general as Maxwell's equations in curved spacetime or non-rectilinear coordinate systems.

In physics, Maxwell's equations in curved spacetime govern the dynamics of the electromagnetic field in curved spacetime or where one uses an arbitrary coordinate system. These equations can be viewed as a generalization of the vacuum Maxwell's equations which are normally formulated in the local coordinates of flat spacetime. But because general relativity dictates that the presence of electromagnetic fields induce curvature in spacetime, Maxwell's equations in flat spacetime should be viewed as a convenient approximation.

In the fields of computer vision and image analysis, the Harris affine region detector belongs to the category of feature detection. Feature detection is a preprocessing step of several algorithms that rely on identifying characteristic points or interest points so to make correspondences between images, recognize textures, categorize objects or build panoramas.

In mathematical physics, the Belinfante–Rosenfeld tensor is a modification of the energy–momentum tensor that is constructed from the canonical energy–momentum tensor and the spin current so as to be symmetric yet still conserved.

Lagrangian field theory is a formalism in classical field theory. It is the field-theoretic analogue of Lagrangian mechanics. Lagrangian mechanics is used to analyze the motion of a system of discrete particles each with a finite number of degrees of freedom. Lagrangian field theory applies to continua and fields, which have an infinite number of degrees of freedom.

In machine learning, the kernel embedding of distributions comprises a class of nonparametric methods in which a probability distribution is represented as an element of a reproducing kernel Hilbert space (RKHS). A generalization of the individual data-point feature mapping done in classical kernel methods, the embedding of distributions into infinite-dimensional feature spaces can preserve all of the statistical features of arbitrary distributions, while allowing one to compare and manipulate distributions using Hilbert space operations such as inner products, distances, projections, linear transformations, and spectral analysis. This learning framework is very general and can be applied to distributions over any space $on which a sensible kernel function may be defined. For example, various kernels have been proposed for learning from data which are: vectors in, discrete classes/categories, strings, graphs/networks, images, time series, manifolds, dynamical systems, and other structured objects. The theory behind kernel embeddings of distributions has been primarily developed by Alex Smola, Le Song, Arthur Gretton, and Bernhard Schölkopf. A review of recent works on kernel embedding of distributions can be found in.$

The streamline upwind Petrov–Galerkin pressure-stabilizing Petrov–Galerkin formulation for incompressible Navier–Stokes equations can be used for finite element computations of high Reynolds number incompressible flow using equal order of finite element space by introducing additional stabilization terms in the Navier–Stokes Galerkin formulation.

References

1 2 3 Chan, Bert Wang-Chak (2019-10-15). "Lenia: Biology of Artificial Life". Complex Systems. 28 (3): 251–286. arXiv: 1812.05433 . doi:10.25088/ComplexSystems.28.3.251.
↑ "Lenia". chakazul.github.io. Retrieved 2021-10-12.
↑ Roberts, Siobhan (2020-12-28). "The Lasting Lessons of John Conway's Game of Life". The New York Times. ISSN 0362-4331 . Retrieved 2021-10-13.
↑ "The virtual creatures competition". virtualcreatures.github.io. Retrieved 2021-10-12.
↑ "ALife Art Award 2018". ALIFE Art Award 2018. Retrieved 2021-10-12.
↑ "2020 ISAL Awards: Winners".
↑ "Lenia". chakazul.github.io. Retrieved 2021-10-13.
1 2 Gilpin, William (2019-09-04). "Cellular automata as convolutional neural networks". Physical Review E. 100 (3): 032402. arXiv: 1809.02942 . doi: 10.1103/PhysRevE.100.032402 . ISSN 2470-0045.
↑ Mordvintsev, Alexander; Randazzo, Ettore; Niklasson, Eyvind; Levin, Michael (2020-02-11). "Growing Neural Cellular Automata". Distill. 5 (2): e23. doi: 10.23915/distill.00023 . ISSN 2476-0757.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[:0-1] 1 2 3 Chan, Bert Wang-Chak (2019-10-15). "Lenia: Biology of Artificial Life". Complex Systems. 28 (3): 251–286. arXiv: 1812.05433 . doi:10.25088/ComplexSystems.28.3.251.

[2] "Lenia". chakazul.github.io. Retrieved 2021-10-12.

[3] Roberts, Siobhan (2020-12-28). "The Lasting Lessons of John Conway's Game of Life". The New York Times. ISSN 0362-4331 . Retrieved 2021-10-13.

[4] "The virtual creatures competition". virtualcreatures.github.io. Retrieved 2021-10-12.

[5] "ALife Art Award 2018". ALIFE Art Award 2018. Retrieved 2021-10-12.

[6] "2020 ISAL Awards: Winners".

[7] "Lenia". chakazul.github.io. Retrieved 2021-10-13.

[Gilpin_032402-8] 1 2 Gilpin, William (2019-09-04). "Cellular automata as convolutional neural networks". Physical Review E. 100 (3): 032402. arXiv: 1809.02942 . doi: 10.1103/PhysRevE.100.032402 . ISSN 2470-0045.

[9] Mordvintsev, Alexander; Randazzo, Ettore; Niklasson, Eyvind; Levin, Michael (2020-02-11). "Growing Neural Cellular Automata". Distill. 5 (2): e23. doi: 10.23915/distill.00023 . ISSN 2476-0757.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[9]

[8]