lili/amodem
1
0
Fork 0
mirror of https://github.com/romanz/amodem.git synced 2026-04-28 03:56:44 +08:00
Code Issues Packages Projects Releases Wiki Activity

levinson: reuse previous vectors

Browse Source
This commit is contained in:
Roman Zeyde
2015-04-08 11:07:35 +03:00
parent 18462289f8
commit 3a33338425
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
amodem/levinson.py
View File

@@ -16,8 +16,8 @@ def solver(t, y):
prev_b = b[-1] prev_b = b[-1]
ef = sum(t[n-i] * prev_f[i] for i in range(n)) 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)) eb = sum(t[i+1] * prev_b[i] for i in range(n))
f_ = np.concatenate([f[-1], [0]]) f_ = np.concatenate([prev_f, [0]])
b_ = np.concatenate([[0], b[-1]]) b_ = np.concatenate([[0], prev_b])
det = 1.0 - ef * eb det = 1.0 - ef * eb
f.append((f_ - ef * b_) / det) f.append((f_ - ef * b_) / det)
b.append((b_ - eb * f_) / det) b.append((b_ - eb * f_) / det)