Ошибка компилятора gfortran при вызове рекурсивной функции log_gamma

Я загрузил гипергеометрическую функцию 2F1, но оказалось, что она не компилируется на моем компьютере. Это из Эта статья. я использую

GNU Fortran (Built by Jeroen for the R-project) 8.3.0

поставляется с RTools 4.0. Можете ли вы понять, почему он не компилируется и как это можно решить? В приведенном ниже коде я просто оставил одну строку, которая генерирует ошибку. Ошибка указана далее.

MODULE HYP_2F1_MODULE
  !--------------------------------------------------------------------
  IMPLICIT NONE
  INTEGER, PARAMETER :: PR=KIND(1.0D0)
  REAL(PR)     ::  ONE=1.0D0
CONTAINS
  !
END MODULE HYP_2F1_MODULE
!
!----------------------------------------------------------------------
RECURSIVE FUNCTION LOG_GAMMA(Z) RESULT(RES)

  USE HYP_2F1_MODULE
  IMPLICIT NONE
  COMPLEX(PR),INTENT(IN) :: Z
  COMPLEX(PR)  :: RES
  !
     RES = LOG_GAMMA( ONE -z);
  !
END FUNCTION LOG_GAMMA

Вот сообщение об ошибке

 testZ.f90:18:22:

      RES = LOG_GAMMA( ONE - Z);
                 1
Error: 'x' argument of 'log_gamma' intrinsic at (1) must be REAL

Сообщение об ошибке, показанное в вашем вопросе, не соответствует коду, который вы разместили.

Bob Jarvis - Слава Україні 02.04.2022 17:01

@bob только что обновил код, затем сообщение об ошибке

Denis Cousineau 02.04.2022 17:03

Ну, это не компилируется, и я не могу понять, почему. Это взято (как есть) из репозитория НАСА. Я пытался преобразовать ONE в сложный CMPLX (ONE, ZERO) перед вычитанием, но все равно не компилируется.

Denis Cousineau 02.04.2022 17:09

RES LOG_GAMMA четко указано как COMPLEX, поэтому я не понимаю, почему в нем говорится, что (1) должно быть REAL.

Denis Cousineau 02.04.2022 17:09

О, я вижу, я не заметил, что сама функция называется LOG_GAMMA. LOG_GAMMA также является встроенной функцией Fortran 2018. Интересно, почему внутреннее не затенено здесь.

Vladimir F Героям слава 02.04.2022 17:11

Разве функция не должна затенять внутреннюю?

veryreverie 02.04.2022 17:12

Проблема будет в рекурсивном атрибуте и связанном предложении результата и может быть ошибкой компилятора.

Vladimir F Героям слава 02.04.2022 17:13

Хороший улов, я не знал, что LOG_GAMMA встроен. Переименую без проблем. Если вы хотите превратить это в ответ? Спасибо.

Denis Cousineau 02.04.2022 17:15

Я предполагаю, что решения включают в себя использование более нового компилятора или использование -std=f95 для отключения встроенного.

veryreverie 02.04.2022 17:15

@veryreverie Очень хороший момент, для старых кодов -std=f95 будет работать. Даже -f2003 достаточно.

Vladimir F Героям слава 02.04.2022 17:19

Кажется, что правильная функция вызывается, если типы совпадают.

Vladimir F Героям слава 02.04.2022 17:22
Стоит ли изучать PHP в 2023-2024 годах?
Стоит ли изучать PHP в 2023-2024 годах?
Привет всем, сегодня я хочу высказать свои соображения по поводу вопроса, который я уже много раз получал в своем сообществе: "Стоит ли изучать PHP в...
Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией
Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией
В JavaScript одним из самых запутанных понятий является поведение ключевого слова "this" в стрелочной и обычной функциях.
Приемы CSS-макетирования - floats и Flexbox
Приемы CSS-макетирования - floats и Flexbox
Здравствуйте, друзья-студенты! Готовы совершенствовать свои навыки веб-дизайна? Сегодня в нашем путешествии мы рассмотрим приемы CSS-верстки - в...
Тестирование функциональных ngrx-эффектов в Angular 16 с помощью Jest
В системе управления состояниями ngrx, совместимой с Angular 16, появились функциональные эффекты. Это здорово и делает код определенно легче для...
Концепция локализации и ее применение в приложениях React ⚡️
Концепция локализации и ее применение в приложениях React ⚡️
Локализация - это процесс адаптации приложения к различным языкам и культурным требованиям. Это позволяет пользователям получить опыт, соответствующий...
Пользовательский скаляр GraphQL
Пользовательский скаляр GraphQL
Листовые узлы системы типов GraphQL называются скалярами. Достигнув скалярного типа, невозможно спуститься дальше по иерархии типов. Скалярный тип...
1
11
75
2
Перейти к ответу Данный вопрос помечен как решенный

Ответы 2

Ответ принят как подходящий

Существует одноименная встроенная функция Fortran 2008+.

16.9.119 LOG_GAMMA (X)
1 Description. Logarithm of the absolute value of the gamma function.
2 Class. Elemental function.
3 Argument. X shall be of type real. Its value shall not be a negative integer or zero.

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

Также подумайте, не может ли встроенная функция с тем же именем делать то, что вам нужно от вашей функции.

Это ошибка. F2018, 15.6.2.2, стр. 319: если появляется RESULT, имя результата функции — имя-результата, и все вхождения имени функции в операторах исполнительной части в его области действия относятся к самой функции.

steve 02.04.2022 18:06

@steve сообщил как gcc.gnu.org/bugzilla/show_bug.cgi?id=105138

Vladimir F Героям слава 02.04.2022 18:22

Спасибо. Я добавил цитируемый текст для F2018 в контрольный журнал.

steve 02.04.2022 18:55

Был отправлен патч для исправления gfortran. См. URL-адрес опубликованного PR @Vladimir. Не знаю, когда, и даже если, он превратится в дерево.

steve 03.04.2022 01:29

В дополнение к предложениям @VladimirF вы можете добавить флаг компилятора -std=f95, чтобы обеспечить соответствие стандарту Fortran 95. Если вы не добавили -fall-intrinsics, это отключит все встроенные функции из стандартов после Fortran 95 и заставит работать устаревший код.

Также -std=f2003 в этом случае, если нужно что-то из Fortran 2003.

Vladimir F Героям слава 02.04.2022 17:26

Другие вопросы по теме