We present a massively parallel implementation of time-dependent density functional theory in real space, aimed at computing optical absorption spectra of realistic systems with hundreds of atoms from first principles. We provide details of the formalism and discuss its implementation, optimization, and efficient parallelization, as well as remaining limitations, in detail. The capabilities of the code are illustrated by calculations of optical properties of hydrogenated silicon quantum dots.