Fortran Coder

标题: 关于调用函数时的报错 [打印本页]

作者: sqs 时间: 2022-5-26 23:08
标题: 关于调用函数时的报错
在module中我分别有两个函数，其中一个函数在调用另外一个的时候，报错，类型如下;

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

 46 |   function get_r(k) result(r)

      |  1

......

   78 |         call get_r(k)

      |                     2

Error: ‘get_r’ at (1) has a type, which is not consistent with the CALL at (2)

其中对应的函数程序为：

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

function get_r(k) result(r)

      implicit none

      real(8) , intent(in)  :: k(2)

      complex(8) :: r(2, 6, 6) 

      ! r is a three dimension array

      r = 0.d0

      r(1, 1, 1) = -a0*r0

      r(1, 4, 4) = a0*r0

      r(1, 5, 5) = -a0*r0

        

  end function get_r

    

    

  function get_v(k) result(v)

        implicit none

        real(8), intent(in)  :: k(2)

        complex(8) :: v(2, 6, 6) 

        complex(8) :: f1(2)

        

        complex(8) :: H(6, 6), r(2, 6, 6)

        integer :: i

        

        v = 0.d0

      

        f1(:) = CI*(rprimd(:, 1)*exp(2.d0*Pi*CI*k(1))+rprimd(:, 2)*exp(2.d0*Pi*CI*k(2)))



        v(:, 1, 2) = t1*f1;        v(:, 2, 1) = dconjg(v(:, 1, 2))

        v(:, 3, 4) = t1*f1;        v(:, 4, 3) = dconjg(v(:, 3, 4))

        v(:, 5, 6) = t1*f1;        v(:, 6, 5) = dconjg(v(:, 5, 6))

        

        call get_H(k)

        

        call get_r(k)

      

        do i = 1, 2, 1

          v(i, :, :) = (v(i, :, :) - CI*(Matmul(r(i, :, :), H) - MatMul(H, r(i, :, :))))/h_bar

        end do

        

  end function get_v

请问这种错误类型是什么意思呢，另外应该如何改正。

作者: necrohan 时间: 2022-5-27 08:41
get_r调用自己？核实一下46行是不是get_r

作者: fcode 时间: 2022-5-27 08:43
get_r 是函数，不能用 call 调用。必须 r = get_r(k)

作者: zjk0112 时间: 2022-5-27 09:43

fcode 发表于 2022-5-27 08:43
get_r 是函数，不能用 call 调用。必须 r = get_r(k)

我想顺便问下，如果是subroutine可以这样调用吗？

作者: fcode 时间: 2022-5-27 09:46
可以

作者: sqs 时间: 2022-5-27 12:15

fcode 发表于 2022-5-27 08:43
get_r 是函数，不能用 call 调用。必须 r = get_r(k)

好的感谢！,另外我还有个问题就是我在主程序当中的子程序subroutine去调用module当中的函数或者程序的时候，出现了函数没有implicit type ,请问这是什么错误呢

作者: sqs 时间: 2022-5-27 12:21

sqs 发表于 2022-5-27 12:15
好的感谢！,另外我还有个问题就是我在主程序当中的子程序subroutine去调用module当中的函数或者程序的时 ...

也就是在module外的子程序调用module 内的子程序时也需要在implicit none前面加use module name 么

作者: fcode 时间: 2022-5-27 12:23
是的。
现代Fortran代码，应该只有 Module 和主程序。没有“外部子程序和外部函数”，即，把所有函数都放入对应的 module 中。
这样，只需在 module A 中 use 其他module B，则module A下的子程序和函数均可调用module B 中的（public的）函数。

作者: sqs 时间: 2022-5-27 12:49

fcode 发表于 2022-5-27 12:23
是的。
现代Fortran代码，应该只有 Module 和主程序。没有“外部子程序和外部函数”，即，把所有函数都放 ...

好的，非常感谢！

作者: sqs 时间: 2022-5-27 13:53

fcode 发表于 2022-5-27 12:23
是的。
现代Fortran代码，应该只有 Module 和主程序。没有“外部子程序和外部函数”，即，把所有函数都放 ...

另外关于子程序参数我还有一点点问题，就是加入子程序当中我定义的输出参量是A的话，在调用时，我将输出参量同样定义为A然后再输出，以及解耦将输出参量定义为其他名称如B，这两种有什么区别呢。以及第一种为什么不会出现问题呢，明明我对一个变量定义了两次，但是还是可以正常地输出

作者: fcode 时间: 2022-5-27 14:42
没有什么区别。程序单元和程序单元之间的命名空间是独立的。
程序单元A（比如主程序）中的变量 x 和程序单元B（比如子程序）中的同名变量 x，是独立的。互不影响。

就像我家的电视机和你家的电视机，都叫电视机，但不是同一个。
注意，虚参归属于子程序的程序单元，而实参归属于调用者的程序单元。

一个鲜活的，生动的，经典的例子：

subroutine 结婚（男方，女方）

call 结婚（张三，李四）
call 结婚（郭靖，黄蓉）

作者: sqs 时间: 2022-5-27 15:13

fcode 发表于 2022-5-27 14:42
没有什么区别。程序单元和程序单元之间的命名空间是独立的。
程序单元A（比如主程序）中的变量 x 和程序单 ...

也就是说，我在主程序当中调用的时候，虽然输出变量定义的名称还是一样的，但是对应到的内存空间其实是不一样了是吧（还是解耦的）

作者: fcode 时间: 2022-5-27 15:17
我没太听明白你的意思，反正主程序里的 x 变量，不一定是子程序里的 x 变量。除非实参和虚参名字一样且对应。

作者: sqs 时间: 2022-5-27 15:30

fcode 发表于 2022-5-27 15:17
我没太听明白你的意思，反正主程序里的 x 变量，不一定是子程序里的 x 变量。除非实参和虚参名字一样且对应 ...

比如说,
module abc内有一subroutine efg(a , b , c ),其中intent(in) 是a, b.intent（out）是 c
在主程序当中调用时能否如下：
call efg(a,b,c)
这样算是实参和虚参名字一样且对应么

作者: sqs 时间: 2022-5-27 15:44

fcode 发表于 2022-5-27 15:17
我没太听明白你的意思，反正主程序里的 x 变量，不一定是子程序里的 x 变量。除非实参和虚参名字一样且对应 ...

或者按照您的说法，就是
subroutine 结婚（郭靖，黄蓉）
调用的时候我能用
call 结婚（郭靖，黄蓉）
当然对应的变量在子程序和主程序当中得分别定义一遍，这样会有问题么

作者: fcode 时间: 2022-5-27 16:09

sqs 发表于 2022-5-27 15:30
比如说,
module abc内有一subroutine efg(a , b , c ),其中intent(in) 是a, b.intent（out）是 c
在主程 ...

subroutine efg(a,b,c) 这里的 a b c 是虚参。
call efg(a,b,c) 这里的 a b c 是实参。

实参和虚参名字可以一样，也可以不一样。
我说的对应，是位置对应。

在我的例子中。男方，女方，是虚参。郭靖黄蓉是实参。
实参和虚参需要对应（否则郭靖成了女方，黄蓉成了男方岂不是笑话了？）
但是名字可以随意，一样，或不一样，都可以。
有一些极端的情况下，名字一样但位置不对应。
subroutine efg(a,b,c)
call efg(c,b,a)
此时，实参的a对应虚参的c，实参的c对应虚参的a。

记住，实参属于调用者，虚参属于被调用者。他们的命名空间是独立的（这样非常有好处，尤其是程序变得很庞大之后）

实参必须在调用者中定义，虚参必须在被调用者中定义。（不管他们变量名字是否相同）编译器会帮你检查实参虚参的类型是否一致。（毕竟，男方女方需要是一个人，而实参不能是条狗）

为了让你理解。我把之前的例子扩展一下：
subroutine 结婚（男方，女方）
subroutine 贷款（主贷人，连带贷款人）

call 结婚（郭靖，黄蓉）
call 贷款（黄蓉，郭靖）!第一次贷款
call 贷款（郭靖，黄蓉）!第二次贷款

作者: sqs 时间: 2022-5-27 16:21

fcode 发表于 2022-5-27 16:09
subroutine efg(a,b,c) 这里的 a b c 是虚参。
call efg(a,b,c) 这里的 a b c 是实参。

我明白啦，感谢感谢！
如果命名空间是独立的，我就可以反复运用同样一个符号去计算，我感觉是很方便的哈哈哈
十分感谢！

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