python2022/hadamard.py at main · ksprski/python2022 · GitHub

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
46
47
48
49
50
import sys, math
size = int(sys.argv[1])

h = [[0]*size for i in range(size)]
h[0][0] = 1
l = int(math.log2(size))

for m in range(l):
	for i in range(2**m):
		for j in range(2**m):
			h[i][j+2**m] = h[i][j]
			h[i+2**m][j] = h[i][j]
			h[i+2**m][j+2**m] = -h[i][j]


t = [[0] * size for i in range(size)]
for i in range(size):
        for j in range(size):
                t[i][j] = h[j][i]


p = [[0] * size for i in range(size)]
for i in range(size):
        for j in range(size):
                for k in range(size):
                        p[i][j] += h[i][k] * t[k][j]


for i in range(size):
	s = ''
	for j in range(size):
		s += str(h[i][j]) + '\t'

	print(s,'\n','\n')

for i in range(size):
	y = ''
	for j in range(size):
		y += str(t[i][j]) + '\t'

	print(y,'\n')

for i in range(size):
	u = ''
	for j in range(size):
		u += str(p[i][j]) + '\t'

	print(u)