scripts/nonlin_symmetry.m

Go to the documentation of this file.

function rmodel = nonlin_symmetry(steps, combined, M_by_N_ratio, noof_ei_extensions, random, explicit, num_cpus)
% function rbmodel = nonlin_symmetry(steps, combined, M_by_N_ratio, noof_ei_extensions, random, explicit, num_cpus)
% small script demonstrating the burgers equation with explicit fv
% discretization and RB model reduction
% example is meant to demonstrate automatic symmetry detection by
% the algorithm.
% These results are presented at HYP2008

% Bernard Haasdonk 3.6.2008

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%% Select here, what is to be performed with the burgers data %%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%step = 1 % single detailed simulation with given data and plot. Run
         % this with varying parameters mu until sure that scheme
         % is stable. Modify dt or the data-functions accordingly,
         % until a nice parameter-domain with uniformly stable
         % detailed scheme is obtained.
%step = 2 % generate colateral reduced basis of L_E operator
%step = 3 % compare single detailed simulation with and without
         % interpolated space operator
%step = 4 % generate dummy reduced basis from single trajectory and check, if
         % ei_interpolation with projection on this space maintains
         % result. A simple reduced simulation can also be
         % performed. All results should be visually identical.
%step = 5 % generate reduced basis
%step = 6 % time measurement of reduced simulation and demo rb gui
%step = 7 % generate error-landscape over varying N and M.
          % This can take several hours!!!

if nargin == 0
  steps = [0,6];
end

%steps = [1,2,3,4,5,7,8];

if nargin <= 1
  combined           = true;
  M_by_N_ratio       = 0;
  noof_ei_extensions = 1;
  random             = false;
  explicit           = true;
  num_cpus           = 4;
end

rmodel = [];

for si=1:length(steps)

step = steps(si);

% output-filenames in rbmatlabresult
detailedfname = 'nonlin_symmetry_detailed.mat';
testdir = 'nonlin_symmetry_test_10';
results_fname_base = 'nonlin_symmetry_step7_result';

plot_params.axis_equal = 1;
plot_params.clim = [0,1];
plot_params.plot = @fv_plot;

params = [];
params.explicit = explicit;
params.model_type = 'nonlin_symmetry';
params.step7_outputfile = fullfile(rbmatlabresult, [params.model_type, '_step7_output']);
params.step8_outputfile = fullfile(rbmatlabresult, [params.model_type, '_step8_output']);

extrapar.Mstrich = 100;

% grid parameters
switch step
  case 0
    if ~isempty(getenv('MACHINEFILE'))
      machinefile = getenv('MACHINEFILE');
      FakeMPI.MPI_Init(machinefile, fullfile(rbmatlabhome, 'fakeMPI', 'run_parCompErrs.sh'));
      FakeMPI.MPI_Run;
    end

    descr  = nonlin_symmetry_descr(params);

    if ~combined
      descr.RB_Mmax_small = 0;
    else
      descr.RB_Mmax_small = 20;
    end

    descr.RB_extension_M_by_N_r = M_by_N_ratio; % 2.5;
    descr.RB_refinement_theta   = 0.15;
    descr.RB_noof_ei_extensions = noof_ei_extensions;
    descr.RB_train_num_intervals = 25;
    descr.RB_val_size            = 4;
    descr.RB_stop_Mmax           = 600;
    descr.RB_stop_Nmax           = 250;
    descr.RB_stop_timeout        = 12*60*60;
    descr.RB_stop_epsilon        = 1e-9;
    descr.RB_epsilon_M_by_N_r    = 1e-1;

    if random
      descr.RB_train_sample_mode        = 'random';
      descr.RB_train_num_intervals      = 400;
      descr.RB_stop_max_val_train_ratio = inf;
    end

    [dmodel, rmodel] = gen_models(descr);

    params.mu_ranges = rmodel.mu_ranges;
    params.mu_names  = rmodel.mu_names;

    dmodel.num_cpus   = num_cpus;

    model_data = gen_model_data(dmodel);
 case 1 % single detailed simulation and plot
  step1_detailed_simulation;
 case 5 % reduced basis
   detailed_data = gen_detailed_data(rmodel, model_data);
   save(fullfile(rbmatlabresult, detailedfname), 'detailed_data', 'rmodel');
   FakeMPI.MPI_Finalize;
 case 6
  load(fullfile(rbmatlabresult, detailedfname), 'detailed_data', 'rmodel');
  step6_demo_rb_gui;
end

end

%| \docupdate