Sequential Guard Rings in Two Dimensions
In analogous 2D sequential case, one could dimension array PHI(?) to PHI(14,14) to hold updated points only. However then points on the edge would need special treatment so that one uses boundary values in update
- Rather dimension PHI(16,16) to include internal and boundary points
- Run loops over x(I) and y(J) from 2 to 15 to only cover internal points
- Preload boundary values in PHI(1, . ), PHI( . , 16), PHI(.,1), PHI( . ,16)
- This is easier and faster as no conditionals (IF statements) in inner loops