Numerical computation --- such as numerical solution of a PDE, or quadrature --- can modelled as a statistical inverse problem in its own right. In particular, we can apply the Bayesian approach to inversion, so that a posterior distribution is induced over the object of interest (e.g. the PDE's solution) by conditioning a prior distribution on the same finite information that would be used in a...