| author | Francesco Abbate <francesco.bbt@gmail.com> | 2011年09月14日 10:02:29 +0200 |
|---|---|---|
| committer | Francesco Abbate <francesco.bbt@gmail.com> | 2011年10月19日 14:43:24 +0200 |
| commit | bba69adf6f3a6303c38f46f4d0ba4b7cc7828c32 (patch) | |
| tree | fc22d8cb1d53e4230fef1b1af5cb82904b4aa9d1 | |
| parent | ce039ed3b092f32a6b8435b9c52d3d2b356a2306 (diff) | |
| download | gsl-shell-bba69adf6f3a6303c38f46f4d0ba4b7cc7828c32.tar.gz | |
| -rw-r--r-- | num/lmfit.lua.in | 18 |
diff --git a/num/lmfit.lua.in b/num/lmfit.lua.in index 5d7db424..b79fc121 100644 --- a/num/lmfit.lua.in +++ b/num/lmfit.lua.in @@ -228,19 +228,19 @@ local function qrsolv(r, p, lambda, diag, qtb, x, sdiag, wa) gsl.gsl_matrix_set(r, k, k, new_rkk) gsl.gsl_vector_set(wa, k, new_wak) end - end - -- Accumulate the transformation in the row of s + -- Accumulate the transformation in the row of s - for i = k + 1, n-1 do - local rik = gsl.gsl_matrix_get (r, i, k) - local sdiagi = gsl.gsl_vector_get (sdiag, i) + for i = k + 1, n-1 do + local rik = gsl.gsl_matrix_get (r, i, k) + local sdiagi = gsl.gsl_vector_get (sdiag, i) - local new_rik = cosine * rik + sine * sdiagi - local new_sdiagi = -sine * rik + cosine * sdiagi + local new_rik = cosine * rik + sine * sdiagi + local new_sdiagi = -sine * rik + cosine * sdiagi - gsl.gsl_matrix_set(r, i, k, new_rik) - gsl.gsl_vector_set(sdiag, i, new_sdiagi) + gsl.gsl_matrix_set(r, i, k, new_rik) + gsl.gsl_vector_set(sdiag, i, new_sdiagi) + end end end |