The rapidity (energy) evolution of  hadronic observables in  scattering of a dilute perturbative projectile on a dense gluonic target is described in QCD by the JIMWLK equation. This is a functional non-linear equation that is consistent with the QCD unitarity. JIMWLK reduces to the linear BFKL equation when the scattering probability is low.