Fix pivot offset calculation for negative incx

2022-01-17 00:08:20 +01:00
parent eca2f50b48
commit afa0cece5c
1 changed files with 2 additions and 3 deletions
--- a/lapack/laswp/generic/laswp_k_8.c
+++ b/lapack/laswp/generic/laswp_k_8.c
@@ -78,10 +78,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;