Fast calculation of boundary crossing probabilities for Poisson processes

The boundary crossing probability of a Poisson process with $n$ jumps is a fundamental quantity with numerous applications. We present a fast $O(n^2 \log n)$ algorithm to calculate this probability for arbitrary upper and lower boundaries.