Question Dummy Argument

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:

---

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

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:

---

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:

---

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:

---

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

---

alors que si tu déclares ta routine

Code:

---

1 2 3 4

subroutine lala(titi, tata, tutu, toto) int :: titi, tata, tutu int, optional :: toto

---

l'appel est simplement

Code:

---

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

Bonjour,

Citation:

---

Envoyé par oliverjack

La routine demande les arguments suivant:

Code:

---

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.

Citation:

---

Envoyé par Ehouarn

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.