[calculs matriciels] rapidité d'inversion
Bonjour,
Je voudrais inverser des matrices de tailles conséquentes (entre 500*500 jusqu'à + de 10 000*10 000), pour ça j'utilise actuellement la fonction vigra::linalg::inverse de la bibliothèque vigra, mais ça rame, forcément ...Et surtout, ça met plus de temps que d'effectuer mes calculs sous Matlab (vexant) où bien sûr de genre de fonction est optimisée. J'aurai juste aimé savoir si quelqu'un avait eu ce genre de problème et savait comment optimiser les temps de calcul avec 'inverse' ?
Merci d'avance !
[calculs matriciels] rapidité d'inversion
Bonjour et merci pour ta réponse !
Et oui, j'ai finalement opté pour un calcul séquentiel.
Mon problème est du type A.x=B. Je connais A, je cherche à estimer une matrice B en inversant l'opérateur x. Dans mon cas, x correspond à une matrice triangulaire inférieure (suis-je bête), je procède donc par approche séquentielle plutôt que d'inverser cet énorme bloc x :roll: , en temps de calcul c'est forcément incomparable ...
Merci encore,
Sophie