Basic HTML version of Foils prepared 20 October 1997

Foil 36 The 9 Fortran Phases in Runge Kutta Update

From Fox Presentation Fall 1995 CPS615 Basic Simulation Track for Computational Science -- Fall Semester 97. by Geoffrey C. Fox


1 1)XdeltaA(i) = V(i)
2 Xuse(i) = X(i) running over all i in Processor
3 2)VdeltaA(i) = MPGrav(i) involves communication
4 3)XdeltaB(i) = V(i) + h*VdeltaA(i)/2
5 Xuse(i) = X(i) +h*XdeltaA(i)/2
6 4)VdeltaB(i) = MPGrav(i) involves communication
7 5)XdeltaC(i) = V(i) + h*VdeltaB(i)/2
8 Xuse(i) = X(i) +h*XdeltaB(i)/2
9 6)VdeltaC(i) = MPGrav(i) involves communication
10 7)XdeltaD(i) = V(i) + h*VdeltaC(i)
11 Xuse(i) = X(i) +h*XdeltaC(i)
12 8)VdeltaD(i) = MPGrav(i) involves communication
13 9)X(i) becomes (X(i) + h*(XdeltaA(i)+2*XdeltaB(i)+2*XdeltaC(i)+XdeltaD(i))/6 )
14 V(i) becomes (V(i) + h*(VdeltaA(i)+2*VdeltaB(i)+2*VdeltaC(i)+VdeltaD(i))/6 )

in Table To:


© 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 Fri Oct 2 1998