demo_rb_steps_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 oop
% version of models. As example a thermalblock model is
% used. 
%
% See also
% demo_rb_steps_struct.m for the same steps using a struct-oriented
% approach 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 and generate model
disp('generation of problem description and model:');
params = [];
params.B1 = 2;
params.B2 = 2;
params.numintervals_per_block = 20;
descr = thermalblock_model(params);
[dmodel, rmodel] = gen_models(descr);
disp('Please type dbcont after inspecting the workspace variables');
keyboard;

% generate model_data 
disp('generation of model_data (grid, fem-info, etc.):');
model_data = gen_model_data(dmodel);
disp('Please type dbcont after inspecting the workspace variables');
keyboard;

% perform detailed simulation and plot
plot_params = [];
plot_params.axis_equal = 1;
plot_params.axis_tight = 1;
sim_data = detailed_simulation(dmodel,model_data);
plot_sim_data(dmodel,model_data,sim_data,plot_params);
disp('Please type dbcont after inspecting the workspace variables');
keyboard;

% generation of high dimensional offline data for RB-scheme
disp('computation of high dimensional offline-data (basis by greedy):');
detailed_data = gen_detailed_data(rmodel,model_data);
disp('Please type dbcont after inspecting the workspace variables');
keyboard;

% generation of low dimensional offline data for RB-scheme
disp('computation of low dimensional offline-data:');
reduced_data = gen_reduced_data(rmodel,detailed_data);
disp('Please type dbcont after inspecting the workspace variables');
keyboard;

% perform online simulation
rmodel.N = 4;
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,1,1,1]);
sim_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:');
sim_data = rb_reconstruction(rmodel,detailed_data,sim_data);
plot_sim_data(dmodel,model_data,sim_data,plot_params);
disp('Please type dbcont after inspecting the workspace variables');
keyboard;

% complete interactive gui:
disp('complete interactive gui:');
demo_rb_gui(rmodel,detailed_data,[],plot_params);