percolation-model/PercolationModelPatterns.py at master · DuncanForgan/percolation-model · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
# Written 7/7/14 by dh4gan

# supplies patterns to be added to a cellular automaton obeying
# Conway's Game of Life Rules

from PercolationModel import PercolationModel2D


def add_block(cell, icentre,jcentre):
    '''
    Adds a 2x2 block into the system, with bottom left corner (icentre, jcentre)
    '''

    extent = 2
    cell.clear(icentre,jcentre,extent)
    cell.grid[icentre:icentre+2,jcentre:jcentre+2] = 1


def add_beehive(cell, icentre,jcentre):
    '''
    Adds a beehive into the system, with (icentre,jcentre) being the inner left blank square
    '''

    extent = 7
    cell.clear(icentre,jcentre,extent)

    cell.grid[icentre-1,jcentre:jcentre+2] = 1
    cell.grid[icentre+1,jcentre:jcentre+2] = 1

    cell.grid[icentre,jcentre-1] = 1
    cell.grid[icentre,jcentre+2] = 1

def add_blinker(cell,icentre,jcentre):
    '''
    Adds a vertical line of 3 blocks (period 1)
    '''

    extent = 4
    cell.clear(icentre,jcentre,extent)
    cell.grid[icentre,jcentre-1:jcentre+2] = 1

def add_loaf(cell,icentre, jcentre):
    '''
    Adds a loaf, with (icentre,jcentre) in the bottom left corner (blank)
    '''

    cell.grid[icentre,jcentre+1:jcentre+3] = 1
    cell.grid[icentre+1:icentre+3,jcentre+3] = 1
    cell.grid[icentre+1,jcentre] = 1
    cell.grid[icentre+2,jcentre+1] = 1
    cell.grid[icentre+3,jcentre+2] = 1

def add_boat(cell,icentre, jcentre):
    '''
    '''
    extent = 4
    cell.clear(icentre,jcentre,extent)

    indices = cell.getVonNeumannNeighbourhood(icentre,jcentre)
    print indices
    for element in indices:
        cell.grid[element[0],element[1]] = 1

    cell.grid[icentre+1,jcentre-1] = 1

def add_toad(cell,icentre,jcentre):
    extent = 3
    cell.clear(icentre,jcentre,extent)

    cell.grid[icentre-1:icentre+2, jcentre]=1
    cell.grid[icentre:icentre+3, jcentre-1]=1


def add_beacon(cell,icentre,jcentre):
    '''
    Adds two 2x2 blocks, which repeat a pattern of period 2
    '''

    extent = 3
    cell.clear(icentre,jcentre,extent)

    add_block(cell,icentre+2,jcentre)
    add_block(cell,icentre,jcentre+2)


def add_pulsar(cell,icentre,jcentre):
    '''
    Adds a pulsar, a period 3 oscillator
    '''
    extent = 8
    cell.clear(icentre,jcentre,extent)

    # Start with inner cross

    # North
    cell.grid[icentre+2:icentre+5,jcentre+1] = 1
    cell.grid[icentre+2:icentre+5,jcentre-1] = 1

    # South
    cell.grid[icentre-4:icentre-1,jcentre+1] = 1
    cell.grid[icentre-4:icentre-1,jcentre-1] = 1

    # East
    cell.grid[icentre+1,jcentre+2:jcentre+5] = 1
    cell.grid[icentre-1,jcentre+2:jcentre+5] = 1

    # West
    cell.grid[icentre+1,jcentre-4:jcentre-1] = 1
    cell.grid[icentre-1,jcentre-4:jcentre-1] = 1


    # Now do surrounding bars - quadrant at a time

    cell.grid[icentre+6,jcentre+2:jcentre+5] = 1
    cell.grid[icentre+2:icentre+5, jcentre+6] = 1

    cell.grid[icentre-4:icentre-1, jcentre+6] = 1
    cell.grid[icentre-6, jcentre+2:jcentre+5] = 1

    cell.grid[icentre-6, jcentre-4:jcentre-1] = 1
    cell.grid[icentre-4:icentre-1, jcentre-6] = 1

    cell.grid[icentre+2: icentre+5, jcentre-6] = 1
    cell.grid[icentre+6, jcentre-4:jcentre-1] = 1

def add_glider(cell,icentre,jcentre):
    '''
    Adds a glider, with (icentre,jcentre) being the bottom left tile (alive)
    '''

    cell.grid[icentre,jcentre:jcentre+3] = 1
    cell.grid[icentre+1,jcentre+2] = 1
    cell.grid[icentre+2,jcentre+1] = 1

def add_spaceship(cell,icentre,jcentre):
    '''
    Adds a lightweight spaceship, with (icentre,jcentre) being the bottom left tile (alive)
    '''
    cell.grid[icentre,jcentre:jcentre+4]=1
    cell.grid[icentre:icentre+3,jcentre]=1
    cell.grid[icentre+3,jcentre+1] = 1
    cell.grid[icentre+3,jcentre+4] = 1
    cell.grid[icentre+1,jcentre+4] = 1

def add_glider_gun(cell,icentre,jcentre):
    '''
    Adds a Gosper glider gun - two blocks with patterns inbetween to set up the gun
    '''
    add_block(cell,icentre+5,jcentre+2)

    # Make first of inner patterns
    cell.grid[icentre+4:icentre+7,jcentre+12] = 1
    cell.grid[icentre+3,jcentre+13] =cell.grid[icentre+7,jcentre+13] = 1
    #cell.grid[icentre+2,jcentre+14] =cell.grid[icentre+8,jcentre+14] = 1
    cell.grid[icentre+2,jcentre+14:jcentre+16] =cell.grid[icentre+8,jcentre+14:jcentre+16] = 1
    cell.grid[icentre+5,jcentre+16] =1
    cell.grid[icentre+3,jcentre+17] =cell.grid[icentre+7,jcentre+17] = 1
    cell.grid[icentre+4:icentre+7,jcentre+18] = 1
    cell.grid[icentre+5,jcentre+19] =1

    # Now second pattern
    cell.grid[icentre+6:icentre+9, jcentre+22:jcentre+24] = 1
    cell.grid[icentre+5,jcentre+24] =cell.grid[icentre+9,jcentre+24] = 1
    cell.grid[icentre+4:icentre+6,jcentre+26] = cell.grid[icentre+9:icentre+11,jcentre+26] = 1


    add_block(cell,icentre+7,jcentre+36)