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];
  }
}