amg_setup.m

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

%AMG_SETUP calls the specified routines for coarse point and intergrid weight selection

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

function amg_setup(level);

amg_globals;

if level ~= COARSEST %if the current level is not the coarsest level
    %call appropriate setup algorithm to select coarse points and get
    %interpolation weights
    
    if SETUP_ALG == RUGE_STUEBEN
        RugeStuebenCoarsen(level);
    elseif SETUP_ALG == AMGm
        AMGmCoarsen(level);
    elseif SETUP_ALG == SMOOTH_AGGREGATE
        SmoothAggregateCoarsen(level);
    end
    
    if SETUP_OPT == AT_ONCE %if we are coarsening more than one level
        amg_setup(level+1); %call setup for next coarse grid
    end
end