Yes, I think we should call this function instead of using our own implementation. Before I thought our QL method is used to solve tridiagonal matrix only so it maybe faster, but now it seems this doesn't matter since the tridiagonal matrix is really small so probably we should use numpy.linalg.eig function instead of our own.