Fortran DLL в Delphi

Я не знаю Delphi, но вот что вы можете сделать, чтобы создать DLL, доступную из языка C. Я надеюсь, что вы найдете это полезным. Вот ваш код F77, измененный, чтобы сделать его совместимым с функциями модуля iso_c_binding Fortran:

      SUBROUTINE MYSUB1(NoEquations,INTARR1) bind(C,name = "MYSUB1")
         !DEC$ ATTRIBUTES DLLEXPORT :: MYSUB1
         use, intrinsic :: iso_c_binding, only: IK => c_int32_t
         integer(IK), intent(in), value :: NoEquations
         integer(IK), intent(out)       :: INTARR1(NoEquations)
         integer                        :: I
         DO 100, I=1,NoEquations
             INTARR1(I) = I
  100    CONTINUE
         RETURN
      END SUBROUTINE MYSUB1

C***************************************************************
C***************************************************************

      SUBROUTINE MYSUB2(NoEquations,INTARR1,INTARR2)
     +bind(C,name = "MYSUB2")
         !DEC$ ATTRIBUTES DLLEXPORT :: MYSUB2
         use, intrinsic :: iso_c_binding, only: IK => c_int32_t
         integer(IK), intent(in), value :: NoEquations
         integer(IK), intent(in)        :: INTARR1(NoEquations)
         integer(IK), intent(out)       :: INTARR2(NoEquations)
         integer                        :: I
         DO 100, I=1,NoEquations
             INTARR2(I) = INTARR1(I)
  100    CONTINUE
         RETURN
      END SUBROUTINE MYSUB2

Обратите внимание на множество тонких, но важных изменений, которые я внес в ваш код, чтобы сделать его совместимым с C:

1. bind(C,name = "MYSUB1") исправляет имя подпрограммы, поэтому вам не нужны дополнительные вторые директивы компилятора (которые я удалил из вашего кода).
2. Также обратите внимание на атрибут value, который указывает компилятору передавать по значению, как это делается в C.
3. Также обратите внимание, что я определяю типы целых чисел в интерфейсах подпрограммы как c_int32_t, чтобы быть совместимыми с типами целых чисел процессора C.
4. Также обратите внимание, что я преобразовал ваши массивы предполагаемого размера INTARR1 и INTARR2 в массивы явной формы INTARR1(NoEquations), INTARR1(NoEquations).

Поскольку в вашем коде есть директивы компилятора Intel !DEC$, я предполагаю, что вы используете компилятор Intel Fortran в Windows. Обратите внимание, что я удалил ваши вторые строки директив в обеих подпрограммах Fortran.

Теперь, предположив, что вы сохранили приведенный выше код в файле с именем mysubs.F, а затем скомпилировали этот файл с помощью компилятора Intel Fortran ifort в командной строке Intel-Windows, как в следующей команде,

ifort mysubs.F /dll /out:libsubs

создаст DLL с именем mysubs.dll в текущей папке и выведет на экран следующее сообщение,

Intel(R) Visual Fortran Intel(R) 64 Compiler for applications running on Intel(R) 64, Version 19.1.1.216 Build 20200306
Copyright (C) 1985-2020 Intel Corporation.  All rights reserved.

Microsoft (R) Incremental Linker Version 14.16.27027.1
Copyright (C) Microsoft Corporation.  All rights reserved.

-out:mysubs.dll
-dll
-implib:mysubs.lib
mysubs.obj
   Creating library mysubs.lib and object mysubs.exp

Чтобы протестировать эту DLL, вы можете попробовать следующий код C, хранящийся в main.c,

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

void MYSUB1(int32_t, int32_t []);
void MYSUB2(int32_t, int32_t [], int32_t []);

int main(int argc, char *argv[])
{
    const int32_t NoEquations = 5;
    int32_t INTARR1[NoEquations];
    int32_t INTARR2[NoEquations];
    int loop;

    // C rules for argument passing apply here

    MYSUB1(NoEquations,INTARR1);
    printf("\nINTARR1:\n"); for(loop = 0; loop < NoEquations; loop++) printf("%d ", INTARR1[loop]);

    MYSUB2(NoEquations,INTARR1,INTARR2);
    printf("\nINTARR2:\n"); for(loop = 0; loop < NoEquations; loop++) printf("%d ", INTARR2[loop]);

    return 0;
}

Компилируя этот код C с помощью компилятора Intel C,

icl main.c -c

печатает следующее на экране,

Intel(R) C++ Intel(R) 64 Compiler for applications running on Intel(R) 64, Version 19.1.1.216 Build 20200306
Copyright (C) 1985-2020 Intel Corporation.  All rights reserved.

main.c

и генерирует основной объектный файл C. Наконец, свяжите объектный файл C с библиотекой DLL Fortran, чтобы сгенерировать исполняемый файл с помощью следующей команды:

icl main.obj mysubs.lib -o main.exe

который печатает следующее на экране,

Intel(R) C++ Intel(R) 64 Compiler for applications running on Intel(R) 64, Version 19.1.1.216 Build 20200306
Copyright (C) 1985-2020 Intel Corporation.  All rights reserved.

Microsoft (R) Incremental Linker Version 14.16.27027.1
Copyright (C) Microsoft Corporation.  All rights reserved.

-out:main.exe
main.obj
mysubs.lib

Проверка вызова DLL. Просто вызовите сгенерированный исполняемый файл main.exe,

main.exe

который печатает на экране,

INTARR1:
1 2 3 4 5
INTARR2:
1 2 3 4 5

Теперь, чтобы вызвать эту DLL из Delphi, просто предположим, что вы вызываете функции C с прототипами, указанными в основном коде C. Вот и все. Больше не нужно иметь дело с Fortran из Delphi.

Заключительный совет:

1. Fortran имеет мощные стандартные функции взаимодействия, подобные тем, которые я добавил в ваш код F77, которые могут легко связать практически любой код Fortran с любым языком (через C).

2. Держитесь подальше от FORTRAN77, которому почти полвека, даже Fortran 90 уже более 30 лет. Последний стандарт Fortran был выпущен в 2018 году, что вместе с Fortran 2008 делает Fortran чрезвычайно мощным, высокоуровневым, быстрым, изначально векторизованным, параллельным, общим и распределенным параллельным языком программирования для численных вычислений.