/* cqrs.c - Part of the trqt package submitted for the 2006 CinC challenge
Copyright (C) 2006 DA Tironi, R Sassi and LT Mainardi
This software is released under the terms of the GNU General
Public License (http://www.gnu.org/copyleft/gpl.html).
*/
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <math.h>
#define LVETT1 10000
#define LRIGAFILE 43
typedef struct{
char s1[4];
char tempo[8];
char s2[2];
char campione[7];
char s3[5];
char anntyp[1];
char s4[4];
char subtyp[1];
char s5[4];
char chan[1];
char s6[4];
char num[1];
}Rigagrezza;
typedef struct{
char tempo[8];
long campione;
char anntyp[1];
short subtyp;
short chan;
short num;
}Riga;
typedef struct{
short seq;
long riga;
long q;
long t;
long freq;
}Misurazioni;
void filtra(Riga *righe, Misurazioni *misurazioni, int length){
int c1=0, c2=0, c3=0, c4=0, c5=0, c6=0;
char sequenza[]="(p)(N)(t)";
while(*(righe[c4].anntyp) != 'N')
c4++;
c5=c4;
for(c2=c4;c2<length;c2++){
if(c1==9){
misurazioni[c3].riga = (c2-9)+1;
misurazioni[c3].seq = (c3==0 ||misurazioni[c3].riga==misurazioni[c3-1].riga+9+c6)?1:0;
misurazioni[c3].q = righe[c2-6].campione;
misurazioni[c3].t = righe[c2-1].campione;
misurazioni[c3].freq = righe[c2-5].campione - righe[c5].campione;
c1=0;
c3++;
c6=0;
}
if(*(righe[c2].anntyp) == 'N'){
c5=c4;
c4=c2;
}
if(*(righe[c2].anntyp) == sequenza[c1])
c1++;
else if(c1==8 && righe[c2].num>3 && *(righe[c2].anntyp)=='t' && *(righe[c2-1].anntyp)=='t')
c6++;
else
c1=0;
}
}
int main(int argc, char *argv[ ]){
int c1=0, c2=0; //contatori
Misurazioni misurazioni[LVETT1/9];
Rigagrezza rigagrezza; //buffer di lettura delle righe del file
Riga righe[LVETT1]; //elenco delle righe
FILE *fin; //puntatore al file
for(c2=0;c2<LVETT1/9;c2++)
misurazioni[c2].riga=0;
if((fin = fopen(argv[1], "r")) == NULL)
return -1;
while(fread(&rigagrezza, LRIGAFILE, 1, fin)){
if(c1<=LVETT1){
strcpy(righe[c1].tempo, rigagrezza.tempo);
righe[c1].campione = atoi(rigagrezza.campione);
strcpy(righe[c1].anntyp, rigagrezza.anntyp);
righe[c1].subtyp = atoi(rigagrezza.subtyp);
righe[c1].chan = atoi(rigagrezza.chan);
righe[c1].num = atoi(rigagrezza.num);
c1++;
}
}
if(c1>LVETT1){
printf("ERROR: a few lines are missing!");
return -1;
}
filtra(righe, misurazioni, c1);
for(c2=0;misurazioni[c2].riga>0;c2++)
printf("%li\t%li\t%li\t%i\n", misurazioni[c2].freq, misurazioni[c2].q,
misurazioni[c2].t, misurazioni[c2].seq);
return 0;
}