Fortran Coder

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

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

[复制链接]

47

帖子

15

主题

0

精华

专家

F 币
311 元
贡献
158 点
跳转到指定楼层
楼主
发表于 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

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

47

帖子

15

主题

0

精华

专家

F 币
311 元
贡献
158 点
沙发
 楼主| 发表于 2022-5-27 12:15:28 | 显示全部楼层
fcode 发表于 2022-5-27 08:43
get_r 是函数,不能用 call 调用。必须 r = get_r(k)

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

47

帖子

15

主题

0

精华

专家

F 币
311 元
贡献
158 点
板凳
 楼主| 发表于 2022-5-27 12:21:22 | 显示全部楼层
sqs 发表于 2022-5-27 12:15
好的感谢!,另外我还有个问题就是我在主程序当中的子程序subroutine去调用module当中的函数或者程序的时 ...

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

47

帖子

15

主题

0

精华

专家

F 币
311 元
贡献
158 点
地板
 楼主| 发表于 2022-5-27 12:49:49 | 显示全部楼层
fcode 发表于 2022-5-27 12:23
是的。
现代Fortran代码,应该只有 Module 和 主程序。没有“外部子程序和外部函数”,即,把所有函数都放 ...

好的,非常感谢!

47

帖子

15

主题

0

精华

专家

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

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

47

帖子

15

主题

0

精华

专家

F 币
311 元
贡献
158 点
6#
 楼主| 发表于 2022-5-27 15:13:39 | 显示全部楼层
fcode 发表于 2022-5-27 14:42
没有什么区别。程序单元和程序单元之间的命名空间是独立的。
程序单元A(比如主程序)中的变量 x 和程序单 ...

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

47

帖子

15

主题

0

精华

专家

F 币
311 元
贡献
158 点
7#
 楼主| 发表于 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)
这样算是实参和虚参名字一样且对应么

47

帖子

15

主题

0

精华

专家

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

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

47

帖子

15

主题

0

精华

专家

F 币
311 元
贡献
158 点
9#
 楼主| 发表于 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 是实参。

我明白啦,感谢感谢!
如果命名空间是独立的,我就可以反复运用同样一个符号去计算,我感觉是很方便的哈哈哈
十分感谢!
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

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

GMT+8, 2024-5-2 19:02

Powered by Tencent X3.4

© 2013-2024 Tencent

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