Exact simulation of the stationary distribution of the FIFO M/G/c queue
We present an exact simulation algorithm for the stationary distribution of the customer delayDfor first-in–first-out (FIFO) M/G/cqueues in which ρ=λ/μ<1. We assume that the service time distributionG(x)=P(S≤x),x≥0 (with mean 0<E(S)=1/μ<∞), and its corresponding equilibrium distributionGe(x)=μ∫0xP(S>y)dyare such that samples of them can be simulated. We further assume thatGhas a finite second moment. Our method involves the general method of dominated coupling from the past (DCFTP) and we use the single-server M/G/1 queue operating under the processor sharing discipline as an upper bound. Our algorithm yields the stationary distribution of the entire Kiefer–Wolfowitz workload process, the first coordinate of which isD. Extensions of the method to handle simulating generalized Jackson networks in stationarity are also remarked upon.