I believe there is a square root missing (sj) in the calculation of the inverse wavelet transform. It should be: iW = (dj * np.sqrt(dt) / wavelet.cdelta * wavelet.psi(0) * (np.real(W) / np.sqrt(sj)).sum(axis=0))
I believe there is a square root missing (sj) in the calculation of the inverse wavelet transform. It should be:
iW = (dj * np.sqrt(dt) / wavelet.cdelta * wavelet.psi(0) *
(np.real(W) / np.sqrt(sj)).sum(axis=0))