Efficient implementation of numerical methods for solving bearing cavitation problems using symmetric system solvers