stencilf90.f
来自「数值计算工具库,C语言编写的,可以直接调用.」· F 代码 · 共 20 行
F
20 行
SUBROUTINE stencilf90(A, B, n, iters)
IMPLICIT NONE
INTEGER, INTENT( IN ) :: n, iters
DOUBLE PRECISION, DIMENSION (n,n,n) :: A, B
DOUBLE PRECISION :: c
INTEGER :: count
c = 1 / 7.
DO count=1,iters
A(2:N-1,2:N-1,2:N-1) = c * (B(2:N-1,2:N-1,2:N-1) + B(3:N,2:N-1,2:N-1) &
+ B(1:N-2,2:N-1,2:N-1) + B(2:N-1,3:N,2:N-1) + B(2:N-1,1:N-2,2:N-1) &
+ B(2:N-1,2:N-1,3:N) + B(2:N-1,2:N-1,1:N-2))
B(2:N-1,2:N-1,2:N-1) = c * (A(2:N-1,2:N-1,2:N-1) + A(3:N,2:N-1,2:N-1) &
+ A(1:N-2,2:N-1,2:N-1) + A(2:N-1,3:N,2:N-1) + A(2:N-1,1:N-2,2:N-1) &
+ A(2:N-1,2:N-1,3:N) + A(2:N-1,2:N-1,1:N-2))
END DO
END SUBROUTINE
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?