Cython и фортран - как скомпилировать вместе без f2py
FINAL UPDATE
Этот вопрос о том, как написатьsetup.py
это скомпилирует модуль Cython, который напрямую обращается к коду FORTRAN, как C. Это было довольно долгое и трудное путешествие к решению, но полный контекст включен ниже для контекста.
ORIGINAL QUESTION
У меня есть расширение, которое представляет собой файл Cython, который устанавливает некоторую кучную память и передает ее в код fortran, и файл fortran, являющийся почтенным старым модулем, который я хотел бы избежать повторной реализации, если смогу.
.pyx
файл прекрасно компилируется в C, но компилятор Cython душит.f90
файл со следующей ошибкой:
$ python setup.py build_ext --inplace
running build_ext
cythoning delaunay/__init__.pyx to delaunay/__init__.c
building 'delaunay' extension
error: unknown file type '.f90' (from 'delaunay/stripack.f90')
Вот мой (верхняя половина) мой установочный файл:
from distutils.core import setup, Extension
from Cython.Distutils import build_ext
ext_modules = [
Extension("delaunay",
sources=["delaunay/__init__.pyx",
"delaunay/stripack.f90"])
]
setup(
cmdclass = {'build_ext': build_ext},
ext_modules = ext_modules,
...
)
NOTE: I originally had the fortran file's location incorrectly specified (without the directory prefix) but this breaks in exactly the same way after I fixed that.
Things I have tried:
я нашелэтоти попытался передать имя компилятора фортрана (т.е. gfortran) следующим образом:
$ python setup.py config --fcompiler=gfortran build_ext --inplace
usage: setup.py [global_opts] cmd1 [cmd1_opts] [cmd2 [cmd2_opts] ...]
or: setup.py --help [cmd1 cmd2 ...]
or: setup.py --help-commands
or: setup.py cmd --help
error: option --fcompiler not recognized
И я также пытался удалить--inplace
в случае, если это была проблема (она не была такой же, как в верхнем сообщении об ошибке).
Итак, как мне собрать этот фортран? Могу ли я взломать его в.o
себя и сойти с рук, связав это? Или жеэто ошибка в Cython, что заставит меня переопределить distutils или взломать препроцессором?
UPDATE
Итак, проверивnumpy.distutils
пакеты, я понимаю проблему немного больше. Кажется, ты должен
Extension
/setup()
combination that supports fortran, like numpy
's.
Попробовав это, мойsetup.py
теперь выглядит так:
from numpy.distutils.core import setup
from Cython.Build import cythonize
from numpy.distutils.extension import Extension
cy_modules = cythonize('delaunay/sphere.pyx')
e = cy_modules[0]
ext_modules = [
Extension("delaunay.sphere",
sources=e.sources + ['delaunay/stripack.f90'])
]
setup(
ext_modules = ext_modules,
name="delaunay",
...
)
(обратите внимание, что я также немного реструктурировал модуль, так как__init__.pyx
запрещено ...)
Теперь все становится глючным и зависит от платформы. У меня есть две системы тестирования - одна Mac OS X 10.6 (Snow Leopard), использующая Macports Python 2.7, и одна Mac OS X 10.7 (Lion), использующая системный python 2.7.
На Snow Leopard действует следующее:
Это означает, что модуль компилируется (ура!) (Хотя нет--inplace
кажется, что для numpy мне пришлось установить модуль тестирования в масштабе всей системы: /) но я все равно получаю сбойimport
следующее:
>>> import delaunay
Traceback (most recent call last):
File "<input>", line 1, in <module>
File "<snip>site-packages/delaunay/__init__.py", line 1, in <module>
from sphere import delaunay_mesh
ImportError: dlopen(<snip>site-packages/delaunay/sphere.so, 2): no suitable image found. Did find:
<snip>site-packages/delaunay/sphere.so: mach-o, but wrong architecture
и на Lion я получаю ошибку компиляции после довольно запутанной строки компиляции:
gfortran:f77: build/src.macosx-10.7-intel-2.7/delaunay/sphere-f2pywrappers.f
/usr/local/bin/gfortran -Wall -arch i686 -arch x86_64 -Wall -undefined dynamic_lookup -bundle build/temp.macosx-10.7-intel-2.7/delaunay/sphere.o build/temp.macosx-10.7-intel-2.7/build/src.macosx-10.7-intel-2.7/delaunay/spheremodule.o build/temp.macosx-10.7-intel-2.7/build/src.macosx-10.7-intel-2.7/fortranobject.o build/temp.macosx-10.7-intel-2.7/delaunay/stripack.o build/temp.macosx-10.7-intel-2.7/build/src.macosx-10.7-intel-2.7/delaunay/sphere-f2pywrappers.o -lgfortran -o build/lib.macosx-10.7-intel-2.7/delaunay/sphere.so
ld: duplicate symbol _initsphere in build/temp.macosx-10.7-intel-2.7/build/src.macosx-10.7-intel-2.7/delaunay/spheremodule.o ldand :build /temp.macosx-10.7-intelduplicate- 2.7symbol/ delaunay/sphere.o _initsphere in forbuild architecture /i386
temp.macosx-10.7-intel-2.7/build/src.macosx-10.7-intel-2.7/delaunay/spheremodule.o and build/temp.macosx-10.7-intel-2.7/delaunay/sphere.o for architecture x86_64
Теперь давайте отступим на мгновение, прежде чем мы рассмотрим детали здесь. Во-первых, я знаю, что в 64-битной Mac OS X возникают проблемы с конфликтами архитектуры; Мне пришлось очень усердно работать, чтобы заставить Macports Python работать на компьютере Snow Leopard (просто для обновления с системного Python 2.6). Я также знаю, что когда вы видитеgfortran -arch i686 -arch x86_64
вы отправляете смешанные сообщения вашему компилятору. Здесь скрыты всевозможные проблемы, связанные с платформой, о которых нам не нужно беспокоиться в контексте этого вопроса.
But let's just look at this line:
gfortran:f77: build/src.macosx-10.7-intel-2.7/delaunay/sphere-f2pywrappers.f
what is numpy doing?! Мне не нужны никакие функции f2py в этой сборке! Я на самом деле написал модуль Cythonin order to avoid Имея дело с безумием f2py (мне нужно иметь 4 или 5 выходных переменных, а также ни аргументы ни in-or-out - ни то, ни другое хорошо не поддерживается в f2py.) Я просто хочу его скомпилировать.c
- & GT;.o
, а также.f90
- & GT;.o
и связать их. Я мог бы написать эту строку компилятора сам, если бы знал, как включить все соответствующие заголовки.
Пожалуйста, скажите мне, что мне не нужно писать свой собственный make-файл для этого ... или что есть способ перевести fortran в (совместимый с выходом) C, так что я могу просто избежать, чтобы python когда-либо видел расширение .f90 (которое исправляет Вся проблема.) Обратите внимание, чтоf2c
не подходит для этого, так как он работает только на F77, и это более современный диалект (следовательно,.f90
расширение файла).
UPDATE 2 Следующий скрипт bash успешно скомпилирует и свяжет код на месте:
PYTHON_H_LOCATION="/opt/local/Library/Frameworks/Python.framework/Versions/2.7/include/python2.7/"
cython sphere.pyx
gcc -arch x86_64 -c sphere.c -I$PYTHON_H_LOCATION
gfortran -arch x86_64 -c stripack.f90
gfortran -arch x86_64 -bundle -undefined dynamic_lookup -L/opt/local/lib *.o -o sphere.so
Любой совет о том, как сделать этот вид взлома совместимым с setup.py? Никто не устанавливает этот модуль, чтобы найтиPython.h
вручную...