This paper discussed the computational accuracy of rough-surface point-contact isothermal elastohydrodynamic lubrication (EHL) analysis by investigating the effects of differential scheme, viscosity-pressure, and shear-thinning models. An EHL experiment with multitransverse ridges was employed as simulated target. Four differential schemes, including the combined and the separate first-order and second-order backward schemes, were investigated. It is found that the separate second-order backward scheme offers the best results based on the comparison with the experimental data, with which two roughness derivatives may be fully or partially canceled each other; thus, the discretization error induced by roughness can be reduced. The consistency of differential schemes is an important issue for the separate schemes. The Yasutomi free-volume viscosity-pressure model and the Eyring rheological model are found to yield the numerical simulations the closest to experimental results.