static void saxpy(int n, sf_double_complex a, const sf_complex *x, sf_complex *y) /* y += a*x */ { ... }