def test_dwt_idwt_inverse():
"""dwt/idwt should invert each other (to machine precision)"""
from electrolib.filters.dwtfile import WaveletFilter as W
x = np.sin(np.arange(0,10,0.1))
w = W.wavelet('db2')
zi = W.zero_state(w)
cA, cD, zf = W.dwt(x, w, zi)
x_rec = W.idwt(cA, cD, w)
assert np.allclose(x, x_rec)
