Forum |  HardWare.fr | News | Articles | PC | S'identifier | S'inscrire | Shop Recherche
636 connectés 

  FORUM HardWare.fr
  Programmation
  C

  Multiplication matricielle

 
 


 Mot :   Pseudo :  
 
Bas de page
Auteur Sujet :

Multiplication matricielle

n°1452925
ngkreator
Posté le 05-10-2006 à 17:13:33  profilanswer
 

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 ;)

mood
Publicité
Posté le 05-10-2006 à 17:13:33  profilanswer
 

n°1452932
_darkalt3_
Proctopathe
Posté le 05-10-2006 à 17:16:15  profilanswer
 

=+ ?
c'est pas += ?

Message cité 1 fois

---------------
Töp of the plöp
n°1452939
ngkreator
Posté le 05-10-2006 à 17:19:21  profilanswer
 

_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
n°1453020
franceso
Posté le 05-10-2006 à 18:47:06  profilanswer
 

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
n°1453028
ngkreator
Posté le 05-10-2006 à 19:04:28  profilanswer
 

Plus rapide?

n°1453037
jesus_chri​st
votre nouveau dieu
Posté le 05-10-2006 à 19:31:02  profilanswer
 

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
n°1453208
bb138
La vie est belle ...
Posté le 06-10-2006 à 09:44:20  profilanswer
 

jesus_christ +1 :)

n°1453244
franceso
Posté le 06-10-2006 à 10:46:04  profilanswer
 

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
n°1453618
jesus_chri​st
votre nouveau dieu
Posté le 06-10-2006 à 19:32:08  profilanswer
 

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

n°1453622
0x90
→
Posté le 06-10-2006 à 19:37:04  profilanswer
 

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 :/ )


Aller à :
Ajouter une réponse
  FORUM HardWare.fr
  Programmation
  C

  Multiplication matricielle

 

Sujets relatifs
multiplication de 2 matricesboucle while: table de multiplication
Multiplication de deux matricesmultiplication et overflow
[Assembleur] Problème avec la multiplication[SQL] Multiplication de colonne
Algo FFT et multiplication rapideinversion matricielle
multiplication de matrices[C] Multiplication de polynômes. Ce code est-il OK ? [Résolu]
Plus de sujets relatifs à : Multiplication matricielle

Forum MesDiscussions.Net, Version 2010.2
(c) 2000-2011 Doctissimo
Page générée en 0.075 secondes

Copyright © 1997-2022 Hardware.fr SARL (Signaler un contenu illicite / Données personnelles) / Groupe LDLC / Shop HFR