*** Welcome to piglix ***

Sobol sequence


Sobol sequences (also called LPτ sequences or (ts) sequences in base 2) are an example of quasi-random low-discrepancy sequences. They were first introduced by the Russian mathematician Ilya M. Sobol (Илья Меерович Соболь) in 1967.

These sequences use a base of two to form successively finer uniform partitions of the unit interval and then reorder the coordinates in each dimension.

Let Is = [0,1]s be the s-dimensional unit hypercube, and f a real integrable function over Is. The original motivation of Sobol was to construct a sequence xn in Is so that

and the convergence be as fast as possible.

It is more or less clear that for the sum to converge towards the integral, the points xn should fill Is minimizing the holes. Another good property would be that the projections of xn on a lower-dimensional face of Is leave very few holes as well. Hence the homogeneous filling of Is does not qualify because in lower dimensions many points will be at the same place, therefore useless for the integral estimation.

These good distributions are called (t,m,s)-nets and (t,s)-sequences in base b. To introduce them, define first an elementary s-interval in base b a subset of Is of the form

where aj and dj are non-negative integers, and for all j in {1, ...,s}.

Given 2 integers , a (t,m,s)-net in base b is a sequence xn of bm points of Is such that for all elementary interval P in base b of hypervolume λ(P) = bt−m.


...
Wikipedia

...