pfqn_gldsingle.m

%{
%{
 % @file pfqn_gldsingle.m
 % @brief Exact normalizing constant for single-class load-dependent models.
%}
%}
 
%{
%{
 % @brief Exact normalizing constant for single-class load-dependent models.
 % @fn pfqn_gldsingle(L, N, mu, options)
 % @param L Service demand vector (Mx1).
 % @param N Population (scalar).
 % @param mu Load-dependent rate matrix (MxN).
 % @param options Solver options.
 % @return lG Logarithm of normalizing constant.
 % @return G Normalizing constant.
%}
%}
function [lG,G]=pfqn_gldsingle(L,N,mu,options)
% G=PFQN_GLDSINGLE(L,N,MU)
 
if nargin<4
    options = [];
end
 
[M,R]=size(L);
if R>1
    line_error(mfilename,'multiclass model detected. pfqn_gldsingle is for single class models.');
end
Nscal = N(1); % codegen: ensure scalar loop bound
g = zeros(M+1, Nscal+1, Nscal+2);
for n=1:Nscal
    g(0 +1,n +1, 1 +1)=0;
end
for m=1:M
    for tm=1:(Nscal+1)
        g(m +1,0 +1,tm +1)=1;
    end
    for n=1:Nscal
        for tm=1:(Nscal-n+1)
            g(m +1, n +1, tm +1)= g(m-1 +1, n +1, 1 +1)+L(m)*g(m +1, n-1 +1, tm+1 +1)/mu(m,tm);
        end
    end
end
G = g(M +1,Nscal +1,1 +1);
lG = log(G);
end