Removed the numeric array bounds

This commit is contained in:
sago007
2018-02-06 20:33:17 +01:00
parent 8c43d3368f
commit 509f7baeed

View File

@@ -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<nh;i++) {
wksp[j++] = C2*a[i]+C1*a[i+nh]+C0*a[i+1]+C3*a[i+nh1];
wksp[j++] = C3*a[i]-C0*a[i+nh]+C1*a[i+1]-C2*a[i+nh1];
wksp[j++] = C2*a[i-1]+C1*a[i+nh-1]+C0*a[i]+C3*a[i+nh1-1];
wksp[j++] = C3*a[i-1]-C0*a[i+nh-1]+C1*a[i]-C2*a[i+nh1-1];
}
}
for (i=1;i<=n;i++) {
a[i]=wksp[i];
a[i-1]=wksp[i];
}
}