demo_rb_steps_struct.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 struct
% version of models. As example a thermalblock model is
% used. 
%
% See also
% demo_rb_steps_oop.m for the same steps using a object-oriented
% approach as model

% Bernard Haasdonk 11.4.2007

help demo_rb_steps_struct;

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;
model = thermalblock_model_struct(params)
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(model)
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(model,model_data)
plot_sim_data(model,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(model,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(model,detailed_data)
disp('Please type dbcont after inspecting the workspace variables');
keyboard;

% perform online simulation
model.N = 4
reduced_data = model.reduced_data_subset(model,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(model, [0.1,1,1,1]);
sim_data = rb_simulation(model, 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(model,detailed_data,sim_data)
plot_sim_data(model,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(model,detailed_data,[],plot_params);