Scalable inference for crossed random effects models