oqn_fourqueues.m

clear node jobclass solver AvgTable
 
%% a large multiclass example with PS and FCFS
model = Network('MyNetwork');
 
%% Block 1: nodes
node{1} = Source(model, 'Source');
node{2} = Queue(model, 'Queue1', SchedStrategy.FCFS);
node{3} = Queue(model, 'Queue2', SchedStrategy.FCFS);
node{4} = Queue(model, 'Queue3', SchedStrategy.PS);
node{5} = Queue(model, 'Queue4', SchedStrategy.FCFS);
node{6} = Sink(model, 'Sink');
 
%% Block 2: classes
jobclass{1} = OpenClass(model, 'Class1', 0);
jobclass{2} = OpenClass(model, 'Class2', 0);
jobclass{3} = OpenClass(model, 'Class3', 0);
 
node{1}.setArrival(jobclass{1}, Exp.fitMean(5.000000)); % (Source of Customers,Class0)
node{1}.setArrival(jobclass{2}, Exp.fitMean(8.000000)); % (Source of Customers,Class1)
node{1}.setArrival(jobclass{3}, Exp.fitMean(7.000000)); % (Source of Customers,Class2)
node{2}.setService(jobclass{1}, Exp.fitMean(0.300000)); % (WebServer,Class0)
node{2}.setService(jobclass{2}, Exp.fitMean(0.500000)); % (WebServer,Class1)
node{2}.setService(jobclass{3}, Exp.fitMean(0.600000)); % (WebServer,Class2)
node{3}.setService(jobclass{1}, Exp.fitMean(1.100000)); % (Storage1,Class0)
node{3}.setService(jobclass{2}, Exp.fitMean(1.300000)); % (Storage1,Class1)
node{3}.setService(jobclass{3}, Exp.fitMean(1.500000)); % (Storage1,Class2)
node{4}.setService(jobclass{1}, Exp.fitMean(2.000000)); % (Storage2,Class0)
node{4}.setService(jobclass{2}, Exp.fitMean(2.100000)); % (Storage2,Class1)
node{4}.setService(jobclass{3}, Exp.fitMean(1.900000)); % (Storage2,Class2)
node{5}.setService(jobclass{1}, Exp.fitMean(1.500000)); % (Storage3,Class0)
node{5}.setService(jobclass{2}, Exp.fitMean(0.900000)); % (Storage3,Class1)
node{5}.setService(jobclass{3}, Exp.fitMean(2.300000)); % (Storage3,Class2)
 
%% Block 3: topology
P = model.initRoutingMatrix(); % initialize routing matrix 
P{1,1}(1,2) = 1; % (Source of Customers,Class0) -> (WebServer,Class0)
P{1,1}(2,3) = 2.500000e-01; % (WebServer,Class0) -> (Storage1,Class0)
P{1,1}(2,4) = 2.500000e-01; % (WebServer,Class0) -> (Storage2,Class0)
P{1,1}(2,5) = 2.500000e-01; % (WebServer,Class0) -> (Storage3,Class0)
P{1,1}(2,6) = 2.500000e-01; % (WebServer,Class0) -> (Out,Class0)
P{1,1}(3,2) = 1; % (Storage1,Class0) -> (WebServer,Class0)
P{1,1}(4,2) = 1; % (Storage2,Class0) -> (WebServer,Class0)
P{1,1}(5,2) = 1; % (Storage3,Class0) -> (WebServer,Class0)
P{2,2}(1,2) = 1; % (Source of Customers,Class1) -> (WebServer,Class1)
P{2,2}(2,3) = 2.500000e-01; % (WebServer,Class1) -> (Storage1,Class1)
P{2,2}(2,4) = 2.500000e-01; % (WebServer,Class1) -> (Storage2,Class1)
P{2,2}(2,5) = 2.500000e-01; % (WebServer,Class1) -> (Storage3,Class1)
P{2,2}(2,6) = 2.500000e-01; % (WebServer,Class1) -> (Out,Class1)
P{2,2}(3,2) = 1; % (Storage1,Class1) -> (WebServer,Class1)
P{2,2}(4,2) = 1; % (Storage2,Class1) -> (WebServer,Class1)
P{2,2}(5,2) = 1; % (Storage3,Class1) -> (WebServer,Class1)
P{3,3}(1,2) = 1; % (Source of Customers,Class2) -> (WebServer,Class2)
P{3,3}(2,3) = 2.500000e-01; % (WebServer,Class2) -> (Storage1,Class2)
P{3,3}(2,4) = 2.500000e-01; % (WebServer,Class2) -> (Storage2,Class2)
P{3,3}(2,5) = 2.500000e-01; % (WebServer,Class2) -> (Storage3,Class2)
P{3,3}(2,6) = 2.500000e-01; % (WebServer,Class2) -> (Out,Class2)
P{3,3}(3,2) = 1; % (Storage1,Class2) -> (WebServer,Class2)
P{3,3}(4,2) = 1; % (Storage2,Class2) -> (WebServer,Class2)
P{3,3}(5,2) = 1; % (Storage3,Class2) -> (WebServer,Class2)
model.link(P);
 
%%
options = lineDefaults;
options.seed = 23000;
options.samples = 1e6;
options.cutoff=1;
solver={};
solver{end+1} = CTMC(model,options); % CTMC is infinite on this model
%solver{end+1} = FLD(model,options);
solver{end+1} = MVA(model,options);
solver{end+1} = MAM(model,options);
solver{end+1} = JMT(model,options);
solver{end+1} = DES(model,options);
%solver{end+1} = SSA(model,options);
%solver{end+1} = NC(model,options);
for s=1:length(solver)
    fprintf(1,'SOLVER: %s\n',solver{s}.getName());
    AvgTable{s} = solver{s}.getAvgTable()
    AvgTable{s}
end