A mixed Monte Carlo and PDE variance reduction method for foreign exchange options under the Heston-CIR model