Basic HTML version of Foils prepared
22 February 1998
Foil 45 Data Parallel Pipeline Algorithm for Grav -- Part II
From
Fox Presentation Fall 1995 CPS615 Basic Simulation Track for Computational Science --
1998 Enhancements
.
by
Geoffrey C. Fox
do k=1, (N-1)/2 ! Loop over shifts of circulating copy
!Shift Circulating Arrays Xc Mc Ac to the right
Xc = cshift (Xc, dim=2, shift = -1)
Mc = cshift (Mc, dim=2, shift = -1)
Ac = cshift (Ac, dim=2, shift = -1)
! calculate R to be distance over 3-D cordinates
D= Xc-X
R = sqrt (spread (sum (D*D, dim=1), dim=1, ncopies=3)))
D = D/R**3
A = A + Mc*D ! fixed acceleration
Ac = Ac - M*D ! circulating acceleration
end do
©
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 Feb 22 1998
Basic HTML version of Foils prepared 22 February 1998
