Basic HTML version of Foils prepared 14 October 1997
Foil 41 Data Parallel Pipeline Algorithm for Grav -- Part II
From Fox Presentation Fall 1995 CPS615 Basic Simulation Track for Computational Science -- Fall Semester 95/96/97. byNancy McCracken and 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)))
