分段错误 - fortran 90 - 对分子例程

Question

program bisect
real(8) ::  output
call bisection(3.d0,4.d0,2.d0,output)
print*, output
end program bisect

subroutine bisection(a,b,error,result)
real(8):: a,b,error,c,result
logical:: proceed
  proceed = .true.
  do while (proceed)
  if (sin(a)*sin(b).lt. 0.d0) then
     c=(a+b)/2
     if (sin(a)*sin(c).lt.0.d0) then
        b=c
     else
        a=c
     end if
  else  
     stop 'cannot be bisected'
  end if 

  if (abs(a-b).lt. error) then
     proceed = .false.
  end if
 end do
 result= a
end subroutine bisection

同一代码的一个版本上传here。

这是我能想到的最小示例。这会在使用 gfortran 以及在网站上运行可执行文件时产生分段错误。

Answer 1

虚拟参数 a 和 b 与作为常量的实际参数相关联。常量是不可定义的——你的程序试图改变“3.0d0”或“4.0d0”的值。如果你的程序成功了，那么整个宇宙都会爆发混乱。

我强烈推荐：

• 使用模块过程。这允许编译器检查实际参数是否与虚拟参数一致。
• 在您的虚拟参数定义中使用 INTENT 规范。这允许编译器检查需要修改的内容是否是可修改的，以及不可修改的内容是否不是。

解决您的问题的方法是在您的主程序中使用适当的变量来保存初始值 3.0d0 和 4.0d0，并将这些可修改的变量传递给您的子程序。或者，您可以在子例程中创建虚拟参数的临时副本。在 F2003 中，VALUE 属性可用于自动执行此操作。

虽然我们这样做了 - 在所有范围内都使用 IMPLICIT NONE。总是。