INTEGER flushcal_(INTEGER *dummy)
{
flushcal();
return (0L);
}
/* Function calopen reads the specified calibration file; if called with
a NULL argument, it takes the value of the environment variable WFDBCAL
to be the name of the calibration file. If the calibration file name
does not begin with "+", calopen empties the calibration list first;
otherwise, the "+" is discarded before attempting to open the file,
which may be in any directory in the WFDB path. If the file can be read,
its contents are appended to the calibration list. */
FINT calopen(char *cfname)
{
WFDB_FILE *cfile;
char buf[128], *p1, *p2, *p3, *p4, *p5, *p6;
/* If no calibration file is specified, return immediately. */
if (cfname == NULL && (cfname = getenv("WFDBCAL")) == NULL &&
(cfname = DEFWFDBCAL) == NULL)
return (0);
if (*cfname == '+') /* don't empty the calibration list */
cfname++; /* discard the '+' prefix */
else flushcal(); /* empty the calibration list */
/* Quit if file can't be found or opened. */
if ((cfile = wfdb_open(cfname, (char *)NULL, WFDB_READ)) == NULL) {
wfdb_error("calopen: can't read calibration file %s\n", cfname);
return (-2);
}
/* Read a line of the calibration file on each iteration. See wfdbcal(5)
for a description of the format. */
while (wfdb_fgets(buf, 127, cfile)) {
/* ignore leading whitespace */
for (p1 = buf; *p1 == ' ' || *p1 == '\t' || *p1 == '\r'; p1++)
;
/* skip comments, empty lines, and improperly formatted lines */
if (*p1 == '#' ||
(p1 = strtok(p1, "\t")) == NULL || /* signal type */
(p2 = strtok((char *)NULL, " \t")) == NULL || /* low, or '-' */
(p3 = strtok((char *)NULL, " \t")) == NULL || /* high, or '-' */
(p4 = strtok((char *)NULL, " \t")) == NULL || /* pulse type */
(p5 = strtok((char *)NULL, " \t")) == NULL || /* scale */
(p6 = strtok((char *)NULL, " \t\r\n")) == NULL) /* units */
continue;
/* This line appears to be a correctly formatted entry. Allocate
memory for a calibration list entry. */
SUALLOC(this_cle, 1, (sizeof(struct cle)));
/* Fill in the fields of the new calibration list entry. */
SSTRCPY(this_cle->sigtype, p1);
if (strcmp(p2, "-") == 0) {
this_cle->caltype = WFDB_AC_COUPLED;
this_cle->low = 0.0;
}
else {
this_cle->caltype = WFDB_DC_COUPLED;
this_cle->low = atof(p2);
}
if (strcmp(p3, "-") == 0)
this_cle->high = this_cle->low = 0.0;
else
this_cle->high = atof(p3);
if (strcmp(p4, "square") == 0)
this_cle->caltype |= WFDB_CAL_SQUARE;
else if (strcmp(p4, "sine") == 0)
this_cle->caltype |= WFDB_CAL_SINE;
else if (strcmp(p4, "sawtooth") == 0)
this_cle->caltype |= WFDB_CAL_SAWTOOTH;
/* otherwise pulse shape is undefined */
this_cle->scale = atof(p5);
SSTRCPY(this_cle->units, p6);
this_cle->next = NULL;
/* Append the new entry to the end of the list. */
if (first_cle) {
prev_cle->next = this_cle;
prev_cle = this_cle;
}
else
first_cle = prev_cle = this_cle;
}
(void)wfdb_fclose(cfile);
return (0);
}
int fetchsignals(void)
{
int first = 1, framelen, i, imax, imin, j, *m, *mp, n;
WFDB_Calinfo cal;
WFDB_Sample **sb, **sp, *sbo, *spo, *v;
WFDB_Time t, ts0, tsf;
/* Do nothing if no samples were requested. */
if (nosig < 1 || t0 >= tf) return (0);
/* Open the signal calibration database. */
(void)calopen("wfdbcal");
if (tfreq != ffreq) {
ts0 = (WFDB_Time)(t0*tfreq/ffreq + 0.5);
tsf = (WFDB_Time)(tf*tfreq/ffreq + 0.5);
}
else {
ts0 = t0;
tsf = tf;
}
/* Allocate buffers and buffer pointers for each selected signal. */
SUALLOC(sb, nsig, sizeof(WFDB_Sample *));
SUALLOC(sp, nsig, sizeof(WFDB_Sample *));
for (n = framelen = 0; n < nsig; framelen += s[n++].spf)
if (sigmap[n] >= 0) {
SUALLOC(sb[n], (int)((tf-t0)*s[n].spf + 0.5), sizeof(WFDB_Sample));
sp[n] = sb[n];
}
/* Allocate a frame buffer and construct the frame map. */
SUALLOC(v, framelen, sizeof(WFDB_Sample)); /* frame buffer */
SUALLOC(m, framelen, sizeof(int)); /* frame map */
for (i = n = 0; n < nsig; n++) {
for (j = 0; j < s[n].spf; j++)
m[i++] = sigmap[n];
}
for (imax = framelen-1; imax > 0 && m[imax] < 0; imax--)
;
for (imin = 0; imin < imax && m[imin] < 0; imin++)
;
/* Fill the buffers. */
isigsettime(t0);
for (t = t0; t < tf && getframe(v) > 0; t++)
for (i = imin, mp = m + imin; i <= imax; i++, mp++)
if ((n = *mp) >= 0) *(sp[n]++) = v[i];
/* Generate output. */
printf(" { \"signal\":\n [\n");
for (n = 0; n < nsig; n++) {
if (sigmap[n] >= 0) {
char *p;
int delta, prev;
if (!first) printf(",\n");
else first = 0;
printf(" { \"name\": %s,\n", p = strjson(sname[n])); SFREE(p);
if (s[n].units) {
printf(" \"units\": %s,\n", p = strjson(s[n].units));
SFREE(p);
}
else
printf(" \"units\": \"mV\",\n");
printf(" \"t0\": %ld,\n", (long)ts0);
printf(" \"tf\": %ld,\n", (long)tsf);
printf(" \"gain\": %g,\n",
s[n].gain ? s[n].gain : WFDB_DEFGAIN);
printf(" \"base\": %d,\n", s[n].baseline);
printf(" \"tps\": %d,\n", (int)(tfreq/(ffreq*s[n].spf)+0.5));
if (getcal(sname[n], s[n].units, &cal) == 0)
printf(" \"scale\": %g,\n", cal.scale);
else
printf(" \"scale\": 1,\n");
printf(" \"samp\": [ ");
for (sbo = sb[n], prev = 0, spo = sp[n]-1; sbo < spo; sbo++) {
delta = *sbo - prev;
printf("%d,", delta);
prev = *sbo;
}
printf("%d ]\n }", *sbo - prev);
}
}
printf("\n ]%s", nann ? ",\n" : "\n }\n");
flushcal();
for (n = 0; n < nsig; n++)
SFREE(sb[n]);
SFREE(sb);
SFREE(sp);
SFREE(v);
SFREE(m);
return (1); /* output was written */
}
