From 90afc1f0572a82bd6eb4532350b368c685b3cc7b Mon Sep 17 00:00:00 2001 From: John Darrington Date: Tue, 15 Nov 2011 15:20:59 +0100 Subject: [PATCH 4/5] sweep.c: Move repeated call out of if ... else This line is identical for all cases, so it only needs to appear once. --- lib/linreg/sweep.c | 4 +--- 1 files changed, 1 insertions(+), 3 deletions(-) diff --git a/lib/linreg/sweep.c b/lib/linreg/sweep.c index c218456..7667e3a 100644 --- 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; - 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; - gsl_matrix_set (B, i, j, tmp); } else { tmp = gsl_matrix_get (A, i, k) / sweep_element; - gsl_matrix_set (B, i, j, tmp); } + gsl_matrix_set (B, i, j, tmp); } } /* -- 1.7.2.5