FORALL
Parallelism must be explicit. Eg in F90 style:
U(2:N-1, 2:N-1) =
& 0.25 * (U(2:N-1, 1:N-2) + U(2:N-1, 3:N)
& + U(1:N-2, 2:N-1) + U(3:N, 2:N-1))
When array syntax unwieldy, can use FORALL statement:
FORALL (I = 2:N-1, J = 2:N-1)
& U(I, J) = 0.25 * (U(I, J-1) + U(I, J+1) +
& U(I-1, J) + U(I+1, J))