Method chebyshev

engine/solver/amg_python.py:51–60

(self, A, x, b, iterations=1)

    # https://github.com/pyamg/pyamg/blob/5a51432782c8f96f796d7ae35ecc48f81b194433/pyamg/relaxation/relaxation.py#L586
    def chebyshev(self, A, x, b, iterations=1):
        coefficients = self.chebyshev_coeff
        x = np.ravel(x)
        b = np.ravel(b)
        for _i in range(iterations):
            residual = b - A*x
            h = coefficients[0]*residual
            for c in coefficients[1:]:
                h = c*residual + A*h
            x += h

    def calc_spectral_radius(self, A):
        t = time.perf_counter()