Discussion :

[feynman]

Bonjour,
j'ai écrit ce programme qui calcule la solution d'un système triangulaire supérieur.
Le problème est que, si j'utilise double précision, je n'obtiens pas des valeurs de A ni de X (A étant la matrice et X le vecteur solution); alors où est le truc ?
Merci pour votre aide
Voilà mon programme :

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
   program trianglee
              parameter(n=3)
              double precision A(n,n),b(n),X(n)
              integer i,j
              do i=1,n
                 do j=1,n
                    if((i.lt.j).or.(i.eq.j)) then
                       A(i,j)=float(i*j)
                     else
                       A(i,j)=0.
                     end if
                   enddo
               enddo
 
              print 10,((A(i,j),j=1,n),i=1,n)
10            format('A=',3(/,1X,3(D4.1,1X)))
              do i=1,n
                 b(i)=1.
              enddo
              Call triangle(n,A,b,X)
              print 20,(X(i),i=1,n)
20            format('X=',3(D4.1,1X))
              end
 
                subroutine triangle(M,C,D,Y)
 
                integer j,i
                double precision C(M,M),D(M),Y(M),S
                Y(M)=D(M)/C(M,M)
                   do i=M-1,1,-1
                      S = 0.
                      do j=i+1,M
                          S=S+C(i,j)*Y(j)
                      enddo
                    Y(i)=(D(i)-S)/C(i,i)
                    enddo
                return
                end

[FR119492]

Salut !

Tu es vraiment incorrigible.

1. Tu ignores toujours l'instruction Implicit None.
2. Tu nous dis que tu ne trouves pas la solution, mais tu ne nous dis pas ce que tu trouves.
3. Tu écris ((i.lt.j).or.(i.eq.j)) alors que (i.le.j) serait beaucoup plus simple.
4. Tu utilises un format D4.1 alors que tu n'as manifestement pas compris la signification du 4 et du 1.
5. Tu ne semble pas avoir encore compris qu'avant d'appeler à l'aide, tu devrais chercher un peu toi-même.
6. La méthode la plus simple pour déboguer un programme consiste à imprimer suffisamment de résultats intermédiaires.

Jean-Marc Blanc

[feynman]

Salut !

Tu es vraiment incorrigible.

1. Tu ignores toujours l'instruction Implicit None.
2. Tu nous dis que tu ne trouves pas la solution, mais tu ne nous dis pas ce que tu trouves.
3. Tu écris ((i.lt.j).or.(i.eq.j)) alors que (i.le.j) serait beaucoup plus simple.
4. Tu utilises un format D4.1 alors que tu n'as manifestement pas compris la signification du 4 et du 1.
5. Tu ne semble pas avoir encore compris qu'avant d'appeler à l'aide, tu devrais chercher un peu toi-même.
6. La méthode la plus simple pour déboguer un programme consiste à imprimer suffisamment de résultats intermédiaires.

Jean-Marc Blanc

Cher Jean Marc,
1.Oui je n utilise pas implicit none je ne vois pas son importance
2.Voila ce que je trouve :

Code :

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

1
2
3
4
5
6
7
   elarwadi@lmarech3:~$ g77 -o triangle triangle.f
elarwadi@lmarech3:~$ ./triangle
A=
 **** **** ****
 **** **** ****
 **** **** ****
X=**** **** ****

3. Si je comprends ce que le nombre 4 et le nombre 1 represente:
si tu as Dw.d alors le w c est le nombre de place desirés et le d est le nombre des nombres apres la virgules
voila

[feynman]

Quand je met real a la place de double precision ca marche, mais quand je le remplace par double precision, je tombe dans le truc

[gregory2526]

Il faut bien être clair sur l'intérêt du IMPLICIT NONE

Il oblige à déclarer le type de chaque variable. Sans cela, toutes les variables sont des réel, sauf celles dont la première lettre est entre I et P (il me semble)

L'un des intérêts :
Si par exemple tu te trompe dans un nom de variable (ex : tu as une variable "epsilon" et a un moment donné tu écris "epslion"), le compilateur va de mettre un message d'erreur et tu n'auras pas de problème. Si tu n'a pas utilisé implicit none, tu n'auras pas d'erreur de compilation mais une erreur à l'utilisation du programme, et tu auras plus de mail à découvrir le problème.

[TheOyoStyledMan]

Bonjour,

je te prie de bien vouloir comprendre l'intérêt du IMPLICIT NONE. Avec moi cela fera trois personnes qui te le diront.

Merci gregory2526 pour ton explication.

PS: on ne te dit pas cela pour t'embêter, mais sincèrement, mon boulot consiste en partie à corriger ces genres d'erreurs et c'est la première chose que l'on regarde (et qui nous permet de "juger" le code auquel on a à faire). Mettre un IMPLICIT NONE permet premièrement à la personne qui programme d'être plus claire et sûre sur ce qu'elle veut faire et deuxièmement à une personne qui relit de mieux comprendre le code.