This article is also available in Japanese.

A source object generates rays with some spatial and angular probability. For example, a Source Point generates rays with exactly the same starting coordinates, but with ray direction cosines uniformly chosen over some allowable range. When we trace rays from this source, we must trace "enough" rays to accurately sample the source.

The probability function for a source is usually sampled with a random number generator. The goal of a random number generator is to produce a series of numbers which have no correlation with each other. Then by tracing lots of rays (generating lots of random numbers) we can accurately sample the source probability function.

{All computer-based random-number generation algorithms are "quasi-random" in thay they are limited (ultimately by the number of bits of the operating system, but generally by algorithm assumptions) to have a period over which the random number sequence repeats. Hence no implemented random number generator is "truly" random. In this article, I refer to ZEMAX' long-period random number generator as "truly" random to distinguish it from the Sobol algorithm.} 

Sobol sampling takes a different approach. Rather than generating random numbers, it generates a uniform distribution in probability space. This is not just a simple grid, but a distribution which appears qualitatively random, but cleverly "fills in" previously unsampled regions of the probability function.

The zip file attached (which can be downloaded from the last page of this article) shows a simple optical system. Two rectangular sources illuminate a detector. These source objects are designed to produce a uniform, rectangular distribution of rays. The bottom one uses ZEMAX' "true" random number generator, whilst the top one uses a Sobol sampling scheme.

If we trace just 10,000 rays per source, it is hard to see the difference between the sampling schemes:



Tracing 100,000 rays shows the differences more clearly



The Sobol sampling scheme produces a more uniform, though still qualitatively random distribution than the truly random distribution. This shows a more "clumpy" distribution, with hot spots and darker  regions. Zooming in on a region in the random-sampled source distribution shows:



The red pixel above indicates where three rays landed in the same pixel, and there are also regions where no rays landed. Hence this region is noisier than the Sobol-illuminated region.