scm_lower_bound.m

Go to the documentation of this file.

function [constant_LB, constant_UB] = scm_lower_bound(model, reduced_data)
%[constant_LB, constant_UB] = scm_lower_bound(model, reduced_data)
%
% function which prepares the inputarguments for the online-phase of the
% SCM using the scm_offline_data via reduced_data.scm_offline_data. Then
% the function performs the online-phase of the scm for the current para-
% meter set in model. Output is the lower and optionally also the upper bound.

% Dominik Garmatter 03.09 2012


% preparing the input for the Online-Phase
mu = get_mu(model);
scm_offline_data = reduced_data.scm_offline_data;
M_alpha = scm_offline_data.M_alpha;
M_plus = scm_offline_data.M_plus;
desired_constant = scm_offline_data.desired_constant;
% preparing the coefficients to the current mu
old_decomp_mode = model.decomp_mode;
model.decomp_mode = 2;
if strcmp(model.rb_problem_type, 'lin_evol') == 1
   [L_I,L_E] = model.operators_ptr(model);
   coeffs = [1;L_E;L_I]; % this equals the form of the components used in the offline stage of the SCM
elseif strcmp(model.rb_problem_type, 'lin_stat') == 1
   coeffs = model.operators(model);
   coeffs = coeffs(:);
else
   error('rb_problem_type unknown!')
end
model.decomp_mode = old_decomp_mode;
% compute \Theta(mu)
if desired_constant == 1
    Theta_mu = coeffs(:);
elseif desired_constant == 2
    Q = size(coeffs,1);
    Theta_mu = zeros(Q^2,1);
    for i = 1:Q
      for j = 1:Q
        Theta_mu((i-1)*Q+j) = coeffs(i)*coeffs(j);
      end
    end
end
% Online-Phase
scm_results = scm_online(mu, Theta_mu, scm_offline_data, M_alpha, M_plus, desired_constant);
% Output
constant_LB = scm_results.Lower_Bound;
if nargout > 1
  constant_UB = scm_results.Upper_Bound;
end