From: N. Brouard <brouard@ined.fr>
Date: Mon, 27 Mar 2017 05:50:23 +0000 (+0000)
Subject: Summary: Temporary
X-Git-Tag: imach-099s7~249
X-Git-Url: https://henry.ined.fr/git/?a=commitdiff_plain;h=7c5f5c7a12cd0f2bf5158aee1208e2dbf58bec60;p=.git

Summary: Temporary
---

diff --git a/src/imach.c b/src/imach.c
index 15cbc1d..43836dc 100644
--- a/src/imach.c
+++ b/src/imach.c
@@ -1,6 +1,9 @@
 /* $Id$
   $State$
   $Log$
+  Revision 1.255  2017/03/08 16:02:28  brouard
+  Summary: IMaCh version 0.99r10 bugs in gnuplot fixed
+
   Revision 1.254  2017/03/08 07:13:00  brouard
   Summary: Fixing data parameter line
 
@@ -4129,7 +4132,16 @@ void ludcmp(double **a, int n, int *indx, double *d)
     big=0.0; 
     for (j=1;j<=n;j++) 
       if ((temp=fabs(a[i][j])) > big) big=temp; 
-    if (big == 0.0) nrerror("Singular matrix in routine ludcmp"); 
+    if (big == 0.0){
+      printf(" Singular Hessian matrix at row %d:\n",i);
+      for (j=1;j<=n;j++) {
+	printf(" a[%d][%d]=%f,",i,j,a[i][j]);
+	fprintf(ficlog," a[%d][%d]=%f,",i,j,a[i][j]);
+      }
+      fflush(ficlog);
+      fclose(ficlog);
+      nrerror("Singular matrix in routine ludcmp"); 
+    }
     vv[i]=1.0/big; 
   } 
   for (j=1;j<=n;j++) {