Fortran Coder

标题: 编译链接问题 [打印本页]

作者: leeyazhou 时间: 2017-11-1 08:14
标题: 编译链接问题

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

        !MAIN PROGRAM

        PROGRAM 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. [A]”
链接出错“hw1.obj : error LNK2001: unresolved external symbol _CDABS@4”
*********************************************************************
而换这样写就没问题

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

      !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:50
本帖最后由 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 ，直接写

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode

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 的源代码。

祝好！

作者: leeyazhou 时间: 2017-11-1 10:27

vvt 发表于 2017-11-1 08:50
第一个问题
complex*8 a , b
等效于

谢谢！理解。
这里的代码对编译器版本（或者fortran版本）有要求吧

作者: vvt 时间: 2017-11-1 11:21
本帖最后由 vvt 于 2017-11-1 11:26 编辑

支持90的编译器就可以

作者: pasuka 时间: 2017-11-1 11:22
首先，上周lz的帖子就有热心版友li913给出本站的函数重载说明
http://fcode.cn/guide-61-1.html
其次，翻翻gfortran的在线手册，譬如sqrt函数
https://gcc.gnu.org/onlinedocs/gcc-7.2.0/gfortran/SQRT.html#SQRT

作者: leeyazhou 时间: 2017-11-1 15:30
谢谢

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)