Я пытаюсь написать функцию, которая дает единственное решение для недостаточно представленной системы уравнений (например, матрица, которая описывает систему, шире, чем высокая). Чтобы сделать это, я искал в документации LAPACK способ уменьшения строки до формы с уменьшенным эшелоном, аналогично функции rref () в Mathematica и TI-калькуляторах. Самым близким, с которым я столкнулся, была http://software.intel.com/en-us/forums/intel-math-kernel-library/topic/53107/ эта крошечная нить. Эта нить, однако, подразумевает, что простое взятие верхней треугольной матрицы «U» (и деление каждой строки по диагонали) такое же, как и у уменьшенной эшелонированной формы матрицы, что, как я полагаю, не имеет место. Я мог бы написать код для rref () сам, но я не верю, что смог бы добиться производительности, которой славится LAPACK.
1) Есть ли лучший способ просто найти какое-то конкретное решение для недостаточно представленной системы?
2) Если нет, то есть ли у LAPACK способ уменьшить матрицу строк?
Спасибо!