plot_convergence.m

function plot_convergence(no_solvers,extras,test,problem, ...
                                fname)
if nargin < 5
    fname = 'xxx';
end


% setup for plot_cc
for i = 1:no_solvers
    res(i).resvec = extras{i,problem}.residual;%residual(:,i);
    res(i).description = test(i).descriptions;
    con_viol(i).resvec = extras{i,problem}.constraint_violation;% constraint_violation(:,i);
    con_viol(i).description = test(i).descriptions;
    Lag(i).resvec = extras{i,problem}.Lagrangian;%Lagrangian(:,i);
    Lag(i).description = test(i).descriptions;
    figure(i)
    semilogy(extras{i,problem}.primalfeas,'r-o')
    hold on, semilogy(extras{i,problem}.dualfeas,'k-x') 
    semilogy(extras{i,problem}.complementray,'b-.')
    %    semilogy(extras{i,problem}.mu.^0.5,'g-*')
    hold off
    xlabel('IP iteration number')    
    legend('Primal feasibility','Dual feasibility', ...
           'Complementarity')%,'\mu^{1/2}')
    saveas(i,['img/',fname,int2str(i),'.png'],'png');
end

%plot_cc('paper',2,'blash_minres_residuals',res)