在Fortran中编写复杂矩阵

Question

如何在Fortran中将复杂（nXn）矩阵写入文件？ 例如：

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

在这个例子中一个送为2N×N写入到文件，即实部和虚元素。 （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'）中被考虑，所以我们取虚部的绝对值。相应地编辑格式描述符，以便用于读取输出文件的程序将能够读取它。