comsol_stat_detailed_simulation.m

function sim_data = comsol_stat_detailed_simulation(model,model_data)
%function sim_data = comsol_stat_detailed_simulation(model,model_data)
%
% function performing a detailed simulation using the comsol_model.
% then sets the result to a rbmatlab structure sim_data
%
% Oliver Zeeb, 01.09.2012

if model.use_comsol %let comsol do all the work!
        comsol_model = model_data.comsol_model;
        comsol_model = comsol_set_mu_from_rbmatlab_model(model,comsol_model); %set the mu in the comsol model

        %activate stationary solver ( if deactivated in COMSOL)
        comsol_model.sol(model.comsol_tags.sol).feature('s1').active(true);

        % execute COMSOL Simulation
        comsol_model.sol(model.comsol_tags.sol).runAll;

        %get solution from comsol:
        sim_data.U = mphgetu(comsol_model);
        sim_data.comsol_model=comsol_model;

        % compute output
        if model.compute_output_functional
          sim_data.s = comsol_model.result().numerical(model.comsol_tags.output).getReal();
        end;
        
else  %do the work by hand --> use only matlab!
    %     1. get operator matrices --> they are in model data!
    %     2. get parameter
    %     3. linear combination of 1. and 2
    %     4. Solve the linear system
    %     5. blow up the solution vector
    %     (6. add particular solution ud)
    %     7. compute output functional with the full solution!
    
    model.comsol_get_eliminated_data = 1;
    % 1.
    model.decomp_mode = 1;
    [A_comp, f_comp] = model.operators(model,model_data);
    
    % 2.
    model.decomp_mode = 2;
    [A_coeff,f_coeff] = model.operators(model,model_data);
    
    % 3.
    Ac = lincomb_sequence(A_comp, A_coeff);
    fc = lincomb_sequence(f_comp, f_coeff);
    
    % 4.
    u = Ac\fc;
    
    % 5. blow up the solution to the full dimension:
    sim_data.U = model_data.ind_vectors.Null*u;
    
    % 6. add particular solution ud
    %sim_data.U = sim_data.U + model_data.operators.ud;
    
    % 7.
    if model.compute_output_functional
        model.comsol_get_eliminated_data = 0;
        model.decomp_mode = 1;
        output_comp = model.operators_output(model,model_data);
        model.decomp_mode = 2;
        output_coeff = model.operators_output(model,model_data);
        output_operator = lincomb_sequence(output_comp, output_coeff);
        sim_data.s = output_operator' * sim_data.U;
        sim_data.s_with_ud = output_operator' * (sim_data.U + model_data.operators.ud);
    end
    
    
end %of model.use_comsol
    


%add some fields to sim_data
sim_data.ind_vectors = model_data.ind_vectors;