function [X_apx, iterations] = smooth(level, b, u_in)

% Ryan McKenzie
% Department of Computational Sciences
% University of Kentucky

amg_globals;

if SMOOTHER == JACOBI
    [X_apx, iterations] = jacobi_smoother(level, b, u_in);
elseif SMOOTHER == SOR
    [X_apx, iterations] = SOR_smoother(level, b, u_in);
else
    X_apx = u_in;
    iterations = 0;
end