projects
/
pspp-builds.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
sweep.c: Move repeated call out of if ... else
[pspp-builds.git]
/
lib
/
linreg
/
sweep.c
diff --git
a/lib/linreg/sweep.c
b/lib/linreg/sweep.c
index c218456418c75ccc450d452eeba21f1fcbbf5cda..7667e3a555db063502eae81dcdd8188921630872 100644
(file)
--- a/
lib/linreg/sweep.c
+++ b/
lib/linreg/sweep.c
@@
-106,20
+106,18
@@
reg_sweep (gsl_matrix * A, int last_col)
tmp = gsl_matrix_get (A, i, j) -
gsl_matrix_get (A, i, k)
* gsl_matrix_get (A, j, k) / sweep_element;
tmp = gsl_matrix_get (A, i, j) -
gsl_matrix_get (A, i, k)
* gsl_matrix_get (A, j, k) / sweep_element;
- gsl_matrix_set (B, i, j, tmp);
}
else if (j > k)
{
tmp = gsl_matrix_get (A, i, j) -
gsl_matrix_get (A, i, k)
* gsl_matrix_get (A, k, j) / sweep_element;
}
else if (j > k)
{
tmp = gsl_matrix_get (A, i, j) -
gsl_matrix_get (A, i, k)
* gsl_matrix_get (A, k, j) / sweep_element;
- gsl_matrix_set (B, i, j, tmp);
}
else
{
tmp = gsl_matrix_get (A, i, k) / sweep_element;
}
else
{
tmp = gsl_matrix_get (A, i, k) / sweep_element;
- gsl_matrix_set (B, i, j, tmp);
}
}
+ gsl_matrix_set (B, i, j, tmp);
}
}
/*
}
}
/*