comsol_get_init_values.m

function init_values = comsol_get_init_values(model,detailed_data)
%function init_values = comsol_get_init_values(model,detailed_data)
%
% function extracting the init_values from a comsol model. After an
% time-dependend simulation, these values have to be added to the solution.
% This is due to the way COMSOL assembles the matrices!
%
% Oliver Zeeb, 2013

%disp('comsol_get_init_values --> my_init_values!!! decomp_mode sauber implementieren!!!');
if model.use_comsol
    switch model.decomp_mode
        case 0
            init_values = mphgetu(detailed_data.comsol_model,'solnum',1);
            init_values = init_values(detailed_data.ind_vectors.com_DOF_ind);
        case 1
            init_values = mphgetu(detailed_data.comsol_model,'solnum',1);
            init_values = {init_values(detailed_data.ind_vectors.com_DOF_ind)};
            %cut off the DBC_DOFs and put everything in a cell
        case 2
            init_values = 1;
    end
else %~model.use_comsol
    switch model.decomp_mode
        case 0
            disp('decomp_mode 0 not yet implemented');
            keyboard;
        case 1
            model.decomp_mode = 1;
            [L_I_comp, ~, b_comp] = model.operators_ptr(model,detailed_data);
            model.decomp_mode = 2;
            [L_I_coeff, ~, b_coeff] = model.operators_ptr(model,detailed_data);
            L_I = lincomb_sequence(L_I_comp, L_I_coeff);  
            b = lincomb_sequence(b_comp, b_coeff);
            init_values = detailed_data.ind_vectors.Null*(L_I \ b) ;
            if model.comsol_get_eliminated_data
                init_values = init_values(detailed_data.ind_vectors.com_DOF_ind);
            end
            init_values = {init_values};
        case 2
            init_values = 0;
    end
end