blob: 10f56d9f6a846fc87eef6ce74e67d325ebf33aa8 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
--- gnuplot-4.0.0/src/matrix.c.old 2004-04-13 18:23:58.000000000 +0100
+++ gnuplot-4.0.0/src/matrix.c 2005-03-01 12:38:07.068232748 +0000
@@ -290,9 +290,12 @@
*d = 1.0;
for (ar = a, lim = &(a[n]); ar < lim; ar++) {
large = 0.0;
- for (ac = *ar, limc = &(ac[n]); ac < limc;)
- if ((temp = fabs(*ac++)) > large)
+ for (ac = *ar, limc = &(ac[n]); ac < limc;){
+ temp = *ac++;
+ temp = fabs(temp);
+ if (temp > large)
large = temp;
+ }
if (large == 0.0)
int_error(NO_CARET, "Singular matrix in LU-DECOMP");
*dp++ = 1 / large;
|
