newton_time_comparison.m

csample = 0.4:0.2:1;

Ms = [];
Ns = [];

crbtimes = [];
rbtimes  = [];
sizes    = [];
maxerrs  = [];

for c=1:length(csample)
  newdir = fullfile(rbmatlabtemp, ['time_run_fixed', num2str(c)])
  mkdir(newdir);
  outfile = newton_steps(0);
  load(outfile);
  model.Mmax         = ceil(csample(c)*model.MMax);
  model.RB_stop_Nmax = ceil(csample(c)*model.RB_stop_Nmax);
  Ms = [Ms, model.Mmax]; Ns = [Ns, model.RB_stop_Nmax];
  save(outfile, 'model', 'model_data');

  newton_steps(2);
  outfile = newton_steps(5);
  tmp = load(outfile);
  crbtimes = [crbtimes, tmp.detailed_data.ei_info{1}.elapsed_time]
  rbtimes  = [rbtimes,  tmp.detailed_data.RB_info.elapsed_time]
  sizes    = [sizes,    [size(tmp.detailed_data.RB,2);...
                         size(tmp.detailed_data.QM{1},2)] ];
  copyfile(outfile, fullfile(newdir, 'detailed'));
  clear 'tmp';
  outfile = newton_steps(7);
  copyfile(outfile, fullfile(newdir, 'errs'));
  tmp = load(outfile);
  maxerrs = [maxerrs, tmp.output{1}.errs(1)];
end

for c=1:length(csample)
  newdir = fullfile(rbmatlabtemp, ['time_run_adapted', num2str(c)])
  mkdir(newdir);
  outfile = newton_steps(0);
  load(outfile);
  model.Mmax         = ceil(csample(c)*model.MMax);
  model.RB_stop_Nmax = ceil(csample(c)*model.RB_stop_Nmax);
  Ms = [Ms, model.Mmax]; Ns = [Ns, model.RB_stop_Nmax];

  model.RB_generation_mode = 'greedy_refined';
  model.RB_numintervals    = ones(size(model.RB_numintervals));
  save(outfile, 'model', 'model_data');

  newton_steps(2);
  outfile = newton_steps(5);
  tmp = load(outfile);
  crbtimes = [crbtimes, tmp.detailed_data.ei_info{1}.elapsed_time]
  rbtimes  = [rbtimes,  tmp.detailed_data.RB_info.elapsed_time]
  sizes    = [sizes,    [size(tmp.detailed_data.RB,2);...
                         size(tmp.detailed_data.QM{1},2)] ];
  copyfile(outfile, fullfile(newdir, 'detailed'));
  clear 'tmp';
  outfile = newton_steps(7);
  copyfile(outfile, fullfile(newdir, 'errs'));
  tmp = load(outfile);
  maxerrs = [maxerrs, tmp.output{1}.errs(1)];
end