Computing the starting state for Gibbs-Duhem integration