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.
real(kind=8), dimension(1:N) :: u,v,w
u = v*w