...1991

Unfortunately, Press' code implementing the DWT for wavelets of varying lengths has a serious programming error in the sorting step.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...coef.ma

This direct implementation works fine for Daubechies' 4-point and 6-point wavelets. For wavelets with more coefficients, solving the nonlinear equation system does not succeed in a reasonable amount of time.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

...c_k.

All coefficients of Daubechies base functions of orders 4 to 20 can be found in Daubechies' book 1992. The coefficients I implemented in cofi.rt are taken from it. Daubechies' book also discusses a more efficient method of numerically computing the coefficients c_k than my implementation coef.ma.

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.