Computational Methods in Coupled Electron-Ion Monte Carlo Simulations