A scalable, numerically stable, high-performance tridiagonal solver using GPUs