Basic HTML version of Foils prepared August 4 1996

Foil 17 PetaFlop Kernel: d). 3-D Jacobi kernel -- Actual Code.

From Summary of Working Groups at PAWS and PetaSoft Meetings PAWS(Mandalay Beach) and PetSoft(Bodega Bay) -- April 23 and June 17-19,96. by Geoffrey C. Fox


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



Northeast Parallel Architectures Center, Syracuse University, npac@npac.syr.edu

If you have any comments about this server, send e-mail to webmaster@npac.syr.edu.

Page produced by wwwfoil on Sun Apr 11 1999