Lapack der Reihe Reduktion

Question

Ich versuche, eine Funktion zu schreiben, das eine einzige Lösung für ein unterrepräsentierte Gleichungssystem erzeugt (zum Beispiel der Matrix, die das System beschreibt, ist breiter als hoch ist). Um dies zu tun, ich habe für eine Art und Weise in der LAPACK Dokumentation gesucht eine Matrix ist es reduziert gestaffelte Form, ähnlich die Funktion Rref () in beiden Mathematica und TI-Rechner rudern-zu verringern. Der nächstgelegene ich kam über war http: / /software.intel.com/en-us/forums/intel-math-kernel-library/topic/53107/ dieser winzige Faden. Dieser Thread ist, scheint jedoch, dass einfach zu implizieren, die „U“ obere Dreiecksmatrix nimmt (und jede Zeile durch die diagonalen Teilungs) ist das gleiche wie die reduzierten Echelon Form einer Matrix, die ich glaube nicht der Fall zu sein. Ich kann Rref () selbst codieren, aber ich glaube nicht, dass ich die Leistung LAPACK berühmt erreichen könnte.

1) Gibt es einen besseren Weg, um jede eine spezifische Lösung für ein unterrepräsentierte System einfach zu bekommen?
2) Wenn nicht, gibt es eine Möglichkeit für LAPACK rudern-zu verringern, um eine Matrix?

Danke!

Answer 1

Eine häufig verwendete Methode für diese die kleinsten Quadrate Lösung finden lapack des sgelsx