syr2: Y uses incy, not incx

Author: Mark Gates

src/library/blas/gens/clTemplates/her2.cl

@@ -138,8 +138,8 @@ __kernel void %PREFIXher2_CL_kernel( __global const %TYPE* _A, __global const %T
 				%CONJUGATE(1, res1);
 			#endif
             %MUL( res5, alpha, res1 );
-            res1 = Y[c * incx];
-            res3 = Y[r * incx];
+            res1 = Y[c * incy];
+            res3 = Y[r * incy];
             #ifndef HER2_ROWMAJOR
 				%CONJUGATE(1, res1);
             #endif
@@ -461,8 +461,8 @@ __kernel void %PREFIXher2_CU_kernel( __global const %TYPE* _A, __global const %T
 				%CONJUGATE(1, res1);
 			#endif
             %MUL( res5, alpha, res1 );
-            res1 = Y[c * incx];
-            res3 = Y[r * incx];
+            res1 = Y[c * incy];
+            res3 = Y[r * incy];
             #ifndef HER2_ROWMAJOR
 				%CONJUGATE(1, res1);
             #endif

src/library/blas/gens/clTemplates/syr2.cl

@@ -128,9 +128,9 @@ __kernel void %PREFIXsyr2_CL_kernel( __global const %TYPE* _A, __global const %T
 
             %TYPE res1, res2;
             res1 = alpha * X[c * incx];
-			res2 = alpha * X[r * incx];
-            res1 = res1 * Y[r * incx];
-			res2 = res2 * Y[c * incx];
+            res2 = alpha * X[r * incx];
+            res1 = res1 * Y[r * incy];
+            res2 = res2 * Y[c * incy];
 
 			A[r + c * lda] += (res1 + res2);
         }

src/library/blas/gens/clTemplates/syr2_her2.cl

@@ -148,8 +148,8 @@ __kernel void %PREFIXsyr2_her2_CL_kernel( __global %TYPE* _A, __global const %TY
 				#endif
 			#endif
             %MUL( res5, alpha, res1 );
-            res1 = Y[c * incx];
-            res3 = Y[r * incx];
+            res1 = Y[c * incy];
+            res3 = Y[r * incy];
             #ifdef HER2_ONLY
 				#ifndef HER2_ROWMAJOR
 					%CONJUGATE(1, res1);
@@ -514,8 +514,8 @@ __kernel void %PREFIXsyr2_her2_CU_kernel( __global %TYPE* _A, __global const %TY
 				#endif
 			#endif
             %MUL( res5, alpha, res1 );
-            res1 = Y[c * incx];
-            res3 = Y[r * incx];
+            res1 = Y[c * incy];
+            res3 = Y[r * incy];
             #ifdef HER2_ONLY
 				#ifndef HER2_ROWMAJOR
 					%CONJUGATE(1, res1);