1.13.10/doc/cmp__rb__lin__nonlin__evol_8m_source.html

% small script performing a reduced basis simulation with the

% linear and nonlinear scheme and plotting the results

% the detailed_data for the nonlinear case can be modified and the

% simulation parameters switched, hereby checking suitability of a

% colateral basis for different parameter settings.


% Bernard Haasdonk 24.5.2007


% prepare and run nonlin_evol simulation


% get_nonlin_evol_params

% params.Mmax = 100;

% detailed_data = rb_detailed_prep(params);

 clear;

% load demo_nonlin_evol_detailed_data_LE_on_RB_wo_offset;

% load demo_nonlin_evol_detailed_data3;

 load demo_nonlin_evol_detailed_data_8_trajectories;

% load demo_nonlin_evol_detailed_data_LE_on_RB;

 offline_data = rb_offline_prep(detailed_data,params);


%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{1});


 % M usually less than N !!


 reduced_data = rb_online_prep(offline_data,params);

 %params.k = 1.25e-8;

 params.k = 0;

 %params.c_init = 1;

 params.c_init = 0.5;

 params.beta = 0.5;

 %params.k = 5e-8;


 disp('performing nonlinear simulation:');

 simulation_data = rb_simulation(reduced_data,params);

 Unonlin = rb_reconstruction(detailed_data,simulation_data);

 params.title = 'Solution from nonlinear rb-simulation';

 plot_element_data_sequence([],Unonlin,params);

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

 %plot_rb_reconstruction(detailed_data,simulation_data,params);


 disp('performing linear simulation:');

 % run linear simulation

 lparams = params;

 lparams.rb_problem_type = 'lin_evol';

 lparams.L_I_inv_norm_bound = 1;

 lparams.L_E_norm_bound = 1;

 %use same detailed-data, i.e. reduced basis

 offline_data = rb_offline_prep(detailed_data,lparams);

 reduced_data = rb_online_prep(offline_data,lparams);

 simulation_data = rb_simulation(reduced_data,lparams);

 Ulin = rb_reconstruction(detailed_data,simulation_data);

 lparams.title = 'Solution from linear rb-simulation';

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

 plot_element_data_sequence([],Ulin,lparams);


 % plot difference

 params.title = 'Difference of linear and nonlinear';

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

 plot_element_data_sequence([],Ulin-Unonlin,params);


 l2_errors = fv_l2_error(Ulin,Unonlin,detailed_data.W);


 figure,plot(l2_errors);


 title('l2-difference');

% TO BE ADJUSTED TO NEW SYNTAX

%| \docupdate