Optimisation of multiple time-step hybrid Monte Carlo Wang–Landau simulations in the isobaric–isothermal ensemble for the determination of phase equilibria