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