用 Fortran 将复数矩阵写入文件

Question

如何将 Fortran 中的复数 (n×n) 矩阵写入文件？ 例如：

DO I=1,N
       write(14,'(100g15.5)') ( M(i,j), j=1,n )
ENDDO  

在此示例中，将 2n×n 个元素写入文件，即实数和虚数。 除了两个元素 Re(a11) Im(a11)，我如何将其写为一个元素 Re(a11)+iIm(a11)？

Answer 1

使用内部函数 REAL 和 AIMAG 写入复数的各个实部和虚部：

CHARACTER(LEN=3),DIMENSION(n,n) :: imag_unit = '+i*'

WHERE(AIMAG(M)<0.)imag_unit = '-i*'

DO I=1,N
  write(14,'(100(g15.5,a,g15.5,2x))') ( REAL(M(i,j)),imag_unit(i,j),&
                                        ABS(AIMAG(M(i,j))), j=1,n )
ENDDO 

说明：此代码定义一个字符串矩阵，当虚部为正时，其值为“+i”，当虚部为正时，值为“-i”消极的。因为负虚部在格式（'-i'）中被考虑，所以我们取虚部的绝对值。相应地编辑格式描述符，以便用于读取输出文件的程序能够读取它。