Permalink
Browse files

include cmath and use std::abs

  • Loading branch information...
1 parent c75d864 commit 1c2397a71e053ed779fa003a9265977e7240f094 @ernfrid ernfrid committed Sep 3, 2016
Showing with 5 additions and 4 deletions.
  1. +5 −4 src/lib/bamrc/BasicStat.cpp
@@ -5,6 +5,7 @@
#include <iomanip>
#include <iostream>
#include <memory>
+#include <cmath>
BasicStat::BasicStat(bool is_indel)
: read_count(0)
@@ -58,15 +59,15 @@ void BasicStat::process_read(bam_pileup1_t const* base) {
if(q2_val > -1) {
//this is in read coordinates. Ignores clipping as q2 may be clipped
- sum_q2_distance += (float) abs(base->qpos - q2_val) / (float) base->b->core.l_qseq;
+ sum_q2_distance += (float) std::abs(base->qpos - q2_val) / (float) base->b->core.l_qseq;
num_q2_reads++;
}
- distances_to_3p.push_back( (float) abs(base->qpos - three_prime_index) / (float) base->b->core.l_qseq);
- sum_3p_distance += (float) abs(base->qpos - three_prime_index) / (float) base->b->core.l_qseq;
+ distances_to_3p.push_back( (float) std::abs(base->qpos - three_prime_index) / (float) base->b->core.l_qseq);
+ sum_3p_distance += (float) std::abs(base->qpos - three_prime_index) / (float) base->b->core.l_qseq;
sum_of_clipped_lengths += clipped_length;
float read_center = (float)clipped_length/2.0;
- sum_event_location += 1.0 - abs((float)(base->qpos - left_clip) - read_center)/read_center;
+ sum_event_location += 1.0 - std::abs((float)(base->qpos - left_clip) - read_center)/read_center;
}
else {

0 comments on commit 1c2397a

Please sign in to comment.