Practical Low-Dimensional Halfspace Range Space Sampling

We develop, analyze, implement, and compare new algorithms for creating varepsilon-samples of range spaces defined by halfspaces which have size sub-quadratic in 1/varepsilon, and have runtime linear in the input size and near-quadratic in 1/varepsilon. The key to our solution is an efficient construction of partition trees. Despite not requiring any techniques developed after the early 1990s, apparently such a result was not ever explicitly described. We demonstrate that our implementations, including new implementations of several variants of partition trees, do indeed run in time linear in the input, appear to run linear in output size, and observe smaller error for the same size sample compared to the ubiquitous random sample (which requires size quadratic in 1/varepsilon). This result has direct applications in speeding up discrepancy evaluation, approximate range counting, and spatial anomaly detection.