lin_evol_rb_derivative_operators_coefficients.m

function [coeff_LL_I, coeff_LL_E, coeff_bb] =  lin_evol_rb_derivative_operators_coefficients(model)
%function [LL_I, LL_E, bb, K_II, K_IE, K_EE, m_I, m_E, m] = ...
%    lin_evol_rb_derivative_operators_coefficients(model, [detailed_data])
%
% Write docu!!
%
% Markus Dihlmann 01.06.2010


  %!!!!!!!!!!!!!!
  %Only valable for advection problem !!!!
  %!!!!!!!!!!!!!!!
    coeff_LL_I = model.rb_derivative_operator_coefficients_Li(model);
    coeff_LL_E = -model.rb_derivative_operator_coefficients_Le(model);
    coeff_bb_E = model.rb_derivative_operator_coefficients_b(model);
    coeff_bb_I = model.rb_derivative_operator_coefficients_bI(model); %b from implicit operator
    if isfield(model, 'rb_derivative_operator_coefficients_B_neu_ex')
        %Neumann boundary conditions exist
        coeff_L_E_neu = -model.rb_derivative_operator_coefficients_L_E_neu(model);
        coeff_bb_E_neu = model.rb_derivative_operator_coefficients_B_neu_ex(model);
        coeff_bb = [coeff_bb_I,coeff_bb_E, coeff_bb_E_neu];
        coeff_LL_E = [coeff_LL_E,coeff_L_E_neu];
    else
        %without neumann
        coeff_bb = [coeff_bb_I,coeff_bb_E];
    end
  % pairs: products of sigmas:
  %K_II = repmatrows(L_I(:), Q_L_I) .*  repmat(L_I(:), Q_L_I, 1); 
  %K_IE = repmatrows(L_I(:), Q_L_E) .*  repmat(L_E(:), Q_L_I, 1); 
  %K_EE = repmatrows(L_E(:), Q_L_E) .*  repmat(L_E(:), Q_L_E, 1); 
  %m_I  = repmatrows(L_I(:), Q_b)   .*  repmat(b(:)  , Q_L_I, 1); 
  %m_E  = repmatrows(L_E(:), Q_b)   .*  repmat(b(:)  , Q_L_E, 1); 
  %m    = repmatrows(b(:)  , Q_b)   .*  repmat(b(:)  , Q_b  , 1); 
  
end

%| \docupdate