Fast Bayesian Coresets via Subsampling and Quasi-Newton Refinement