hammer_marlowe_stroud.py 文件源码

def __init__(self, index):
        if index == 1:
            self.degree = 2
            data = [
                (fr(1, 4), _r(1/sqrt(5))),
                ]
        elif index == 2:
            self.degree = 2
            data = [
                (fr(1, 4), _r(-1/sqrt(5))),
                ]
        else:
            assert index == 3
            self.degree = 3
            data = [
                (-fr(4, 5), _s4()),
                (fr(9, 20), _r(fr(1, 3))),
                ]

        self.bary, self.weights = untangle(data)
        self.points = self.bary[:, 1:]
        return