Recently, Quasi-Monte Carlo (QMC) methods have been implemented in a reliable adaptive algorithm. This raises the possibility of combining... Show moreRecently, Quasi-Monte Carlo (QMC) methods have been implemented in a reliable adaptive algorithm. This raises the possibility of combining adaptive QMC with efficiency improvement techniques for independent and identically distributed (IID) Monte Carlo (MC) such as control variates (CV). The challenge for adding CV to QMC is that the optimal CV coefficient for QMC is generally not the same as that for MC. Here we propose a method for imple- menting CV in a reliable adaptive QMC algorithm. One merit of using CV with MC is that theoretically the efficiency is always no worse than vanilla MC. Our method is implemented in an efficient way so that the extra cost for CV is tolerable, and the overall time savings can be substantial. We test our algorithm on various problems including option pricing and mul- tivariate normal probability estimation for dimensions from 4 to 64. The same tests are performed on adaptive QMC algorithm without CV as a comparison. Our results show that with good CV, the cost of adaptive QMC is greatly reduced compared to vanilla QMC. M.S. in Applied Mathematics, July 2016 Show less
