RError.com

RError.com Logo RError.com Logo

RError.com Navigation

  • 主页

Mobile menu

Close
  • 主页
  • 系统&网络
    • 热门问题
    • 最新问题
    • 标签
  • Ubuntu
    • 热门问题
    • 最新问题
    • 标签
  • 帮助
主页 / 问题

问题[fortran]

Martin Hope
ΝNL993
Asked: 2022-07-12 00:02:20 +0000 UTC

Fortran 书籍和学习资源

  • 1

目前,我开始对Fortran编程语言感兴趣,它不是一种非常流行的语言,因此没有足够的教育材料,Fortran 有几个版本也是一个问题:66、77、95、 98、05 和 08(以及其他)。如果对您来说不难,请注明不同版本的材料。

如果有任何俄语材料,那么毫不犹豫地发送(这是优先事项),如果是英语,则添加到已经存在的列表中。

鉴于答案包含大量材料(按 Fortran 标准),那么作为专家,您可以指出哪些材料值得阅读,哪些材料最好绕过。(此选项是为那些找不到更多材料的人提供的)

用书回答格式

书名(年份) - 描述,来自作者 - 链接

对于翻译版本:

书名(书名) - 描述(描述),来自作者(来自翻译作者的翻译) - 链接

文章或文件/文档的答案格式(如 doc、docx、pdf、djvu)

姓名

描述(如果有),年份,

链接)

在批评问题之前,请阅读此评论(由主持人撰写)。

fortran книги
  • 1 个回答
  • 133 Views
Martin Hope
Wissen
Asked: 2020-02-10 00:25:06 +0000 UTC

为什么计算结果正常输出到终端,输出到.txt文件却出错?

  • 1

我正在使用 gfortran,一个 CodeBlocks 框架

program main
  integer, parameter :: n = 3, nsteps = 100
  real, parameter :: a = 0.0, b =10e-14
  complex(8) ::  x(0:n),h
  x = (/1.0, 0.0, 0.0,0.0/)
  h = (b - a)/nsteps
  print *,x
  call rk4sys(n,h,x,nsteps)
end program main

subroutine xpsys(n,x,f)
  complex(8), dimension (0:n) ::  x
  complex(8), dimension (0:n) ::  f
  integer n
  complex(8) :: E,d, hbar,W,omega0,omega,gamma1, delta

  d=3.335641E-30
  hbar=(1.0546E-34,0.0)
  omega = 2.482/(6.24E18*hbar)
  W = 1.25E+10
  omega0= 2.4743/(6.24E18*hbar)
  delta=0.0
  gamma1= 1E15
  E=(1/hbar)*W;


  f(0)= (0.,-1.)*d*(x(1)*E - x(2)*conjg(E))+gamma1*x(3)
  f(1)= conjg((0.,1.)*(x(2)*delta+ d*E*(x(0)- x(3)))-gamma1*x(2)/2);
  f(2)= (0., 1.)*(x(2)*delta + d*E*(x(0)- x(3)))-gamma1*x(2)/2
  f(3)= (0., 1.)*d*(x(1)*E- x(2)*conjg(E))-gamma1*x(3)

  !DrhoDt(1)=-1i*d.*(rho(2).*(E) - rho(3).*conj(E))+gamma.*rho(4);+
  !DrhoDt(2)=conj(1i.*(rho(3).*(delta)+ d.*E.*(rho(1)- rho(4)))-gamma.*rho(3)/2);
  !DrhoDt(3)=1i.*(rho(3).*(delta)+ d.*E.*(rho(1)- rho(4)))-gamma.*rho(3)/2;
  !DrhoDt(4)=1i.*d.*(rho(2).*(E)- rho(3).*conj(E))-gamma.*rho(4);
end subroutine xpsys

subroutine rk4sys(n,h,x,nsteps)
  complex(8), dimension (0:nsteps,0:3) ::  Xmass
  complex(8) ::  x(0:n)
  complex(8), allocatable :: y(:)
  complex(8), allocatable :: f(:,:)
  integer :: i, k, n
  complex(8) :: h
  allocate (y(0:n), f(0:n,4))

  Do i=0,nsteps
    Do j=0,3
        Xmass(i,j)=0
    End Do

  End Do


  Xmass(0,0)=x(0)
  Xmass(0,1)=x(1)
  Xmass(0,2)=x(2)
  Xmass(0,3)=x(3)

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))
      Xmass(k,i) = x(i)
    end do in4
    end do out
  k=0

  OPEN(10,FILE='res.txt')
  !WRITE(10,'(5A25)') 'T','X1','X2','X3',' X4'
  DO k=0,nsteps
    print *,k, real (Xmass(k,0)),AIMAG(Xmass(k,0)), real (Xmass(k,1)), &
         AIMAG(Xmass(k,1)),real (Xmass(k,2)),AIMAG(Xmass(k,2)), &
         real (Xmass(k,3)),AIMAG(Xmass(k,3))

    WRITE(10,'(I3, 8Z19.18)')k, real (Xmass(k,0)),AIMAG(Xmass(k,0)), real (Xmass(k,1)), &
         AIMAG(Xmass(k,1)),real (Xmass(k,2)),AIMAG(Xmass(k,2)), &
         real (Xmass(k,3)),AIMAG(Xmass(k,3))


  END DO
  CLOSE(10)
  Xmass=0
 end subroutine rk4sys

这显示在终端中 - 一切都很好

1 0.90719088395362735 0.0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000Y来-0.290397082089698 0.00000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000来

.txt 中的输出不按顺序

1 003FED07B52D39ECD6 000000000000000000 000000000000000000 00BFD294ECFFDF317C 000000000000000000

fortran
  • 1 个回答
  • 10 Views
Martin Hope
Wissen
Asked: 2020-02-06 23:18:37 +0000 UTC

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

  • 3

我正在尝试使用 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 个回答
  • 10 Views

Sidebar

Stats

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

    我看不懂措辞

    • 1 个回答
  • Marko Smith

    请求的模块“del”不提供名为“default”的导出

    • 3 个回答
  • Marko Smith

    "!+tab" 在 HTML 的 vs 代码中不起作用

    • 5 个回答
  • Marko Smith

    我正在尝试解决“猜词”的问题。Python

    • 2 个回答
  • Marko Smith

    可以使用哪些命令将当前指针移动到指定的提交而不更改工作目录中的文件?

    • 1 个回答
  • Marko Smith

    Python解析野莓

    • 1 个回答
  • Marko Smith

    问题:“警告:检查最新版本的 pip 时出错。”

    • 2 个回答
  • Marko Smith

    帮助编写一个用值填充变量的循环。解决这个问题

    • 2 个回答
  • Marko Smith

    尽管依赖数组为空,但在渲染上调用了 2 次 useEffect

    • 2 个回答
  • Marko Smith

    数据不通过 Telegram.WebApp.sendData 发送

    • 1 个回答
  • Martin Hope
    Alexandr_TT 2020年新年大赛! 2020-12-20 18:20:21 +0000 UTC
  • Martin Hope
    Alexandr_TT 圣诞树动画 2020-12-23 00:38:08 +0000 UTC
  • Martin Hope
    Air 究竟是什么标识了网站访问者? 2020-11-03 15:49:20 +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
    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