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

sqs

在下面的程序当中，我使用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就会报错呢

唐汉

你在主程序里面声明了一个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 01:35
你在主程序里面声明了一个f，在子程序里面又声明了一个f，而且，在子程序里面的f没有输出。像这种情况你可 ...

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

gps99

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

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

gps99

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

另外，还有个老办法加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

布衣龙共

楼主这个过程，是非常适合用函数的。
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-3 21:24
楼主这个过程，是非常适合用函数的。
PS：任何程序，无论大小，都不推荐使用common，这东西就让它逝去吧。
...

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

sqs

布衣龙共 发表于 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

这两者又有什么区别呢

唐汉

sqs 发表于 2022-4-4 16:05
而且我不增加输出端，而是直接对主程序当中的变量进行赋值，这样可以么[mw_shl_code=fortran,true]progra ...

这样做是可以的，那是因为contains之后的所有子程序依然在一个“scope”里面。你在一个子程序里面也可以用contains来包含一个或多个子程序，他们也会分享上一层子程序的中的变量。但是我不太建议这么做。这样非常容易混淆。当你的程序越写越长了，我还是建议用module把这些都分开。

fcode

sqs 发表于 2022-4-4 16:05
而且我不增加输出端，而是直接对主程序当中的变量进行赋值，这样可以么[mw_shl_code=fortran,true]progra ...

如楼上所言，这样写的话，func函数就和主程序耦合在一起了。
当你需要把 func函数拷贝到其他地方使用时，就会有点“难受”。就不如单纯的解耦之后的函数那样，直接拷贝后直接调用。

对于小型的程序，当然怎么写差别都不大。对于大型程序，为了提高代码的复用性，推荐各程序单元间尽可能解耦。