Fix pivot index for negative increments

2022-01-17 00:05:33 +01:00
parent 5188aede5d
commit 1b49ef8dcf
1 changed files with 2 additions and 3 deletions
--- a/lapack/laswp/generic/laswp_k_1.c
+++ b/lapack/laswp/generic/laswp_k_1.c
@@ -57,10 +57,9 @@ int CNAME(BLASLONG n, BLASLONG k1, BLASLONG k2, FLOAT dummy1, FLOAT *a, BLASLONG
  a--;
  k1 --;

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

  if (n  <= 0) return 0;