интероперабельность c-fortran - производные типы с указателями

Question

У меня есть длинный код fortran, который должен быть использован для использования с python. Я решил использовать интерфейс Fortran-> C-> Python.

У меня проблема: я вывел типы в Fortran-модулей, которые содержат

double precision, allocatable 

типа в качестве членов.

При попытке компиляции с ifort я получаю (с gfortran что-то похожее):

Each component of a derived type with the BIND attribute shall be a nonpointer,   
nonallocatable data component with interoperable type and type parameters 

Это на самом деле с соглашением с компилятором документации Intel и стандарт Fortran 2003: пункт 15.2.5.

Есть ли способ получить доступ к «типу» Фортрана с выделенным или указателем внутри?

Answer 1

Не совсем. Даже совершенно новый TS 29113 не предлагает решения. Он имеет TYPE(*), но вы должны знать, как декодировать структуру в C. Вы можете сделать это путем обратного проектирования.

На самом деле существует макрос для размера структуры в упомянутых TS. См. CFI_attribute_allocatable. Но проблема с разыменованием остается. Существует стандартный дескриптор массива, но этот TS еще не поддерживается некоторыми компиляторами (в частности, gfortran).

Возможно, просто передать указатель на структуру общей функции, такой как qsort. Затем вы используете только часть памяти, которая, как оказалось, содержит также некоторый указатель или выделяемый дескриптор, но C не должен знать о них.

Другая возможность состоит в том, чтобы построить производный тип, который содержит только type(c_ptr) вместо указателей Фортрана и использует c_loc(), чтобы заполнить их. Тогда обязательно должны быть только смежные массивы.