sawine@3
|
1 |
import copy
|
sawine@3
|
2 |
from caugrid import Grid
|
sawine@3
|
3 |
|
sawine@3
|
4 |
class Rule(object):
|
sawine@3
|
5 |
|
sawine@3
|
6 |
def iterate(self, oldgrid):
|
sawine@3
|
7 |
grid = copy.deepcopy(oldgrid)
|
sawine@3
|
8 |
for x in xrange(oldgrid.minx - 1, oldgrid.maxx + 2):
|
sawine@3
|
9 |
for y in xrange(oldgrid.miny - 1, oldgrid.maxy + 2):
|
sawine@3
|
10 |
#print x, y, self.neighbours(oldgrid, x, y)
|
sawine@3
|
11 |
n = self.neighbours(oldgrid, (x, y))
|
sawine@3
|
12 |
if n > 2:
|
sawine@3
|
13 |
grid.set((x, y))
|
sawine@3
|
14 |
return grid
|
sawine@3
|
15 |
|
sawine@3
|
16 |
def neighbours(self, grid, (testx, testy)):
|
sawine@3
|
17 |
n = 0
|
sawine@3
|
18 |
#print "testing ", testx, testy
|
sawine@3
|
19 |
for x in range(testx - 1, testx + 2):
|
sawine@3
|
20 |
for y in range(testy - 1, testy + 2):
|
sawine@3
|
21 |
if (x, y) in grid.cells:
|
sawine@3
|
22 |
n += 1
|
sawine@3
|
23 |
#print x, y, n
|
sawine@3
|
24 |
return n
|
sawine@3
|
25 |
|
sawine@3
|
26 |
|
sawine@3
|
27 |
class Rule2(object):
|
sawine@3
|
28 |
|
sawine@3
|
29 |
def iterate(self, oldgrid):
|
sawine@3
|
30 |
grid = copy.deepcopy(oldgrid)
|
sawine@3
|
31 |
for x in xrange(oldgrid.minx - 1, oldgrid.maxx + 2):
|
sawine@3
|
32 |
for y in xrange(oldgrid.miny, oldgrid.maxy + 2):
|
sawine@3
|
33 |
#print "testing ", x, y,
|
sawine@3
|
34 |
if (x+1, y) in oldgrid.cells or (x, y-1) in oldgrid.cells:
|
sawine@3
|
35 |
grid.set((x, y), 1)
|
sawine@3
|
36 |
return grid
|
sawine@3
|
37 |
|
sawine@3
|
38 |
import random
|
sawine@3
|
39 |
import caumarshal as marshal
|
sawine@3
|
40 |
|
sawine@3
|
41 |
class PotentialGrowth(object):
|
sawine@3
|
42 |
|
sawine@3
|
43 |
def __init__(self):
|
sawine@3
|
44 |
random.seed()
|
sawine@3
|
45 |
|
sawine@3
|
46 |
def iterate(self, oldgrid):
|
sawine@3
|
47 |
grid = copy.deepcopy(oldgrid)
|
sawine@3
|
48 |
for cell in oldgrid.cells.iteritems():
|
sawine@3
|
49 |
pos = cell[0]
|
sawine@3
|
50 |
value = cell[1]
|
sawine@3
|
51 |
if value > 1:
|
sawine@3
|
52 |
new_pos, new_value = self.grow(grid, pos)
|
sawine@3
|
53 |
grid.set(new_pos, new_value)
|
sawine@3
|
54 |
grid.set(pos, value - 1)
|
sawine@3
|
55 |
return grid
|
sawine@3
|
56 |
|
sawine@3
|
57 |
def grow(self, grid, (x, y)):
|
sawine@3
|
58 |
neighbours = [(x-1, y), (x, y+1), (x+1, y), (x, y-1)]
|
sawine@3
|
59 |
neighbours = [n for n in neighbours if n not in grid.cells]
|
sawine@3
|
60 |
try:
|
sawine@3
|
61 |
pos = random.choice(neighbours)
|
sawine@3
|
62 |
value = 4
|
sawine@3
|
63 |
except IndexError:
|
sawine@3
|
64 |
pos = (x, y)
|
sawine@3
|
65 |
value = grid.cells[pos]
|
sawine@3
|
66 |
return pos, value
|
sawine@3
|
67 |
|
sawine@3
|
68 |
import sys
|
sawine@3
|
69 |
import math
|
sawine@3
|
70 |
|
sawine@3
|
71 |
def ruleTest():
|
sawine@3
|
72 |
rule = Rule()
|
sawine@3
|
73 |
grid = Grid()
|
sawine@3
|
74 |
grid.set((0, 0))
|
sawine@3
|
75 |
grid.set((0, 1)).set((1, 0)).set((0, -1)).set((-1, 0))
|
sawine@3
|
76 |
|
sawine@3
|
77 |
|
sawine@3
|
78 |
print "iteration radius(diff) area(diff) pi"
|
sawine@3
|
79 |
|
sawine@3
|
80 |
iterations = int(sys.argv[1])
|
sawine@3
|
81 |
|
sawine@3
|
82 |
for i in range(iterations):
|
sawine@3
|
83 |
A = len(grid)
|
sawine@3
|
84 |
r = grid.width() / 2.0
|
sawine@3
|
85 |
r_ideal = math.sqrt(A / math.pi)
|
sawine@3
|
86 |
A_ideal = r**2 * math.pi
|
sawine@3
|
87 |
pi = A / r**2
|
sawine@3
|
88 |
print "%i %f(%f) %i(%i) %f" % (i, r, r - r_ideal,
|
sawine@3
|
89 |
A, A - A_ideal, pi)
|
sawine@3
|
90 |
|
sawine@3
|
91 |
grid = rule.iterate(grid)
|
sawine@3
|
92 |
marshal.dumpGrid(grid, "grid.cfg")
|
sawine@3
|
93 |
#print
|
sawine@3
|
94 |
|
sawine@3
|
95 |
def rule2Test():
|
sawine@3
|
96 |
rule = Rule2()
|
sawine@3
|
97 |
grid = Grid()
|
sawine@3
|
98 |
grid.set((0, 0)).set((-1, 0)).set((0, 1))
|
sawine@3
|
99 |
olda = 1
|
sawine@3
|
100 |
|
sawine@3
|
101 |
iterations = int(sys.argv[1])
|
sawine@3
|
102 |
|
sawine@3
|
103 |
for i in range(iterations):
|
sawine@3
|
104 |
a = (float(len(grid)) - grid.width()) * 4.0
|
sawine@3
|
105 |
pi = 1.0 / ((grid.width()-1)**2 / ((a + olda) / 2.0))
|
sawine@3
|
106 |
print i, grid.width(), len(grid), pi
|
sawine@3
|
107 |
grid = rule.iterate(grid)
|
sawine@3
|
108 |
olda = a
|
sawine@3
|
109 |
marshal.dumpGrid(grid, "grid.cfg")
|
sawine@3
|
110 |
|
sawine@3
|
111 |
def potentialTest():
|
sawine@3
|
112 |
rule = PotentialGrowth()
|
sawine@3
|
113 |
grid = Grid()
|
sawine@3
|
114 |
grid.set((0, 0), 4)
|
sawine@3
|
115 |
|
sawine@3
|
116 |
iterations = int(sys.argv[1])
|
sawine@3
|
117 |
|
sawine@3
|
118 |
for i in range(iterations):
|
sawine@3
|
119 |
A = len(grid)
|
sawine@3
|
120 |
r = grid.width() / 4.0 + grid.height() / 4.0
|
sawine@3
|
121 |
r_ideal = math.sqrt(A / math.pi)
|
sawine@3
|
122 |
A_ideal = r**2 * math.pi
|
sawine@3
|
123 |
pi = A / r**2
|
sawine@3
|
124 |
print "%i %f(%f) %i(%i) %f" % (i, r, r - r_ideal,
|
sawine@3
|
125 |
A, A - A_ideal, pi)
|
sawine@3
|
126 |
|
sawine@3
|
127 |
grid = rule.iterate(grid)
|
sawine@3
|
128 |
marshal.dumpGrid(grid, "grid.cfg")
|
sawine@3
|
129 |
|
sawine@3
|
130 |
def main():
|
sawine@3
|
131 |
#ruleTest()
|
sawine@3
|
132 |
#rule2Test()
|
sawine@3
|
133 |
potentialTest()
|
sawine@3
|
134 |
|
sawine@3
|
135 |
if __name__ == "__main__":
|
sawine@3
|
136 |
main()
|