  Saturation physics is expected to be relevant at sufficiently small parton
momentum fractions $x$ in high-energy proton- (or deuteron-)ion collisions at
RHIC and the LHC. Accordingly, these collisions provide the best available
testing ground for the saturation model. However, producing precise numerical
predictions from the model is a complicated task; the state of the art in this
area involves next-to-leading order QCD calculations, which are difficult to do
numerically. Here I'll review recent progress in extracting numerical
predictions from saturation models and matching them to experimental results.
