We present a blue noise sampling algorithm that is parallel, runs fast on a GPU, exhibits excellent spectrum, works in arbitrary dimensions, requires no pre-computed dataset, and allows adaptive sampling.