CELLULAR_AUTOMATON is a FORTRAN90 program which demonstrates the behavior of the 1D Cellular Automaton rule #30.
This program carries out iterations of the 1D cellular automaton known as rule 30.
Given an initial linear array of 0's and 1's, rule 30 produces a new array using transformations based on each value and the value of its left and right neighbors, as follows:
111 110 101 100 011 010 001 000 V V V V V V V V 0 0 0 1 1 1 1 0Note that there are 256 = 2^8 possible ways to fill in this output chart, and that rule 30 gets its index by the fact that (0,0,0,1,1,1,1,0) can be interpreted as the binary representation of 30.
For instance, if the current values of X(4), X(5) and X(6) are 0, 1 and 1, respectively, then the new value of X(5) will be 1.
The first and last entries of the array must be treated specially, since they don't have a left or right neighbor. One simple treatment is to assume that there are phantom neighbors whose values are both 0. Another is to enforce periodic boundary conditions.
CELLULAR_AUTOMATON is available in a C version and a C++ version and a FORTRAN77 version and a FORTRAN90 version and a MATLAB version.
You can go up one level to the FORTRAN90 source codes.