redução fileira de Lapack

Question

Eu estou tentando escrever uma função que produz uma única solução para um sistema de sub-representadas de equações (por exemplo, a matriz que descreve o sistema é mais largo do que alto). A fim de fazer isso, eu tenho procurado na documentação LAPACK uma maneira de linha-reduzir uma matriz para a sua forma reduzida escalão, semelhante ao rref function () em ambas as calculadoras Mathematica e TI. O mais próximo que eu deparei foi http: / /software.intel.com/en-us/forums/intel-math-kernel-library/topic/53107/ este minúsculo fio. Esta discussão, no entanto, parece implicar que simplesmente tomar o "U" matriz triangular superior (e dividindo cada linha pela diagonal) é o mesmo que a forma escalonada reduzida de uma matriz, o que eu não acredito que seja o caso. Eu poderia codificar-se rref () mesmo, mas eu não acredito que eu poderia alcançar o LAPACK desempenho é famosa.

1) Existe uma melhor maneira de simplesmente obter qualquer solução de um específico para um sistema de sub-representados?
2) Se não, há uma maneira para LAPACK a linha-reduzir uma matriz?

Obrigado!

Answer 1

Um método frequentemente utilizado para isso é a solução menos quadrado, ver LAPACK é sgelsx