From abdf6cdd73ec1b83865566e67ff43ce9d50adf2b Mon Sep 17 00:00:00 2001
From: Jiri Pittner <jiri@pittnerovi.com>
Date: Mon, 15 May 2023 16:45:22 +0200
Subject: [PATCH] sanity check of ipiv from cblas_dgesv

---
 nonclass.cc | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

diff --git a/nonclass.cc b/nonclass.cc
index 0a38db7..9f0dcb7 100644
--- a/nonclass.cc
+++ b/nonclass.cc
@@ -26,7 +26,7 @@
 #include "qsort.h"
 #include "fortran.h"
 
-#define IPIV_DEBUG
+#undef IPIV_DEBUG
 
 
 namespace LA {
@@ -155,7 +155,11 @@ static void linear_solve_do(NRMat<double> &A, double *B, const int nrhs, const i
 		for (int i=0; i<n; ++i)  {double t=A[i][i]; if(!finite(t) || std::abs(t) < EPSDET ) {*det=0.; break;} else *det *=t;}
 		//find out whether ipiv are numbered from 0 or from 1
 		int shift=1;
-	        for (int i=0; i<n; ++i) if(ipiv[i]==0) shift=0;
+	        for (int i=0; i<n; ++i)
+			{
+			if(ipiv[i]==0) shift=0;
+			if(ipiv[i]<0 || ipiv[i]>n) laerror("problem with ipiv in clapack_dgesv");
+			}
 #ifdef IPIV_DEBUG
        std::cout <<"shift = "<<shift<<std::endl;
 #endif