thacher.py 文件源码

def __init__(self, n):
        self.degree = 2
        r = sqrt(3) / 6
        data = [
            (1.0, [n * [2*r]]),
            (+r, _s(n, -1, r)),
            (-r, _s(n, +1, r)),
            ]

        self.points, self.weights = untangle(data)
        reference_volume = 2**n
        self.weights *= reference_volume
        return