Multiplication matricielle

Multiplication matricielle - C - Programmation

Marsh Posté le 05-10-2006 à 17:13:33    

J'essaye de faire un programme pour multiplier des matrices carrées.  
 
Pour ce faire j'utilise la ligne suivante dans une boucle:
 

Code :
  1. MatriceResultat[ligne][colonne]=+(matriceA[ligne][curseur]*matriceB[curseur][colonne]);


Cela équivaut à la formule: http://gersoo.free.fr/Images/docs/stras/stra02.gif
 
J'utilise cette ligne dans une boucle (le curseur allant de 0 à la taille de la matrice-1). Je répète tout ça pour chaque "case" de la matrice du résultat biensûr.
 
Ca donne:
 

Code :
  1. void multiplier(int matA[M][M],int matB[M][M],int resultat[M][M]){
  2.      int ligne,colonne,curs;
  3.      zero(resultat); //mise à zéro de tous les coefficients de la matrice resultat
  4.    
  5.      //On parcours chaque ligne:
  6.      for(ligne=0;ligne<M;ligne++)
  7.         //On parcours chaque coefficient de la ligne:
  8.         for(colonne=0;colonne<M;colonne++)
  9.            //Calcul d'un coefficient (ligne,colonne) de la matrice resultat:
  10.            for(curseur=0;curseur<M;curseur++){
  11.            resultat[ligne][colonne]=+(matA[ligne][curseur]*matB[curseur][colonne]);
  12.            }
  13. }


 
Le problème est que ça ne fonctionne pas comme vous vous en doutez. Mais je voudrai juste savoir si vous auriez fait la même chose que moi pour la formule du dessus?
 
Merci pour votre aide ;)

Reply

Marsh Posté le 05-10-2006 à 17:13:33   

Reply

Marsh Posté le 05-10-2006 à 17:16:15    

=+ ?
c'est pas += ?

Message cité 1 fois

---------------
Töp of the plöp
Reply

Marsh Posté le 05-10-2006 à 17:19:21    

_darkalt3_ a écrit :

=+ ?
c'est pas += ?


 
 [:cupra]  Je suis vraiment tête en l'air!
 
Je viens de me remotiver à refaire du C après 6 mois d'inactivité [:totoz]  C'est impressionant tous les automatismes que j'ai perdu!
 
Merci beaucoup sinon je vais voir ça.
 
Edit: Oui c'est bon c'est ça. Tu  m'a évité de perdre beaucoup de temps.
 

Code :
  1. /*Programme d'opérations sur des matrices carrées d'ordre M*/
  3. Voilà mon programme on sait jamais ça peut servir pour ceux qui tomberont sur le sujet par google ou autre:
  5. #include <stdio.h>
  6. #define M 4
  8. //Définition des fonctions
  10. void afficher(int tab[M][M]){
  11.      int ligne,colonne;
  12.      printf("\n" );
  13.      for(ligne=0;ligne<M;ligne++){
  14.         for(colonne=0;colonne<M;colonne++)
  15.            printf("\t%d",tab[ligne][colonne]);
  16.         printf("\n" );
  17.         }
  18.      }
  19.    
  20. void unite(int tab[M][M]){
  21.      int ligne,colonne;
  22.    
  23.      for(ligne=0;ligne<M;ligne++){
  24.         for(colonne=0;colonne<M;colonne++)
  25.            if(colonne==ligne)
  26.               tab[ligne][colonne]=1;
  27.            else
  28.               tab[ligne][colonne]=0;
  29.      }
  30. }
  32. void zero(int tab[M][M]){
  33.      int ligne,colonne;
  34.    
  35.      for(ligne=0;ligne<M;ligne++)
  36.         for(colonne=0;colonne<M;colonne++)
  37.         tab[ligne][colonne]=0;
  38. }
  40. void initialiser(int tab[M][M]){
  41.      int ligne,colonne;
  42.    
  43.      for(ligne=0;ligne<M;ligne++)
  44.         for(colonne=0;colonne<M;colonne++)
  45.         tab[ligne][colonne]=ligne+colonne;
  46. }
  48. void multiplier(int matA[M][M],int matB[M][M],int res[M][M]){
  49.      int ligne,colonne,curs;
  50.      zero(res);
  51.    
  52.      for(ligne=0;ligne<M;ligne++)
  53.         for(colonne=0;colonne<M;colonne++)
  54.            for(curs=0;curs<M;curs++){
  55.            res[ligne][colonne]+=(matA[ligne][curs]*matB[curs][colonne]);
  56.            }
  57. }
  59. void additionner(int matA[M][M],int matB[M][M],int res[M][M]){
  60.      int ligne,colonne,curs;
  61.      zero(res);
  62.    
  63.      for(ligne=0;ligne<M;ligne++)
  64.         for(colonne=0;colonne<M;colonne++)
  65.         res[ligne][colonne]=matA[ligne][colonne]+matB[ligne][colonne];
  66.          
  67. }
  69. //Début du programmme!
  71. void main()
  72. {
  73.      int matriceA[M][M],matriceB[M][M],resultat[M][M];
  74.    
  75.      //Les matrices A et B sont converties en matrices unité:
  76.      unite(matriceA); unite(matriceB);
  77.      afficher(matriceA); afficher(matriceB);
  78.    
  79.      //La matrice qui contiendra le résultat est remplie de zéros:
  80.      zero(resultat); afficher(resultat);
  81.    
  82.      //On multiplie les 2 matrices A et B: on insère le résultat dans la matrice resultat:
  83.      multiplier(matriceA,matriceB,resultat); afficher(resultat);
  84.    
  85.      system("PAUSE" );
  86. }


Message édité par ngkreator le 05-10-2006 à 17:21:55
Reply

Marsh Posté le 05-10-2006 à 18:47:06    

juste une remarque en passant : en pratique il est souvent plus rapide d'effectuer la sommation sur une variable temporaire qu'on recopie après dans la case de la matrice (ne me demande pas pourquoi, je sais pas vraiment)
 

Code :
  1. void multiplier(int matA[M][M],int matB[M][M],int resultat[M][M]){
  2.      int ligne,colonne,curs;
  3.      zero(resultat); //mise à zéro de tous les coefficients de la matrice resultat
  4.    
  5.      //On parcours chaque ligne:
  6.      for(ligne=0;ligne<M;ligne++)
  7.         //On parcours chaque coefficient de la ligne:
  8.         for(colonne=0;colonne<M;colonne++){
  9.            //Calcul d'un coefficient (ligne,colonne) de la matrice resultat:
  10.            int tmp = 0;
  11.            for(curseur=0;curseur<M;curseur++)
  12.               tmp += (matA[ligne][curseur]*matB[curseur][colonne]);
  13.            }
  14.            resultat[ligne][colonne] = tmp;
  15. }


 
je sais pas si ça t'intéresse, mais au cas où...


Message édité par franceso le 05-10-2006 à 18:48:07

---------------
TriScale innov
Reply

Marsh Posté le 05-10-2006 à 19:04:28    

Plus rapide?

Reply

Marsh Posté le 05-10-2006 à 19:31:02    

oui, car la somme se fait dans une variable temporaire qui a toutes les chances d'être dans un registre.
au passage, une déclaration
 
register int tmp = 0;
 
serait + explicite

Message cité 1 fois
Reply

Marsh Posté le 06-10-2006 à 09:44:20    

jesus_christ +1 :)

Reply

Marsh Posté le 06-10-2006 à 10:46:04    

jesus_christ a écrit :

oui, car la somme se fait dans une variable temporaire qui a toutes les chances d'être dans un registre.

Il y a aussi il me semble le fait que si resultat est un pointeur, ton compilo fait extrêment gaffe aux problèmes d'aliasing et du coup certaines optimisations ne peuvent pas être faites.


---------------
TriScale innov
Reply

Marsh Posté le 06-10-2006 à 19:32:08    

exact, il vaut mieux déréférencer, calculer puis remettre en mémoire ;)

Reply

Marsh Posté le 06-10-2006 à 19:37:04    

Ca risquerais pas un ptit peut de faire mal si les entiers sont un poil trop gros ?
 
( ce qui va vite quand on fait la somme de plusieurs produits :/ )

Reply

Sujets relatifs:

Leave a Replay

Make sure you enter the(*)required information where indicate.HTML code is not allowed