(lvl, iterations=DEFAULT_NITER, omega=1.0, withrho=True)
return smoother
Set up Jacobi NE.
Set up Jacobi NE.
Set up Jacobi NE .
def setup_jacobi_ne(lvl, iterations=DEFAULT_NITER, omega=1.0, withrho=True):
    """Set up Jacobi NE."""
    matrix_asformat(lvl, 'A', 'csr')
    if withrho:
        omega = omega/rho_D_inv_A(lvl.Acsr)**2

    def smoother(A, x, b):
        relaxation.jacobi_ne(lvl.Acsr, x, b, iterations=iterations,
                             omega=omega)
    return smoother
def setup_jacobi_ne ( lvl , iterations = DEFAULT_NITER , omega = 1.0 , withrho = True ) : matrix_asformat ( lvl , 'A' , 'csr' ) if withrho : omega = omega / rho_D_inv_A ( lvl . Acsr ) ** 2 def smoother ( A , x , b ) : relaxation . jacobi_ne ( lvl . Acsr , x , b , iterations = iterations , omega = omega ) return smoother