HPC Magazine February 2014 - An Introduction to Performance Programming - part I.
Listing 1: A trivial C version of the BLAS saxpy function.
void saxpy(int n, double alpha, double *x, double *y) { int i; for (i = 0 ; i < n ; i++) { y[i] = alpha * x[i] + y[i]; } }