(lvl, iterations=DEFAULT_NITER, omega=1.0, withrho=True)
return smoother
Set up weighted-Jacobi.
Set up weighted-Jacobi.
Set up weighted - Jacobi .
def setup_jacobi(lvl, iterations=DEFAULT_NITER, omega=1.0, withrho=True):
    """Set up weighted-Jacobi."""
    if withrho:
        omega = omega/rho_D_inv_A(lvl.A)

    def smoother(A, x, b):
        relaxation.jacobi(A, x, b, iterations=iterations, omega=omega)
    return smoother
def setup_jacobi ( lvl , iterations = DEFAULT_NITER , omega = 1.0 , withrho = True ) : if withrho : omega = omega / rho_D_inv_A ( lvl . A ) def smoother ( A , x , b ) : relaxation . jacobi ( A , x , b , iterations = iterations , omega = omega ) return smoother