coarse_solve.m

五点差分型多重网格方法：各种插值算子的比较）

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