La réduction de la ligne de Lapack
Question
Je suis en train d'écrire une fonction qui produit une solution unique à un système d'équations sous-représentés (par exemple la matrice qui décrit le système est plus large que haut). Pour ce faire, je l'ai cherché dans la documentation LAPACK un moyen de réduire les maisons en rangée d'une matrice à sa forme échelonnée réduite, similaire à la Rref () fonction dans les deux calculateurs Mathematica et TI. Le plus proche je suis tombé était http: / /software.intel.com/en-us/forums/intel-math-kernel-library/topic/53107/ ce petit fil. Ce fil, cependant, semble impliquer que prenant simplement la matrice triangulaire supérieure « U » (et en divisant chaque ligne par la diagonale) est la même que la forme échelonnée réduite d'une matrice, que je ne pense pas que ce soit le cas. Je pourrais coder Rref () moi-même, mais je ne crois pas que je pourrais atteindre les performances LAPACK est célèbre pour.
1) Y at-il une meilleure façon d'obtenir simplement une solution spécifique à un système sous-représentés?
2) Dans le cas contraire, est-il un moyen pour LAPACK à la ligne-réduire une matrice?
Merci!
La solution
Une méthode souvent utilisée pour c'est la solution la moins carrée, voir la sgelsx de lapack