dif8.dsp

载dsp上实现16dct和8dct的程序

.module/ram do_DIF8;        {*** 2 8-way DIFs ***}
.external tmp;
.entry DIF8;

DIF8:
        i0=^tmp;
        i1=^tmp+4;
        i2=^tmp;
        i3=^tmp+4;
        m0=5;

        ax1=dm(i3,m1);
        af=pass ax1, ax0=dm(i2,m1);
        ar=ax0+af, ax1=dm(i3,m1), my0=pm(i6,m6);
        ar=ax0-af, dm(i0,m1)=ar;
        mr=ar*my0(ss);
        af=pass ax1, ax0=dm(i2,m1);
        ar=ax0+af, ax1=dm(i3,m1), my0=pm(i6,m6);
        ar=ax0-af, dm(i0,m1)=ar;
        mr=ar*my0(ss), dm(i1,m1)=mr1;
        af=pass ax1, ax0=dm(i2,m1);
        ar=ax0+af, ax1=dm(i3,m0), my0=pm(i6,m6);
        ar=ax0-af, dm(i0,m1)=ar;
        mr=ar*my0(ss), dm(i1,m1)=mr1;
        af=pass ax1, ax0=dm(i2,m0);
        ar=ax0+af,                my0=pm(i6,m7);
        ar=ax0-af, dm(i0,m0)=ar;
        mr=ar*my0(ss), dm(i1,m1)=mr1;
        dm(i1,m0)=mr1;

        ax1=dm(i3,m1);
        af=pass ax1, ax0=dm(i2,m1);
        ar=ax0+af, ax1=dm(i3,m1), my0=pm(i6,m6);
        ar=ax0-af, dm(i0,m1)=ar;
        mr=ar*my0(ss);
        af=pass ax1, ax0=dm(i2,m1);
        ar=ax0+af, ax1=dm(i3,m1), my0=pm(i6,m6);
        ar=ax0-af, dm(i0,m1)=ar;
        mr=ar*my0(ss), dm(i1,m1)=mr1;
        af=pass ax1, ax0=dm(i2,m1);
        ar=ax0+af, ax1=dm(i3,m1), my0=pm(i6,m6);
        ar=ax0-af, dm(i0,m1)=ar;
        mr=ar*my0(ss), dm(i1,m1)=mr1;
        af=pass ax1, ax0=dm(i2,m1);
        ar=ax0+af,                my0=pm(i6,m6);
        ar=ax0-af, dm(i0,m1)=ar;
        mr=ar*my0(ss), dm(i1,m1)=mr1;
        dm(i1,m1)=mr1;

        rts;                               {*** end 2 8-way DIFs ***}
.endmod;