编译链接问题

leeyazhou

[Fortran] 纯文本查看 复制代码

        !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] 纯文本查看 复制代码

      !MAIN PROGRAM
        PROGRAM HOMEWORK0
         complex*8 a,b
        a=cdabs(cdsqrt((-4.d0,0.d0)))
        write(*,*)a        
                stop
        END PROGRAM HOMEWORK0

********************************************************************
请问这是什么问题？能介绍一下fortran是如何编译链接的。

vvt

第一个问题
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] 纯文本查看 复制代码

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

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

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

vvt

支持90的编译器就可以

pasuka

首先，上周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

谢谢