a scalable, numerically stable tridiagonal solver implemented in NVidia CUDA