我有一个最初针对英特尔编译器的大型 Fortran 代码库。我现在正准备用 gfortran 进行编译。不幸的是，代码中充斥着 Intel 风格的预处理指令，例如：
!DEC$ IF DEFINED (MYDIRECTIVE)
REAL, DIMENSION(:,:,:), ALLOCATABLE :: my_real_var
!DEC$ ENDIF
从我通过谷歌搜索和 gfortran 文档可以看出，除了 C 风格的预处理之外，没有内部 gfortran 支持，例如：
#if defined MYDIRECTIVE
REAL, DIMENSION(:,:,:), ALLOCATABLE :: my_real_var
#endif
有没有其他人遇到过这个问题并提出了一个优雅的解决方案？显然，我可以编写一个 shell 脚本，在将代码传递给 gfortran 进行编译之前调用外部预处理器，但这对我来说似乎不是一个很好的解决方案。
有什么想法吗？谢谢大师！
有没有办法捕获整数异常，gfortran或者ifort像捕获浮点异常一样？
考虑这个简单的程序来计算阶乘：
program factorial
use, intrinsic :: iso_fortran_env
implicit none
integer(8) :: fac
real(REAL64) :: facR
integer,parameter :: maxOrder = 30
integer :: i
fac = 1 ; facR = 1.e0_REAL64
do i=2,maxOrder
fac=fac*i ; facR=facR*real(i,REAL64)
write(*,*) i, fac, facR
enddo ! i
end program
在某些时候会出现溢出 -integer(8)如此处所示，它将发生在 21 左右。但是如果没有使用浮点数作为参考的计算，我无法确定......
我试图通过 Intel Fortran 编译器在 Windows 8、x64 上使用 f2py 来实现。
我尝试过两件事：
1) 通过 Windows 64x 安装程序安装。http://www.lfd.uci.edu/~gohlke/pythonlibs/ NumPy 已正确安装，但执行 f2py -c --help-fcompiler 时未找到编译器 我得到：
Fortran compilers found:
Compilers available for this platform, but not found:
--fcompiler=absoft Absoft Corp Fortran Compiler
--fcompiler=compaqv DIGITAL or Compaq Visual Fortran Compil
--fcompiler=g95 G95 Fortran Compiler
--fcompiler=gnu GNU Fortran 77 compiler
--fcompiler=gnu95 GNU Fortran 95 compiler
--fcompiler=intelem Intel Fortran Compiler for 64-bit apps
--fcompiler=intelev Intel Visual Fortran Compiler for Itani
--fcompiler=intelv Intel Visual Fortran Compiler for …我经常看到人们使用该OPEN语句而不明确指定STATUS. 在 Fortran 90 和 2008 标准中，这样说STATUS：
如果指定 UNKNOWN，则状态取决于处理器。如果省略此说明符，则默认值为 UNKNOWN。
我将此解释为，如果STATUS省略，任何事情都可能发生，具体取决于您使用的机器。
然而，从进行一些测试来看，默认行为（当STATUS省略时）似乎是REPLACE. 但我找不到 gfortran 编译器手册（来自https://gcc.gnu.org/onlinedocs/）中记录的这种行为。
问题：这REPLACE确实是 gfortran 和 ifort 等流行编译器的默认行为吗？如果是这样，这实际上有记录吗（但我只是碰巧没有找到它）？
当我最初无法预测数组的确切大小时，我需要在 Fortran 90 中使用动态数组。所以我编写了一段代码，每次将新元素添加到数组末尾时，该代码都应该扩展可分配数组：
subroutine DArray()
double precision, dimension(:), allocatable :: list
allocate(list(1))
list(1) = 1.1
call AddToList(list, 2.2)
call AddToList(list, 3.2)
call AddToList(list, 4.2)
call AddToList(list, 5.2)
print *, list(1)
print *, list(2)
print *, list(3)
print *, list(4)
print *, list(5)
end
subroutine AddToList(list, element)
double precision :: element
double precision, dimension(:), allocatable :: list
double precision, dimension(:), allocatable :: clist
if(allocated(list)) then
isize = size(list)
allocate(clist(isize+1))
do i=1,isize
clist(i) = list(i)
end do
clist(i+1) = element
deallocate(list) …我正在尝试将一个结构从 Fortran 传递到 C，其中 Fortran 中的结构有一个可分配的。我想在 fortran 中的结构内分配数组并在 C 中读取它。但是，当我尝试在 CI 中打印分配的数组时，收到错误消息分段错误（核心转储）。我使用的是英特尔编译器版本 19.0.1.144。
我对 C 比较陌生，所以为了测试这个想法，我在 fortran 中分配了一个一维数组并将其传递给 C。这非常有效。您可以在 fortran 代码中看到 ff 是可分配数组，并且为它分配了维度 ii（等于 5）。我初始化数组ff并使用类型(C_ptr)的指针来存储ff数组的地址。我在 C 中打印 ff ，它打印出正确的值 ff ，即 14.7,15.7,16.7,17.7,18.7。使用相同的想法，我声明了一个带有整数 ival、实数 cval 和可分配变量 dval 的结构体 data1。我正在为结构类型 data1 分配一个指针 pdata。我分配了指针和数组 dval。我在 C Main 中调用 fortran 来打印结构数量，当我打印数组 d （或 dval）时出现错误。另外，C中的ival和cval的值也不正确。
SUBROUTINE Simulation(ii,ffp,cdata) BIND(C)
use, intrinsic :: iso_c_binding
integer(kind=4), intent(in) :: ii
type(C_PTR), intent(out) :: ffp
real (C_double), dimension(:), allocatable, target, save :: ff
integer(kind=4) :: err, i
Type :: data1 …当 gfortran 对方程（例如 x = y*z ）进行向量化时，编译器将仅使用向量寄存器（例如 YMM）和向量操作码，当且仅当编译器知道所有数据 (x,y,z) 都可以是加载到所述向量寄存器上。然而，如果编译器不知道是否所有数据都可以加载到寄存器中，它将生成冗余代码，以便可以操作剩余的数据。
此冗余代码通常不会使用最高可用的 SIMD 向量寄存器，并且会使可执行文件膨胀。我还没有进行广泛的测试，但常见的推理使我相信它可能会导致 ICACHE 缺失或阻止函数/子例程内联优化。
我想删除这个冗余代码，因为我知道我的所有数据 (x,y,z) 将完美地适合 YMM 向量寄存器。
我之前曾在这里更详细地描述过我的观察结果：https: //www.cfd-online.com/Forums/main/231759-fortran- assembly-how-remove-redundant-non-vectorized-code.html
不过，我想在这里举一些简单的例子来展示我的观点。
对于以下操作：x(:) = y*z其中x，y，z都是可分配数组，编译器不知道该操作是否可以完全矢量化，因为它不知道数组的长度是多少。如果 x,y,z 是整数数组并且长度是 8 的倍数，我们可以推断整个操作可以简单地使用 AVX2 向量化操作来完成。然而，由于数组长度在编译时未知，编译器必须生成可用于任意长度数组的冗余代码。
上帝螺栓链接：https://gcc.godbolt.org/z/1fdzK8
这是为该x(:) = y*z部分生成的程序集：
.L7:
vmovdqu ymm1, YMMWORD PTR [r13+0+rax]
vpmulld ymm0, ymm1, YMMWORD PTR [r12+rax]
vmovdqu YMMWORD PTR [r14+rax], ymm0
add rax, 32
cmp rdx, rax
jne .L7
mov rdx, r15
and rdx, -8
lea rax, [rdx+1]
cmp rdx, …有几个关于 SO 的精彩讨论已经涵盖了如何在 Linux 上生成可执行共享库：
在 C/C++ 中，这似乎相对简单；基本上有两个部分：
.interp内容，向 ELF添加一个部分（因为ld不包括共享库的部分）：const char interp_section[] __attribute__((section(".interp"))) = "/path/to/dynamic/linker";-Wl,-e,entry_point有谁知道如何使用 Fortran 编写的库来实现这一目标？具体来说，如何将一个.interp部分添加到使用编译的共享库中ifort？
我正在处理的 Fortran 代码有几行类似于
WRITE(filename, '(A16,"_",I4,".dat")') filename, indx
该代码已在许多不同的平台和几乎所有主要编译器上成功编译并运行了数百次。但突然间，最新的（或者，无论如何，新的）英特尔编译器不喜欢它了。它给出警告信息"forrtl: .... Internal file write-to-self; undefined results"。该行执行后，"filename"原本合理的字符数组变为空白。
我想问题在于它既filename是写入的输入又是内部写入的目的地。修复很容易。它可以将filename目的地替换为类似的东西filename_tmp。但正如我所说，到目前为止，这从来都没有必要。
所以我想知道，文件名作为输入和目标是否违反了 Fortran 标准，但这些年来所有这些编译器都对此视而不见，而现在英特尔变得越来越严格？还是英特尔太“势利”了？还是彻头彻尾的越野车？
以下代码可在 GNU gfortran 和 Intel ifort 中编译。但只有gfortran编译版本才能成功运行。
program fort_tst
use iso_c_binding
INTEGER, POINTER :: a(:)
TYPE(C_PTR) :: ptr
INTEGER, POINTER :: b(:)
ALLOCATE(a(5))
ptr = c_loc(a)
CALL c_f_pointer(ptr,b,[5])
DEALLOCATE(b)
end program fort_tst
英特尔编译代码中的错误是：
forrtl: severe (173): A pointer passed to DEALLOCATE points to an object that cannot be deallocated
Image PC Routine Line Source
fort_tst 000000000040C5A1 Unknown Unknown Unknown
fort_tst 0000000000403A17 Unknown Unknown Unknown
fort_tst 0000000000403812 Unknown Unknown Unknown
libc-2.17.so 00002AAAAB20F555 __libc_start_main Unknown Unknown
fort_tst 0000000000403729 Unknown Unknown Unknown
gfortran …
fortran pointers gfortran fortran-iso-c-binding intel-fortran