testing_script_rb_operators_decomp.m

Go to the documentation of this file.

% Testskript für lin_evol_opt_rb_operators
% zum Testen der einzelnen decomp_modes

clear

load PCA_basis_coarse16_3_50

model.t=30;
model.compute_derivative_info=0;

if model.compute_derivative_info == 1
model.decomp_mode = 0;
[LL_I_0, LL_E_0, bb_0, M_E_0, M_b_0, M_EE_0, M_Eb_0, M_bb_0, M_I_0, ...
    M_II_0, M_IE_0, M_Ib_0, M_EdEd_0, M_IdId_0, M_bdbd_0, M_Ed_0, ...
    M_Id_0, M_bd_0, M_IEd_0, M_IId_0, M_Ibd_0, M_EEd_0, ...
    M_EId_0, M_Ebd_0, M_EdId_0, M_Edbd_0, M_Idbd_0] = ...
    lin_evol_opt_rb_operators(model, detailed_data);

model.decomp_mode = 1;
[LL_I_1, LL_E_1, bb_1, M_E_1, M_b_1, M_EE_1, M_Eb_1, M_bb_1, M_I_1, ...
    M_II_1, M_IE_1, M_Ib_1, M_EdEd_1, M_IdId_1, M_bdbd_1, M_Ed_1, ...
    M_Id_1, M_bd_1, M_IEd_1, M_IId_1, M_Ibd_1, M_EEd_1, ...
    M_EId_1, M_Ebd_1, M_EdId_1, M_Edbd_1, M_Idbd_1] = ...
    lin_evol_opt_rb_operators(model, detailed_data);

model.decomp_mode = 2;
[LL_I_2, LL_E_2, bb_2, M_E_2, M_b_2, M_EE_2, M_Eb_2, M_bb_2, M_I_2, ...
    M_II_2, M_IE_2, M_Ib_2, M_EdEd_2, M_IdId_2, M_bdbd_2, M_Ed_2, ...
    M_Id_2, M_bd_2, M_IEd_2, M_IId_2, M_Ibd_2, M_EEd_2, ...
    M_EId_2, M_Ebd_2, M_EdId_2, M_Edbd_2, M_Idbd_2] = ...
    lin_evol_opt_rb_operators(model, detailed_data);

else %model.compute_derivative_info=1
model.decomp_mode = 0;
[LL_I_0, LL_E_0, bb_0, M_E_0, M_b_0, M_EE_0, M_Eb_0, M_bb_0, M_I_0, ...
    M_II_0, M_IE_0, M_Ib_0] = ...
    lin_evol_opt_rb_operators(model, detailed_data);

model.decomp_mode = 1;
[LL_I_1, LL_E_1, bb_1, M_E_1, M_b_1, M_EE_1, M_Eb_1, M_bb_1, M_I_1, ...
    M_II_1, M_IE_1, M_Ib_1] = ...
    lin_evol_opt_rb_operators(model, detailed_data);

model.decomp_mode = 2;
[LL_I_2, LL_E_2, bb_2, M_E_2, M_b_2, M_EE_2, M_Eb_2, M_bb_2, M_I_2, ...
    M_II_2, M_IE_2, M_Ib_2] = ...
    lin_evol_opt_rb_operators(model, detailed_data);
end