pssage d'un tableau de string entre Fortran DLL et python avec ctypes

**lakhmesi** · 08/09/2023, 15h58

Bonjour à tous,

je débute avec Ctypes et j'ai un problème de passage d'un tableau de string d'une ddl Fortran vers python.

j'ai un fichier à lire avec la dll avec la structure suivante :

aaaaa
bbbbb
ccccc
4.0 5.0 6.5
8.0 5.1 6.5
3.0 5.2 6.5

j'essaie de récupérer un tableau de string et un autre tableau de réel avec les codes suivants:

mon code fortran

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
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
 
module read_file
           USE iso_fortran_env
           USE,INTRINSIC :: ISO_C_BINDING
          implicit none
 
 
      contains
 
          subroutine OPEN_FILE(FicWxx,path_cptr,lenpath_cptr)  BIND(C)
			  !DEC$ ATTRIBUTES DLLEXPORT :: OPEN_FILE
              integer (KIND=C_INT)         , intent(in) ,value     :: FicWxx
              type(c_ptr), value :: path_cptr
              integer(c_int), value :: lenpath_cptr
              character(len=lenpath_cptr,kind=c_char), pointer :: PATH
			  logical  :: lexist
 
              call c_f_pointer(path_cptr, PATH)                                   			
 
 
              inquire(file=trim(PATH),exist=lexist)
              if(.not.lexist) then
                  stop
              else
                  open(FicWxx,file=trim(PATH))  
              end if
 
          end subroutine OPEN_FILE
 
!--------------------------- read file   ----------------------------------------
 
 
		  subroutine GET_Char_COL(FicWxx,numel,Py_Elem ) BIND(C) 
			  !DEC$ ATTRIBUTES DLLEXPORT :: GET_Char_COL
 
			  integer (KIND=C_INT)       ,value            , intent(in)    :: FicWxx
			  integer (KIND=C_INT)                         , intent(in)    :: numel
			  type(c_ptr)                , value             :: Py_Elem
              character(KIND=C_char,len=16)      ,pointer   :: Elem  
 
			  integer                                                        :: n
 
			  call c_f_pointer(Py_Elem,Elem)
 
              do n = 1, numel
                 read(FicWxx,*)Elem
             end do 
		  end subroutine GET_Char_COL
 
       !     
          subroutine GET_FLOAT(FicWxx,numel,param) BIND(C)
              !DEC$ ATTRIBUTES DLLEXPORT :: GET_FLOAT
 
              integer, parameter                                       :: Ncol=3
			  integer (KIND=C_INT)   ,value              , intent(in)  :: FicWxx
              integer (KIND=C_INT)   ,value              , intent(in)  :: numel
              real  (KIND=C_FLOAT),dimension(numel,Ncol), intent(out)  :: param
 
              integer                                                  :: n
 
              do n = 1, numel
	            read(FicWxx,*) param(n,:)
              end do
 
          end subroutine GET_FLOAT
 
      end module read_file

mon code python:

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
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
import ctypes
import numpy as np
from numpy.ctypeslib import ndpointer
import os
 
 
def Import_DLL():         
 
    PATH = os.path.dirname(os.path.realpath(__file__))
    Lib = ctypes.CDLL(os.path.join(PATH,"./DLL_read.dll"))
    return Lib
 
 
def Decla_Arg():           
 
    Lib.open_file.argtypes = [ctypes.c_int, ctypes.c_char_p,ctypes.c_int]
    Lib.open_file.restype = None
 
    Lib.get_char_col.argtypes = [ctypes.c_int, ctypes.POINTER(ctypes.c_int), ndpointer(dtype=np.dtype('a16'))]
    Lib.get_char_col.restype = None
 
    Lib.get_float.argtypes = [ctypes.c_int, ctypes.c_int, ndpointer(dtype=ctypes.c_float)]
    Lib.get_float.restype = None
 
 
def open_file(FicWxx,PATH):
    Lib.open_file(FicWxx,PATH.encode("utf-8"),len(PATH))
    return
 
def get_char_col(FicWxx,nliais):
    z=[' ' for i in range(nliais)]
    #z=[]
    elem=np.asarray(z, dtype = np.dtype('a16'))
    print(elem)
    nliais = ctypes.c_int(nliais)
    Lib.get_char_col(FicWxx, nliais, elem)
    return elem
 
def get_res_type(FicWxx, nliais):
 
    param=np.zeros(shape=(3 , nliais),dtype=ctypes.c_float)
    nliais = ctypes.c_int(nliais)
    Lib.get_float(FicWxx, nliais,param)
    return param
 
# --------------  ___main___  --------------------------------
if __name__ == "__main__":
 
    Lib=Import_DLL()
    Decla_Arg()
    nliais=3
    FicWxx = 15
    PATH =".......\\test.txt"
    open_file(FicWxx,PATH)
 
    elem = get_char_col(FicWxx, nliais)
    print("elem =", elem)
 
    param = get_res_type(FicWxx, nliais)
    print("param =", param)

j'ai une erreur : in get_char_col
Lib.get_char_col(FicWxx, nliais, elem) -> OSError: exception: stack overflow.

PS: si j'ai un fichier sans les string , le code marche bien pour récupérer que le tableau de réel.

Merci pour votre aide.

**wiztricks** · 08/09/2023, 16h17

Envoyé par lakhmesi

Merci pour votre aide.

Pour imaginer pouvoir vous aider, il faut installer le compilateur fortran, compiler votre code, reproduire l'erreur,... et "tester" des idées de solutions avant de la partager... c'est beaucoup de travail pour un intérêt discutable.

Depuis 15ans que je traîne dans le coin, j'ai jamais vu ce genre de question: j'imagine que ceux qui se lancent là dedans ont les bases pour s'en sortir et/ou savent chercher tout seul dans les documentations et Internet.

Mais bon, une bouteille à la mer n'a jamais tué personne.

- W

**lakhmesi** · 08/09/2023, 16h29

Envoyé par wiztricks

Pour imaginer pouvoir vous aider, il faut installer le compilateur fortran, compiler votre code, reproduire l'erreur,... et "tester" des idées de solutions avant de la partager... c'est beaucoup de travail pour un intérêt discutable.

Depuis 15ans que je traîne dans le coin, j'ai jamais vu ce genre de question: j'imagine que ceux qui se lancent là dedans ont les bases pour s'en sortir et/ou savent chercher tout seul dans les documentations et Internet.

Mais bon, une bouteille à la mer n'a jamais tué personne.

- W

cher monsieur ,

j'ai crée ma dll et j'ai compilé mon code. c'est comme ça que je sais que ça ne marche pas et j'ai une erreur.

j'ai déjà essayé de résoudre le problème sans sucés alors je l'ai posté ici.

si ça vous dérange , alors quitter le forum et voila .

**wiztricks** · 08/09/2023, 16h42

Envoyé par lakhmesi

si ça vous dérange , alors quitter le forum et voila .

Un forum de discussion n'est pas une hotline pour vous servir... on peut raconter ce qu'on veut sur votre prose comme dire qu'il va falloir être patient avant de voir arriver une solution.

Lorsque vous en aurez marre d'attendre, vous pourrez toujours acheter un support payant chez Intel.

- W

**lakhmesi** · 08/09/2023, 16h53

Envoyé par wiztricks

Un forum de discussion n'est pas une hotline pour vous servir... on peut raconter ce qu'on veut sur votre prose comme dire qu'il va falloir être patient avant de voir arriver une solution.

Lorsque vous en aurez marre d'attendre, vous pourrez toujours acheter un support payant chez Intel.

- W

je vous ai pas sollicité personnellement pour m'aider et me trouver une solution.

je poste dans l'optique que quelqu'un d'autre a déja rencontré le problème avant moi et me donner des pistes à explorer .

**wiztricks** · 08/09/2023, 17h20

Envoyé par lakhmesi

je poste dans l'optique que quelqu'un d'autre a déja rencontré le problème avant moi et me donner des pistes à explorer .

Un peu de recherche sur Internet montre que des tas de gens ont eu ce genre de problème. On en parle ici par exemple.

Si vous ignorez l'existence de f2py, difficile de faire croire que vous avez beaucoup cherché.

- W

**lakhmesi** · 08/09/2023, 18h37

Envoyé par wiztricks

Un peu de recherche sur Internet montre que des tas de gens ont eu ce genre de problème. On en parle ici par exemple.

Si vous ignorez l'existence de f2py, difficile de faire croire que vous avez beaucoup cherché.

- W

j'ai dèja fait un tour dans Stack Overflow avant de poster ici .

je connais bien f2py et déjà utilisé mais dans mon cas je suis obligé d'utiliser Ctypes.

Merci pour ton aide.

**lakhmesi** · 11/09/2023, 12h02

Bonjour à tous ,

j'ai trouvé la solution au problème. Je la partage , ça pourrait servir un jour

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
subroutine GET_Char_COL(FicWxx,numel,Py_Elem ) BIND(C) 
          !DEC$ ATTRIBUTES DLLEXPORT :: GET_Char_COL
 
          integer (KIND=C_INT)       ,value            , intent(in)    :: FicWxx
          integer (KIND=C_INT)                         , intent(in)    :: numel
          type(c_ptr)                , value               :: Py_Elem
          character(KIND=C_char,len=16) ,dimension(:)     ,pointer   :: Elem  
 
          integer                                                        :: n
 
          call c_f_pointer(Py_Elem,Elem,[numel])  ! numel is Rank-one array .The size must be equal to the rank of Py_Elem 
 
          do n = 1, numel
             read(FicWxx,*)Elem(n)
         end do 
      end subroutine GET_Char_COL

pssage d'un tableau de string entre Fortran DLL et python avec ctypes

Interfaçage autre langage Python

Discussions similaires

Partager

Partager