test_derivative_lin_evol_opt.m

function test_derivative_lin_evol_opt(step)

switch step
% rb_derivative_simulation
    case 1
params = [];
params.coarse_factor =  16;
model = advection_fv_output_opt_model(params);

model.RB_generation_mode=  'PCA_trajectories';
M=[0,0,0;0,0,1;0,1,0;1,0,0;0,1,1;1,0,1;1,1,1;1,1,0];

for i = 1:size(M,1)
model.RB_mu_list{i} = M(i,:);
end

model_data = gen_model_data(model);
detailed_data = gen_detailed_data(model, model_data);
save('PCA_basis_coarse16.mat', 'model', 'detailed_data', 'model_data');

%% load basis
    case 2
load('PCA_basis_coarse16.mat');


params = [];
params.coarse_factor =  16;
model = advection_fv_output_opt_model(params);

reduced_data = gen_reduced_data(model, detailed_data);

sim_data_der = rb_derivative_simulation(model, reduced_data);


model_data = gen_model_data(model);
model.decomp_mode=0;
sim_data = model.rb_reconstruction(model,model_data, detailed_data, sim_data_der);

keyboard;

end