Несоответствие интерфейса в фиктивной процедуре 'f' при передаче функции в подпрограмму

Как насчет:

interface
  real pure function f(y)
    real, dimension(:), intent(in) :: y
  end function
end interface

Когда вы передаете аргумент solve в функцию, размер массива будет передан автоматически. Вам не нужно делать эту часть интерфейса.

Если вы хотите обеспечить безопасность, как указано в вашем комментарии к решению MSB, вам следует использовать -fcheck=bounds, и компилятор сгенерирует проверки во время выполнения для предполагаемых и отложенных массивов форм. См. справочную страницу gfortran для получения дополнительной информации о -fcheck. Однако вы немного потеряете в скорости.

У вас есть решение, но зачем его объяснять... если фиктивный аргумент имеет явный интерфейс (что он и делает здесь), то есть требование, чтобы характеристики процедуры передавались как фактический аргумент должен совпадать с аргументом фиктивного аргумента, с некоторыми исключениями, связанными с чистотой и элементарными внутренними свойствами. Среди прочего, характеристики процедуры включают в себя характеристики ее фиктивных аргументов.

К характеристикам фиктивного аргумента, среди прочего, относится его форма. Если эта форма не является постоянным выражением, характеристики включают «точную зависимость [формы] от сущностей в выражении».

Блок интерфейса для фиктивного аргумента f объявляет размер массива SIZE(x). x — это предполагаемая переменная формы, связанная с хостом — ее размер может меняться во время выполнения, поэтому SIZE(x) не является константой. Следовательно, это выражение и сущности в нем становятся характеристикой фиктивного аргумента.

Процедура модуля g объявляет массив размером 2. Это явно константа.

Независимо от значения непостоянного выражения для размера фиктивного аргумента f, эти характеристики размера массива (какое-то выражение против константы) не совпадают - отсюда и ошибка.

Когда вы заменяете SIZE(x) константой 2, характеристики, очевидно, совпадают. Когда вы изменяете предполагаемую форму x на постоянный размер 2, тогда SIZE(x) становится постоянным выражением со значением 2, поскольку это постоянное выражение, все, что имеет значение, это его значение, поэтому характеристики двух аргументов затем совпадают. Когда вы изменяете как фиктивный аргумент f, так и фиктивный аргумент g, чтобы принять форму (:), характеристики совпадают.

Вот демонстрация, показывающая, как передавать выделяемый массив.

Некоторые советы:

• Используйте модули, чтобы избежать громоздкого интерфейса.
• Добавьте дополнительную информацию о размере матрицы при передаче массива в фактическую функцию. Например, f(y, sizeinfo), чтобы внутри вашей фактической функции вы могли правильно объявить размер входной матрицы. Распределяемый массив можно передать подпрограмме solve, поэтому размер можно получить с помощью size(mat) в вашей подпрограмме solve.

Итак, исправленная версия выглядит так:

module foo

contains

  subroutine solve(x, f)
    real, dimension(:), intent(inout) :: x
    real,external::f
        integer::sizeinfo

    print *,'x=', x
    sizeinfo = size(x)
    print *, 'f(x)=',f(x,sizeinfo)
  end subroutine

    real function g(y,sizeinfo)
    integer::sizeinfo
    real, dimension(sizeinfo) :: y

    g = sum(y)
  end function
end module

Вот основная программа:

program main
use foo

real, dimension(2) :: x = (/1.0, 2.0/)

call solve(x, g)

end program

И результат:

x=   1.000000       2.000000    
f(x)=   3.000000