for(int k=0; k<N; k++) {
float c = cos(TWOPI * (float)k / (float)N);
float s = sin(TWOPI * (float)(-dir*k) / (float)N);
#if !(defined(__LCC__) && __LCC__ <= 121)
simd_vector v = {c,s,c,-s};
cs[k] = v;
#else
cs[k] = _mm_set_ps(c,s,c,-s);
#endif
}
};