résoudre

Je porte un code existant de MATLAB vers C ++ et j'ai un système linéaire pour résoudre (plutôt que la forme la plus typique A x = b )$xA=b$$Ax=b$

La matrice est dense et de forme générale, mais ne dépasse pas 1000x1000. Donc, dans MATLAB, la solution est trouvée par la fonction ou la notation barre oblique$A$mrdivide(b,A)x = b/A;

Comment résoudre ce problème dans mon code C ++ à l'aide des routines BLAS et LAPACK?

Je connais la routine LAPACK DGESVqui résout pour x .$Ax=b$$x$

Donc, une pensée que j'avais était de faire quelques manipulations en utilisant des identités de transposition matricielle:

$(xA)^T=b^T$

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

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

Résolvez ensuite la forme finale en utilisant l' DGESVopération sur l' transposé . (donc coût pour transposer A et coût pour résoudre le système)$A^T$$A$

Existe-t-il une approche plus efficace ou autrement meilleure ?

Je travaille avec les classes matricielles et vectorielles ainsi que l'implémentation BLAS de la bibliothèque BOOST uBLAS ainsi que les liaisons aux routines de la bibliothèque LAPACK. J'ai utilisé cette configuration avec succès pour d'autres opérations et j'espère trouver une solution limitée à ces bibliothèques.

De plus, je dois noter que je n'effectue ce type d'opération que quelques fois pendant la configuration du code, donc les performances ne sont pas une préoccupation critique.

Peut-être que cette documentation MATLAB sur mrdivideest utile pour les autres.

Réponse triviale pour le carré : qui résout également pour A T x = b quand .$A$dgesvx$A^T x = b$TRANS = 'T'

Veuillez noter qu'avec BLAS ou LAPACK, vous n'avez pratiquement pas à transposer (permutation d'éléments en mémoire) une matrice: la plupart des sous-programmes ont un TRANSargument pour s'adapter au fonctionnement sur la matrice de transposition ou sur une matrice stockée avec une disposition de mémoire différente. (La transposition équivaut à changer la disposition de la mémoire contiguë de Fortran à une configuration contiguë C et vice versa.)

$A$

$A$