https://hgpu.org/?p=8309
Orthogononalization on a general purpose graphics processing unit with double double and quad double arithmetic