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

If you are reusing a particular b frequently you could define:

 b(3, 3)=reshape([1, 1, 1, 2, 2, 2, 3, 3, 3], [3, 3])

then you just can do:

 a=a*b

..