2019-10-26 19:55:06 +00:00
|
|
|
#!/usr/bin/env python3
|
|
|
|
|
|
|
|
|
|
from math import sqrt, sin, cos, acos, pi
|
2019-11-29 00:25:00 +00:00
|
|
|
import util
|
2019-10-26 19:55:06 +00:00
|
|
|
|
|
|
|
|
max_x = 280//2
|
|
|
|
|
max_y = 192//2
|
|
|
|
|
|
|
|
|
|
def f(t, k):
|
|
|
|
|
t = float(t)
|
|
|
|
|
r = k/cos(0.4*acos(sin(2.5*(t+pi/2))))
|
|
|
|
|
return r*cos(t),r*sin(t)
|
|
|
|
|
|
|
|
|
|
coords = []
|
|
|
|
|
for k_mul in range(500):
|
|
|
|
|
for t_mul in range(int(pi*1000+1)):
|
|
|
|
|
a, b = f(float(t_mul/100), float(k_mul)/10.0)
|
|
|
|
|
x = round(max_x//2+a*1.2)
|
|
|
|
|
y = round(max_y//2+b)
|
|
|
|
|
if y % 2 != 0:
|
|
|
|
|
continue
|
|
|
|
|
if x < 0 or x >= max_x or y < 0 or y >= max_y:
|
|
|
|
|
continue
|
|
|
|
|
coords.append((x,y))
|
|
|
|
|
|
2019-11-29 00:25:00 +00:00
|
|
|
unique_coords = util.unique(coords)
|
|
|
|
|
unique_vals = util.vals_1bit(unique_coords)
|
|
|
|
|
|
|
|
|
|
util.write("../../../src/fx/fx.hgr.slow.star.data.a", unique_vals)
|
2019-10-26 20:10:35 +00:00
|
|
|
|
|
|
|
|
unique_vals.reverse()
|
2019-11-29 00:25:00 +00:00
|
|
|
util.write("../../../src/fx/fx.hgr.slow.star.in.data.a", unique_vals)
|
|
|
|
|
|
|
|
|
|
mult_coords = util.ripple(util.radial_multiply(unique_coords))
|
|
|
|
|
util.write("../../../src/fx/fx.hgr.full.of.stars.data.a", util.vals_1bit(mult_coords))
|
