Summary: Temporary

diff --git a/src/imach.c b/src/imach.c
index 15cbc1dab607cd40331c5c8743d7954bda6479a0..43836dcdf30b48485f029d9b1500df5a0367d4bb 100644 (file)
--- 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++) {