Почему метод наименьших квадратов не работает для матрицы высокого, а не полного ранга?

У меня есть матрица A формы (224, 45). Это ранг 44. Когда я пробую этот код:

solution = np.linalg.lstsq(A, rhs)[0]

Я получаю решение с очень высокими значениями, что-то вроде 1e14.

Когда я пытаюсь

solution = np.linalg.lstsq(A.T.dot(A), A.T.dot(rhs))[0]

Я получаю хорошее решение, как и ожидалось.

Вопрос в том, нужно ли мне использовать np.linalg.lstsq только для квадратных симметричных матриц или что-то в этом роде? В документах об этом ничего нет.