Sortir d'une subroutine

Version imprimable

07/08/2007, 14h31
Butterfly83

Sortir d'une subroutine

Bonjour,
alors j aimerais avoir un peu d aide, sachant que je suis novice en Fortran. Alors je raisous un système et suivant que ma méthode converge ou pas je continu ou je veux diviser mon pas de temps pour adapter, en gros, mon maillage.
J ai pensé à faire un "goto" mais a priori ca ne marche pas, sachant que mon pas de temps et toutes les interpolation qui en découle est définit dans un autre fichier.

Si quelqu un a une petite idée ce serrait sympa.
Merci.

fait montrer comment sont imbriqué (en gros) tes boucles, sinon tu a la boucle DO avec condition final et la commande EXIT:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
 
Boucletemporelle:do while(t<TempsFinal)
  Boucle_pas_de_temps:Do
    <.................>
    <tout ton calcul + calcul de l'erreur (servant à adapter le pas)
    <.................>
    if (erreur<seuil)then 
      t=t+dt
      c=(Constante/erreur)**(1.d0/5.d0)   !facteur d"adaptation du pas          
      if (c>2.d0) c=2.d0
      if (c<5.d-1) c=5.d-1 
      dt=c*dt
      EXIT
    else
      dt=dt/2.d0
    endif
  end do Boucle_pas_de_temps
end do boucle temporelle

07/08/2007, 16h47
Butterfly83

En gros ce que j ai c est:

Boucle de temps:
subroutine(1) qui fait tout les calculs
subroutine(2) qui dit si ca converge
Fin boucle de temps

et je voudrais que quand (2) me dit que ca converge pas retourner au début de la boucle temps.
07/08/2007, 16h48
Butterfly83

?le facteur d adaptation c est quoi?
07/08/2007, 16h49
Butterfly83

la subroutine (2) est dans la (1), c est très simplifié mais en gros c est ce qui ce passe.
07/08/2007, 16h52
Butterfly83

Mais le schéma que tu m as donné ressemble à ce que j ai dans la tête et ce que je veux programmer.
08/08/2007, 11h25
genteur slayer

typiquement, cela ne change rien, à condition de bien laissé le 'EXIT' dans la boucle (et pas dans une subroutine) mais dans ton cas, plutot que de dire si convergence ou pas, ta subroutine2 devrais te renvoyé une estimation de la majoration de l'erreur et tu fais des conditions de sortie sur cette erreur.... en gros le facteur c permet d'agrandir le pas te temps lorsque l'erreur est petite et de le diminuer si l'erreur et un peu trop grande mais sans recalculer le pas que l'on viend de faire, en fait, on adapte le par de temps de l'itération d'après.... (cela permet d'améliaurer les convergences) cependant il faut que 1/2<c<2 sinon l'écart entre les différents pas de temps est un peu trop grand...
08/08/2007, 11h38
Butterfly83
Je dit peut etre une connerie et je m en excuse si c est le cas mais de la où je suis j ai aucun moyen de modifié le pas de temps ce que mon programme fait c est qu a chaque pas de temps il calcul la valeur de mon inconnue a tout les pas d'espace et après il fait t=t+dt.
D après ce que tu me dit je suis tentée de faire:
Code:

1 2 3 4 if (not.converge) exit dt=dt/2 endif
un truc dans le genre....

Merci pour tes réponse et encore une fois désolée si je suis pas très réactive, je trouve ce langage un peu rustique..
08/08/2007, 11h45
genteur slayer

déjà, oublie pas les balises [CODE], ensuite, oublie pas le THEN, ensuite(bis), le EXIT met le en fin d'instruction sinon ce qu'il y a parès n'est pas calculer, ENFIN, si tu peux faire "dt=dt/2.d0" alors tu peux faire "dt=c*dt"...

au passage oublie pas le ".D0" signifiant que le "2" est un double precision et non un entier sinon tu te retrouverai avec la division euclidienne....
08/08/2007, 12h48
Butterfly83

Tout d abord je te remercie pour ton aide c est gentil.

Ensuite j ai bien fait ce que tu m as dit mais à la compilation y a une erreur:
"unclassificable statement (starts "exit")"

je comprend pas.
08/08/2007, 13h11
Butterfly83

Je crois savoir pourquoi ça me donne cette erreur et c est très con ce que j ai fait en fait je suis pas dans une bouble...
08/08/2007, 13h38
genteur slayer

quand tu écris sur ce forum, y a un icone en forme de # en haut et cela t'ouvre les balise code:

Code:

là je suis dans les balises code

cela permet de mieux voir où sont tes morceaux de codes... c tout!!!

ensuite pour EXIT, certain compilateur comme g95 comprennent une fonction exit(errcode) qui quitte complètement le programme et revoie errcode au système d'exploitation. lorsque que tu es dans une boucle DO, (qu'elle soit infinie ou non) la commande EXIT te sort de cette boucle (et cette boucle seulement)

sinon je te renvoie à un super cours :http://www.idris.fr/data/cours/lang/...tran_cours.pdf tu va au chapitre 5.2 et y a tout d'expliqué sur les boucles...

le chapitre commence par les GOTO mais si l'on veux programmer de façon propre, à mon avis, cette instruction n'est à utiliser que si on ne peux VRAIMENT pas faire autrement.... il y a toujours au moins deux solution à un problème....
08/08/2007, 13h44
Butterfly83

Merci...
Je vais jeter un oeil au cours.

D'après ce que je comprend en fait "exit" va pas me servir a grand chose moi je veux que quand il se rend compte qu il ne converge pas il revienne au debut du programme et qu il divise le pas de temps d avant en fait...
08/08/2007, 14h07
genteur slayer
il faut que l'avancement en temps se fasse à la fin de ta boucle temporelle:
Code:

1 2 3 4 5 6 7 8 if (converge)then t=t+dt n=n+1 !ça c'est pour tes stockages en tableau (par exemple) else dt=dt/2 cycle endif
si il a convergé, alors il avance dans le temp et dans les indices de tableau, sinon divise le pas de temps par deux et revien au début de la boucle
08/08/2007, 15h13
Butterfly83

Excuse moi encore. Mais je comprend^pas ca me met encore cette erreur
"unclassifiable statement (starts "cycle")"
08/08/2007, 17h30
genteur slayer

t'es pas dans une boucle....

ou alors peut-être est-ce là des instruction Fortran95 mais il me semble pas....

quel est ton compilo/environement?
09/08/2007, 11h28
Butterfly83

Finalement je m en suis sortie avec un "goto" c est pas très élégant mais a priori ca à l air de marcher.

Je te remercie pour ton aide et à bientot pour de prochaines aventures..:D