krylov.py 文件源码

def __init__(self, n):
        self.weights = numpy.full(n, 2*sympy.pi/n)
        self.points = numpy.column_stack([
            [sympy.cos(sympy.pi * sympy.Rational(2*k, n)) for k in range(n)],
            [sympy.sin(sympy.pi * sympy.Rational(2*k, n)) for k in range(n)],
            ])
        self.degree = n - 1
        return