F2PY - Zugriff auf Modulparameter aus Unterprogramm

Question

Ich kann f2py nicht dazu bringen, einen Parameter von einem Modul in einer separaten Unterroutine zu referenzieren, wo es verwendet wird, um eine Eingabearraydimension zu definieren. I.e.

! File: testmod.f90 
MODULE testmod 
INTEGER, PARAMETER :: dimsize = 20 
END MODULE testmod 

und der Parameter dimsize in einem Unterprogramm referenziert werden muss (nicht im Modul enthielt) in einer anderen Datei, die der Einstiegspunkt für meinen Python-Modul sein wird: die Paramter werden in ein Modul defeind

! File testsub.f90 
SUBROUTINE testsub(arg) 
USE testmod 
REAL, INTENT(IN) :: arg(dimsize) 
END SUBROUTINE testsub 

ich kompilieren wie folgt aus:

f2py -m testmod -h testmod.pyf testsub.f90 
pgf90 -g -Mbounds -Mchkptr -c -fPIC testmod.f90 -o testmod.o 
pgf90 -g -Mbounds -Mchkptr -c -fPIC testsub.f90 -o testsub.o 
f2py -c testmod.pyf testmod.o testsub.o 

aber diese Fehlermeldung erhalten:

testmodmodule.c: In function 'f2py_rout_testmod_testsub': 
testmodmodule.c:180: error: 'dimsize' undeclared (first use in this function) 

Ich habe versucht, testsub.g90 enthalten die folgende Anweisung, wie vorgeschlagen ni andere Beiträge ändern:

SUBROUTINE testsub(arg) 
USE testmod 
!f2py integer, parameter :: dimsize 
REAL, INTENT(IN) :: arg(dimsize) 
END SUBROUTINE testsub 

, aber ohne Erfolg. Ich muss das Unterprogramm vom Modul getrennt halten.

Wie bekomme ich f2py, um die Variable dimsize korrekt aufzulösen?

TIA

Answer 1

Obwohl ich es nicht getestet habe, Ich denke Sie haben es fast mit Ihrem ursprünglichen Code. Wir machen etwas ähnliches für einige unserer Code, aber mit gfortran.

Sie sollten die Datei testmod.f90 nicht benötigen f2py. Sie sollten kompilieren es nur zu einer Objektdatei wie würden Sie, wenn das normal Fortran waren:

pgf90 -g -Mbounds -Mchkptr -c -fPIC testmod.f90 -o testmod.o 

Dann sollten Sie in der Lage sein, Ihre testsub.f90 in ein Python-verwendbares Modul zu kompilieren mit:

f2py --fcompiler=pgf90 --f90flags="-g -Mbounds -Mchkptr" -c testsub.f90 -m testsub testmod.o 

Das sollte eine testsub.so, oder gleichwertiges, erstellen lassen Sie import testsub und dann testsub.testsub(my_arg) in Python.

Answer 2

Es ist lange her, seit diese Frage irgendeine Aktivität bekommen hat, aber ich habe das Problem herausgefunden, also dachte ich, ich würde es für jeden anderen veröffentlichen, der dieses Problem hat. Das Problem ist, dass, während der Fortran-Code vollkommen gültig ist, der von F2PY erzeugte C-Wrapper nicht weiß, was Dimsize ist, wenn er überprüft, dass die Arrays die richtige Größe haben (innerhalb des C-Wrappers ist er in irgendeine Art von Modulstruktur eingekapselt).

einfach diesen Code ändern:

SUBROUTINE testsub(arg) 
USE testmod 
!f2py integer, parameter :: dimsize 
REAL, INTENT(IN) :: arg(dimsize) 
END SUBROUTINE testsub 

dazu:

SUBROUTINE testsub(arg) 
USE testmod 
!f2py integer, intent(aux) :: dimsize 
REAL, INTENT(IN) :: arg(dimsize) 
END SUBROUTINE testsub 

Dies wird "Define Hilfs C Variable in f2py Wrapperfunktion erzeugt". Wie die Scipy docs sagen, ist dies "nützlich, um Parameterwerte zu speichern, so dass sie in Initialisierungsausdruck anderer Variablen zugegriffen werden können".

Ich habe dies nicht mit Ihrem Code getestet, aber ich habe es mit einer ähnlichen Situation getestet, die ich hatte.