Stan Tomov wrote:Hello,
I think the bug is actually in the example file, in particular, you have used expressions like sizeof(x1), which in your case is 32 (4 doubles), instead of sizeof(x) or just sizeof(double).
mgates3 wrote:incx = 0 is an error condition; you are not passing a vector. Most of the BLAS routines such as DGEMV explicitly state that INCX must not be zero. The CPU implementation actually DOES swaps of the first element, but for even n, it does an even number of swaps, so the end result is the same. Try with n=3 to see that the first element is swapped. Nonetheless, the user should not call dswap with incx=0 or incy=0.
Users browsing this forum: Google [Bot], hsahasra and 2 guests