Fix pivot offset calculation for negative incx

2022-01-17 00:10:21 +01:00
parent 3b6293f5a0
commit 57e2a72f40
1 changed files with 2 additions and 3 deletions
--- a/lapack/laswp/generic/zlaswp_k_2.c
+++ b/lapack/laswp/generic/zlaswp_k_2.c
@@ -60,10 +60,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;