y = y + a * xwhere x and y are vectors, and a is a scalar. DAXPY is short for Double precision A times X Plus Y.

Timings were taken on a 100 MHz IBM RS/6000 Model 43P, using KAI's optimizing C++ compiler (KCC). The sharp drop-off around N=2000 is due to cache effects. Switch settings: KCC +K3 -O3, f77 -O3.