I'm creating a program to compute the LU factorization of a matrix, and I was wondering if there was a way to compute an LU factorization without using Gaussian Elimination.
I'm mainly worried about floating point errors creeping in during Gaussian elimination, as an error in computing one row would then affect all subsequent operations with that row.
Thanks for your time!