dom_dec_operators.m

Go to the documentation of this file.

function [A, r, A_gamma, r_gamma] = dom_dec_operators(model, ...
                                                  model_data)
% function [A, r, A_gamma, r_gamma] =
% dom_dec_operators(model,model_data)
%
% function computing all required matrizes  for the simulation
%
% Required fields of model:
%  decomp_mode: if 0 or 1, matrizes of the detailed systems are
%  assembled, if 2 , coefficients for the matrix of the reduced
%  system are computed

% I. Maier, 19.07.2011

model.base_model.decomp_mode = model.decomp_mode;

A = cell(1,2);
r = cell(1,2);
for i = 1:2

  if isfield(model_data,'operators')
    [A{i}, r{i}] = model_data.operators(i);
  else
    
    if model.decomp_mode < 2
      model_data.df_info = model_data.df_infos{i};
    end;
    [A{i}, r{i}] = model.base_model.operators(model.base_model, ...
                                              model_data);
  end;

end;

A_gamma = cell(1,2);
r_gamma = cell(1,2);

% modification on the interface
if model.decomp_mode == 0
  
  for i = 1:2
    A_gamma{i} = A{i};
    r_gamma{i} = r{i};
    % delete and insert dirichlet rows:
    gids = model_data.gamma_dofs{i};
    ndofs = model_data.df_infos{i}.ndofs;
    
    A_tmp = sparse(gids,gids,ones(size(gids)),ndofs,ndofs);
    
    A_gamma{i}(gids,:) = 0;
    r_gamma{i}(gids) = 0;
    
    A_gamma{i} = A_gamma{i} + A_tmp;
  end;
  
end;