lin_evol_opt_gen_detailed_data_separate_bases_parallel.m

Go to the documentation of this file.

function detailed_data =  lin_evol_opt_gen_detailed_data_separate_bases_parallel(model, model_data)
%function detailed_data =
%lin_evol_opt_gen_detailed_data_separate_bases(model, model_data)
%
% function generating the detailed data for models with sensitivity
% derivatives. Separate bases for the approximation of the solution and the
% derivatives are generated.
%   
%   
% Optional fields of model:
%   
%
% Generated fields of detailed_data:
%   RB:             The reduced bases for the solution
%   RB_der:         cell array with thr reduced bases for the derivative solution
%   RB_info:        Usual info field
%   RB_der_info:    cell array containing the RBinfo for the derivative
%                   reduced bases
%

% Markus Dihlmann 31.03.2011

%get number of parameters to know how many derivative_rb have to be
%generated
nr_parameters = length(get_mu(model));

Nmax = model.RB_stop_Nmax;
epsilon_tol = model.RB_stop_epsilon;
% check some fields in model
if isfield(model,'RB_mixed_sensitivity_basis_generation')
    if model.RB_mixed_sensitivity_basis_generation
        disp('Warning: The field model.RB_mixed_sensitivity_basis_generation is probably set to the wrong value');
        model.RB_mixed_sensitivity_basis_generation=0;
    end
end


%1. start generating a normal reduced basis
%-----------------------------------------
compute_derivative_info_old = model.compute_derivative_info;
model.compute_derivative_info =0;
if isfield(model,'base_model')&&isfield(model.base_model,'compute_derivative_info')
    model.base_model.compute_derivative_info = 0;
end

%set fields, so that the error and the dofs of the solution is returned
model.get_estimator_from_sim_data = @(sim_data)get_estimator_from_sim_data_mixed(sim_data,0);
model.get_dofs_from_sim_data = @(sim_data)get_dofs_from_sim_data_for_sensitivity_extension(sim_data,0);
model.get_rb_size = @(model, detailed_data)lin_evol_opt_get_rb_size(model, detailed_data,0);

if isfield(model,'RB_solution_Nmax')
    model.RB_stop_Nmax=model.RB_solution_Nmax;
else
    model.RB_stop_Nmax = 2*Nmax;
end
if isfield(model,'RB_solution_epsilon')
    model.RB_stop_epsilon = model.RB_solution_epsilon;
else
    model.RB_stop_epsilon = epsilon_tol/10;
end

detailed_data = lin_evol_gen_detailed_data(model, model_data);
RB_finished = detailed_data.RB;
RB_info = detailed_data.RB_info;
model_data.RB = RB_finished;
filecache_clear;

%2. generate a reduced basis for every derivative in a loop
%-----------------------------------------------------------
model.RB_stop_Nmax = Nmax;
model.RB_stop_epsilon = epsilon_tol;
model.compute_derivative_info =1;
if isfield(model,'base_model')&&isfield(model.base_model,'compute_derivative_info')
    model.base_model.compute_derivative_info = 1;
end
RB_der_info = cell(1,nr_parameters);
RB_der = cell(1,nr_parameters);
for i=1:length(RB_der)
    %empty initialization with correct size for lin_evol_opt_rb_operators
    %to work a priori for all p
    RB_der{i} = zeros(size(RB_finished,1),1);
end
model_data.RB_der = RB_der;
%model_data.RB = RB_finished;
%compute_derivative_indices_old=model.compute_derivative_indices;
%indices_to_compute = find(compute_derivative_indices_old);
matlabpool open
parfor p=1:nr_parameters
    disp('-----------------------------------------------');
    disp(['Generating RB of derivative ',num2str(p),'...'])
    disp('-----------------------------------------------');
    %set fields, so that the right error is produced
    model_help = model;
    model_help = setfield_in_model_and_base_model(model_help,'set_rb_in_detailed_data',@(detailed_data,RB)lin_evol_opt_set_rb_in_detailed_data_separate_bases(detailed_data,RB,p));
    model_help = setfield_in_model_and_base_model(model_help,'get_estimator_from_sim_data',@(sim_data)get_estimator_from_sim_data_mixed(sim_data,1));
    model_help = setfield_in_model_and_base_model(model_help,'get_dofs_from_sim_data', @(sim_data)get_dofs_from_sim_data_for_sensitivity_extension(sim_data,1));
    model_help = setfield_in_model_and_base_model(model_help,'get_rb_size', @(model, detailed_data)lin_evol_opt_get_rb_size(model, detailed_data,p));
    model_help = setfield_in_model_and_base_model(model_help,'get_rb_from_detailed_data',@(detailed_data)lin_evol_opt_get_rb_from_detailed_data_separate_bases(detailed_data,p));
    
%     model_help.set_rb_in_detailed_data = @(detailed_data,RB)lin_evol_opt_set_rb_in_detailed_data_separate_bases(detailed_data,RB,p);
%     model_help.get_estimator_from_sim_data = @(sim_data)get_estimator_from_sim_data_mixed(sim_data,1);
%     model_help.get_dofs_from_sim_data = @(sim_data)get_dofs_from_sim_data_for_sensitivity_extension(sim_data,1);
%     model_help.get_rb_size = @(model, detailed_data)lin_evol_opt_get_rb_size(model, detailed_data,p);
%     model_help.get_rb_from_detailed_data=@(detailed_data)lin_evol_opt_get_rb_from_detailed_data_separate_bases(detailed_data,p);
%     model_help.separate_bases_derivative_nr = p;
    %model.optimization.params_to_optimize = [model.optimization.params_to_optimize,1];
    %model.compute_derivative_indices = zeros(1,size(compute_derivative_indices_old,2));
    %model.compute_derivative_indices(indices_to_compute(p)) = 1;
    model_help.optimization.params_to_optimize = zeros(1,nr_parameters);
    model_help.optimization.params_to_optimize(p) = 1;
    if isfield(model_help,'base_model')&&isfield(model_help.base_model,'optimization')
        model_help.base_model.optimization.params_to_optimize = zeros(1,nr_parameters);
        model_help.base_model.optimization.params_to_optimize(p) = 1;
    end
    
    detailed_data_buf = lin_evol_gen_detailed_data(model_help, model_data);
    RB_der_info{p} = detailed_data_buf.RB_info;
    RB_der{p} = detailed_data_buf.RB_der{p};
    %model_data.RB_der{p} = RB_der{p};
    filecache_clear;
end
matlabpool close
detailed_data.RB_info = RB_info;
detailed_data.RB_der_info = RB_der_info;
detailed_data.RB = RB_finished;
detailed_data.RB_der = RB_der;

model.compute_derivative_info = compute_derivative_info_old;
%model.compute_derivative_indices = compute_derivative_indices_old;

function model = setfield_in_model_and_base_model(model,field, entry)
model = setfield(model,field,entry);
if isfield(model,'base_model')&&isfield(model.base_model,field)
    model.base_model = setfield(model.base_model,field,entry);
end