我在Fortran的常数数字精度方面遇到了麻烦.
我是否需要将每一个都0.1写成0.1d0双精度？我知道编译器有一个标志,例如-fdefault-real-8在gfortran中解决了这种问题.它是一种便携可靠的方式吗？我怎么能检查标志选项是否真的适用于我的代码？
我使用F2py在我的Python代码中调用Fortran代码,即使我给出了一个未指定的标志,它也不会报告错误,这就是让我担心的问题.
我在F90中有类似的代码:
real(8), dimension(10,10,10) :: A
do i = 1, 1000
print*,A(i,1,1)
enddo
我很惊讶这个工作,它比简单地循环3个维度更快i,j,k.
有人可以解释为什么这有效吗？
fortran ×2