solver_nc_cacheqn_retrieval_analyzer.m

function [QN,UN,RN,TN,CN,XN,lG,hitprob,missprob,delayedprob,hitproblist,itemprob,latency,runtime,method] = solver_nc_cacheqn_retrieval_analyzer(sn, options)
% [...] = SOLVER_NC_CACHEQN_RETRIEVAL_ANALYZER(SN, OPTIONS)
%
% NC analyzer for a CLOSED integrated cache-queueing model whose Cache node has a
% delayed-hit retrieval system. Delegates to da_cacheqn_retrieval with an NC
% network solve (load-dependent when the coalescing fetch station is present).
% Returns the true cache hit/miss probabilities (hit = P(item cached),
% miss = 1 - hit); the delayed-hit fraction folds into miss (delayedprob = 0).
%
% Copyright (c) 2012-2026, Imperial College London
% All rights reserved.
 
T0 = tic;
[res, hitprob, missprob, delayedprob, ~, ~] = da_cacheqn_retrieval(sn, @netsolve, options);
 
QN = res.QN; UN = res.UN; RN = res.RN; TN = res.TN; XN = res.XN;
if isfield(res,'CN'), CN = res.CN; else, CN = NaN(1, sn.nclasses); end
if isfield(res,'lG'), lG = res.lG; else, lG = NaN; end
 
K = sn.nclasses;
hitproblist = NaN(K, numel(sn.nodeparam{find(sn.nodetype==NodeType.Cache,1)}.itemcap));
itemprob = [];
latency = NaN(1, K);
method = 'fpi';
runtime = toc(T0);
 
    function res = netsolve(snit)
        res = struct();
        if ~isempty(snit.lldscaling) || ~isempty(snit.cdscaling)
            [res.Q,res.U,res.R,res.T,res.C,res.X,res.lG,res.runtime] = solver_ncld_analyzer(snit, options);
        else
            [res.Q,res.U,res.R,res.T,res.C,res.X,res.lG,res.runtime] = solver_nc_analyzer(snit, options);
        end
        res.QN = res.Q; res.UN = res.U; res.RN = res.R; res.TN = res.T; res.XN = res.X; res.CN = res.C;
    end
end