Fortran Do 循环的上限

Question

在以下代码中，DO 循环的上限在循环内修改为

    integer :: i1, i2, n1, n2
    n1 = 4 ; n2 = 1

    do i1 = 1, n1
    do i2 = 1, n2
        print *, "i1 =", i1, "i2 =", i2
        n1 = 2
        n2 = 2
    enddo
    enddo

其中 gfortran4.8 和 ifort14.0 给出以下结果：

i1 =           1 i2 =           1
i1 =           2 i2 =           1
i1 =           2 i2 =           2
i1 =           3 i2 =           1
i1 =           3 i2 =           2
i1 =           4 i2 =           1
i1 =           4 i2 =           2

这表明在进入每个 DO 循环时界限是固定的（即，尽管 n1 在循环内修改为 2，但 i1 的上限固定为 4）。此行为与 C/C++ 的行为相反，其中相应的代码

int n1 = 4, n2 = 1;

for( int i1 = 1; i1 <= n1; i1++ )
for( int i2 = 1; i2 <= n2; i2++ )
{
    printf( "i1 = %d i2 = %d\n", i1, i2 );
    n1 = 2;
    n2 = 2;
}

给出

i1 = 1 i2 = 1
i1 = 1 i2 = 2
i1 = 2 i2 = 1
i1 = 2 i2 = 2

这是合理的，因为 for 循环的定义（即开始、结束和增量条件）。所以我的问题是：是否可以假设 DO 循环的上述行为是通用的（= 由 Fortran 标准定义）或编译器相关的（即 ifort 或 gfortran 的不同编译器可能表现不同）？

Answer 1

Fortran 标准指定循环的迭代次数DO在循环开始时固定，使用公式（至少直到 F95，此时删除了使用非整数循环控制表达式的能力）

iterations = MAX(0,INT((final-value - initial-value + step-size)/step-size)

这样做的结果之一就是您所看到的行为：更改循环体内用于控制循环的任何值都不会更改将执行的迭代次数。