Вызов функции C из Fortran с помощью указателей

Question

Я пытаюсь создать функцию в C, которая вызывается из Fortran. Прежде всего. Код fortran скомпилируется с помощью f77, а c скомпилирован с gcc. Оба они скомпилированы в .so библиотеки. Функция c будет считывать память на устройстве по адресу и размеру, указанному fortran. Я могу видеть адрес, и размер передается в c, но у меня возникают проблемы с заполнением данных в c-функции и возвратом его в fortran. См. Соответствующий код ниже. Мое предположение заключается в том, что в распределении памяти или в синтаксисе указателя для переменной данных есть что-то неправильное.

C

void copymemory_(uint32_t *addr, int *size, uint8_t *data) 
{ 

    int i; 
    printf("addr %x \n", *addr); 
    printf("size %i \n", *size); 

    for (i = 0; i<*size; i++) 
    { 
     *data[i] = i; 
     printf("memory %i \n",*data[i]); 
    } 

} 

Fortran

integer memory(4) 
call copymemory(z'dead', 4, memory) 

DO i = 1,memsize 
    call printf(memory(i)) 
END DO 

Answer 1

У меня есть несколько точек в код.

Пожалуйста, поставьте компилируемый код и опишите вывод этого точного кода!

Это касается #include для стандартных заголовков, если вы хотите, чтобы люди отлаживали ваш код, упростили их, чтобы им не пришлось искать, какие строки вы пропустили, потому что они казались вам «очевидными». Просто вставьте все. Код должен компилироваться при копировании отсюда!

Даже исполняемая часть вашего кода не компилируется в моем компиляторе. Мне пришлось изменить *data на data. Вы действительно скопировали свой фактический код?

cfun.c: In function ‘copymemory_’: 
cfun.c:13:9: error: invalid type argument of unary ‘*’ (have ‘int’) 
     *data[i] = i; 
     ^
cfun.c:14:31: error: invalid type argument of unary ‘*’ (have ‘int’) 
     printf("memory %i \n",*data[i]); 
         ^

Ваш Fortran код содержит вызов некоторой printf подпрограммы. Где это определено? Он присутствует в вашем фактическом коде? Я сомневаюсь. Скопируйте в StackOverflow всегда полный и компилируемый код.

Таким образом, после фиксации очевидного вашего настоящего кода является:

#include <stdio.h> 
#include <stdint.h> 

void copymemory_(uint32_t *addr, int *size, uint8_t *data) 
{ 

    int i; 
    printf("addr %x \n", *addr); 
    printf("size %i \n", *size); 

    for (i = 0; i<*size; i++) 
    { 
     data[i] = i; 
     printf("memory %i \n",data[i]); 
    } 

} 


implicit none 

integer :: memsize = 4 

integer memory(4) 
integer i 

call copymemory(z'dead', 4, memory) 

DO i = 1,memsize 
    print *, (memory(i)) 
END DO 

END 

Он не откажет, но memory в Fortran содержит мусор. Это необходимо, потому что это integer, и вы рассматриваете его как int8_t. C.

Так что либо рассматривайте его как массив из четырех целых чисел в C, либо скопируйте его байтом по байту, но затем вы должны передать правильное количество байтов для копирования. Из вашего описания не ясно, какой из них является ваше намерение, так что я покажу только одну возможность:

void copymemory_(uint32_t *addr, int *size, uint32_t *data) 

Выход правильно, то:

> gfortran-4.10 -fsanitize=address cfun.c ffun.f90 
> ./a.out 
addr dead 
size 4 
memory 0 
memory 1 
memory 2 
memory 3 
      0 
      1 
      2 
      3