stokes_detailed_simulation.m

Go to the documentation of this file.

function sim_data = stokes_detailed_simulation(model, model_data)
%function sim_data = stokes_detailed_simulation(model, model_data)
% Performs a detailed simulation of a (Navier-)stokes problem
%
% In general this can also be used for every linear stationary problem.
% Uses a fixed point defect correction algorithm.

% I. Maier, 22.01.2014


% initialization
start_total = tic;

model.decomp_mode = 0;

% assemble linear terms
[A, f] = model_data.operators(model, model_data);

% solution variable
uh = Fem.DiscFunc([], model_data.df_info);

% prepare algorithm
params = cell(1, 4);
if model.has_nonlinearity    
    S(1).type = '.';
    S(1).subs = 'uh';
    S(2).type = '.';
    S(2).subs = 'dofs';
    model.uh = uh;
    params{1} = @(x) model.operators(subsasgn(model, S, x), model_data);
    if isfield(model, 'fp_tolerance')
        params{2} = model.fp_tolerance;
    end
    if isfield(model, 'fp_maxiter')
        params{3} = model.fp_maxiter;
    end
    if isfield(model, 'fp_damping')
        params{4} = model.fp_damping;
    end
end

% execute algorithm
[uh.dofs, defect, niterations] = ...
    VecMat.fixed_point_algorithm(A, f, params{:}, max(model.verbose-2, 0));

% return results:
sim_data.uh = uh;
sim_data.defect = sqrt(defect' * model.get_inner_product_matrix(model_data) * defect);
sim_data.niterations = niterations;
sim_data.total_time = toc(start_total);

end