Reproducible Triangular Solvers for High-Performance Computing