aristov

.SECTION/DM vars;

.var/circ signal[80] =
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;

.SECTION/PM varsp;

.var/circ repr[80] = "g_re.txt";
.var/circ impr[80] = "g_im.txt";

.SECTION/PM program;

jump start; rti; rti; rti;
rti; rti; rti; rti;
rti; rti; rti; rti;
rti; rti; rti; rti;
rti; rti; rti; rti;
rti; rti; rti; rti;
rti; rti; rti; rti;
rti; rti; rti; rti;
rti; rti; rti; rti;
rti; rti; rti; rti;
jump filter; rti; rti; rti;
rti; rti; rti; rti;

start: ar = 4;
    ax0=0x0800;
dm(0x3ffb) = ar;//TSCALE
ar = 999;
dm(0x3ffc) = ar;//TCOUNT
dm(0x3ffd) = ar;//TPERIOD

i4 = repr; l4 = length(repr); m4 = 1; m5 = 0;
i0 = signal; l0 = length(signal); m0 = 0; m1 = 1;
my0 = pm(i4,m4);

i5 = impr; l5 = length(impr);
my1 = pm(i5,m4);

ay1 = 0x00ff;
sr1 = 0;

ifc = 0x00ff;
nop;

ena timer;
imask = b#0000000001;
nop;
zacycl: jump zacycl;

filter: mx0 = io(0);
dm(i0,m1) = mx0, mr = 0;
cntr = 78;
do filt until ce;
    mr = mr + mx0 * my0 (ss), mx0 = dm(i0,m1), my0 = pm(i4,m4);
filt:
    mr = mr - mx0 * my1 (ss), my1 = pm(i5,m4);
mr = mr + mx0 * my0 (ss), mx0 = dm(i0,m0), my0 = pm(i4,m4);
mr = mr - mx0 * my1 (ss), my1 = pm(i5,m4);
mr = mr + mx0 * my0 (ss), my0 = pm(i4,m4);
mr = mr - mx0 * my1 (ss), my1 = pm(i5,m4);
mr = mr(rnd);
if mv sat mr;
io(1) = mr1;
rti;