 
     for (; index < n; index += stride) {
         auto x_ =  x[incx * index];
-        auto y_ =  y[incx * index];
+        auto y_ =  y[incy * index];
 
         float c_abs = hypot(x_.x, x_.y);
         float c_arg = atan2f(x_.y, x_.x);


    for (; index < n; index += stride) {
        auto x_ =  x[incx * index];
        auto y_ =  y[incx * index];

        float c_abs = hypot(x_.x, x_.y);
        float c_arg = atan2f(x_.y, x_.x);


    for (; index < n; index += stride) {
        auto x_ =  x[incx * index];
        auto y_ =  y[incy * index];

        double c_abs = hypot(x_.x, x_.y);
        double c_arg = atan2(x_.y, x_.x);
