Un’operazione vettoriale di base è l’operazione saxpy, si tratta in realtà di una combinazione tra una moltiplicazione vettoriale e un’addizione tra vettori. L’espressione che ci permette di eseguire una tale operazione è la seguente:
y: = alfa * x + y
dove alfa è uno scalare, ma x e y sono vettori.
Ecco una subroutine semplice per questa operazione:
subroutine saxpy (n, alpha, x, y) integer n real alpha, x(*), y(*) c c Saxpy: Compute y := alpha*x + y, c where x and y are vectors of length n (at least). c c Local variables integer i c do 10 i = 1, n y(i) = alpha*x(i) + y(i) 10 continue c return end