Fix offset calculation for negative incx

2022-01-17 00:09:14 +01:00
parent afa0cece5c
commit 3b6293f5a0
1 changed files with 2 additions and 3 deletions
--- a/lapack/laswp/generic/zlaswp_k_1.c
+++ b/lapack/laswp/generic/zlaswp_k_1.c
@@ -59,10 +59,9 @@ int CNAME(BLASLONG n, BLASLONG k1, BLASLONG k2, FLOAT dummy1, FLOAT dummy4,
  lda *= 2;
  k1 --;

-#ifndef MINUS
 ipiv += k1;
-#else
-  ipiv -= (k2 - 1) * incx;
+#ifdef MINUS
+  ipiv -= (k2 - k1 - 1) * incx;
 #endif

  if (n  <= 0) return 0;