Fix pivot offset calculation for negative incx

2022-01-17 00:07:33 +01:00
parent 0e9e951306
commit eca2f50b48
1 changed files with 2 additions and 3 deletions
--- a/lapack/laswp/generic/laswp_k_4.c
+++ b/lapack/laswp/generic/laswp_k_4.c
@@ -65,10 +65,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;