We present a new numerical program, HTurbo, which provides fast and numerically precise predictions for Higgs boson production cross sections. The present version of the code implements the perturbative QCD expansion up to the next-to-next-to-leading order also combined with the resummation of the large logarithmic corrections at small transverse momenta up to next-to-next-to-leading logarithmic accuracy and it includes the Higgs boson production through gluon fusion and decay in two photons with the full dependence on the final-state kinematics. Arbitrary kinematical cuts can be applied to the final states in order to obtain fiducial cross sections and associated kinematical distributions. We present a benchmark comparison with the predictions obtained with the numerical programs HRes and HNNLO programs for which HTurbo represents an improved reimplementation.