关于调用函数时的报错
在module中我分别有两个函数,其中一个函数在调用另外一个的时候,报错,类型如下;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)
其中对应的函数程序为:
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
请问这种错误类型是什么意思呢,另外应该如何改正。
get_r调用自己?核实一下46行是不是get_r get_r 是函数,不能用 call 调用。必须 r = get_r(k) fcode 发表于 2022-5-27 08:43
get_r 是函数,不能用 call 调用。必须 r = get_r(k)
我想顺便问下,如果是subroutine可以这样调用吗? 可以 fcode 发表于 2022-5-27 08:43
get_r 是函数,不能用 call 调用。必须 r = get_r(k)
好的感谢!,另外我还有个问题就是我在主程序当中的子程序subroutine去调用module当中的函数或者程序的时候,出现了函数没有implicit type ,请问这是什么错误呢 sqs 发表于 2022-5-27 12:15
好的感谢!,另外我还有个问题就是我在主程序当中的子程序subroutine去调用module当中的函数或者程序的时 ...
也就是在module外的子程序调用module 内的子程序时也需要在implicit none前面加use module name 么
是的。
现代Fortran代码,应该只有 Module 和 主程序。没有“外部子程序和外部函数”,即,把所有函数都放入对应的 module 中。
这样,只需在 module A 中 use 其他module B,则module A下的子程序和函数均可调用module B 中的(public的)函数。
fcode 发表于 2022-5-27 12:23
是的。
现代Fortran代码,应该只有 Module 和 主程序。没有“外部子程序和外部函数”,即,把所有函数都放 ...
好的,非常感谢! fcode 发表于 2022-5-27 12:23
是的。
现代Fortran代码,应该只有 Module 和 主程序。没有“外部子程序和外部函数”,即,把所有函数都放 ...
另外关于子程序参数我还有一点点问题,就是加入子程序当中我定义的输出参量是A的话,在调用时,我将输出参量同样定义为A然后再输出,以及解耦将输出参量定义为其他名称如B,这两种有什么区别呢。以及第一种为什么不会出现问题呢,明明我对一个变量定义了两次,但是还是可以正常地输出
页:
[1]
2