调用子程序为什么每一次得出的结果都不一样呢

sqs · 发表于 2022-4-3 01:02:05

在下面的程序当中，我使用gfortran编译器进行编译，但是为什么每次我得到的结果都不一样呢？

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

program para
        implicit none
        real(8)::f
        call func(1.0d0,1.0d0)
        print*,f
        end program
        subroutine func(x,y)
                implicit none
                real(8),intent(in) :: x, y
                real(8) :: f
                f = 2*x+3*y
        end subroutine
~                                                                                                                                                                                    
~

得出的结果如下：
另外为什么我这个地方用real,intent(out)::f就会报错呢

唐汉 · 发表于 2022-4-3 01:35:53

你在主程序里面声明了一个f，在子程序里面又声明了一个f，而且，在子程序里面的f没有输出。像这种情况你可以直接用function，或者subroutine但是加上输出intent(out)

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

program main

    use iso_fortran_env, only: real64
    implicit none

    real(real64) :: f
    call func(1._real64, 1._real64, f)
    print *, f

contains

    subroutine func(x, y, ret)
        real(real64), intent(in) :: x, y
        real(real64), intent(out) :: ret

        ret = 2._real64*x + 3._real64*y
    end subroutine func

end program main

sqs · 发表于 2022-4-3 13:47:08

唐汉发表于 2022-4-3 01:35
你在主程序里面声明了一个f，在子程序里面又声明了一个f，而且，在子程序里面的f没有输出。像这种情况你可 ...

好的，那为什么我每次执行a.out的时候输出的结果都不一致呢

gps99 · 发表于 2022-4-3 20:52:29

sqs 发表于 2022-4-3 13:47
好的，那为什么我每次执行a.out的时候输出的结果都不一致呢

子程序没有输出位。主程序中的 f等于没有赋值，每次运行 print的是个随机地址的初始值

gps99 · 发表于 2022-4-3 21:02:11

本帖最后由 gps99 于 2022-4-3 21:03 编辑

增加输出位是最直接的办法，上面的大咖已经给出代码。

另外，还有个老办法加common公共块。数据块内地址传递，和变量名无关。小程序可以用，数据多的话用module数据块代替

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

program para
        implicit none
        real(8)::f
        common /d1/ f
        call func(1.0d0,1.0d0)
        print*,f
end program para

subroutine func(x,y)
        implicit none
        real(8),intent(in) :: x, y
        real(8) :: f
        common /d1/ f
        f = 2*x+3*y
end subroutine func

布衣龙共 · 发表于 2022-4-3 21:24:25

本帖最后由布衣龙共于 2022-4-3 21:29 编辑

楼主这个过程，是非常适合用函数的。
PS：任何程序，无论大小，都不推荐使用common，这东西就让它逝去吧。

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

program para
  use , intrinsic :: iso_fortran_env, only: R8 => real64
  implicit none
  write(*,*) func(1.0_R8,1.0_R8)
contains
  Pure real(R8) Function func(x,y) result(f)
    real(R8),intent(in) :: x, y
    f = 2._R8*x+3._R8*y
  End Function func
end program para

sqs · 发表于 2022-4-4 15:55:21

布衣龙共发表于 2022-4-3 21:24
楼主这个过程，是非常适合用函数的。
PS：任何程序，无论大小，都不推荐使用common，这东西就让它逝去吧。
...

啊哈，但是一旦我需要有三四个以上的输出结果，那这个时候就只能使用subroutine而不能使用function了呀，是这样吧

sqs · 发表于 2022-4-4 16:05:30

布衣龙共发表于 2022-4-3 21:24
楼主这个过程，是非常适合用函数的。
PS：任何程序，无论大小，都不推荐使用common，这东西就让它逝去吧。
...

而且我不增加输出端，而是直接对主程序当中的变量进行赋值，这样可以么

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

program main
    use iso_fortran_env, only: real64
    implicit none
    real(real64) :: f
    call func(1._real64, 1._real64)
    print *, f
contains
    subroutine func(x, y)
        real(real64), intent(in) :: x, y
        f = 2._real64*x + 3._real64*y
    end subroutine func
end program main

这两者又有什么区别呢