ctypesを使用してFortranサブルーチンにPython文字列を渡す

Question

ctypesを使用して共有ライブラリのFortranサブルーチンに引数を渡そうとしています。今ここに私の簡単なFortranコードです：私はこれが正常に動作するようです

gcc -shared -fPIC -o libtest.o fstr_test.f90 

でこれをコンパイル

MODULE test_module 
INCLUDES 
SUBROUTINE fstr_test(file_or_extension, ierr, iopen) 
    IMPLICIT NONE 
    INTEGER, INTENT(out) :: ierr 
    INTEGER, OPTIONAL :: iopen 
    CHARACTER(LEN=*), INTENT(in) :: file_or_extension 
    WRITE(6,*) file_or_extension,ierr,iopen 
    RETURN 
END SUBROUTINE fstr_test 
END MODULE test_module 

、私は実際にはPython 3での共有ライブラリをロードすることができますここでは、コード非クラッシュです：

は、

from ctypes import * 
libtest = cdll.LoadLibrary("libtest.so") 
f = libtest.test_module_MOD__fstr_test 
file = b'wout_li383.nc' 
ierr = c_int(0) 
iopen = c_int(0) 
cfile=c_char_p(b"test_str") 
f(cfile,byref(ierr),byref(iopen)) 

ただし、文字列を正しく渡すことができません。上の形式では、文字列は空であり、2つの整数はコンソールに正しく出力されます。基本的に私が行った唯一の進歩は、Pythonをクラッシュさせることです。私はここで解決策を試してみました：

Passing string to Fortran DLL using ctypes and Python

が、それはPythonの2話して、多分ので、私のために動作するようには思えないのですか？例

f.argtype = [c_char_p,c_int_c_int] 
c = c_char_p(b"test.txt") 
f(c,byref(ierr),byref(iopen)) 

については

は、カーネルをクラッシュしませんが、文字列は空です。 byrefに変更する：

同じ効果。長さ引数を追加して長さを渡して渡すと、 "Dead Kernel"というメッセージが表示され、カーネルを再起動する必要があります。 seのエラーメッセージはありません。

f.argtype = [c_char_p,c_int,c_int_c_int] 
file = b'text.txt' 
c = c_char_p(file) 
c_len = c_int(len(file)) 
f(c,byref(c_len),byref(ierr),byref(iopen)) 

これは明らかです。私はエラーメッセージを受け取っていない、カーネルはちょうど死んでいる。

Answer 1

いくつかの問題は、引数リストの最後に隠れた長さ（少なくともgfortranが他のコンパイラで何をしているか不明）が必要です。あなたはどちらかそれにByRefのを必要としないと

SUBROUTINE fstr_test(file_or_extension, ierr, iopen) 
    IMPLICIT NONE 
    INTEGER, INTENT(out) :: ierr 
    INTEGER, OPTIONAL :: iopen 
    CHARACTER(LEN=*), INTENT(in) :: file_or_extension 
    WRITE(6,*) file_or_extension,ierr,iopen 
    ierr=1 
    RETURN 
    END SUBROUTINE fstr_test 

と

gfortran -shared -fPIC -o libtest.so fstr_test.f90 

次にサブルーチンでIERR変数を変更c_char_p

でのpythonをバイト文字列をラップしていませんコードは次のようにする必要があります。

import ctypes 

lib=ctypes.CDLL("./libtest.so") 
f=getattr(lib,'__test_module_MOD_fstr_test') 

f.argtypes=[ctypes.c_char_p,ctypes.c_int,ctypes.c_int,ctypes.c_int] 
f.restype=None 

x=b'abcdef' 
ierr=ctypes.c_int(0) 
iopen=0 
f(x,ctypes.byref(ierr),iopen,len(x)) 
print(ierr.value) 

私はオプションの引数letを取得していませんが、だから私は空のintを追加しました、これはまだ文字列を渡す必要がありますし、ierrを取得します。