 // (See accompanying file LICENSE or copy at
 //  http://opensource.org/licenses/MIT)
 //=======================================================================
-
 #include "egblas/axpy.hpp"
+#include <iostream>
 
 #include "complex.hpp"
 

    const auto stride = blockDim.x * gridDim.x;

    for (; index < n; index += stride) {
        y[incy * index] = alpha * x[incx * index] + y[incy * index];
    }
}


    } else if (alpha == 0.0f) {
        axpy_kernel0_run(n, y, incy);
    } else {
        if (incx == 1 && incy == 1) {
            axpy_kernel_run_flat(n, alpha, x, y);
        } else {
            axpy_kernel_run(n, alpha, x, incx, y, incy);
        }
    }
}

