HPC Magazine octobre 2013 - Programmer pour la performance (1ère partie).

Listing 1 - Version C basique de la fonction saxpy de BLAS.


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