adaptive_basisgen_admos2007.m

Go to the documentation of this file.

% script for generating the MCMDS basis generation images

% Bernard Haasdonk 22.2.2010

% a adaptive training set extension:
step = 1; % generate different bases      
% step = 2; % generate images overfitting*_small.eps 
% step = 3; % generate diagrams of training runtimes and errors
% step = 4; % generate diagrams of test errors and error ratios
% step = 5; % generate diagrams of max test errors

switch step
 case 1
  disp('to be implemented!!');
 case 2 
  figure;
  load(['C:\Dokumente und Einstellungen\Bernard\Eigene Dateien\',...
        'TempMatlab\basisgen\RB_uniform_fixed_64.mat']);
  params.num_plot_intervals = [40,40];
  offline_data = rb_offline_prep(detailed_data,params);
  params.N = size(detailed_data.RB,2);
  online_data = rb_online_prep(offline_data,params);
  params.no_lines = 1;
  plot_error_estimator(online_data,params);
  c = get(gca,'Children');
  cdata = get(c(end),'Cdata');
  cdata = log10(cdata);
  set(c(end),'Cdata',cdata);
  set(gca,'Clim',[-8,-7]);
  title('log_{10}(\Delta(\mu))');
  colorbar;
  %colorbar;
  hold on;
  par.numintervals = params.RB_numintervals;
  par.range = params.mu_ranges;
  par.color = [0,0,0];
  pgrid = cubegrid(par);
  plot(pgrid,par);
  %
  figure;
  load(['C:\Dokumente und Einstellungen\Bernard\Eigene Dateien\',...
        'TempMatlab\basisgen\RB_uniform_refined_2_2_r1.mat']);
  params.num_plot_intervals = [40,40];
  offline_data = rb_offline_prep(detailed_data,params);
  params.N = size(detailed_data.RB,2);
  online_data = rb_online_prep(offline_data,params);
  params.no_lines = 1;
  plot_error_estimator(online_data,params);
  c = get(gca,'Children');
  cdata = get(c(end),'Cdata');
  cdata = log10(cdata);
  set(c(end),'Cdata',cdata);
  set(gca,'Clim',[-8,-7]);
  title('log_{10}(\Delta(\mu))');
  colorbar;
  hold on;
  params.no_lines = 0;
  params.color = [0,0,0];
  plot(detailed_data.RB_info.MMesh_list{end},params);
  %CLIM : 2.07876e-008 7.30869e-008
  %
  figure;
  load(['C:\Dokumente und Einstellungen\Bernard\Eigene Dateien\',...
        'TempMatlab\basisgen\RB_adaptive_refined_2_2_r1_smax1000.mat']);
  params.num_plot_intervals = [40,40];
  offline_data = rb_offline_prep(detailed_data,params);
  params.N = size(detailed_data.RB,2);
  online_data = rb_online_prep(offline_data,params);
  params.no_lines = 1;
  plot_error_estimator(online_data,params);
  c = get(gca,'Children');
  cdata = get(c(end),'Cdata');
  cdata = log10(cdata);
  set(c(end),'Cdata',cdata);
  set(gca,'Clim',[-8,-7]);
  title('log_{10}(\Delta(\mu))');
  colorbar;
  hold on;
  params.no_lines = 0;
  params.color = [0,0,0];
  plot(detailed_data.RB_info.MMesh_list{end},params);
  
 case 3 % runtime diagrams
    
  %%%%% run basisgen_main with the following settings in the header:
  %
  %test_RB_methods = {...
  %    'RB_uniform_refined_2_2_r1',...
  %    'RB_uniform_refined_2_2_r5',...
  %    'RB_uniform_fixed_64',...
  %    'RB_uniform_fixed_256',...
  %    'RB_uniform_refined_2_2_r1',...
  %    'RB_uniform_refined_2_2_r5',...
  %    'RB_adaptive_refined_2_2_r1_smax1000_theta5e-2',...
  %    'RB_adaptive_refined_2_2_r5_smax1000_theta5e-2',...
  %    'RB_adaptive_refined_2_2_r1_smax1000_theta1e-1',...
  %    'RB_adaptive_refined_2_2_r5_smax1000_theta1e-1',...
  %                };
  % test_RB_indicators = {'estimator'}; % only estimators
  % compute_RB_methods = test_RB_methods;
  
  %%%%%% then perform the following on the generated files for getting 
  %%%%%% the figures
  
  %fns = {...
  %    'RB_uniform_refined_2_2_r1_tested',...
  %    'RB_uniform_fixed_256_tested',...
  %    'RB_adaptive_refined_2_2_r1_smax1000_theta5e-2_tested'};
  
  
  % uniform-refined approaches nur bis t=3600:
  fns = {'RB_uniform_refined_p3_i2_r1_tested',...
         'RB_uniform_refined_p3_i3_r1_tested',...
         'RB_adaptive_refined_p3_i2_r1_smax1000_theta5e-2_tested',...
         'RB_adaptive_refined_p3_i3_r1_smax1000_theta5e-2_tested',...
         'RB_uniform_fixed_p3_64_tested',...
         'RB_uniform_fixed_p3_125_tested',...
        };
  
  params.plot_train_times = 1;
  params.plot_train_estimators = 0;
  params.plot_test_estimators = 0;
  params.plot_test_errors = 0;
  params.plot_linestyles = {'-','-.','-','-.','-','-.'};
  params.plot_linecolors = {'r','r',[0,0.6,0],[0,0.6,0],'b','b'};
  params.plot_linewidths = [2,2,2,2,2,2];
  params.plot_legends_location = 'Northwest';
  params.plot_legends = {'uniform-refined 2^3',...
                    'uniform-refined 3^3',...
                    'adaptive-refined 2^3',...
                    'adaptive-refined 3^3',...
                    'uniform-fixed 4^3',...
                    'uniform-fixed 5^3'};
  
  plot_basisgen_results(fns,params);
  
  % uniform-refined approaches jetzt ohne timeout:
  fns = {'RB_uniform_refined_p3_i2_r1_notimeout_tested',...
         'RB_uniform_refined_p3_i3_r1_notimeout_tested',...
         'RB_adaptive_refined_p3_i2_r1_smax1000_theta5e-2_tested',...
         'RB_adaptive_refined_p3_i3_r1_smax1000_theta5e-2_notimeout_tested',...
         'RB_uniform_fixed_p3_64_tested',...
         'RB_uniform_fixed_p3_125_tested'};
  %'RB_uniform_fixed_p3_216_notimeout_tested',...
  
  params.plot_train_times = 0;
  params.plot_train_estimators = 1;
  params.plot_test_estimators = 0;
  params.plot_test_errors = 0;
  params.plot_legends_location = 'NorthEast';
  plot_basisgen_results(fns,params);
  
  % => Speichern der Bilder als 
  %    computation_times.eps und
  %    train_errors.eps
  
 case 4:
  
  %%%%% run basisgen_main with the following settings in the header:
  %
  %test_RB_methods = {...
  %    'RB_uniform_refined_2_2_r1',...
  %    'RB_uniform_refined_2_2_r5',...
  %    'RB_uniform_fixed_64',...
  %    'RB_uniform_fixed_256',...
  %    'RB_uniform_refined_2_2_r1',...
  %    'RB_uniform_refined_2_2_r5',...
  %    'RB_adaptive_refined_2_2_r1_smax1000_theta5e-2',...
  %    'RB_adaptive_refined_2_2_r5_smax1000_theta5e-2',...
  %    'RB_adaptive_refined_2_2_r1_smax1000_theta1e-1',...
  %    'RB_adaptive_refined_2_2_r5_smax1000_theta1e-1',...
  %                };
  % test_RB_indicators = {'estimator'}; % only estimators
  % compute_RB_methods = test_RB_methods;
  
  %%%%%% then perform the following on the generated files for getting 
  %%%%%% the figures
  
  %fns = {...
  %    'RB_uniform_refined_2_2_r1_tested',...
  %    'RB_uniform_fixed_256_tested',...
  %    'RB_adaptive_refined_2_2_r1_smax1000_theta5e-2_tested'};
  
  % uniform-refined approaches jetzt ohne timeout:
  fns = {'RB_uniform_refined_p3_i2_r1_notimeout_tested',...
         'RB_uniform_refined_p3_i3_r1_notimeout_tested',...
         'RB_adaptive_refined_p3_i2_r1_smax1000_theta5e-2_tested',...
         'RB_adaptive_refined_p3_i3_r1_smax1000_theta5e-2_notimeout_tested',...
         'RB_uniform_fixed_p3_64_tested',...
         'RB_uniform_fixed_p3_125_tested'};
  %'RB_uniform_fixed_p3_216_notimeout_tested',...
  
  params.plot_train_times = 0;
  params.plot_train_estimators = 0;
  params.plot_test_estimators = 1;
  params.plot_test_errors = 1;
  params.plot_linestyles = {'-','-.','-','-.','-','-.'};
  params.plot_linecolors = {'r','r',[0,0.6,0],[0,0.6,0],'b','b'};
  params.plot_linewidths = [2,2,2,2,2,2];
  params.plot_legends_location = 'North';
  params.plot_legends = {'uniform-refined 2^3',...
                    'uniform-refined 3^3',...
                    'adaptive-refined 2^3',...
                    'adaptive-refined 3^3',...
                    'uniform-fixed 4^3',...
                    'uniform-fixed 5^3'};
  
  plot_basisgen_results(fns,params);
  close(gcf-1);
  close(gcf);
  
  params.plot_legends_location = 'NorthWest';
  plot_basisgen_results(fns,params);
  close(gcf-2);
  close(gcf-1);
  
  % => Manuelles Speichern der Bilder als 
  %    test_errors.eps und
  %    test_error_ratio.eps
  
 case 5:
  
  fns = {'RB_uniform_refined_p3_i2_r1_notimeout_tested',...
         'RB_uniform_refined_p3_i3_r1_notimeout_tested',...
         'RB_adaptive_refined_p3_i2_r1_smax1000_theta5e-2_tested',...
         'RB_adaptive_refined_p3_i3_r1_smax1000_theta5e-2_notimeout_tested',...
         'RB_uniform_fixed_p3_64_tested',...
         'RB_uniform_fixed_p3_125_tested'};
  %'RB_uniform_fixed_p3_216_notimeout_tested',...
  
  params.plot_train_times = 0;
  params.plot_train_estimators = 0;
  params.plot_test_estimators = 1;
  params.plot_test_errors = 0;
  params.plot_max_test_vs_time = 1;
  params.plot_linestyles = {'-','-.','-','-.','-','-.'};
  params.plot_linecolors = {'r','r',[0,0.6,0],[0,0.6,0],'b','b'};
  params.plot_linewidths = [2,2,2,2,2,2];
  params.plot_legends_location = 'NorthEast';
  params.plot_legends = {'uniform-refined 2^3',...
                    'uniform-refined 3^3',...
                    'adaptive-refined 2^3',...
                    'adaptive-refined 3^3',...
                    'uniform-fixed 4^3',...
                    'uniform-fixed 5^3'};
  
  plot_basisgen_results(fns,params);
  close(gcf-3);
  close(gcf-2);
  close(gcf-1);
  set(gca,'Xlim',[0,3600]);
  set(gca,'Ylim',[1e-7,2e-4]);
  
end;