Multi-dimensional Transform Inversion with Applications to the Transient M/g/1 Queue

We develop an algorithm for numerically inverting multi-dimensional transforms. Our algorithm applies to any number of continuous variables (Laplace transforms) and discrete variables (generating functions). We use the Fourier-series method; i.e., the inversion formula is the Fourier series of a periodic function constructed by aliasing. This amounts to an application of the Poisson summation formula. By appropriately exponentially damping the given function, we control the aliasing error. We choose the periods of the multi-dimensional periodic function so that each infinite series is a finite sum of nearly alternating infinite series; then we apply the Euler transformation to compute the infinite series from finitely many terms. The multidimensional inversion algorithm enables us, evidently for the first time, to quickly and accurately calculate probability distributions from several classical transforms in queueing theory. For example, we apply our algorithm to invert the two-dimensional transforms of the joint distribution of the duration of a busy period and the number served in that busy period, and the time-dependent of the transient queue-length and workload distributions, in the M/G/1 queue. In other related work, we have applied the inversion algorithms here to calculate time-dependent distributions in the transient BMAP/G/1 queue (with a batch Markovian arrival process) and the piecewise-stationary M t / G t /1 queue.