https://hgpu.org/?p=1665
Triangular matrix inversion on Graphics Processing Unit