Efficient numerical solver of the thermohydrodynamic Reynolds equation for journal bearings