Compute an `LU` factorization of a real tridiagonal matrix `A` using elimination with partial pivoting and row interchanges - View it on GitHub