https://hgpu.org/?p=3639
Parallel calculation of the median and order statistics on GPUs with application to robust regression