Discussion :

[salseropom]

Bonjour,
je suis en train d'écrire un code d'optimisation et j'utilise pour cela la libraire Harwell. La fonction que j'appelle est vb01ad. Cette fonction appelle la fonction deriv.f (calculant une estimation des dérivée de la fonction coût) que l'utilisateur implémente. Mais deriv n'est pas un argument de vb01ad : l'utilisateur doit appeler une fonction deriv.f (j'espère être suffisamment clair).

Il se trouve que dans mon code, j'ai une seconde fonction coût, donc j'appelle 2 fois la fonction vb01ad. Je dois donc faire une seconde fonction deriv.f ...

j'ai donc 2 problèmes :
1) 2 fonctions différentes ne peuvent pas avoir le même nom
2) comment dire prend la fonction deriv1.f ou deriv2.f (qui doivent d'ailleurs s'appeler deriv.f)

En C, je pourrais faire qqch du genre ("pointeur de fonction")

deriv=deriv1
ici j'appelle une 1e fois vb01ad

deriv=deriv2
ici j'appelle une 2e fois vb01ad

Pouvez-vous m'aider s'il vous plait ?

Le code que je reprends est en fortran 77 mais je peux le réécrire en fortran 90 s'il le faut

Merci d'avance

[Sylvain Bergeron]

Moins élégant et moins efficace qu'en C, mais voici ma suggestion :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 
      function deriv(x)
      common/SelecteurPourDeriv/ iSelect
      if (iSelect == 1)then
         code pour le cas #1
      else
         code pour le cas #2
      endif
      end
 
      program 
      ...
      common/SelecteurPourDeriv/ iSelect
      ...
      iSelect = 1
      appel à vb01ad pour le cas #1
      ...
      iSelect = 2
      appel à vb01ad pour le cas #2
      ...

F77 ou F95 ?

C'est certain que j'utiliserais le format libre et un module plutôt qu'un common pour passer le paramètre si j'avais le choix.

[salseropom]

OK merci !
je n'avais pas pensé au common !

[Ladgalen]

Bonjour

Je fais ressortir ce post pour savoir avec un module s'il y a une autre façon de faire ?

Par exemple est ce qu'on peut créer un type dérivé dans lequel on a le nom de la fonction a utiliser ?