problème de Subroutine en fortran 90

[invite6236f8da]

Bonjour,

J'essaie de faire tourner un code écrit en fortran90. Le code en question contient une subroutine, dès que je la met comme internal subroutine avec ( contains avant le end program) le code tourne correctement, mais une fois j'opte pour une subroutine externe (après le end program) le code perd conscience et fournit des résultats erronés. sachant que je ne change rien dans la subroutine, je ne fais que copier/coller après : end program.


quelqu'un pourrait élucider ce problème ?
D'avance merci.
-----

[JPL]

Ce message s'était égaré à un endroit où personne ne pouvait répondre.

[Calvert]

Salut,

peut-être que mettre le code en question (ou un petit exemple montrant le même comportement) serait utile. En l'état, c'est difficile d'aider.

[invite6236f8da]

  contains   		 
		  subroutine thomas(k,m,a,b,c,F,V)
		  integer m,k
		  real, dimension(200):: a,b,c,alpha,beta,y,F,V

		 alpha(k)=a(k)
		 do i=k+1,m
		 beta(i)=b(i)/alpha(i-1)
		 alpha(i)=a(i)-beta(i)*c(i-1)
		 enddo

		 
		 y(k)= F(k)
		 do i=k+1,m
		 y(i)= F(i)-(beta(i)*y(i-1))
		 enddo
		
		V(m)=y(m)/alpha(m)
		do i=m-1,k,-1
		V(i)= ((y(i)-c(i)*V(i+1))/alpha(i))
		enddo

		  return 
		 end subroutine thomas

[A voir en vidéo sur Futura]

[Calvert]

Re,

alors, le problème vient de ta variable 'i' dans la sous-routine.

Si tu mets ta routine dans le corps du programme avec un "contains", fortran comprends que le "i" de la sous-routine est le même que le "i" du corps principal du programme (une variable déclarée dans le corps principal est visible jusqu'au "end program"). Chez moi, avec gfortran, ça compile, mais ça se comporte mal: comme tu appelles ta sous-routine dans une boucle sur "i", et que tu changes "i" dans la sous-routine, c'est la panique à bord). Si je remplace ton "i" par un "i2", ça marche bien.

Si tu mets ta routine hors du programme principal, en l'état, ça marche parfaitement bien. Mais seulement parce que tu n'as pas mis de "implicit none" dans ta sous-routine. Sinon, la variable "i" n'est pas déclarée.

Bref :
- Je recommanderais de mettre des "implicit none" systématiquement pour chaque sous-routine.
- Quand tu utilises des sous-routines internes, il faut être très prudent avec les noms des variables : éviter d'employer les mêmes variables que dans le programme principal si elles ne sont pas sensées avoir le même usage.

[invite6236f8da]

Effectivement, c'est exactement le problème !
Merci infiniment Calvert