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
|
This is easier and faster as no conditionals (IF statements) in inner loops |