TD-CARMA: Painless, accurate, and scalable estimates of gravitational-lens time delays with flexible CARMA processes