Fonction : produit matriciel [Résolu]

[Wenneguen]

Bonjour,

je dois créer une fonction effectuant le produit de 2 matrices mais j'avoue que j'ai un peu de mal. Voilà ce que j'ai fait :

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
let produit_matriciel A B =

 let resultat = make_vect (vect_length A) (make_vect (vect_length B.(0)) 0) 
 in
	let r = ref 0
	in
 	
	for i=0 to vect_length A 
	for j=0 to vect_length A.(0)
	for k=0 to vect_length A.(0)
	do
	  r:= !r + A.(i).(k)*B.(k).(j)
	done;
	resultat.(i).(j) <- !r

	resultat
;;

Évidemment ça ne marche pas, caml me met l'erreur suivant :

Code :

1
2
3
4
5
Entrée interactive:
>	for j=0 to vect_length A.(0)
>	^^^
Erreur de syntaxe.
#

Quel est le problème ?

Merci.

[gorgonite]

la syntaxe du for... il faut des do ?

[Wenneguen]

Oui en effet merci ! J'ai changé quelques trucs, mais le produit calculé n'est pas le bon ( de plus toutes les lignes de la matrice rendue sont les mêmes... )

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
let produit_matriciel a b =

 let resultat = make_vect (vect_length a) (make_vect (vect_length b.(0)) 0) 
 in
 	
	for i=0 to vect_length a -1 do
	 for j=0 to vect_length b.(0) -1 do
	  for k=0 to vect_length a.(0) -1 do
	    resultat.(i).(j) <- resultat.(i).(j) + a.(i).(k)*b.(k).(j)
	  done; 
	 done;
	done;
	
	resultat
;;

[Ptival]

Code :

let resultat = make_vect (vect_length a) (make_vect (vect_length b.(0)) 0)

Le problème, c'est qu'ici tu crées un vecteur de taille "taille de a", dont chaque élément est le même vecteur. Du coup quand tu modifies une ligne, tu modifies toutes les lignes.

Le problème et la solution sont expliqués dans la FAQ du débutant.

[Wenneguen]

En effet c'est exactement ça, merci beaucoup !