关于调用函数时的报错

sqs

在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

get_r调用自己？核实一下46行是不是get_r

fcode

get_r 是函数，不能用 call 调用。必须 r = get_r(k)

zjk0112

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

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

fcode

可以

sqs

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

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

sqs

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

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

fcode

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

sqs

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

好的，非常感谢！

sqs

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

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