https://hgpu.org/?p=25736
Least Squares on GPUs in Multiple Double Precision