comsol_evol_gen_model_data.m

function model_data = comsol_evol_gen_model_data(model)
%function model_data = comsol_stat_gen_model_data(model)
% 
% This function which extracts and computes the model_data from the comsol_model
% - comsol_model
% - the grid
% - the operators needed for detailed_simulation
% - inner-product matrices
% - Index Vectors
%
% Oliver Zeeb, 04.06.13

if model.use_comsol %let comsol do all the work!
    comsol_model = mphload(model.comsol_model_name); 
    comsol_model = comsol_set_mu_from_rbmatlab_model(model,comsol_model);
    
    %disable scaling so that alle the extracted matrices can be combined.
    comsol_model.sol(model.comsol_tags.sol).feature('v1').set('scalemethod', 'none'); 

    model_data.comsol_model = comsol_model;
    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);
    model_data.operators.init_values =  mphgetu(comsol_model,'solnum',1);
    model_data.comsol_tags = model.comsol_tags;

    %extract the operators from COMSOL and store them in model_data
    model.decomp_mode = 1;
    [L_I_comp, L_E_comp, b_comp] = model.operators_ptr(model,model_data);
    model_data.operators.L_I_comp = L_I_comp;
    model_data.operators.L_E_comp = L_E_comp;
    model_data.operators.b_comp = b_comp;
    if model.compute_output_functional
        [v, s_l2norm] = model.operators_output(model, model_data);
        model_data.operators.operators_output = v;
        model_data.operators.output_operator_l2norm = s_l2norm;
    end

    %compute index vectors
    str = mphmatrix(comsol_model,'sol1','Out',{'Null'},'initmethod','init','initsol','zero');
    model_data.ind_vectors.Null = str.Null;
    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);
    
    % inner product matrices:
    model_data.inner_product_matrices = comsol_calculate_inner_product_matrices(model,model_data);
    

    
else %do the work by hand --> use only matlab! --> load model data from file
    load([model.comsol_model_name' '_model_data']);
    
end