Fortran Coder

查看: 4836|回复: 16
打印 上一主题 下一主题

[子程序] 关于调用函数时的报错

[复制链接]

49

帖子

17

主题

0

精华

熟手

F 币
270 元
贡献
164 点
跳转到指定楼层
楼主
发表于 2022-5-26 23:08:43 | 只看该作者 回帖奖励 |正序浏览 |阅读模式
在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

请问这种错误类型是什么意思呢,另外应该如何改正。
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

49

帖子

17

主题

0

精华

熟手

F 币
270 元
贡献
164 点
17#
 楼主| 发表于 2022-5-27 16:21:38 | 只看该作者
fcode 发表于 2022-5-27 16:09
subroutine efg(a,b,c)  这里的 a b c 是虚参。
call efg(a,b,c) 这里的 a b c 是实参。

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

2022

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1598 元
贡献
689 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

16#
发表于 2022-5-27 16:09:46 | 只看该作者
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 贷款(郭靖,黄蓉)!第二次贷款


49

帖子

17

主题

0

精华

熟手

F 币
270 元
贡献
164 点
15#
 楼主| 发表于 2022-5-27 15:44:20 | 只看该作者
fcode 发表于 2022-5-27 15:17
我没太听明白你的意思,反正主程序里的 x 变量,不一定是子程序里的 x 变量。除非实参和虚参名字一样且对应 ...

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

49

帖子

17

主题

0

精华

熟手

F 币
270 元
贡献
164 点
14#
 楼主| 发表于 2022-5-27 15:30:41 | 只看该作者
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)
这样算是实参和虚参名字一样且对应么

2022

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1598 元
贡献
689 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

13#
发表于 2022-5-27 15:17:20 | 只看该作者
我没太听明白你的意思,反正主程序里的 x 变量,不一定是子程序里的 x 变量。除非实参和虚参名字一样且对应。

49

帖子

17

主题

0

精华

熟手

F 币
270 元
贡献
164 点
12#
 楼主| 发表于 2022-5-27 15:13:39 | 只看该作者
fcode 发表于 2022-5-27 14:42
没有什么区别。程序单元和程序单元之间的命名空间是独立的。
程序单元A(比如主程序)中的变量 x 和程序单 ...

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

2022

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1598 元
贡献
689 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

11#
发表于 2022-5-27 14:42:56 | 只看该作者
没有什么区别。程序单元和程序单元之间的命名空间是独立的。
程序单元A(比如主程序)中的变量 x 和程序单元B(比如子程序)中的同名变量 x,是独立的。互不影响。

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

一个鲜活的,生动的,经典的例子:

subroutine 结婚(男方,女方)

call 结婚(张三,李四)
call 结婚(郭靖,黄蓉)

49

帖子

17

主题

0

精华

熟手

F 币
270 元
贡献
164 点
10#
 楼主| 发表于 2022-5-27 13:53:42 | 只看该作者
fcode 发表于 2022-5-27 12:23
是的。
现代Fortran代码,应该只有 Module 和 主程序。没有“外部子程序和外部函数”,即,把所有函数都放 ...

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

49

帖子

17

主题

0

精华

熟手

F 币
270 元
贡献
164 点
9#
 楼主| 发表于 2022-5-27 12:49:49 | 只看该作者
fcode 发表于 2022-5-27 12:23
是的。
现代Fortran代码,应该只有 Module 和 主程序。没有“外部子程序和外部函数”,即,把所有函数都放 ...

好的,非常感谢!
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-11-23 16:55

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表