scripts/thermalblock.m

Go to the documentation of this file.

function thermalblock(step)
% thermalblock example
%  step 1: detailed simulation with different parameters
%  step 2: reduced basis simulation
%  step 3: demo_rb_gui

if nargin<1
  step = 1;
end;

switch step
 case 1 % different detailed simulations

  disp('detailed simulation:');
  
  
  descr      = thermalblock_descr;
  dmodel     = gen_detailed_model(descr);

  model_data = gen_model_data(dmodel);
  dmodel = set_mu(dmodel,[1,1,1,1,1,1,1,1,1]);

  sim_data = detailed_simulation(dmodel,model_data);
  plot_params.title = ['detailed solution , \mu_i=1, output s=',...
                       num2str(sim_data.s)];
  figure, plot_sim_data(dmodel,model_data,sim_data,plot_params);
  dmodel = set_mu(dmodel,[0.1,0.1,1,1,0.1,0.1,1,0.1,1]);

  sim_data = detailed_simulation(dmodel,model_data);
  plot_params.title = ['detailed solution , \mu_i=0.1/1, output s=',...
                       num2str(sim_data.s)];
  figure, plot_sim_data(dmodel,model_data,sim_data,plot_params);

 case 2 % rb treatment and comparison to detailed

  % for non-snapshot parameter nice error behaviour is observed
  % for snapshot parameter error is numerically zero, but
  % estimators due to numerical noise not zero.

  disp('detailed simulation:');
  descr            = thermalblock_descr;
  [dmodel, rmodel] = gen_models(descr);
  model_data       = gen_model_data(dmodel);
  % non-snapshot parameter, error <= Delta
  dmodel = set_mu(dmodel,[0.5,1,1,1,1,1,1,1,1]);
  % non-snapshot parameter with excellent Delta==error!!!
  %model = model.set_mu(model,[1,1,1,1,1,1,1,1,1]);
  % snapshot parameter with error, Delta approximately 0.
  %model = model.set_mu(model,[0.1,1,1,1,1,1,1,1,1]);
  sim_data = detailed_simulation(dmodel,model_data);
  sim_data.mu = get_mu(dmodel);
  plot_params.title = 'detailed solution';
  disp(['output s(mu) = ',num2str(sim_data.s)]);
  figure, plot_sim_data(dmodel,model_data,sim_data,plot_params);
  disp('offline computations I (gen_detailed_data: basis):');
  detailed_data = gen_detailed_data(rmodel,model_data);
  disp('offline computations II (gen_reduced_data: operator components):');
  reduced_data = gen_reduced_data(rmodel,detailed_data);
%  model.N = 1;
%  reduced_data = model.reduced_data_subset(model,reduced_data);

  disp('online reduced simulation:');
  rmodel = set_mu(rmodel, sim_data.mu);
  rb_sim_data = rb_simulation(rmodel,reduced_data);
  rb_sim_data = rb_reconstruction(rmodel, detailed_data, rb_sim_data);
  plot_params.title = 'reduced solution';
  figure, plot_sim_data(rmodel,model_data,rb_sim_data,plot_params);
  disp(['output sN(mu) = ',num2str(rb_sim_data.s)]);

  disp('-------------------------------------------');
  eh = sim_data.uh - rb_sim_data.uh;
  l2_err = fem_l2_norm(eh);
  h10_err = fem_h10_norm(eh);
  disp(['L2-error: ',num2str(l2_err)]);
  disp(['H10-error: ',num2str(h10_err)]);
  disp(['error s(mu)-sN(mu) = ',num2str(sim_data.s-rb_sim_data.s)]);
  disp(['error estimator Delta_u(mu) = ',num2str(rb_sim_data.Delta)]);
  disp(['error estimator Delta_s(mu) = ',num2str(rb_sim_data.Delta_s)]);

  % save detailed_data for demo_rb_gui:
  save(fullfile(rbmatlabtemp,'thermal_block_offline_data.mat'),...
       'dmodel','detailed_data');

 case 3 % demo_rb_gui

  fn = fullfile(rbmatlabtemp,'thermal_block_offline_data.mat');
  if ~exist(fn) % call step 2 first
    disp('calling thermalblock(2) for basis generation')
    thermalblock(2);
  end;

  disp('demo_rb_gui:')
  descr            = thermalblock_descr;
  [~, rmodel] = gen_models(descr);
  load(fn);
  rmodel.N = size(detailed_data.RB,2);
  %model_data       = gen_model_data(dmodel);
  %detailed_data    = gen_detailed_data(dmodel,model_data);
  plot_params.axis_tight = 1;
  plot_params.yscale_uicontrols = 0.5;
  demo_rb_gui(rmodel,detailed_data,[],plot_params);

end;