https://hgpu.org/?p=8844
Orthogonalization on a General Purpose Graphics Processing Unit with Double Double and Quad Double Arithmetic