fontsize_labels = 14;
fontsize_grid   = 12;
fontname = 'Times';

h = zeros(5,1);
e = zeros(5,1);
hold off
clf
set(gca,'FontName',fontname)
set(gca,'FontSize',fontsize_grid)
set(gca,'FontSize',fontsize_labels)
for i=1:5
  E = 8*2^i;
  syscommand = compose('./spectraladjointassimilation  -ts_adapt_dt_max 3.e-3 -E %d -N 2 -ncoeff 5  -a .1 -tao_grtol 1.e-12  -tao_gatol 1.e-12  -tao_max_it 7',[E])
  [status,result] = system(syscommand{1});
  eval(result);
[m,n] = size(history);
  h(i) = history(m,2);
  e(i) = E;
  if (mod(i,2) == 1)
    yyaxis left
    semilogy(history(:,1),history(:,2),'Markersize',6,'LineWidth',2);
    ylabel('Analytic Error');
    yyaxis right
    semilogy(history(:,1),history(:,3),'Markersize',6,'LineWidth',2);
    ylabel('Objective function');
    hold on
  end
end
legend('16 elements','64 elements','256 elements','16 elements','64 elements','256 elements')
xlabel('Iteration')
print('convergencestudy-h','-depsc');

hold off
clf
set(gca,'FontName',fontname)
set(gca,'FontSize',fontsize_grid)
set(gca,'FontSize',fontsize_labels)
loglog(e,h,'-+','Markersize',6,'LineWidth',2);
xlabel('1/h');
ylabel('Analytic Error');
print('convergencestudy-h-2','-depsc');


clf
set(gca,'FontName',fontname)
set(gca,'FontSize',fontsize_grid)
set(gca,'FontSize',fontsize_labels)
h = zeros(4,1);
e = zeros(4,1);
for i=1:4
  N = 1 + 3*i;
  syscommand = compose('./spectraladjointassimilation  -ts_adapt_dt_max 3.e-3 -E 8 -N %d -ncoeff 12  -a .1 -tao_grtol 1.e-12  -tao_gatol 1.e-12  -tao_max_it 20',[N])
  [status,result] = system(syscommand{1});
  eval(result);
history
  e(i) = N;
  [m,n] = size(history);
  h(i) = history(m,2)
  yyaxis left
  semilogy(history(:,1),history(:,2),'Markersize',6,'LineWidth',2);
  ylabel('Analytic Error');
  yyaxis right
  semilogy(history(:,1),history(:,3),'Markersize',6,'LineWidth',2);
  ylabel('Objective function');
  hold on
end
legend('Order 4','Order 7','Order 10','Order 13','Order 4','Order 7','Order 10','Order 13')
xlabel('Iteration')
print('convergencestudy-p','-depsc');

hold off
clf
set(gca,'FontName',fontname)
set(gca,'FontSize',fontsize_grid)
set(gca,'FontSize',fontsize_labels)
loglog(e,h,'-+','Markersize',6,'LineWidth',2);
xlabel('Polynomial order');
ylabel('Analytic Error');
print('convergencestudy-p-2','-depsc');