INVERSION DE MATRICE

Commentaire de maupertuit le 08/10/2008 13:57:46

Ton code à l'air fonctionnel, par contre juste une modif ou deux:
* il manque les delete
* si tu alloues une matrice, alloue les données dans un vecteur continu, ca sera plus rapide à l'excution.

void InvdiagMat(int NbElement, float *MatVect)
{
   int sizeBuffer=NbElement*NbElement;

   float **Mat=new float*[NbElement];
   for(int i=0;i<NbElement;++i)
      Mat[i]=&MatVect[i*NbElement];

   bool testDiag=true;
   for(int i=0;i<NbElement;++i)
      if(Mat[i][i]==0)
      {
         testDiag=false;
         break;
      }

      if(testDiag)
      {
         float *  tmpVect=new float[sizeBuffer];
         float ** temp=new float*[NbElement];
         for(int i=0;i<NbElement;i++)
            temp[i]=&tmpVect[i*NbElement];

         for(int i=0;i<sizeBuffer;i++)
            tmpVect[i]=0.0f;



         for(int i=0;i<NbElement;i++)
         {

            for(int j=0;j<NbElement;j++)
            {
               temp[i][i]=1/Mat[i][i];
               if(j!=i) temp[i][j]=-Mat[i][j]/Mat[i][i];

               for(int k=0;k<NbElement;k++)
               {

                  if(k!=i) temp[k][i]=Mat[k][i]/Mat[i][i];

                  if(j!=i &&k!=i)
                     temp[k][j]=Mat[k][j]-Mat[i][j]*Mat[k][i]/Mat[i][i];
               }

            }

            for(int i=0;i<NbElement;i++)
               for(int j=0;j<NbElement;j++)
                  Mat[i][j]=temp[i][j];


         }

         delete temp;
         delete tmpVect;
      }
      delete Mat;
}