Commit 59e7783e authored by Mathieu Giraud's avatar Mathieu Giraud

core/segment.cpp, core/tools.{h,cpp}: output evalue in >>

parent c8b4cb3b
......@@ -156,6 +156,10 @@ string Segmenter::getInfoLine() const
s += " " + info_extra ;
s += " " + segmented_germline->code ;
s += " " + string(segmented_mesg[because]) ;
if (evalue > NO_LIMIT_VALUE)
s += " " + scientific_string_of_double(evalue);
return s ;
}
......@@ -581,6 +585,7 @@ FineSegmenter::FineSegmenter(Sequence seq, Germline *germline, Cost segment_c)
sequence = seq.sequence ;
Dend=0;
segment_cost=segment_c;
evalue = NO_LIMIT_VALUE;
CDR3start = -1;
CDR3end = -1;
......
......@@ -54,6 +54,13 @@ string string_of_int(int number)
return ss.str();
}
string scientific_string_of_double(double number)
{
stringstream ss;
ss << scientific << number ;
return ss.str();
}
bool is_extended_nucleotide(char nuc) {
switch(nuc) {
......
......@@ -72,7 +72,7 @@ bool pair_occurrence_sort(pair<T, int> a, pair<T, int> b);
string string_of_int(int number);
string scientific_string_of_double(double number);
/**
* @param nuc is A, C, G, T or any extended nucleotide (or lowercase)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment