1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
from sympy import *
q = 19
d = 3
c = [ Matrix(x) for x in [
( 0, 1, 1), (-1, 0, 1), ( 0, 0, 1), ( 1, 0, 1), ( 0, -1, 1),
(-1, 1, 0), ( 0, 1, 0), ( 1, 1, 0), (-1, 0, 0), ( 0, 0, 0), ( 1, 0, 0), (-1,-1, 0), ( 0, -1, 0), ( 1, -1, 0),
( 0, 1,-1), (-1, 0,-1), ( 0, 0,-1), ( 1, 0,-1), ( 0, -1,-1)
]]
w = [Rational(*x) for x in [
(1,36), (1,36), (1,18), (1,36), (1,36),
(1,36), (1,18), (1,36), (1,18), (1,3), (1,18), (1,36), (1,18), (1,36),
(1,36), (1,36), (1,18), (1,36), (1,36)
]]
c_s = sqrt(Rational(1,3))
|
