编译链接问题
!MAIN PROGRAMPROGRAM HOMEWORK0
complex*8 a,b
a=cdsqrt((-4.d0,0.d0))
write(*,*)a
b = cdabs(a)
write(*,*)b
stop
END PROGRAM HOMEWORK0
此程序编译有警告“D:\hw1.f(9) : Warning: This argument's data type is incompatible with this intrinsic procedure; procedure assumed EXTERNAL. ”
链接出错“hw1.obj : error LNK2001: unresolved external symbol _CDABS@4”
*********************************************************************
而换这样写就没问题
!MAIN PROGRAM
PROGRAM HOMEWORK0
complex*8 a,b
a=cdabs(cdsqrt((-4.d0,0.d0)))
write(*,*)a
stop
END PROGRAM HOMEWORK0
********************************************************************
请问这是什么问题？能介绍一下fortran是如何编译链接的。
本帖最后由 vvt 于 2017-11-1 08:51 编辑
第一个问题
complex*8 a , b
等效于
complex(4) a , b
在大多数编译器上，等效于 complex a , b ，是单精度的。
您要使用双精度，请用
complex(8) a , b
或 complex(kind=8) a , b !推荐这种
或 complex*16 a , b !不推荐这种
第二个问题
b = cdabs(a)要求 a 是双精度的，而因为你的第一个问题，导致 a 是单精度，所以会出错。
第三个问题
其实没必要写 cdabs ，直接写
PROGRAM HOMEWORK0
Implicit None
complex(kind=8) :: a , b
a=sqrt((-4.d0,0.d0))
write(*,*)a
b = abs(a)
write(*,*)b
END PROGRAM HOMEWORK0
就可以了，编译器会根据你的参数类型，自动决定用单精度还是双精度，real还是complex
详解 http://whyd.w.fcode.cn
最后一个问题fortran是如何编译链接的。
如果您是想知道如何操作编译器，令其编译链接。那么不同编译器的操作方式不一样。
如果您是想了解编译器内部是如何实现编译链接的。您可以看《编译原理》，还有 gfortran 的源代码。
祝好！
vvt 发表于 2017-11-1 08:50
第一个问题
complex*8 a , b
等效于
谢谢！理解。
这里的代码对编译器版本（或者fortran版本）有要求吧 本帖最后由 vvt 于 2017-11-1 11:26 编辑
支持90的编译器就可以
首先，上周lz的帖子就有热心版友li913给出本站的函数重载说明
http://fcode.cn/guide-61-1.html
其次，翻翻gfortran的在线手册，譬如sqrt函数
https://gcc.gnu.org/onlinedocs/gcc-7.2.0/gfortran/SQRT.html#SQRT
谢谢:-handshake
页:
[1]