4 def buildmats(qref, qweight, mat_dtype="float64"):  argument
9     qref[0] = 0.2
10     qref[1] = 0.6
11     qref[2] = 1. / 3.
12     qref[3] = 0.2
13     qref[4] = 0.2
14     qref[5] = 0.2
15     qref[6] = 1. / 3.
16     qref[7] = 0.6
24         x1 = qref[0 * Q + i]
25         x2 = qref[1 * Q + i]
52 def buildmatshdiv(qref, qweight, mat_dtype="float64"):  argument
57     qref[0] = -1. / np.sqrt(3.)
58     qref[1] = qref[0]
59     qref[2] = qref[0]
60     qref[3] = -qref[0]
61     qref[4] = -qref[0]
62     qref[5] = -qref[0]
63     qref[6] = qref[0]
64     qref[7] = qref[0]
72         x1 = qref[0 * Q + i]
73         x2 = qref[1 * Q + i]
92 def buildmatshcurl(qref, qweight, mat_dtype="float64"):  argument
97     qref[0] = 0.2
98     qref[1] = 0.6
99     qref[2] = 1. / 3.
100     qref[3] = 0.2
101     qref[4] = 0.2
102     qref[5] = 0.2
103     qref[6] = 1. / 3.
104     qref[7] = 0.6
112         x1 = qref[0 * Q + i]
113         x2 = qref[1 * Q + i]