Scalable variational inference for low-rank spatio-temporal receptive fields