test_rb_richards_fv.m

function OK = test_rb_richards_fv
%function OK = test_rb_richards_fv;
%
% function generating a reduced basis for the richards_fv problem for 15
% parameters. This tests the generation of a reduced basis that depends on an
% empirical interpolation of the explicit operator. The test does a complete
% reduced basis simulation for a parameter that is not in the training
% parameter set and assures that the error is kept resonably small.
%
% OK == 1 if test is OK, otherwise 0
%
% Martin Drohmann 16.1.2009

OK = 1;

% we make it somewhat faster by truncating the time interval and restricting
% the collateral and reduced basis to 20 resp. 10 elements.

params.model_type = 'test';
params.model_size = 'small';
params.verbose    = 0;

descr    = richards_fv_descr(params);
bg_descr = richards_fv_bg_descr(params);

[dmodel, rbmodel] = gen_models(descr, bg_descr);

model_data = gen_model_data(dmodel);

rbmodel.Mstrich = 1;

% constructing the collateral and the reduced basis
try
  evalc('detailed_data = gen_detailed_data(rbmodel, model_data);');
catch exception
  OK = 0
  disp('Basis generation failed');
  disp(getReport(exception));
end
disp('basis generation terminated');

%model.M = cellfun(@(x)(size(x,2) - model.Mstrich), detailed_data.BM, 'UniformOutput', true) - model.Mstrich;

rbdd = get_by_description(detailed_data.datatree, 'rb', rbmodel);
reduced_data = gen_reduced_data(rbmodel, rbdd);
disp('offline preparation terminated');

reduced_data = extract_reduced_data_subset(rbmodel, reduced_data);
disp('online preparation terminated');

set_mu(rbmodel, [0.5, 0.4]);
set_mu(dmodel, [0.5, 0.4]);

simulation_data = rb_simulation(rbmodel, reduced_data);
disp('rb simulaton terminated');

simulation_data = rb_reconstruction(rbmodel, rbdd, simulation_data);
disp('reconstruction finished');

sim_data = detailed_simulation(dmodel, model_data);
disp('detailed simulation finished');

err = dmodel.l2_error_sequence_algorithm(sim_data.U, simulation_data.U, model_data);
disp('error computation finished');

save(fullfile(rbmatlabhome, 'test', 'test_rb_richards_fv_data.mat'), 'rbmodel', 'detailed_data');

if err(end) > 1e-3;
  disp('error is too large!!!');
  OK = 0;
end;

%| \docupdate