Separate HTML for Basic Foil 32 Simple Data Parallel Version of N Body Force Computation -- Grav -- I

Basic HTML version of Foils prepared 14 October 1997

Foil 32 Simple Data Parallel Version of N Body Force Computation -- Grav -- I

From Fox Presentation Fall 1995 CPS615 Basic Simulation Track for Computational Science -- Fall Semester 95/96/97. by Nancy McCracken and Geoffrey C. Fox

function Grav(X,M)

C accepts positions of particles X and masses of particles M

C returns accelerations in Grav

C Uses completely parallel calculation, ignoring anti-symmetry of force

integer, parameter :: N =number of particles
real, parameter :: G = gravitational constant
real array (1:3, 1:N) :: X,M,Grav

! calculates acceleration on body i due to body j in entries ( :,i,j )

real array (1:3, 1:N,1:N) :: A, Xi, Xj, Ms, D, R
logical, array (1:3, 1:N,1:N) ::diag
integer i,j,k

! diag is true for diagonal of N by N slices

forall ( k=1:3, i=1:N, j=1:N ) diag = (i.eq.j)

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

Page produced by wwwfoil on Fri Oct 2 1998