A mixed Monte Carlo/PDE variance reduction method under the Heston-CIR model