We compute the renormalized vacuum polarization of a massive scalar field in the Hartle-Hawking state on (2+1)-dimensional rotating, spacelike stretched black hole solutions to topologically massive gravity, surrounded by a Dirichlet mirror that makes the state well defined. The Feynman propagator is written as a mode sum on the complex Riemannian section of the spacetime, and a Hadamard renormalization procedure is implemented by matching to a mode sum on the complex Riemannian section of a rotating Minkowski spacetime. No analytic continuation in the angular momentum parameter is invoked. Selected numerical results are given, demonstrating the numerical efficacy of the method. We anticipate that this method can be extended to wider classes of rotating black hole spacetimes, in particular to the Kerr spacetime in four dimensions.