I was wondering if there is a fast vector multiply

routine analogous to Goto's fast DAXPY routine. I

did a few tests and found that DAXPY is about

3 times as fast

as an F90 vector multiply or vector add.

i.e.

real(kind=8), dimension(1:N) :: u,v,w

u = v*w

thanks