Discussion :

[jeremux]

Bonjour à tous et à toutes,

je dois créer une fonction creer qui va créer un type matrice. Mais j'ai des problèmes, surement au niveau des boucles.

Mon travail:

type info = {mutable c:float;i:int;j:int};;

(** le type info sert a stocker la valeur d’un coefficient et sa
position dans la matrice*)

type matrice ={ nb_lig:int; nb_col:int; mutable tab: info array};;

(*---------------------------------------------creer-------------------------------------*)
let creer n p =
let remplissage n p =
let mat = Array.create (n*p) ({c=0.0;i=1;j=1}) in
for l=0 to (n*p)-1 do
mat.(l)<-{c=0.0;i=1;j=l+1}
done;
for k=0 to (n*p)-1 do
mat.(k)<-{c=0.0;i=k+1;j=mat.(k).j}
done;
mat in
{nb_lig=n; nb_col=p; tab=remplissage n p};;


Mais j'ai un problème, par exemple

creer 2 3 retourne

[|{c = 0.; i = 1; j = 1}; {c = 0.; i = 2; j = 2}; {c = 0.; i = 3; j = 3};
{c = 0.; i = 4; j = 4}; {c = 0.; i = 5; j = 5}; {c = 0.; i = 6; j = 6}|]}

Et le résultat attendu serait plutôt:

[|{c = 0.; i = 1; j = 1}; {c = 0.; i = 1; j = 2}; {c = 0.; i = 1; j = 3};
{c = 0.; i = 2; j = 1}; {c = 0.; i = 2; j = 2}; {c = 0.; i = 2; j = 3}|]}.


Merci.

PS: Un algorithme me suffira vu que je n'ai aucun problème avec le langage lui même.


Jeremux.

[gasche]

Tu devrais faire deux boucles *imbriquées*, au lieu de deux boucles l'une à la suite de l'autre.

Pour chaque i, tu itères (pour chaque j, définir la case (i, j)).

Par ailleurs, attention à l'utilisation du "mutable" devant le coefficient c. Ça peut provoquer des bugs subtils (let test = Array.make 5 {c = 0.; i = 0; j = 0} in test.(0).c <- 1.; test), et tu peux sans doute t'en passer.
Par ailleurs, le fait de représenter des matrices "creuses" par des tableaux à deux dimensions plein me paraît un peu curieux. Si tu ne connais pas à l'avance la taille de la structure de donnée utiliser pour stocker les éléments, une liste paraît plus adaptée qu'un tableau; et pourquoi avoir deux dimensions ?

Edit : je n'avais pas fait attention, ce n'est pas une matrice de dimensions n et p, mais un tableau (à une dimension) de dimension n*p.

[jeremux]

J'aurais bien aimé faire avec les listes, mais dans le cours on nous impose de faire ça avec des tableaux.

J'ai compris pour la boucle, mais comment stocker dans un tableau ? de taille n*p.

[jeremux]

Ok, j'ai trouvé

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
let creer n p =
let remplissage n p = 
 let mat = Array.create (n*p) ({c=0.0;i=1;j=1}) in 
 let k = ref 0 in
 	for t=1 to n do
 		for d=1 to p do
 			mat.(!k)<-{c=0.0;i=t;j=d};
 			k:=!k+1
 		done;
 	done;
 
 
 mat in
{nb_lig=n; nb_col=p; tab=remplissage n p};;

A bientôt et merci!

[gasche]

Tu n'as pas besoin de déclarer explicitement k, puisqu'il est égal à (t-1)*p+d-1. Si tu indexais de 0 à N-1 comme tout le monde, ce serait même directement t*p+d.

[jeremux]

Pour l'indexation j'ai respecté la définition des coefficient pour une matrice en algèbre linéaire (définis par l'indice de ligne et de colonne). Merci, beaucoup pour l'idée sur k, ça me permet de mieux travailler le produit de deux matrices!