Matrix 如何在Fortran中读/写矩阵?
下面我有一个简单的Fortran代码,我从一个文本文件中读取一个矩阵,然后将它传递给一个子例程。但是当在子例程中读取它时,它会给出一个垃圾值。我似乎不知道代码出了什么问题Matrix 如何在Fortran中读/写矩阵?,matrix,fortran,Matrix,Fortran,下面我有一个简单的Fortran代码,我从一个文本文件中读取一个矩阵,然后将它传递给一个子例程。但是当在子例程中读取它时,它会给出一个垃圾值。我似乎不知道代码出了什么问题 program main implicit none real, dimension(:,:), allocatable :: mat open (unit=99, file='inp.txt', status='old', action='read') allocate(mat(2,2)) read(99, *), mat
program main
implicit none
real, dimension(:,:), allocatable :: mat
open (unit=99, file='inp.txt', status='old', action='read')
allocate(mat(2,2))
read(99, *), mat
CALL fun(mat)
end
subroutine fun(mat)
dimension mat(2,2)
write(*,*) mat
return
end
矩阵
mat
作为[2 4;2 10]
读入,但当在子例程中显示时,它给出[1073741824 1082130432;1073741824 1092616192]
如果在Fortran中使用隐式变量声明,变量以I
,J
,K
,L
,M
,或N
表示整数(默认情况下)
由于主程序的作用域中只有implicit none
,而子例程没有,因此real
arraymat
在子例程中被解释为整数。因此,垃圾
只需在子例程中将
mat
声明为real
,即可解决问题 您的mat有一个隐式类型的整数,但您正在传递它reals。。。更好地使用模块和隐式无无处不在…谢谢。在子例程中包括隐式无
,并将mat声明为真实
工作!