  We review recent progress on the calculations on the inclusive forward hadron
production within the saturation formalism. After introducing the concept of
perturbative parton saturation and nonlinear evolution we discuss the formalism
for the forward hadron production at high energy in the leading and
next-to-leading order. Numerical results are presented and compared with the
experimental data on forward hadron production in $dA$ and $pA$. We discuss the
problem of the negativity of the NLO cross section at high transverse momenta,
study its origin in detail and present possible improvements which include the
corrected kinematics and the suitable choice of the rapidity cutoff.
