termcrit.m

来自「基于matlab的约束非线性规划算法库」· M 代码 · 共 19 行

function term = TermCrit(p,betak,f,epsrel,epsabs,cTc,epscon);
%Check if termination
%norm(c)^2 = cTc should be < epscon in order to check the termination criteria
%term is set = 1 if absolute termination is reached: 0.5*norm(f)<epsabs
%term is set = 10 if relative termination is reached: betak/norm(f)<releps
%
term = 0;
if isempty(p)
    term = -1003;
elseif cTc < epscon
    lenf = norm(f);
    if 0.5*lenf^2 <= epsabs
        term = 1;
    elseif betak <= lenf*epsrel
        term = term+10;
    else
        term = 0;
    end
end