Homec4science
Diffusion Nonlinear (master) Browse

R7172/codes/helpers/plot_results_lbfgs.m
0b40aa10684emaster

Last Change
Raw File

plot_results_lbfgs.m
View Options

% Plotting function for lbfgs
plot.grad_mapping = false;
font_size = 20;
lw = 2;
if plot.grad_mapping
figure('Position',[100,500,1140,270]);
subplot(1,3,1)
else
figure('Position',[100,500,760,270]);
subplot(1,2,1)
end
loglog(output_lbfgs.iter, output_lbfgs.feas,'--r', 'Linewidth', lw),
hold on,
loglog(output_apgm.iter, output_apgm.feas,'--g', 'Linewidth', lw),
loglog(output.iter,output.feas,'--k', 'Linewidth', lw),
grid on
title('Feasibility Gap','Interpreter', 'latex', 'FontSize', font_size);
xlabel('iters', 'Interpreter', 'latex', 'FontSize', font_size);
ylabel('$||A(u)-b||$', 'Interpreter', 'latex', 'FontSize', font_size);
ax = gca;
ax.XTick = 10.^(-20:10);
legend({'AL lBFGS', 'AL APGM', 'LinAL'}, 'Interpreter', 'Latex', 'FontSize', 16, 'Location', 'best')
if plot.grad_mapping
subplot(1,3,2)
loglog(output_lbfgs.iter,output_lbfgs.grad_mapping, '--r', 'Linewidth', lw)
hold on,
loglog(output_apgm.iter,output_apgm.grad_mapping, '--g', 'Linewidth', lw)
loglog(output.iter,output.grad_mapping, '--k', 'Linewidth', lw)
grid on
title('Gradient Mapping','Interpreter', 'latex', 'FontSize', font_size);
xlabel('iters', 'Interpreter', 'latex', 'FontSize', font_size);
ylabel('$|G(\gamma_k, u_k)|^2$', 'Interpreter', 'latex', 'FontSize', font_size);
ax = gca;
ax.XTick = 10.^(-20:10);
end
if plot.grad_mapping
subplot(1,3,3)
else
subplot(1,2,2)
end
% loglog(output_lbfgs.iter, abs(output_lbfgs.obj-opt_val)/abs(opt_val), '--r', 'Linewidth', lw)
% hold on,
% loglog(abs(output.obj-opt_val)/abs(opt_val), '--k', 'Linewidth', lw)
loglog(output_lbfgs.iter, abs(output_lbfgs.obj-opt_val), '--r', 'Linewidth', lw);%/abs(opt_val)
hold on,
loglog(output_apgm.iter, abs(output_apgm.obj-opt_val), '--g', 'Linewidth', lw);%/abs(opt_val)
loglog(output.iter, abs(output.obj-opt_val), '--k', 'Linewidth', lw);%/abs(opt_val)
grid on
title('Relative Objective Residual','Interpreter', 'latex', 'FontSize', font_size);
xlabel('iters', 'Interpreter', 'latex', 'FontSize', font_size);
ylabel('${ |f(u)-f^\star| }$', 'Interpreter', 'latex', 'FontSize', font_size);
ax = gca;
ax.XTick = 10.^(-20:10);
%%
% figure('Position',[100,100,494,351]),
% loglog(output_lbfgs.iter, output_lbfgs.beta,'r','Linewidth', lw);grid on;hold on;loglog(output_lbfgs.iter, output_lbfgs.sigma, 'k','Linewidth', lw)
% legend({'$\beta_{k}$', '$\sigma_{k}$'}, 'Interpreter', 'Latex', 'FontSize', font_size, 'Location', 'best')
% xlabel('iters', 'Interpreter', 'latex', 'FontSize', font_size);
% title('Parameter Updates','Interpreter', 'latex', 'FontSize', font_size);
% in case we want to color different update regions. Use function jbfill
% from mathworks file exchange
% roi = y0< y1;
% jbfill(x,y0,lower,roi,'r','none',1,.5);
% roi = (y0 > y1) .* (y0 < y2);
% jbfill(x,y0,lower,roi,'k','none',1,.5);
% roi = y0> y2;
% jbfill(x,y0,lower,roi,'m','none',1,.5);
c4science · Help