Löse

Ich portiere einen vorhandenen Code von MATLAB nach C ++ und habe ein lineares System zum Lösen von (anstelle der typischeren Form A x = b )$xA=b$$Ax=b$

Die Matrix ist dicht und von allgemeiner Form, aber nicht größer als 1000 × 1000. In MATLAB wird die Lösung also durch die Funktion oder die Schrägstrichnotation gefunden$A$mrdivide(b,A)x = b/A;

Wie soll ich das in meinem C ++ - Code mit BLAS- und LAPACK-Routinen lösen?

Ich kenne die LAPACK-Routine, DGESVdie für x auflöst .$Ax=b$$x$

Ein Gedanke, den ich hatte, ist, einige Manipulationen mit Matrix-Transponierten-Identitäten durchzuführen:

$(xA)^T=b^T$

$A^T x^T = b^T$

$x^T = (A^T)^{-1} b^T$

Lösen Sie dann das endgültige Formular, indem Sie DGESVdas transponierte . (also kosten um A zu transponieren und kosten um System zu lösen)$A^T$$A$

Gibt es einen effizienteren oder auf andere Weise besseren Ansatz ?

Ich arbeite mit Matrix- und Vektorklassen sowie der BLAS-Implementierung aus der BOOST-uBLAS-Bibliothek sowie Bindungen an die LAPACK-Bibliotheksroutinen. Ich habe dieses Setup erfolgreich für andere Vorgänge verwendet und hoffe, eine Lösung zu finden, die auf diese Bibliotheken beschränkt ist.

Außerdem sollte ich beachten, dass ich diese Art von Operation nur einige Male während der Code-Einrichtung durchführe, sodass die Leistung kein kritisches Problem darstellt.

Vielleicht MATLAB Dokumentation auf mrdividefür andere hilfreich.

Triviale Antwort für Quadrat : Benutze, was auch für A T x = b gilt, wenn .$A$dgesvx$A^T x = b$TRANS = 'T'

Bitte beachten Sie, dass Sie mit BLAS oder LAPACK eine Matrix kaum transponieren (Elemente im Speicher vertauschen) müssen: Die meisten Unterprogramme haben ein TRANSArgument, das für den Betrieb mit der transponierten Matrix oder einer Matrix mit einem anderen Speicherlayout geeignet ist. (Die Transposition entspricht der Änderung des Fortran-zusammenhängenden Speicherlayouts in einen C-zusammenhängenden Speicher und umgekehrt.)

$A$

$A$