-
Notifications
You must be signed in to change notification settings - Fork 1
/
macwilliams.py
45 lines (29 loc) · 980 Bytes
/
macwilliams.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.colors as colors
import functools
import matplotlib.colors as colors
@functools.cache
def MW(n, p):
if n == 0:
return np.array([1])
mw = np.array(
[[1, 1],
[1, -1]]
)
for i in range(1, n):
mw = np.pad(mw, [(0, 1), (0, 1)]) + np.pad(mw, [(1, 0), (0, 1)]) + np.pad(mw, [(0, 1), (1, 0)]) - np.pad(mw, [(1, 0), (1, 0)])
mw[:, 1:-1] = mw[:, 1:-1] * (p+1)//2
mw = np.mod(mw, p)
return np.mod(mw, p)
def hrbr(m, p):
hr = lambda x: 0.77 if 0 <= x < 0.5 else 0.414
br = lambda x: 5 * x + 0.5 if 0 <= x < 0.1 else -5 * x + 5.5 if 0.9 < x <= 1 else 1
return np.round(255 * colors.hsv_to_rgb((hr(m/p) * (1 - m/p), 1, br(m / p))))
def MWplot(n, p):
mw = MW(n, p)
hsv_mw = np.zeros((n+1, n+1, 3))
for i in range(n+1):
for j in range(n+1):
hsv_mw[i][j] = hrbr(mw[i][j], p)
return hsv_mw