demo_rb_lin_evol_oop.m

Go to the documentation of this file.

% small collection of RBmatlab's rb-simulation abilities
% please inspect the source-file during execution, as the main
% purpose is to demonstrate the possibilities of RBmatlab's
% rb-philosophy
%
% this script demonstrates use for object oriented programming
% version of models. As example a convection-diffusion model is
% used. 
%
% See also
% demo_rb_steps_struct.m for the same steps using a struct as model

% Bernard Haasdonk 11.4.2007

help demo_rb_steps_oop;

disp('Small demonstration of RBmatlab reduced basis commands');
disp(['Please type dbcont after inspecting the workspace variables' ...
      ' at different halt points.']);
disp('Opening figures can be closed.');

% detailled simulation:
clear;
% get problem description
descr = convdiff_descr;

% add some settings to perform output estimation
descr.compute_output_functional = 1;
descr.operators_output = @fv_operators_output;
descr.output_function_ptr = @output_function_box_mean;

% define complete domain for output estimation
descr.sbox_xmin = min(descr.xrange);
descr.sbox_xmax = max(descr.xrange);
descr.sbox_ymin = min(descr.yrange);
descr.sbox_ymax = max(descr.yrange);

disp(descr);
% load models
disp('initializing model:');
[dmodel, rmodel] = gen_models(descr);


%model_data = gen_model_data(dmodel);
%detailed_data = gen_detailed_data(rmodel, model_data);
%save('demos/datafiles/demos_rb_lin_evol_oop.mat', 'detailed_data');

disp(dmodel);
disp('loading precomputed detailed_data:');
load(fullfile('demos', 'datafiles', 'demos_rb_lin_evol_oop.mat'), ...
    'detailed_data');
disp(detailed_data);

% 
% disp('detailed linear evolution simulation:');
% model_data = gen_model_data(dmodel);
% sim_data = detailed_simulation(dmodel, model_data);
% params.plot_title = 'detailed simulation';
% params.plot = @fv_plot;
% plot_sim_data(dmodel, model_data, sim_data, params);
% disp('Please type dbcont after inspecting the workspace variables');
% keyboard;
% % 
% % reduced basis loading and plotting
% 
% clear('sim_data');
% disp('reduced basis plotting:');
% params.title = 'reduced basis vectors';
% params.plot = @fv_plot;
% dd_leaf = get_leaf(detailed_data, rmodel);
% plot_sequence(dd_leaf.RB(:,1:10),detailed_data.model_data.grid,params);
% %plot_mu_frequency(detailed_data.mu_values,params);
% disp('Please type dbcont after inspecting the workspace variables');
% keyboard;


% computation of mu and N-independent matrices 
disp('computation of mu- and N-independent offline-data:');
reduced_data = gen_reduced_data(rmodel, detailed_data);
disp('Please type dbcont after inspecting the workspace variables');
keyboard;

% computation of mu and N-independent matrices 
disp(['Selection of N and computation of mu-independent but',...
      ' N-dependent online-data:']);
rmodel.N = detailed_data.Nmax;
reduced_data = extract_reduced_data_subset(rmodel, reduced_data);
disp('Please type dbcont after inspecting the workspace variables');
keyboard;

% mu-parameter variation: sets the fields of params:
%    params.c_init = 0; params.beta = 1; params.k = 5e-8; 
disp('Selection of mu and online-simulation');
set_mu(rmodel, [0,1,5e-8]);

disp('reduced basis simulation:');
simulation_data = rb_simulation(rmodel, reduced_data);
disp('Please type dbcont after inspecting the workspace variables');
keyboard;

% reconstruction and plot of reduced simulation
disp('reconstruction of reduced simulation:');
simulation_data = rb_reconstruction(rmodel,detailed_data,simulation_data);
params.title = 'reduced simulation';

dmodel.plot_sim_data(detailed_data.model_data, simulation_data, params);
figure;
rb_plot_output_estimation(simulation_data,dmodel);
disp('Please type dbcont after inspecting the workspace variables');