Sampling-based inference for large linear models, with application to linearised Laplace