import sys

 import math

 import rules

 import gridmarshal as marshal

 from grid import Grid

 import argparse

 configs = {rules.Rule1: (((0, 0), 1), ((-1, 0), 1), ((1, 0), 1), ((0, 1), 1), 

 ((0, -1), 1)), 

 rules.Rule2: (((0, 0), 1), ((-1, 0), 1), ((1, 0), 1), ((0, 1), 1), ((0, -1), 1)),

 rules.PotentialGrowth: (((0, 0), 4),)}

 rulemap = {"1": rules.Rule1, "2": rules.Rule2, "3": rules.PotentialGrowth}

 def main(rule, iterations, dumpfile):

 	grid = Grid()

 	for pos, value in configs[type(rule)]:

 		grid.set(pos, value)

 	last_pi = 3.0

 	for i in range(iterations):

 		A = len(grid)

 		r = grid.width() / 4.0 + grid.height() / 4.0

 		r_ideal = math.sqrt(A / math.pi)

 		A_ideal = r**2 * math.pi

 		pi = A / r**2

 		pi = (pi + last_pi) / 2.0

 		last_pi = pi

 		print "%i %f(%f) %i(%i) %f" % (i, r, r - r_ideal, 

 										A, A - A_ideal, pi)

 		grid = rule.iterate(grid)

 		if dumpfile:

 			marshal.dumpGrid(grid, dumpfile)

 if __name__ == "__main__":

 	parser = argparse.ArgumentParser(description="cellular automaton.")

 	parser.add_argument("iterations", type=int, help="number of iterations")

 	parser.add_argument("--rule", "-r", default="3")

 	parser.add_argument("--persistent", "-p")

 	args = parser.parse_args()

 	main(rulemap[args.rule](), args.iterations, args.persistent)