runAnalyzer.m

function [runtime, analyzer] = runAnalyzer(self, options)
% RUNANALYZER Execute the bound-analysis solver.
%
% Dispatches noniterative bound families to the shared bound handler
% solver_ba_analyzer, and hierarchical/iterative families to their
% SolverBA-native analyzers. Each method returns a single (upper or lower)
% bound; use getBounds() to obtain the {lower,upper} bracket for a family.
%
% Copyright (c) 2012-2026, Imperial College London
% All rights reserved.
 
if nargin < 2
    options = self.getOptions;
end
 
sn = self.getStruct;
 
if sn.nclosedjobs <= 0
    line_error(mfilename, 'SolverBA supports closed queueing networks only.');
end
 
% 'default' selects the tightest available noniterative upper throughput
% bound (geometric bound), matching the historical SolverMVA convention.
% 'qr' is a friendly alias for the QRF quadratic reduction bound. 'lr' is the
% LP-based Linear Reduction bound (mapqn_bnd_lr_pf, solved by linprog) and is
% NOT an alias of 'qrf.mmi.linear': that method's name refers only to its
% explicit Aeq/beq constraint representation, while its objective is the
% nonlinear MEM entropy solved by fmincon. Bare 'lr' means 'lr.upper'.
method = options.method;
if strcmp(method,'default')
    method = 'gb.upper';
elseif strcmp(method,'lr')
    method = 'lr.upper';
elseif strcmp(method,'qr')
    method = 'qrf.mmi';
end
 
% Show library attribution (QRF uses the Optimization Toolbox) once.
if options.verbose ~= VerboseLevel.SILENT && ~GlobalConstants.isLibraryAttributionShown()
    libs = SolverBA.getLibrariesUsed([], setfield(options,'method',method)); %#ok<SFLD>
    if ~isempty(libs)
        line_printf('The solver will leverage %s.\n', strjoin(libs, ', '));
        GlobalConstants.setLibraryAttributionShown(true);
    end
end
 
iter = 1;
if startsWith(method,'qrf')
    % QRF (Quadratic Reduction Framework) LP-based bounds for single-class
    % closed networks with PH service. The adapter enforces its own gating.
    analyzer = @(qn) solver_ba_qrf_analyzer(qn, setfield(options,'method',method)); %#ok<SFLD>
    bopts = options; bopts.method = method;
    [QN,UN,RN,TN,CN,XN,runtime] = solver_ba_qrf_analyzer(sn, bopts);
elseif any(strcmp(method, self.listValidMethods()))
    analyzer = @(qn) solver_ba_analyzer(qn, setfield(options,'method',method)); %#ok<SFLD>
    bopts = options; bopts.method = method;
    [QN,UN,RN,TN,CN,XN,lG,runtime,iter] = solver_ba_analyzer(sn, bopts);
else
    line_error(mfilename, ['Unknown bound method ''%s''. Valid methods: %s'], ...
        method, strjoin(self.listValidMethods(), ', '));
end
 
M = sn.nstations;
R = sn.nclasses;
AN = zeros(M,R);
WN = zeros(M,R);
self.setAvgResults(QN,UN,RN,TN,AN,WN,CN,XN,runtime,method,iter);
end