Bonjour,

J'utilise la bibliotheque Eigen pour faire un calcul matriciel.
Je mets a jour Un certain nombre de matrices de maniere iterative. Probelem, au bout de la second boucle, la matrice ne s'inverse plus et j'obtiens des NAN partout. J'ai implementer le meme calcul sous matlab, pourtant la matrice s'inverse.

Un aperçu de mon code (en gros ce qui se passe):
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
 
for(int i=0; i<50; i++) //iterations
{
//mise a jour de toutes les autres matrices
.....
//equations normales
Eigen::MatrixXd S(24,12);
S=(Qll*G.transpose())*(G*Qll*G.transpose()).inverse(); //endroit de l'inversion
}
Toutes les autres matrices se mettent correctement a jour, sauf l'inversion de S...

Merci d'avance pour votre aide