A subdivision algorithm for phase equilibrium calculations at high pressures