Computes the solution of a linear equation, given its LU decomposition.
x = linalg_solvelu ( L , U , b , verbose )
a n-by-n matrix of doubles, lower triangular.
a n-by-n matrix of doubles, upper triangular.
a n-by-1 matrix of doubles
a 1-by-1 matrix of boolean, set to true to display intermediate messages
a n-by-1 matrix of doubles
Solve Ax=b with A=L*U and L lower triangular and U upper triangular. Effective algorithm with vectorization.