Discussion :

[oliverjack]

Bonjour à tous,

Dans un programme je dois résoudre une équation differentielle de type
y=f(y,t).

Pour ce faire je me suis tourné vers les différents solver disponible en fortran.
J'ai essayé rkf45 qui ne semble pas tt a fait approprié pour mon problème. (trop d'évaluations)
Une alternative qui semble intéressante et efficace est lsoda.

La routine demande les arguments suivant:

Code :

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

1
2
3
4
5
6
SUBROUTINE DLSODA (F, NEQ, Y, T, TOUT, ITOL, RTOL, ATOL, ITASK,
     1            ISTATE, IOPT, RWORK, LRW, IWORK, LIW, JAC, JT)
      EXTERNAL F, JAC
      INTEGER NEQ, ITOL, ITASK, ISTATE, IOPT, LRW, IWORK, LIW, JT
      DOUBLE PRECISION Y, T, TOUT, RTOL, ATOL, RWORK
      DIMENSION NEQ(*), Y(*), RTOL(*), ATOL(*), RWORK(LRW), IWORK(LIW)

Mon problème vient de JAC qui est dit optionnel si JT vaut 2
JAC est une matrice jacobienne et je ne veux pas la fournir au programme, si JT vaut 2, la matrice jacobienne est calculée par une subroutine de DLSODA.

Dans la documentation il est bien écrit que JAC est optionnel et que si on ne souhaite pas fournir de matrice jacobienne il suffit d'envoyer " a dummy argument".

C'est la que je suis perdu, pour moi un dummy est un argument qu'on passe a une fonction, ou subroutine. Alors la je veux bien donner un argument mais quoi?

Voici un lien vers la subroutine qui contient la documentation.

http://www.shocksolution.com/math_to...cs/dlsoda.html
Si quelqu'un pouvait m'aider je sèche...

OJ

[rogue-spectre]

Salut,

Je ne connais pas très bien l'utilisation du mot clé "external". Quand je veux un argument optionnel j'utilise le mot clé "optional" par exemple JAC serait optionnel avec un déclaration

Code :

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

1
2
 
     REAL, DIMENSION(:,:), OPTIONAL :: JAC

de plus les arguments optionnels sont de préférence placés à la fin de la liste des arguments car sinon si tu as par exemple

Code :

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

1
2
3
4
 
subroutine lala(titi, toto, tata, tutu)
   int :: titi, tata, tutu
   int, optional :: toto

et que tu veux appeler ta routine sans l'argument toto tu vas devoir faire un truc du genre

Code :

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

1
2
 
    call lala( a, tata=b, tutu=truc)

alors que si tu déclares ta routine

Code :

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

1
2
3
4
 
subroutine lala(titi,  tata, tutu, toto)
   int :: titi, tata, tutu
   int, optional :: toto

l'appel est simplement

Code :

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

1
2
 
 call lala( a, b, truc)

Donc a mon avis dans ton code il faut quand même que tu passes une matrice JAC pour JT=2, même si ta matrice est bidon.

En espérant t'avoir aidé
Bonne journée

[Ehouarn]

Bonjour,

La routine demande les arguments suivant:

Code :

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

1
2
3
4
5
6
SUBROUTINE DLSODA (F, NEQ, Y, T, TOUT, ITOL, RTOL, ATOL, ITASK,
     1            ISTATE, IOPT, RWORK, LRW, IWORK, LIW, JAC, JT)
      EXTERNAL F, JAC
      INTEGER NEQ, ITOL, ITASK, ISTATE, IOPT, LRW, IWORK, LIW, JT
      DOUBLE PRECISION Y, T, TOUT, RTOL, ATOL, RWORK
      DIMENSION NEQ(*), Y(*), RTOL(*), ATOL(*), RWORK(LRW), IWORK(LIW)

Mon problème vient de JAC qui est dit optionnel si JT vaut 2
JAC est une matrice jacobienne et je ne veux pas la fournir au programme, si JT vaut 2, la matrice jacobienne est calculée par une subroutine de DLSODA.

Dans la documentation il est bien écrit que JAC est optionnel et que si on ne souhaite pas fournir de matrice jacobienne il suffit d'envoyer " a dummy argument".

C'est la que je suis perdu, pour moi un dummy est un argument qu'on passe a une fonction, ou subroutine. Alors la je veux bien donner un argument mais quoi?

Voici un lien vers la subroutine qui contient la documentation.

http://www.shocksolution.com/math_to...cs/dlsoda.html

Vu les syntaxes employées, ce n'est pas du Fortran "moderne" (I.e. 90 ou +)qui utiliserait le mot clé "optional" pour préciser qu'un argument est optionel.
Je pencherai plus pour interpréter "a dummy argument can be used" comme: cet argument ne sera pas utilisé, donc mettre là un tableau "bidon", simplement pour être conforme aux arguments d'appel de la routine et sachant que son contenu pourra être n'importe quoi puisqu'il ne sera pas utilisé.

Bonne continuation.

[Sylvain Bergeron]

Vu les syntaxes employées, ce n'est pas du Fortran "moderne" (I.e. 90 ou +)qui utiliserait le mot clé "optional" pour préciser qu'un argument est optionel.
Je pencherai plus pour interpréter "a dummy argument can be used" comme: cet argument ne sera pas utilisé, donc mettre là un tableau "bidon", simplement pour être conforme aux arguments d'appel de la routine et sachant que son contenu pourra être n'importe quoi puisqu'il ne sera pas utilisé.

J'en arrive à la même interprétation. Je pense que tu n'as même pas à passer un tableau bidon, un simple scalaire devrait fonctionner si l'argument n'est pas utilisé par la routine.