Krylov subspace time domain computations of monochromatic sources for multi-frequency optical response