imach/src/imach.c - diff

Return to imach.c CVS log

Up to [Local Repository] / imach / src

version 1.255, 2017/03/08 16:02:28	version 1.256, 2017/03/27 05:50:23
Line 1	Line 1
/* $Id$	/* $Id$
$State$	$State$
$Log$	$Log$
	Revision 1.256 2017/03/27 05:50:23 brouard
	Summary: Temporary

Revision 1.255 2017/03/08 16:02:28 brouard	Revision 1.255 2017/03/08 16:02:28 brouard
Summary: IMaCh version 0.99r10 bugs in gnuplot fixed	Summary: IMaCh version 0.99r10 bugs in gnuplot fixed

Line 4132 void ludcmp(double *a, int n, int indx	Line 4135 void ludcmp(double *a, int n, int indx
big=0.0;	big=0.0;
for (j=1;j<=n;j++)	for (j=1;j<=n;j++)
if ((temp=fabs(a[i][j])) > big) big=temp;	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;	vv[i]=1.0/big;
}	}
for (j=1;j<=n;j++) {	for (j=1;j<=n;j++) {