1.13.10/doc/cmp__rb__detailed__simulation_8m_source.html

% small script performing a reduced basis simulation and a detailed

% simulation for specified parameters. Hereby, visual difference

% can checked.


% Bernard Haasdonk 16.9.2007


clear all;

echo on;

%load demo_nonlin_evol_detailed_data_8_trajectories;

%load(fullfile(rbmatlabresult,'chemnitz_detailed_data'));

load(fullfile(rbmatlabresult,'burgers_fv_detailed_interpol'));

echo off;


offline_data = rb_offline_prep(detailed_data,params);


params.M = 60;

params.N = 40;

%params.N = 25; params.M = length(detailed_data.TM);

%params.N = 123; params.M = length(detailed_data.TM);

%params.N = size(detailed_data.RB,2); params.M = length(detailed_data.TM);


 % M usually less than N !!


 reduced_data = rb_online_prep(offline_data,params);

 %params.c_init = 1;

 %params.k = 5e-8;

 %params.c_init = 0.4;

 %params.beta = 0.4;

 params.c_init_lo = 0.1;

 params.vrot_angle = -pi/5;


 disp('performing reduced simulation:');

 simulation_data = rb_simulation(reduced_data,params);

 Urb = rb_reconstruction(detailed_data,simulation_data);

 params.title = 'Solution from RB-simulation';

 plot_element_data_sequence([],Urb,params);

% plot_element_data_sequence([],Unonlin(:,1:100),params);

 %plot_rb_reconstruction(detailed_data,simulation_data,params);


 disp('performing detailed simulation:');

 % run linear simulation

 Udetailed = detailed_simulation(detailed_data.grid,params);

 params.title = 'Solution from detailed simulation';

 % plot_element_data_sequence([],Ulin(:,1:100),lparams);

 plot_element_data_sequence([],Udetailed,params);


 % plot difference

 params.title = 'Difference of detailed and reduced';

% plot_element_data_sequence([],Ulin(:,1:100)-Unonlin(:,1:100),params);

 plot_element_data_sequence([],Urb-Udetailed,params);


 l2_errors = fv_l2_error(Urb,Udetailed,detailed_data.W);


 figure,plot(l2_errors);


 title('l2-difference');

% TO BE ADJUSTED TO NEW SYNTAX

%| \docupdate