A simple enough question: to do a multidimensional integral, given that one has decided that some sort of Monte Carlo method is appropriate, is there any advantage that a regular MC integration using pseudorandom numbers has over a quasi-Monte Carlo integration using a quasirandom sequence? If so, how would I recognize situations where this advantage would come into play? (And if not, why does anyone ever use plain old Monte Carlo integration?)
Tell me more
×
Computational Science Stack Exchange is a question and answer site for
scientists using computers to solve scientific problems. It's 100% free, no registration required.
|
|
|
Monte Carlo simulations are the method of choice for computation of electron scattering. Tricks like importance sampling are used sometimes, so you might say it's not plain old Monte Carlo. But the main point is probably that an inherently stochastic process is simulated here, while you are asking only about using Monte Carlo for integration. Because nobody else tried to offer an answer, let me try to expand my answer a bit. Assume we have an electron scattering simulation, where only a single number, like a backscattering coefficient, is computed. If we would reformulate this as a multidimensional integral, it would probably be an infinite dimensional integral. On the other hand, during the simulation of a single trajectory, only a finite number of random numbers is required (this number can become quite large, if secondary electron generation is taken into account). If we would use a quasirandom sequence like latin hypercube sampling, we would have to use an approximation with a fixed number of dimensions, and generate a random number for every dimension for each sample point. So I think the difference is whether some kind of high dimensional unit-hypercube is sampled, versus an infinite dimensional probability cloud around the origin. |
||||
|
|
|
Some of my research involves solving large scale stochastic partial differential equations. In which case, traditional monte carlo approximation of the integral of interest converges too slowly for it to be worthwhile in a practical sense... i.e. I don't want to have to run 10 times more simulations just to get a decimal point more accuracy to the integral. Instead, I tend to use other methods like sparse smolyak grids because they offer better accuracy in fewer function evaluations. This is only possible because I can assume a certain degree of smoothness in the function. It is reasonable to conjecture that if you expect that the function you are integrating to have certain structure (like smoothness), it would be best to use the quasi-monte carlo scheme that exploits it. If you really can't make very many assumptions about the function, then monte carlo is the only tool I can think of to handle it. |
|||
|
|
