comsol_stat_gen_model_data.m

function model_data = comsol_stat_gen_model_data(model)
%function model_data = comsol_stat_gen_model_data(model)
%
% Function which extracts the model data from the comsol_model
% - comsol_model
% - the grid
% - operators needed for the detailed solution
% - Index Vectors
%
% Oliver Zeeb, 04.06.13

if model.use_comsol %let comsol do all the work!
        %load comsol_model
        comsol_model=mphload([model.comsol_model_name '.mph']);
        comsol_model = comsol_set_mu_from_rbmatlab_model(model,comsol_model);

        
        % set comsol_model and grid_info in the model_data
        model_data.comsol_model = comsol_model;
        model_data.comsol_model.sol('sol1').feature('s1').active(false);
        comsol_model.mesh(model.comsol_tags.mesh).run;
        comsol_model.sol(model.comsol_tags.sol).runAll;
        model_data.grid_info = mphxmeshinfo(comsol_model);%,'solname', model.comsol_tags.sol);

        %set comsol_model_tags to comsol_data
        model_data.comsol_tags = model.comsol_tags;

        % inner product matrices:
        model_data.comsol_model.sol('sol1').feature('s1').active(false);
        model_data.inner_product_matrices = comsol_calculate_inner_product_matrices(model,model_data);

        %compute index vectors
        model_data.comsol_model.sol('sol1').feature('s1').active(false);
        str = mphmatrix(comsol_model,'sol1','Out',{'Null','ud'},'initmethod','init', 'initsol', 'zero');
        model_data.ind_vectors.Null = str.Null;
        model_data.operators.ud = str.ud;
        model_data.ind_vectors.com_DBC_ind = find(sum(str.Null,2)==0);
        model_data.ind_vectors.com_DOF_ind = find(sum(str.Null,2)~=0);
        model_data.ind_vectors.com_DBC_value = str.ud(model_data.ind_vectors.com_DBC_ind); %ud: particular solution --> DBC_values!
        
        % compute and save operators
        model.decomp_mode = 1;
        [A_comp_full, f_comp_full] = model.operators(model,model_data);
        model_data.operators.A_comp_full = A_comp_full;
        model_data.operators.f_comp_full = f_comp_full;
        for k=1:length(A_comp_full)
            model_data.operators.A_comp_eliminated{k} = A_comp_full{k}(model_data.ind_vectors.com_DOF_ind,model_data.ind_vectors.com_DOF_ind);
        end
        for k=1:length(f_comp_full)
            model_data.operators.f_comp_eliminated{k} = f_comp_full{k}(model_data.ind_vectors.com_DOF_ind);
        end
        
        if model.compute_output_functional
            model_data.operators.output_comp_full = model.operators_output(model,model_data);
            model_data.operators.output_comp_eliminated = model_data.operators.output_comp_full;
            for k=1:length(model_data.operators.output_comp_eliminated) %eliminated version: cut the full operators
                model_data.operators.output_comp_eliminated{k} = model_data.operators.output_comp_eliminated{k}(model_data.ind_vectors.com_DOF_ind);
            end
        end
        
else %do the work by hand --> use only matlab! --> load model data from file
    load([model.comsol_model_name' '_model_data']);
    
end