Fast calculation of boundary crossing probabilities for Poisson processes

We present a fast $O(n^2 \log n)$ algorithm for calculating the probability that a one-dimensional Poisson process will stay within arbitrary boundaries that are bounded by $n$. This algorithm is faster than previous $O(n^3)$ methods, and can be used to compute $p$-values for continuous goodness-of-fit statistics.