(lvl, iterations=DEFAULT_NITER, sweep=DEFAULT_SWEEP,
                          omega=1.0)
return smoother
Set up Gauss-Seidel NR.
Set up Gauss-Seidel NR.
Set up Gauss - Seidel NR .
def setup_gauss_seidel_nr(lvl, iterations=DEFAULT_NITER, sweep=DEFAULT_SWEEP,
                          omega=1.0):
    """Set up Gauss-Seidel NR."""
    matrix_asformat(lvl, 'A', 'csc')

    def smoother(A, x, b):
        relaxation.gauss_seidel_nr(lvl.Acsc, x, b, iterations=iterations,
                                   sweep=sweep, omega=omega)
    return smoother
def setup_gauss_seidel_nr ( lvl , iterations = DEFAULT_NITER , sweep = DEFAULT_SWEEP , omega = 1.0 ) : matrix_asformat ( lvl , 'A' , 'csc' ) def smoother ( A , x , b ) : relaxation . gauss_seidel_nr ( lvl . Acsc , x , b , iterations = iterations , sweep = sweep , omega = omega ) return smoother