Go to the source code of this file.
|
double | det (double **aa, int n) |
|
double det |
( |
double ** |
aa, |
|
|
int |
n |
|
) |
| |
Definition at line 6 of file qdet.c.
9 double d, big, dum, sum, tem;
17 if ((tem=fabs(aa[i][j])) > big) big=tem;
18 if (big == 0.0)
nrerror (
"Singular matrix in routine dim");
24 for (k=1;k<i;k++) sum -= aa[i][k]*aa[k][j];
31 sum -= aa[i][k]*aa[k][j];
33 if ( (dum=vv[i]*fabs(sum)) >= big) {
47 if (aa[j][j] == 0.0) aa[j][j]=
TINY;
50 for (i=j+1;i<=n;i++) aa[i][j] *= dum;
void free_vector(double *v, int nl, int nh)
void nrerror(char *error_text)
double * vector(int nl, int nh)