Various (0, m, 2)-nets with different spectral profiles produced by our (a) Netshuffle and (b-e) Netoptimize algorithms. In each of the five columns, the plot in the top shows a 1024-point set. The middle-right image shows the corresponding frequency power-spectrum, and the bottom plot shows its radial average, both obtained by averaging 100 periodograms of 4096-point sets, using the standard PSA tool with the default settings. The eight subplots show a 64-points net of the same profile, demonstrating the dyadic net property: that the points are stratified over all the elementary intervals. The top set is similarly stratified over all the 11 possible stratifications. The Netshuffle version in (a) is close to Owen's scrambling, but is free of the grid-like artifacts typically seen in its spectrum. The blue noise profiles are obtained by adjusting the ς parameter of a Gaussian optimization process.
We explore the space of (0, m, 2)-nets in base 2 commonly used for sampling. We present a novel constructive algorithm that can exhaustively generate all nets — up to m-bit resolution — and thereby compute the exact number of distinct nets. We observe that the construction algorithm holds the key to defining a transformation operation that lets us transform one valid net into another one. This enables the optimization of digital nets using arbitrary objective functions. For example, we define an analytic energy function for blue noise, and use it to generate nets with high-quality blue-noise frequency power spectra. We also show that the space of (0, 2)-sequences is significantly smaller than nets with the same number of points, which drastically limits the optimizability of sequences.
In this demo, you may interact with 16-point (0, 4, 2)-nets,
and browse the 232 set of distinct ones, identified by
the first four bits of point coordinates.
We show the points as square blocks to indicate the
degree of freedom in defining the trailing bits while
preserving the identity of the net.
The flow diagram visualizes the Netshuffle algorithm.
The nodes are colored to make it easier to follow the flow, and to
identity the dyadic distribution of each coordinate over the other.
The following interactions are available:
In this demo you can compare rendering results of our Blue Nets sampler to state-of-the-art LD samplers; namely, (0, 2)-Sequence, Sobol, and Halton.
The following interactions are available:
Get for free at ACM digital library
|
|
Source Code ZIP (11 KB) |
|
Pre-Optimized Blue Nets ZIP (70 KB) |
TBA
Thanks to the anonymous reviewers for the valuable comments. Thanks to Mohanad Ahmed for his insightful discussions.