gauss_elim.m

来自「五点差分型多重网格方法：各种插值算子的比较）」· M 代码 · 共 14 行

function X_soln = gauss_elim(level, b)

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

amg_globals;

%Solves the system A(level) x_soln = b
%This command will use Gauss Elimination when appropriate 
%When not appropriate i.e. A(level) is near singular or badly scaled
%this command will use QR Decomposition with Pivoting

X_soln = A(level).matrix \ b;