High-precision cross-section data of the deuteron-proton breakup reaction at 130 MeV are presented for 72 kinematically complete configurations. The data cover a large region of the available phase space, divided into a systematic grid of kinematical variables. They are compared with theoretical predictions, in which the full dynamics of the three-nucleon (3N) system is obtained in three different ways: realistic nucleon-nucleon (NN) potentials are combined with model 3N forces (3NF's) or with an effective 3NF resulting from explicit treatment of the Δ-isobar excitation. Alternatively, the chiral perturbation theory approach is used at the next-to-next-to-leading order with all relevant NN and 3N contributions taken into account. The generated dynamics is then applied to calculate cross-section values by rigorous solution of the 3N Faddeev equations. The comparison of the calculated cross sections with the experimental data shows a clear preference for the predictions in which the 3NF's are included. The majority of the experimental data points are well reproduced by the theoretical predictions. The remaining discrepancies are investigated by inspecting cross sections integrated over certain kinematical variables. The procedure of global comparisons leads to establishing regularities in disagreements between the experimental data and the theoretically predicted values of the cross sections. They indicate deficiencies still present in the assumed models of the 3N system dynamics.