Broadcast array multiplication in Fortran 90/95

Question

I was wondering that would there be a better (succinct) way to code this in Fortran? I am trying to multiply each column of a(3, 3) by each value in b(3)<

Answer 1

The do loop can be replaced by a one-liner using FORALL:

forall (i=1:3) a(:, i) = a(:, i) * b(i)