FDTD with auxiliary bath fields for condensed-phase polaritonics: Fundamentals and implementation (APL Computational Physics)