test_rb_basisgen.m

Go to the documentation of this file.

function OK = test_rb_basisgen
%function OK = test_rb_basisgen
% test showing a simple greedy basis generation
%
% function generating a basis as stored in
% demo_lin_evol_detailed_data.mat by performing a greedy basis extension
% on a very small 2x2x2 parameter grid without parameter sampling adaptation.
%
% OK == 1 if test is OK, otherwise 0

% Bernard Haasdonk 20.8.2007

OK = 1;

descr = convdiff_descr;
descr.RB_Nmax                = 20;
descr.RB_train_num_intervals = 1;
descr.RB_val_size            = 0;
descr.RB_refinement_mode     = 'none';
[dmodel, rbmodel]  = gen_models(descr);

model_data = gen_model_data(dmodel);

try
  % only generate 20 basis functions, otherwise too long.
  evalc('detailed_data = gen_detailed_data(rbmodel, model_data);');
catch exception
  OK = 0;
  disp('Execution of greedy algorithm failed with error:');
  disp(getReport(exception));
end

if get_rb_size(detailed_data.datatree) > 20
  disp('RB basis is too large!')
  OK = 0;
end;

max_err_sequence = get_field(detailed_data.datatree, 'max_err_sequence');

if max_err_sequence(20) > 1.2370e-04 * 1.1;
  disp('error(20) is too large !')
  OK = 0;
end;

if max_err_sequence(20) < 1.2370e-04 * 0.1;
  disp('error(20) is too small !')
  OK = 0;
end;