Automatic massively parallel Markov chain Monte Carlo with quantifiable error