My goal is to transform a matrix into upper triangular form in Python. I know the function
scipy.linalg.lu will do LU decomposition and get both upper and lower triangular ones. I need to repeat this procedure several times. If there is functions from
lapack that compute only upper triangular matrix, the performance of my algorithm will be improved significantly.
Could you please suggest such functions? Thank you so much!