Basic HTML version of Foils prepared
2 October 1995
Foil 8 SEQUENTIAL LAPLACE PROGRAMMING
JACOBI ITERATION IN ONE DIMENSION
(constant in y direction)
From
New CPS615 Foils-- D 23 September 95 CPS615 Basic Simulation Track for Computational Science --
Fall Semester 95
.
by
Geoffrey C. Fox
DIMENSION PHINEW(NTOT), PHIOLD(NTOT)
NTOT1 = NTOT-1
BEGIN TEST=0
DO 1 I=2, NTOT1
TEMP = 0.5 * (PHIOLD(I-1) + PHIOLD(I+1))
TEST = AMAX1 (TEST, ABS (TEMP-PHIOLD(I)))
1 PHINEW (I) = TEMP
DO 2 I=2, NTOT1
2 PHIOLD (I) = PHINEW(I)
IF (TEST > CONVG) GO TO BEGIN
ELSE STOP
©
on Tue Oct 7 1997
