Problèmematrice et do

Bonsoir
j'ai un problème avec do...
voici le code

Code:

---

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

! lecture des données DO I=NB,1,-1 ! evaluation de N CALL CALC(Y,X,I,N,...) DO J=1,N WRITE(*,'(40X,F10.2,F10.4)') X(J),Y(I,J) END DO END DO   ! LA PARTIE OU SE TROUVE LE PROBLEME DO I=1,NB DO J=1,N WRITE(*,*) Y(I,J) END DO END DO   END   SUBROUTINE CALL(Y,X,I,N,...) DO J=N-1,1,-1 .   .   .   . Y(I,J)=Y1 END DO return END

---

Les valeurs qui s'affichent:
-----Y(1,N)-------
99.87
99.91
99.95
99.99
100.03
100.07
100.11
100.15
100.19
100.23
100.27
-----------
Les autres valeurs de la matrice prennent la valeur de Y(2,1)
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30
53.30

Bonjour,

Citation:

---

Code:

---

1 2 3 4 5 6 7

! LA PARTIE OU SE TROUVE LE PROBLEME DO J=1,NB DO J=1,N WRITE(*,*) Y(I,J) END DO END DO

---

---

Utiliser deux fois le même indice dans des boucles imbriquées n'est jamais une bonne idée...

c'est d'ailleurs étonnant que le compilateur est autorisé ça!

non c'est juste une faute de frappe, dans le code j'ai mis ça

Code:

---

1 2 3 4 5 6 7

! LA PARTIE OU SE TROUVE LE PROBLEME DO I=1,NB DO J=1,N WRITE(*,*) Y(I,J) END DO END DO

---

Bonjour,
Différentes petites choses:
1/ je suppose que ta subroutine s'appelle CALC et non call comme tu l'as définie...
2/ return ne prend pas de e
Mais je pense plus que tu as une erreur d'indexation dans ta subroutine (partie cachée)...

bon dans la subroutine CALC y a aucun problème, parce qu'elle me donne toutes les valeurs de la matrice Y(NB,N) dans la première boucle,une fois je sors de cette boucle et j'essaie de réafficher cette matrice je ne trouve que Y(1,N) sachant que NB=1 est la dernière valeur envoyer à CALC pour calculer Y(1,N)

dans ce cas, tu écrase la valeur de NB dans la première boucle et il ne t'affiche que la derniere calculer ensuite parce que tu as écrasé tous les autres... mais ton erreur n'est pas dans ce que tu as écris dans ce post...

déjà comme tu as définit Y dans ta routine CALC? est-tu en F77 ou 90 (savoir si on peut simplifier ta syntaxe)? pkoi ne pas faire la double boucle directement dans calc?

je suis en F77
enfaite mon but c pas d'afficher Y(NB,N) pour la 2eme fois mais je veut stocker cette matrice pour la réutiliser dans une autre subroutine pour faire d'autres calculs, NB est 1 constante évaluer avant de rentrer a la 1ere boucle

j'avais bien compris que tu ne cherchais pas juste à l'affiché sinon la première boucle était suffisante....

mais j'en revien à ma question: si ta routine CALC considère son paramètre Y comme un vecteur, à chaque fois que tu l'appelle, il écrasera les précédentes valeurs et donc il ne te restera que la dernière valeur calculée
d'où: comment sont déclarés les paramètre à l'intérieur de CALC?

voici le code

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21

SUBROUTINE CALC(Y,X,I,Q,Y2,N,XL,B,FM,S0,FN)   COMMON/CONST/G,NB DIMENSION Q(NB),FN(NB),S0(NB),FM(NB),B(NB),XL(NB),Y(NB,N),X(N) Y(I,N)=Y2 Y1=Y(I,N) X(N)=XL(I) DX=-XL(I)/(N-1) DO J=N-1,1,-1 DY1=DX*DYX(Y1,Q,S0,FN,B,FM,I) DY2=DX*DYX(Y1+.5*DY1,Q,S0,FN,B,FM,I) DY3=DX*DYX(Y1+.5*DY2,Q,S0,FN,B,FM,I) DY4=DX*DYX(Y1+DY3,Q,S0,FN,B,FM,I) Y1=Y1+(DY1+DY4)/6.+(DY2+DY3)/3. X(J)=X(J+1)+DX Y(I,J)=Y1 END DO   RETURN END

---

Y est un vecteur

j'ai essai avec un exemple simple et ça donne

Code:

---

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

DIMENSION Y(4,42) DO I=1,4 CALL MATRICE(Y,I) DO k=1,4 WRITE(*,*)Y(I,k) END DO WRITE(*,*)'-----------' END DO call affichage(Y) END     SUBROUTINE MATRICE(Y,I) DIMENSION Y(4,42) DO J=1,4 Y(I,J)=J+I END DO RETURN END   SUBROUTINE affichage(Y) DIMENSION Y(4,42) DO I=1,4 DO J=1,4 WRITE(*,*)Y(I,J) END DO END DO RETURN END

---

mais avec mon programme :calim2: