diff --git a/Makefile b/Makefile index 71514357..4324e966 100644 --- a/Makefile +++ b/Makefile @@ -1,4 +1,4 @@ -CC= gcc +CC= cc #CC= clang --analyze CFLAGS= -g -Wall -Wno-unused-function -O2 WRAP_MALLOC=-DUSE_MALLOC_WRAPPERS diff --git a/bwa.c b/bwa.c index c78ebb67..cd564313 100644 --- a/bwa.c +++ b/bwa.c @@ -144,9 +144,9 @@ uint32_t *bwa_gen_cigar2(const int8_t mat[25], int o_del, int e_del, int o_ins, max_del = (int)((double)(((l_query+1)>>1) * mat[0] - o_del) / e_del + 1.); max_gap = max_ins > max_del? max_ins : max_del; max_gap = max_gap > 1? max_gap : 1; - w = (max_gap + abs(rlen - l_query) + 1) >> 1; + w = (max_gap + llabs(rlen - l_query) + 1) >> 1; w = w < w_? w : w_; - min_w = abs(rlen - l_query) + 3; + min_w = llabs(rlen - l_query) + 3; w = w > min_w? w : min_w; // NW alignment if (bwa_verbose >= 4) { diff --git a/bwase.c b/bwase.c index 77c50db9..23d7d253 100644 --- a/bwase.c +++ b/bwase.c @@ -179,7 +179,7 @@ bwa_cigar_t *bwa_refine_gapped_core(bwtint_t l_pac, const ubyte_t *pacseq, int l assert(re <= l_pac); rseq = bns_get_seq(l_pac, pacseq, rb, re, &rlen); assert(re - rb == rlen); - ksw_global(len, seq, rlen, rseq, 5, mat, 5, 1, SW_BW > abs(rlen - len) * 1.5? SW_BW : abs(rlen - len) * 1.5, n_cigar, &cigar32); + ksw_global(len, seq, rlen, rseq, 5, mat, 5, 1, SW_BW > llabs(rlen - len) * 1.5? SW_BW : llabs(rlen - len) * 1.5, n_cigar, &cigar32); assert(*n_cigar > 0); if ((cigar32[*n_cigar - 1]&0xf) == 1) cigar32[*n_cigar - 1] = (cigar32[*n_cigar - 1]>>4<<4) | 3; // change endding ins to soft clipping if ((cigar32[0]&0xf) == 1) cigar32[0] = (cigar32[0]>>4<<4) | 3; // change beginning ins to soft clipping diff --git a/bwt_lite.h b/bwt_lite.h index 4fadcced..30a148f0 100644 --- a/bwt_lite.h +++ b/bwt_lite.h @@ -1,6 +1,7 @@ #ifndef BWT_LITE_H_ #define BWT_LITE_H_ +#include #include typedef struct { diff --git a/bwtgap.h b/bwtgap.h index 7dd61659..68fceb0f 100644 --- a/bwtgap.h +++ b/bwtgap.h @@ -1,6 +1,7 @@ #ifndef BWTGAP_H_ #define BWTGAP_H_ +#include #include "bwt.h" #include "bwtaln.h" diff --git a/bwtsw2_pair.c b/bwtsw2_pair.c index b945128b..c054e90b 100644 --- a/bwtsw2_pair.c +++ b/bwtsw2_pair.c @@ -236,7 +236,7 @@ void bsw2_pair(const bsw2opt_t *opt, int64_t l_pac, const uint8_t *pac, int n, b double diff; G[0] = hits[i]->hits[0].G + a[1].G; G[1] = hits[i+1]->hits[0].G + a[0].G; - diff = fabs(G[0] - G[1]) / (opt->a + opt->b) / ((hits[i]->hits[0].len + a[1].len + hits[i+1]->hits[0].len + a[0].len) / 2.); + diff = abs(G[0] - G[1]) / (opt->a + opt->b) / ((hits[i]->hits[0].len + a[1].len + hits[i+1]->hits[0].len + a[0].len) / 2.); if (diff > 0.05) a[G[0] > G[1]? 0 : 1].G = 0; } if (a[0].G == 0 || a[1].G == 0) { // one proper pair only diff --git a/kthread.c b/kthread.c index 1b13494f..a2bf6810 100644 --- a/kthread.c +++ b/kthread.c @@ -1,6 +1,7 @@ #include #include #include +#include /************ * kt_for() * diff --git a/maxk.c b/maxk.c index fee5765f..3b6c021a 100644 --- a/maxk.c +++ b/maxk.c @@ -3,6 +3,7 @@ #include #include #include +#include #include #include "bwa.h" #include "bwamem.h"