24 #ifndef OTHER_SEQ_HPP_
25 #define OTHER_SEQ_HPP_
27 #include "SequenceSet.hpp"
30 template<
typename SequenceSetType>
32 avg_size(
const SequenceSetType &seqSet)
34 size_t n_seqs = seqSet.n_seqs();
36 for (
size_t i=0; i<n_seqs; ++i)
37 length+= seqSet->size();
42 template<
typename SequenceSetType>
44 max_size(
const SequenceSetType &seqSet)
47 size_t n_seqs = seqSet.n_seqs();
48 for (
size_t i=0; i<n_seqs; ++i)
49 if (seqSet[i].size() > max_len)
50 max_len=seqSet[i].size();
55 template<
typename SequenceType,
typename MemoryType>
57 check_set(
const SequenceSet<SequenceType, MemoryType> &
set)
59 size_t n_seqs =
set.n_seqs();
61 size_t *val_counting =
new size_t[256];
66 for (
size_t i=0; i<n_seqs; ++i)
68 const SequenceType &seq =
set[i];
71 ++val_counting[static_cast<int>(seq[j])];
76 if (val_counting[j] != 0)
81 if (val_counting[j] != 0)
86 if (val_counting[j] != 0)
89 for (j=123; j<256; ++j)
91 if (val_counting[j] != 0)