PetaFlop Kernel: d). 3-D Jacobi kernel -- Actual Code.
PetaFlop Kernel: d). 3-D Jacobi kernel -- Actual Code.
- do k = 1, nz
- do j = 1, ny
- do i = 1, nx
- do m = 1, nc
- do mp = 1, nc
- a(i,j,k,m) = u(mp,m) * (b(i+1,j,k,mp) + b(i-1,j,k,mp) &
- b(i,j+1,k,mp) + b(i,j-1,k,mp) + b(i,j,k+1,mp) + &
- b(i,j,k-1,mp)) + a(i,j,k,m)
- enddo
- enddo
- enddo
- enddo
- enddo
Notes:
17