Блог

Вот программа, которую я написал для тестирования многопоточного чтения файла независимо:

 program he
    use omp_lib

    implicit none
    integer :: i, a(20), funit

    open(unit=2, file='dat/he', status='replace')
    write(2, '(i2/, i2/, i2/, i2)') 1, 3, 5, 7
    close(2)

    !$omp parallel do num_threads(20) private(funit)
    do i = 1, 20
        funit=100 omp_get_thread_num()
        open(funit, file='dat/he', status='old')
        read(funit,'(i2)') a(i)
    end do
    !$omp end parallel do

    print *, a
end program he
  

Я скомпилировал его с помощью компилятора gfortran без каких-либо ошибок: gfortran -fopenmp -o exe a.f90 , и запустил ./exe , но при каждом запуске получал разные ошибки ./exe . Тем не менее, я также скомпилировал его с помощью ifort: ifort -qopenmp -o exe a.f90 , и успешно запустил его ./exe .

Одна из ошибок gfortran во время выполнения

 At line 15 of file fortran_learning_book.f90 (unit = 102)
Fortran runtime error: File already opened in another unit

Program aborted. Backtrace:

Program aborted. Backtrace:
#0  0x10c30aff2
#0  0x10c30aff2

Program aborted.
  

Ваш код не является допустимым стандартным Fortran и поэтому демонстрирует неопределенное поведение с различными компиляторами. Стандарт Fortran довольно ясен, что:

Если файл уже подключен к модулю, оператор OPEN для этого файла с другим модулем не должен выполняться.

(Раздел 9.5.6.1 из ISO / IEC 1539-1: 2010)

В отличие от некоторых других языков, стандартный Fortran не позволяет открывать один и тот же файл несколько раз.

Редактировать: Как отметил IanH, предстоящий стандарт Fortran (неофициально известный как Fortran 2015) официально снимет ограничение, но все равно оставит на усмотрение реализации решение о том, разрешено ли несколько подключений к одному и тому же файлу или нет.

Как GNU Fortran, так и Intel Fortran способны компилировать правильный код Fortran с расширениями OpenMP.

Я подозреваю, что причина вашей проблемы заключается в том, что вы не указали никаких action параметров в своем open заявлении. Действительно, действие по умолчанию для разрешения (или read , write или readwrite ) для вновь открытого файла определяется реализацией. Таким образом, вы потенциально открываете несколько раз параллельно один и тот же файл в readwrite режиме, например, позволяя, в частности write , режим из нескольких независимых модулей, что приводит к потенциальным несоответствиям файлов. Тогда компилятор и / или библиотека времени выполнения могут проверять (или не проверять) согласованность такого рода и сообщать об этом как об ошибке… Судя по тому, как это выглядит, gfrotran и ifort имеют разный подход по крайней мере в одной из двух частей (действие по умолчанию и / или проверка RTL).

Единственный случай, когда такое параллельное открытие одного и того же файла является «законным» и безопасным, — это если вы ограничиваете свой доступ режимом чтения. Поэтому добавление action='read' в ваш open оператор может просто устранить проблему. Действительно, он делает это для меня с помощью gfortran версии 5.4.0.

Теперь мы находимся в серой области (по крайней мере, для меня), и я действительно не знаю, насколько все это совместимо с переносимым / стандартным. Поэтому, хотя я считаю, что это решает вашу непосредственную проблему, я не могу прокомментировать общую переносимость кода между разными компиляторами.