What is the most computationally efficient method to solve large systems of linear equations of the form Ax = b?
In my current approach, I'm using scipy.linalg.lu_factor and scipy.linalg.lu_solve, but I'm not sure if there's a more efficient way to go about this.