Variance reduction for Monte Carlo methods to evaluate option prices under multi-factor stochastic volatility models