#!/usr/bin/env python3 import math def f(x,r=36.3): x = float(x) try: return math.sqrt(r*r*(1.0-(x*x/(r*r*0.6)))) except: return -1 coords = [(255,255)] for i in range(300, 0, -1): a = float(i)/10.0 any = False b = f(a) for x in range(20, 0, -1): y = round(float(x)*b/a) if y < 1 or y > 24: continue for m in range(1, y+1): if (x-1,m-1) not in coords: coords.append((x-1,m-1)) any = True if any: coords.append((255,255)) q2 = [(-x-1,y) for x,y in coords] q2.reverse() q3 = [(-x-1,-y-1) for x,y in coords] q4 = [(x,-y-1) for x,y in coords] q4.reverse() coords = q4 + coords + q2 + q3 with open("../../src/fx/fx.hgr.radial.data.a", "w") as f: for x, y in coords: if x not in range(-40,40): f.write(" !byte 255,255\n") else: f.write(" !byte %s,%s\n" % ((23-y,19-x)))