#include <stdio.h>
#include <ecg/db.h>
#include <ecg/ecgcodes.h>

#define abs(A)	((A) >= 0 ? (A) : -(A))

main(argc, argv)
int argc;
char *argv[];
{
    int filter, time=0, slopecrit, sign, maxslope=0, nslope=0,
        qtime, maxtime, t0, t1, t2, t3, t4, t5, t6, t7, t8, t9,
        ms160, ms200, s2, scmax, scmin = 0, v[DB_MAXSIG];
    DB_Anninfo a;
    DB_Annotation annot;
    static DB_Siginfo s[DB_MAXSIG];

    if (argc < 2) {
        fprintf(stderr, "usage: %s record [threshold]\n", argv[0]);
        exit(1);
    }
    a.name = argv[0]; a.stat = WRITE;
    if (dbinit(argv[1], &a, 1, s, DB_MAXSIG) < 1) exit(2);
    if (sampfreq(NULL) != 250.)
        fprintf(stderr, "warning: %s is designed for 250 Hz input\n",
                argv[0]);
    if (argc > 2) scmin = muvadu(0, atoi(argv[2]));
    if (scmin < 1) scmin = muvadu(0, 1000);
    slopecrit = scmax = 10 * scmin;
    ms160 = strtim("0.16"); ms200 = strtim("0.2"); s2 = strtim("2");
    annot.subtyp = annot.chan = annot.num = 0; annot.aux = NULL;
    (void)getvec(v);
    t9 = t8 = t7 = t6 = t5 = t4 = t3 = t2 = t1 = v[0];

    do {
        filter = (t0 = v[0]) + 4*t1 + 6*t2 + 4*t3 + t4
                - t5         - 4*t6 - 6*t7 - 4*t8 - t9;
        if (time % s2 == 0) {
            if (nslope == 0) {
                slopecrit -= slopecrit >> 4;
                if (slopecrit < scmin) slopecrit = scmin;
            }
            else if (nslope >= 5) {
                slopecrit += slopecrit >> 4;
                if (slopecrit > scmax) slopecrit = scmax;
            }
        }
        if (nslope == 0 && abs(filter) > slopecrit) {
            nslope = 1; maxtime = ms160;
            sign = (filter > 0) ? 1 : -1;
            qtime = time;
        }
        if (nslope != 0) {
            if (filter * sign < -slopecrit) {
                sign = -sign;
                maxtime = (++nslope > 4) ? ms200 : ms160;
            }
            else if (filter * sign > slopecrit &&
                     abs(filter) > maxslope)
                maxslope = abs(filter);
            if (maxtime-- < 0) {
                if (2 <= nslope && nslope <= 4) {
                    slopecrit += ((maxslope>>2) - slopecrit) >> 3;
                    if (slopecrit < scmin) slopecrit = scmin;
                    else if (slopecrit > scmax) slopecrit = scmax;
                    annot.time = strtim("i") - (time - qtime) - 4;
                    annot.anntyp = NORMAL; (void)putann(0, &annot);
                    time = 0;
                }
                else if (nslope >= 5) {
                    annot.time = strtim("i") - (time - qtime) - 4;
                    annot.anntyp = ARFCT; (void)putann(0, &annot);
                }
                nslope = 0;
            }
        }
        t9 = t8; t8 = t7; t7 = t6; t6 = t5; t5 = t4;
        t4 = t3; t3 = t2; t2 = t1; t1 = t0; time++;
    } while (getvec(v) > 0);

    dbquit();
    exit(0);
}