Discussion :

[mounstein]

Bonjour,

alors voilà j'ai un problème pour traduire mon algorithme de décomposition de Cholesky en pascal !!!

Le problème c'est comment traduire les indices de lignes et de colonnes de la matrice !!
Dois-je utiliser un tableau avec 2 dimensions ou y a t-il une autre astuce ??

Merci de répondre

[M.Dlb]

dois je utiliser un tableau avec 2 dimension ou y a t-il une autre astuce??

Non y'a pas d'autre astuce

[mounstein]

est ce que tu peux m'orienter stp??

[mounstein]

est ce qu'il ya quelqu'un pour m'aider svp??

j'ai mis mon algorithme sous langage pascal!! la compilation est faite mais le resultat est etrange!! j'arrive pas a voire la matrice L que je devais voire en final!

[Eric Sigoillot]

Oui et ? On a oublié notre boule de cristal, donc sans code, sans résultat, sans rien en fait, on ne fera pas grand chose...

@++

[mounstein]

oui c'est vrai desolée voici le code!

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
program cholesky;
 
const
   n=3;
   m=3;
 
var
 
   a:array [1..n,1..m] of real;
   r:array [1..n,1..m] of real;
   i,j,k:integer;
   s1:real;
   s2:real;
 
begin
      s1:=0;
      s2:=0;
      writeln('entrez votre matrice A');
       for i:=1 to n do
        for j:=1 to m do
         readln(a[i,j]);
 
       repeat
         begin
 
           for k:=1 to i-1 do
           s1:=s1+sqrt(r[i,i]);
 
           r[i,i]:=sqrt(a[i,i]-s1);
 
           if r[i,i]=0 then
             writeln('matrice non reguliere')
 
           while (j>1) and (j<=n) do
             for k:=1 to i-1 do
             s2:=s2+(r[i,k]*r[j,k]);
 
             r[i,i]:=(a[i,j]-s2)/r[i,i];
 
           end;
 
      until i=n
end.

[Loceka]

Faire des boucles, c'est bien, mais incrémenter/décrémenter les indices, c'est mieux !

En gros ton repeat ne sert à rien si tu ne fait pas varier "i", de même ton while tourne dans le vide vu que "j" ne varie pas...

[mounstein]

j'ai pas bien compris ou dois je incrementer et decrementer?

[j.p.mignot]

votre extrait de code

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
 repeat
         begin
 
           for k:=1 to i-1 do
           s1:=s1+sqrt(r[i,i]);
 
           r[i,i]:=sqrt(a[i,i]-s1);
 
           if r[i,i]=0 then
             writeln('matrice non reguliere')
 
           while (j>1) and (j<=n) do
             for k:=1 to i-1 do
             s2:=s2+(r[i,k]*r[j,k]);
 
             r[i,i]:=(a[i,j]-s2)/r[i,i];
 
           end;
 
      until i=n

pose les problèmes:
1- ( sans importance mais tout de même)
avec repeat begin end until, la paire begin end est suppreflue
2- la boucle se résume à
repeat .... until i=n
mais i n'est pas incrementé dans la boucle et son initialisation est plus que douteuse car elle résulte de la fin d'un index for i=1 to n

Je n'ai pas regardéle code mais à prioris cela doit être quelquechose du type

Code :

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

1
2
3
4
5
6
 
i:=0;
   repeat
    inc(i);
    ..... { liste des choses à faire }
   until i=n;