A minor release --- mostly fine-tuning of Fortran 90:
DO i = 1, n !...required sequential nature an a(i,i) = x(i) ! impediment to optimisation; ENDDO
forall(i = 1:n) (a(i,i) = x(i)