RError.com

RError.com Logo RError.com Logo

RError.com Navigation

  • 主页

Mobile menu

Close
  • 主页
  • 系统&网络
    • 热门问题
    • 最新问题
    • 标签
  • Ubuntu
    • 热门问题
    • 最新问题
    • 标签
  • 帮助
主页 / 问题 / 624178
Accepted
Wissen
Wissen
Asked:2020-02-06 23:18:37 +0000 UTC2020-02-06 23:18:37 +0000 UTC 2020-02-06 23:18:37 +0000 UTC

如何修复 Array 'f' cannot have a deferred shape 错误?

  • 772

我正在尝试使用 GNU Fortran 编译器使用四阶 Runge-Kutta 方法求解具有复杂变量的 SDE。如何解决“Array 'f' cannot have a deferred shape”错误?

program main
      integer, parameter :: n = 3, nsteps = 100
      real, parameter :: a = 0.0, b =1.0
      real ::  x(0:n)

      real :: d, hbar,W,omega0,omega,gamma1, E, delta
      d=3.335641E-30
      hbar=1.0546e-34
      omega = 2.482/(6.24e18*hbar)
      W = 1.25e10
      omega0= 2.4743/(6.24e18*hbar)
      delta=0.0
      gamma1= 1E15
      E=(1/hbar)*W;
      x = (/1.0, 0.0, 0.0,0.0/)
      h = (b - a)/nsteps
      call rk4sys(n,h,x,nsteps)
end program main



subroutine xpsys(n,x,f)
      complex, dimension (0:n) ::  x, f
      integer n

      f(0)= (-(0., 1.))*d*x(1)*E - x(2)*E+gamma1*x(4) ! 00+
      f(1)= ((0., -1.)*(x(2)*(delta)+ d*E*(x(0)- x(3)))-gamma1*x(2)/2)!11
      f(2)= (0., 1.)*(x(2)*(delta)+ d*E*(x(0)- x(3)))-gamma1*x(2)/2!10
      f(3)= (0., 1.)*d*(x(1)*E- x(2)*E)-gamma1*x(3)!01

end subroutine xpsys

subroutine rk4sys(n,h,x,nsteps)
      complex ::  x(0:n), f(:,:)
      real, allocatable :: y(:)
      integer :: i, k, n
      real :: h
      print *,0,x
      allocate (y(0:n), f(0:n,4))
    out:  do k = 1,nsteps
            call xpsys(n,x,f(0,1))
    in1:    do i = 0,n
              y(i) = x(i) + 0.5*h*f(i,1)
            end do in1
            call xpsys(n,y,f(0,2))
    in2:    do i = 0,n
              y(i) = x(i) + 0.5*h*f(i,2)
            end do in2
            call xpsys(n,y,f(0,3))
    in3:    do i = 0,n
              y(i) = x(i) + h*f(i,3)
            end do in3
            call xpsys(n,y,f(0,4))
    in4:    do i = 0,n
              x(i) = x(i) + (h/6.0)* (f(i,1) + 2.0*(f(i,2) + f(i,3)) + f(i,4))
            end do in4
            print *, k, x
          end do out
end subroutine rk4sys
fortran
  • 1 1 个回答
  • 10 Views

1 个回答

  • Voted
  1. Best Answer
    user176262
    2020-02-06T23:52:18Z2020-02-06T23:52:18Z

    在过程中,rk4sys动态分配的数组f必须声明为allocatable:

    subroutine rk4sys(n,h,x,nsteps)
      complex ::  x(0:n)
      complex, allocatable :: f(:,:)
      real, allocatable :: y(:)
      ...
    
    • 3

相关问题

Sidebar

Stats

  • 问题 10021
  • Answers 30001
  • 最佳答案 8000
  • 用户 6900
  • 常问
  • 回答
  • Marko Smith

    Python 3.6 - 安装 MySQL (Windows)

    • 1 个回答
  • Marko Smith

    C++ 编写程序“计算单个岛屿”。填充一个二维数组 12x12 0 和 1

    • 2 个回答
  • Marko Smith

    返回指针的函数

    • 1 个回答
  • Marko Smith

    我使用 django 管理面板添加图像,但它没有显示

    • 1 个回答
  • Marko Smith

    这些条目是什么意思,它们的完整等效项是什么样的

    • 2 个回答
  • Marko Smith

    浏览器仍然缓存文件数据

    • 1 个回答
  • Marko Smith

    在 Excel VBA 中激活工作表的问题

    • 3 个回答
  • Marko Smith

    为什么内置类型中包含复数而小数不包含?

    • 2 个回答
  • Marko Smith

    获得唯一途径

    • 3 个回答
  • Marko Smith

    告诉我一个像幻灯片一样创建滚动的库

    • 1 个回答
  • Martin Hope
    Air 究竟是什么标识了网站访问者? 2020-11-03 15:49:20 +0000 UTC
  • Martin Hope
    Алексей Шиманский 如何以及通过什么方式来查找 Javascript 代码中的错误? 2020-08-03 00:21:37 +0000 UTC
  • Martin Hope
    Qwertiy 号码显示 9223372036854775807 2020-07-11 18:16:49 +0000 UTC
  • Martin Hope
    user216109 如何为黑客设下陷阱,或充分击退攻击? 2020-05-10 02:22:52 +0000 UTC
  • Martin Hope
    Qwertiy 并变成3个无穷大 2020-11-06 07:15:57 +0000 UTC
  • Martin Hope
    koks_rs 什么是样板代码? 2020-10-27 15:43:19 +0000 UTC
  • Martin Hope
    user207618 Codegolf——组合选择算法的实现 2020-10-23 18:46:29 +0000 UTC
  • Martin Hope
    Sirop4ik 向 git 提交发布的正确方法是什么? 2020-10-05 00:02:00 +0000 UTC
  • Martin Hope
    faoxis 为什么在这么多示例中函数都称为 foo? 2020-08-15 04:42:49 +0000 UTC
  • Martin Hope
    Pavel Mayorov 如何从事件或回调函数中返回值?或者至少等他们完成。 2020-08-11 16:49:28 +0000 UTC

热门标签

javascript python java php c# c++ html android jquery mysql

Explore

  • 主页
  • 问题
    • 热门问题
    • 最新问题
  • 标签
  • 帮助

Footer

RError.com

关于我们

  • 关于我们
  • 联系我们

Legal Stuff

  • Privacy Policy

帮助

© 2023 RError.com All Rights Reserve   沪ICP备12040472号-5