From 3a33338425c6206f5be5e2d8e4ed4d883bed994d Mon Sep 17 00:00:00 2001 From: Roman Zeyde Date: Wed, 8 Apr 2015 11:07:35 +0300 Subject: [PATCH] levinson: reuse previous vectors --- amodem/levinson.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/amodem/levinson.py b/amodem/levinson.py index d2f2cd8..5280067 100644 --- a/amodem/levinson.py +++ b/amodem/levinson.py @@ -16,8 +16,8 @@ def solver(t, y): prev_b = b[-1] ef = sum(t[n-i] * prev_f[i] for i in range(n)) eb = sum(t[i+1] * prev_b[i] for i in range(n)) - f_ = np.concatenate([f[-1], [0]]) - b_ = np.concatenate([[0], b[-1]]) + f_ = np.concatenate([prev_f, [0]]) + b_ = np.concatenate([[0], prev_b]) det = 1.0 - ef * eb f.append((f_ - ef * b_) / det) b.append((b_ - eb * f_) / det)