Discussion :

[Julien N]

Bonjour,

Admettons le code Fortran suivant (add_f.f90) :

Code :

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

1
2
3
4
5
6
7
8
9
10
function add(x, y)
 
    implicit none
 
    double precision             :: add
    double precision, intent(in) :: x, y
 
    add = x + y
 
end function add

Je peux le compiler en une library comme ceci :

Code :

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

1
2
gfortran -c -fno_underscoring add_f.90
gfortran -shared -fPIC -o add_f.so add_f.o

Je peux vérifier le contenu de la library :

Code :

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

1
2
3
4
5
6
nm -D add_f.so
00000000000010f9 T add
                 w __cxa_finalize
                 w __gmon_start__
                 w _ITM_deregisterTMCloneTable
                 w _ITM_registerTMCloneTable

Je souhaite maintenant appeler la fonction add depuis un script Python. Je me suis inspiré pour cela du sujet suivant : https://www.developpez.net/forums/d2.../#post11404304

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
import ctypes
 
# Définition de la fonction add depuis la library
libf = ctypes.CDLL('add_f.so')
func = libf.add
 
# Définition du type des arguments et de la sortie
func.argtypes = [ctypes.c_double, ctypes.c_double]
func.restype = ctypes.c_double
 
# test
res = func(1.2, 3)
 
print(res)

Qui me renvoie un messagre d'erreur :

Code :

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

/tmp/geany_run_script_SE8FJ0.sh: line 7: 54151 Segmentation fault

Si je ne spécifie pas le type des arguments j'ai un message d'erreur de type TypeError.

Qu'est-ce qui est incorrect ici ?

J

[wiztricks]

Salut,

Quand je me lance dans des opérations inhabituelles, je cherche un peu sur Internet d'abord. Ce que je n'ai pas manqué de faire en lisant votre sujet. Ah ben ouais, on peut le faire sans trop de problèmes pourvu qu'on passe ses arguments par référence et qu'on fasse un peu attention aux variations dans les types.

- W

[Julien N]

Bonjour Wiztricks,

Je ne saisi pas en quoi mon exemple diffère des exemples que j'ai pu voir sur la doc de ctypes. Il me manque une pièce au puzzle.

Je me suis penché sur cette notion de "par référence" et "par valeur". Pour moi les arguments x et y ne doivent pas être passés par référence. Peut être la valeur de retour, add, doit l'être elle.

Mais si je teste :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
from ctypes import CDLL, POINTER, c_double
 
libf = CDLL('add_f.so')
func = libf.add
func.restype = POINTER(c_double)
 
# Test
res = func(c_double(1.2), c_double(3))
 
print(res)

J'obtiens un pointer <__main__.LP_c_double object at 0x7fac006a4ae8>. Et dès que j'essaie d'accéder à sa valeur par l'attribut .contents, je me retrouve avec un segfault à nouveau.

J

[wiztricks]

Salut,

Je ne saisi pas en quoi mon exemple diffère des exemples que j'ai pu voir sur la doc de ctypes. Il me manque une pièce au puzzle.

La documentation de ctypes vous explique juste comment charger une DLL, fabriquer les arguments des fonctions, les appeler,...

Elle ne vous dira pas que telle fonction de la DLL devra être appelée avec des arguments par référence ou par valeur. Et comme ctypes est plutôt orienté C, çà ne racontera pas grand chose non plus sur Fortran.

Après, si vous voulez réinventer la roue plutôt que de voir comment çà peut se faire en cherchant un peu sur Internet...

- W

[Julien N]

Croyez bien que si je pose une question c'est que je n'ai pas trouvé la réponse par moi-même. Il est plus facile de trouver une réponse lorsqu'on sait ce que l'on recherche. C'est pour cette raison que ma question initiale est :

Qu'est-ce qui est incorrect ici ?

Je ne vois pas en quoi je réinvente la roue en essayant de faire marcher un exemple simple destiné à comprendre comment cela fonctionne.

[wiztricks]

Je ne vois pas en quoi je réinvente la roue en essayant de faire marcher un exemple simple destiné à comprendre comment cela fonctionne.

Tous les exemples que j'ai trouvé disent (et montrent) qu'avec fortran, les arguments sont appelés par référence...

- W