infer_quick_model.m

function model = infer_quick_model(open, stations, classes, servers, jobs, routing)
%QUICKMODEL Generate simple closed queueing network based on given parameters
    numStations = size(stations,2);
    numClasses = size(classes, 1);
   
    if ~exist('servers','var')
        % servers parameter does not exist, so default it to no concurrency
        servers = ones(numStations);
    end
 
    if ~exist('jobs','var')
        % jobs parameter does not exist, so default it to one
        jobs = ones(numClasses);
    end
 
    model = Network('quickModel');
    if open
        node{1} = Source(model, 'mySource');
        node{numStations+2} = Sink(model, 'mySink');
    end
    
    for i=1:numStations
        node{i} = Queue(model, strcat('QueueStation',num2str(i)), stations{i});
        node{i}.setNumberOfServers(servers(i));
        node{i}.setCapacity(Inf);
    end
    
    P = model.initRoutingMatrix;
    for c=1:numClasses
        if ~open
            jobclass{c} = ClosedClass(model, strcat('Class', num2str(c)), jobs(c), node{1});
 
            if ~exist('routing','var')
                simpleRoutes = eye(numStations);
                P{c} = simpleRoutes([setdiff(1:size(simpleRoutes,1), [1]), [1]], :);
            else
                P{c} = routing(c, :);
            end
        else
            jobclass{c} = OpenClass(model, strcat('Class', num2str(c)));
        end
        
        for i=1:numStations
            node{i}.setService(jobclass{c}, Exp.fitMean(classes(c,i)));
        end
    end
 
    if open
        for c=1:numClasses
            P{c} = Network.serialRouting(node);
        end
    end
 
    model.link(P);
    model.sn.classcap = inf;
end