diff --git a/code/client/snd_wavelet.c b/code/client/snd_wavelet.c index 273862ca..c569b43c 100644 --- a/code/client/snd_wavelet.c +++ b/code/client/snd_wavelet.c @@ -30,7 +30,7 @@ Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA static void daub4(float b[], unsigned long n, int isign) { float wksp[4097] = { 0.0f }; - float *a=b-1; // numerical recipies so a[1] = b[0] + float *a=b; unsigned long nh,nh1,i,j; @@ -39,21 +39,21 @@ static void daub4(float b[], unsigned long n, int isign) nh1=(nh=n >> 1)+1; if (isign >= 0) { for (i=1,j=1;j<=n-3;j+=2,i++) { - wksp[i] = C0*a[j]+C1*a[j+1]+C2*a[j+2]+C3*a[j+3]; - wksp[i+nh] = C3*a[j]-C2*a[j+1]+C1*a[j+2]-C0*a[j+3]; + wksp[i] = C0*a[j-1]+C1*a[j]+C2*a[j+1]+C3*a[j+2]; + wksp[i+nh] = C3*a[j-1]-C2*a[j]+C1*a[j+1]-C0*a[j+2]; } - wksp[i ] = C0*a[n-1]+C1*a[n]+C2*a[1]+C3*a[2]; - wksp[i+nh] = C3*a[n-1]-C2*a[n]+C1*a[1]-C0*a[2]; + wksp[i ] = C0*a[n-2]+C1*a[n-1]+C2*a[0]+C3*a[1]; + wksp[i+nh] = C3*a[n-2]-C2*a[n-1]+C1*a[0]-C0*a[1]; } else { - wksp[1] = C2*a[nh]+C1*a[n]+C0*a[1]+C3*a[nh1]; - wksp[2] = C3*a[nh]-C0*a[n]+C1*a[1]-C2*a[nh1]; + wksp[1] = C2*a[nh-1]+C1*a[n-1]+C0*a[0]+C3*a[nh1-1]; + wksp[2] = C3*a[nh-1]-C0*a[n-1]+C1*a[0]-C2*a[nh1-1]; for (i=1,j=3;i