使用隐式循环打印 Fortran 数组

Question

该程序未成功地尝试使用格式规范中的隐式循环在一行上打印一个数组。它通过显式循环成功地完成了这项工作。

program cycle
    implicit none
    integer, dimension(5) :: a=(/1,2,3,4,5/), b=(/11,12,13,14,15/)
    integer :: n, i

    Print *, "Implicit loop"
    print "(i0, 1x)", (a(i)*b(i), i=1,n)
    Print *, "Explicit loop"
    do i=1,n-1
        write(*, '(i0, 1x)', advance = "no") a(i)*b(i)
    end do
    write(*, '(i0)') a(n)*b(n)
end program cycle

结果如下：

 Implicit loop
11
24
39
 Explicit loop
11 24 39

如何使隐式循环在一行上打印所有内容？甚至可能吗？我尝试插入advance="no"，并用各种逗号和括号括起来，但没有奏效。

Answer 1

你的格式字符串只支持两个项目，一个整数和一个空格。之后开始一条新记录并从头开始解释格式字符串。

(i0, 1x)

您必须将其转换为多个项目

(999(i0, 1x))

或（Fortran 2008）

(*(i0, 1x))

然后括号的使用次数与前面的数字一样多。 * 表示无限期。当然，只要 i/o 列表中还有项目需要处理。

Answer 2

目前，我无法发布 cmets，因为我的声誉较低，因此将其写为答案。弗拉基米尔的回答是完美的。但要添加到该答案，您还可以执行以下操作：

         integer, parameter :: n=5    ! this line was missing in the question
         write(fmt, 5)n , "(i0, 1x)"
  5      format('(', I2, A9, ')')
         write(*,*)"Implicit loop, formatted"
         write(*, fmt) (a(i)*b(i), i=1,n)

对于显式循环，您可以在 Fortran 2003 或更高版本中使用 new_line 执行以下操作：

    Print *, "Explicit loop 2"
    do i=1,n  ! notice the loop ends at 'n'
        write(*, '(i0, 1x)', advance = "no") a(i)*b(i)
    end do
    write(*, *) new_line('')

这适用于 intel 和 gfortran。

Answer 3

首先，非高级 I/O 只能在外部文件上执行，不能用于名称列表或列表导向 I/O。

其次，可以将重复计数应用于一组编辑描述符， 括在括号中：

print '(4(i5,f10.3))', (bar(j), foo(j), j=1,4)

（对于整数 bar 和实数 foo）。这相当于写

print '(i5,f10.3,i5,f10.3,i5,f10.3,i5,f10.3)', (bar(j), foo(j), j=1,4)

这样的重复计数可能是嵌套的：

print '(2(2i5,2f10.3))', bar(1),bar(2),foo(1),foo(2),bar(3),bar(4),foo(3),foo(4)

第三，可以通过三种方式给出format 规范。它们如下：

(1) 一个默认的字符表达式（我有偏见的最佳方法）：

print '(f10.3)', foo
write(*, '(f10.3)') foo

或

character(len=*), parameter :: FMT='(f10.3)'

print FMT, foo
write(*, FMT) foo

或

character, parameter :: FMT_ARRAY(*)=[ '(','f','1','0','.','3',')' ]

! Elements of an array expression are concatenated.
print FMT_ARRAY, foo
write(*, FMT_ARRAY) foo 

or

character(4) :: fmt_1(10)
character(3) :: fmt_2(10)
integer      :: i, j

fmt_1(10) = '(f10'
fmt_2(3) = '.3)'

i = 10
j = 3

! Format is built dynamically at execution time from various components.
print fmt_11(i)//fmt_2(j), foo 
write(*, fmt_11(i)//fmt_2(j)) foo 

(2) 星号（仅在代码的测试/开发阶段使用）：

print *, foo
write(*,*) foo

(3) 声明标签（今天不是撒旦！2016 年没有这个位置）：

    print 42, foo
    write(*, 42) foo

42 format(f10.3)