3 %
function initializing some fields of a lin-evol model
5 % generated fields of descr:
6 % rb_problem_type : ''LinEvol''
7 % L_I_inv_norm_bound : '1'
9 % operators_ptr : fv_operators_implicit_explicit()
10 % operators_diff_implicit : fv_operators_zero()
11 % operators_diff_explicit : fv_operators_zero()
12 % operators_conv_implicit : fv_operators_zero()
13 % operators_conv_explicit : fv_operators_zero()
14 % operators_neumann_explicit : fv_operators_zero()
15 % operators_neumann_implicit : fv_operators_zero()
17 % get_inner_product_matrix : extracts inner product matrix from model data.
18 % inner_product : inner product algorithm
19 % mass_matrix : fv_mass_matrix()
21 % l2_error_sequence_algorithm : fv_l2_error()
22 % linfty_error_sequence_algorithm : fv_linfty_error()
23 % error_algorithm : fv_error()
25 % error_estimation : '1'
26 % affinely_decomposed : 'true'
27 % compute_output_functional : 'false'
28 % data_const_in_time : '1'
29 % lxf_lambda : '1.0194e+003'
30 % enable_error_estimator : '1'
32 descr.rb_problem_type = 'LinEvol';
34 descr.L_I_inv_norm_bound = 1;
35 descr.L_E_norm_bound = 1;
37 descr.operators_ptr = @fv_operators_implicit_explicit;
38 descr.operators_diff_implicit = @fv_operators_zero;
39 descr.operators_diff_explicit = @fv_operators_zero;
40 descr.operators_conv_implicit = @fv_operators_zero;
41 descr.operators_conv_explicit = @fv_operators_zero;
42 descr.operators_neumann_explicit = @fv_operators_zero;
43 descr.operators_neumann_implicit = @fv_operators_zero;
46 descr.get_inner_product_matrix = @(detailed_data) detailed_data.W;
47 descr.inner_product = @(descr,descr_data,vecs1,vecs2) ...
48 vecs1
' * descr_data.W * vecs2;
49 descr.mass_matrix = @fv_mass_matrix;
50 descr.divclean_mode = 0;
52 descr.l2_error_sequence_algorithm = @fv_l2_error;
53 descr.linfty_error_sequence_algorithm = @fv_linfty_error;
54 descr.error_algorithm = @fv_error;
55 descr.error_norm = 'l2
';
56 descr.error_estimation = 1; % turn on rb error estimation as default
58 descr.affinely_decomposed = true;
60 descr.compute_output_functional = false;
62 descr.data_const_in_time = 1;
64 % descr.lxf_lambda = 0.1172; % lambda value for Lax-Friedrichs diffusivity
65 % descr.lxf_lambda = 0.001172; % lambda value for Lax-Friedrichs diffusivity
66 descr.lxf_lambda = 1.0194e+003;
68 descr.enable_error_estimator = 1;