function X_soln = coarse_solve(level, b)

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

amg_globals;

initGuess = zeros(size(b));
if C_SOLVER == THE_SMOOTHER
    X_soln = smooth(level, b, initGuess);
elseif C_SOLVER == DIRECT_ELIM
    X_soln = gauss_elim(level, b);
else
    X_soln = initGuess;
end