4 % Function computing all required matrizes
for the simulation
6 % Required fields of model:
7 % decomp_mode:
if 0 or 1, matrizes of the detailed systems are
8 % assembled,
if 2 , coefficients
for the matrix of the reduced
11 % I. Maier, 19.07.2011
13 model.base_model.decomp_mode = model.decomp_mode;
19 if isfield(model_data,
'operators')
20 [A{i}, r{i}] = model_data.operators(i);
23 if model.decomp_mode < 2
24 model_data.df_info = model_data.df_infos{i};
26 [A{i}, r{i}] = model.base_model.operators(model.base_model, ...
35 % modification on the interface
36 if model.decomp_mode == 0
41 %
delete and insert dirichlet rows:
42 gids = model_data.gamma_dofs{i};
43 ndofs = model_data.df_infos{i}.ndofs;
45 A_tmp = sparse(gids,gids,ones(size(gids)),ndofs,ndofs);
47 A_gamma{i}(gids,:) = 0;
50 A_gamma{i} = A_gamma{i} + A_tmp;
function [ A , r , A_gamma , r_gamma ] = dom_dec_operators(model, model_data)
Function computing all required matrizes for the simulation.