pub fn lu_solve<T: ComplexField>( a: &Array2<T>, b: &Array1<T>, ) -> Result<Array1<T>, LuError>
Solve Ax = b using LU decomposition
This is a convenience function that combines factorization and solve.