64 lines
2.4 KiB
C
64 lines
2.4 KiB
C
|
#ifndef __STREAM_H
|
||
|
|
#define __STREAM_H
|
||
|
|
|
||
|
|
#include <vector>
|
||
|
|
#include "config.h"
|
||
|
|
#include "flexitree.h"
|
||
|
|
|
||
|
|
using std::vector;
|
||
|
|
|
||
|
|
class Stream {
|
||
|
|
public:
|
||
|
|
Stream() {};
|
||
|
|
void Init(const Config &cfg, const int intern_id, const int
|
||
|
|
extern_id);
|
||
|
|
void Append(const int next_value);
|
||
|
|
int AddToDB(SeqForest &normal, int &db_size, int total_pairs_read,
|
||
|
|
const Config &cfg) const;
|
||
|
|
void CompareSeq(const Config &cfg, const SeqForest &normal, const
|
||
|
|
int total_pairs_read);
|
||
|
|
int GetMaxHDist(void) const {return max_hdist;}
|
||
|
|
int GetMaxLFC(void) const {return max_lfc;}
|
||
|
|
int Ready(void) const {return ready;}
|
||
|
|
int GetNumAnoms(void) const {return num_anoms;}
|
||
|
|
int GetNumPairsRead(void) const {return num_pairs_read;}
|
||
|
|
int GetNumSeqsFnd(void) const {return num_seqs_fnd;}
|
||
|
|
private:
|
||
|
|
vector<int> current_seq; // current sequence being filled or
|
||
|
|
// processed
|
||
|
|
int num_in_seq; // current_seq is full up through
|
||
|
|
// num_in_seq
|
||
|
|
int num_pairs_read; // the number of input pairs belonging to
|
||
|
|
// this stream that have been read so far
|
||
|
|
int num_anoms; // the number of anomalies found so far
|
||
|
|
int num_seqs_fnd; // the number of (not necessarily unique)
|
||
|
|
// sequences belonging to this stream
|
||
|
|
// found so far
|
||
|
|
int ext_sid; // the external stream id
|
||
|
|
int int_sid; // the internal stream id
|
||
|
|
int max_hdist; // the largest minimum Hamming distance
|
||
|
|
// found in this stream
|
||
|
|
int seq_hdist; // the minimum Hamming distance for
|
||
|
|
// current_seq
|
||
|
|
vector<int> lf; // array for locality frame
|
||
|
|
int seq_lfc; // the locality frame count for this
|
||
|
|
// sequence
|
||
|
|
int max_lfc; // the largest locality frame count
|
||
|
|
// encountered so far
|
||
|
|
int ready; // a flag to indicate whether this stream
|
||
|
|
// has a full sequence ready to be
|
||
|
|
// processed. 0 = no, 1 = yes.
|
||
|
|
int seq_len; // sequence length
|
||
|
|
int ComputeMisses(const SeqForest &normal);
|
||
|
|
void ComputeHDist(const SeqForest &normal);
|
||
|
|
void ComputeLF(const int is_anom, const int lf_size);
|
||
|
|
void ReportSeq(const Config &cfg, const int total_pairs_read,
|
||
|
|
const int is_anom) const;
|
||
|
|
void ReportNewSeq(const Config &cfg, const int total_pairs_read,
|
||
|
|
const int db_size) const;
|
||
|
|
};
|
||
|
|
|
||
|
|
#endif
|
||
|
|
|
||
|
|
|