test_etaqa_vs_ldes_fcfs.m

% Test ETAQA FCFS departure process: SolverMAM dec.mmap vs SolverLDES
%
% Erl(5)/Erl(2)/1 -> Erl(3)/1 FCFS tandem
% Compares MAM (dec.mmap with ETAQA departures) against LDES simulation
 
model = Network('ETAQA-FCFS-Tandem');
 
source = Source(model, 'Source');
queue1 = Queue(model, 'Queue1', SchedStrategy.FCFS);
queue2 = Queue(model, 'Queue2', SchedStrategy.FCFS);
sink   = Sink(model, 'Sink');
 
oclass = OpenClass(model, 'Class1');
source.setArrival(oclass, Erlang.fitMeanAndOrder(2, 5));   % mean=2, SCV=0.2
queue1.setService(oclass, Erlang.fitMeanAndOrder(0.8, 2)); % rho1=0.4
queue2.setService(oclass, Erlang.fitMeanAndOrder(1.2, 3)); % rho2=0.6
 
model.link(Network.serialRouting(source, queue1, queue2, sink));
 
%% Solve
solverMMAP = SolverMAM(model, 'method', 'dec.mmap');
T_mmap = solverMMAP.getAvgTable();
 
solverSRC = SolverMAM(model, 'method', 'dec.source');
T_src = solverSRC.getAvgTable();
 
solverMVA = SolverMVA(model);
T_mva = solverMVA.getAvgTable();
 
solverLDES = SolverLDES(model, 'seed', 23000, 'samples', 5e5);
T_ldes = solverLDES.getAvgTable();
 
%% Display
fprintf('\n=== ETAQA FCFS: MAM vs LDES ===\n');
fprintf('%-10s %10s %10s %10s %10s\n', 'Station', 'LDES', 'MVA', 'dec.source', 'dec.mmap');
 
fprintf('--- Queue Lengths ---\n');
for i = 1:height(T_mva)
    stn = string(T_mva.Station(i));
    if contains(stn, 'Queue')
        fprintf('%-10s %10.4f %10.4f %10.4f %10.4f\n', stn, ...
            T_ldes.QLen(i), T_mva.QLen(i), T_src.QLen(i), T_mmap.QLen(i));
    end
end
 
fprintf('--- Response Times ---\n');
for i = 1:height(T_mva)
    stn = string(T_mva.Station(i));
    if contains(stn, 'Queue')
        fprintf('%-10s %10.4f %10.4f %10.4f %10.4f\n', stn, ...
            T_ldes.RespT(i), T_mva.RespT(i), T_src.RespT(i), T_mmap.RespT(i));
    end
end
 
fprintf('--- Utilizations ---\n');
for i = 1:height(T_mva)
    stn = string(T_mva.Station(i));
    if contains(stn, 'Queue')
        fprintf('%-10s %10.4f %10.4f %10.4f %10.4f\n', stn, ...
            T_ldes.Util(i), T_mva.Util(i), T_src.Util(i), T_mmap.Util(i));
    end
end