# NOT RUN {
## create a random series, calculate a forward
## transform, then invert
x <- rnorm(1024)
x.dwt <- reconstruct(wavDWT(x))
x.modwt <- reconstruct(wavMODWT(x))
all(c(splus2R::vecnorm(x.dwt-x), splus2R::vecnorm(x.modwt-x)) < .Machine$single.eps)
# }
Run the code above in your browser using DataLab