hp_demo_elliptic.m

Go to the documentation of this file.

%demo script for hp-approach routine using a elliptic problem
clear;

model = thermalblock_model_hp();

model.rb_init_data_basis = @RB_init_data_basis_stationary_default;
model.RB_error_indicator = 'estimator';
model.RB_extension_algorithm = @RB_extension_lin_stat_default;
model.get_rb_from_detailed_data = @(detailed_data)detailed_data.RB;



model.RB_train_size = 50;
model.RB_train_rand_seed = 1234;
model.RB_stop_Nmax = 8;
model.RB_stop_epsilon=1e-4;


% boundary box parameter:
hp_params.boundarybox = 0.25;
hp_params.boundarybox_minimum = 0.15;
hp_model=hp_gen_model_elliptic(model,hp_params);
%generate model data
model_data = gen_model_data(model);




%set the first anchor
hp_model = set_mu(hp_model,[1;1]);
%compute a h_refinment and RB generation
hp_model = hp_elliptic(hp_model,model_data);

hp_model = set_mu(hp_model,[0.13;0.13]);

sim_data = hp_simulation(hp_model);