написать сложную матрицу в fortran

Question

Как написать комплекс (nXn) матрица в Fortran для файла? Например:

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

В этом примере один прибудет 2nXn элементы, записываемые в файл то есть реальный и воображаемый. Вместо двух элементов Re (a11) Im (a11), Как я могу записать его как один элемент Re (a11) + iIm (a11)?

Answer 1

Используйте встроенные функции РЕАЛЬНЫЕ и 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'), мы принимаем абсолютное значение мнимой части. Отредактируйте дескриптор формата соответственно, чтобы программа, которую вы используете для чтения выходного файла, сможет ее прочитать.