demo_rb_error_gui.m

Go to the documentation of this file.

function demo_rb_error_gui(varargin)
%function demo_rb_error_gui(model,detailed_data,[reduced_data],...
%                          plot_params[,title, callbackfigure, ...
%                   cbhandle])
% demo gui for detailed simulations and reduced simulation and
% plotting the error
%
% (calling demo_rb_gui after switching some pointers)

% B. Haasdonk 2.3.2012

tlist = varargin;
if nargin < 1
  disp('no model and data specified... generating thermalblock data')
    params = [];
  params.B1 = 2;
  params.B2 = 2;
  params.numintervals_per_block = 50;
  params.mu_range = [0.1;10];
  params.numintervals_per_block = 5;
  model = thermalblock_model_struct(params);
  model_data = gen_model_data(model);
  detailed_data = gen_detailed_data(model,model_data);
  reduced_data = gen_reduced_data(model,detailed_data);
  my_reduced_data.reduced_data = reduced_data;
  my_reduced_data.detailed_data = detailed_data;
  plot_params = [];
  plot_params.axis_equal = 1;
  plot_params.axis_tight = 1;
  plot_params.yscale_uicontrols = 0.7;
  tlist{1} = model;
  tlist{2} = detailed_data;
  tlist{3} = my_reduced_data;
  tlist{4} = plot_params;
end;

model = tlist{1};

model.original_rb_simulation = model.rb_simulation;
model.rb_simulation = @my_rb_simulation;

model.original_gen_reduced_data = model.gen_reduced_data;
model.gen_reduced_data = @my_gen_reduced_data;

model.original_rb_reconstruction = model.rb_reconstruction;
model.rb_reconstruction = @my_rb_reconstruction;

model.original_reduced_data_subset = model.reduced_data_subset;
model.reduced_data_subset = @my_reduced_data_subset;

model.original_get_dofs_from_sim_data = model.get_dofs_from_sim_data;
model.get_dofs_from_sim_data = @(sim_data) ...
    my_get_dofs_from_sim_data(sim_data,model);

model.original_plot_sim_data = model.plot_sim_data;
model.plot_sim_data = @my_plot_sim_data;

model.enable_error_estimator = 0;

tlist{1} = model;
tlist{5} = 'reduced basis error';
demo_rb_gui(tlist{:});

function my_data = my_gen_reduced_data(model,detailed_data)
my_data.detailed_data = detailed_data;
my_data.reduced_data = model.original_gen_reduced_data(model,detailed_data);

function my_sim_data = my_rb_simulation(model,my_data)
my_sim_data.sim_data = detailed_simulation(model,my_data.detailed_data);
if my_data.reduced_data.N ~= model.N
  my_data = extract_reduced_data_subset(model,my_data);  
end;
my_sim_data.rb_sim_data = model.original_rb_simulation(model, ...
                                           my_data.reduced_data);

function my_sim_data = my_rb_reconstruction(model,detailed_data,my_sim_data)
my_sim_data.rb_sim_data = model.original_rb_reconstruction(...
    model,detailed_data,my_sim_data.rb_sim_data);

function reduced_data_subset = my_reduced_data_subset(model,my_data)
reduced_data_subset = my_data;
reduced_data_subset.reduced_data = ...
    model.original_reduced_data_subset(model,my_data.reduced_data);

function dofs = my_get_dofs_from_sim_data(my_sim_data,model);
dofs = model.original_get_dofs_from_sim_data(my_sim_data.sim_data)-...
       model.original_get_dofs_from_sim_data(my_sim_data.rb_sim_data);

function p = my_plot_sim_data(varargin);
%function p = my_plot_sim_data(model,model_data,my_sim_data,plot_params);
model = varargin{1};
sim_data = varargin{3}.sim_data;
dofs_detailed = model.original_get_dofs_from_sim_data(sim_data);
dofs_reduced = model.original_get_dofs_from_sim_data(varargin{3}.rb_sim_data);
sim_data = model.set_dofs_in_sim_data(sim_data,...
                                           dofs_detailed-dofs_reduced);
tlist = varargin;
tlist{3} = sim_data;
p = model.original_plot_sim_data(tlist{:});