2018-11-03 13:23:04 +00:00
|
|
|
#!/usr/bin/env python3
|
|
|
|
|
|
|
|
|
|
import math
|
2019-11-23 18:45:21 +00:00
|
|
|
import util
|
2018-11-03 13:23:04 +00:00
|
|
|
|
2019-10-29 14:11:35 +00:00
|
|
|
radius_x = 280//2
|
|
|
|
|
radius_y = 192//2
|
|
|
|
|
|
2018-11-03 13:23:04 +00:00
|
|
|
def f(x,r=36.3):
|
|
|
|
|
try:
|
|
|
|
|
return math.sqrt(r*r*(1.0-(x*x/(r*r*0.6))))
|
|
|
|
|
except:
|
|
|
|
|
return -1
|
|
|
|
|
|
2019-10-29 14:11:35 +00:00
|
|
|
coords = []
|
|
|
|
|
for i in range(30000, 0, -1):
|
|
|
|
|
a = float(i)/1000.0
|
2018-11-03 13:23:04 +00:00
|
|
|
b = f(a)
|
2019-10-29 14:11:35 +00:00
|
|
|
for x in range(140, 0, -1):
|
2018-11-03 13:23:04 +00:00
|
|
|
y = round(float(x)*b/a)
|
2019-10-29 14:11:35 +00:00
|
|
|
if x < 1 or x > radius_x or y < 1 or y > radius_y:
|
2018-11-03 13:23:04 +00:00
|
|
|
continue
|
|
|
|
|
for m in range(1, y+1):
|
2019-10-29 14:11:35 +00:00
|
|
|
if m % 2 != 0:
|
|
|
|
|
continue
|
|
|
|
|
coords.append((radius_x - x,radius_y - m))
|
|
|
|
|
|
2019-11-28 22:48:56 +00:00
|
|
|
unique_coords = util.unique(coords)
|
|
|
|
|
unique_vals = util.vals_1bit(unique_coords)
|
2019-10-29 14:11:35 +00:00
|
|
|
|
2019-11-29 00:25:00 +00:00
|
|
|
util.write("../../../src/fx/fx.hgr.radial.data.a", unique_vals)
|
|
|
|
|
util.write("../../../src/fx/fx.hgr.radial2.data.a", util.halfripple(unique_vals))
|
|
|
|
|
util.write("../../../src/fx/fx.hgr.radial4.data.a", util.ripple(unique_vals))
|
2019-11-29 03:10:25 +00:00
|
|
|
util.write("../../../src/fx/fx.hgr.radial5.data.a", util.ripple(util.halfripple(unique_vals)))
|
2019-11-28 22:48:56 +00:00
|
|
|
|
|
|
|
|
mult_coords = util.ripple(util.radial_multiply(util.ripple(util.radial_multiply(util.ripple(util.radial_multiply(unique_coords))))))
|
|
|
|
|
util.write("../../../src/fx/fx.hgr.pinwheels.data.a", util.vals_1bit(mult_coords))
|
