From 62ed33e80cd0dbf32988cd225e50306e0f1593fd Mon Sep 17 00:00:00 2001 From: Calvin Morrison Date: Thu, 10 Oct 2013 19:42:00 -0400 Subject: don't depend on math.h's POW by using a bitshift instead: --- kmer_frequency_per_sequence.c | 3 +-- kmer_total_count.c | 2 +- 2 files changed, 2 insertions(+), 3 deletions(-) diff --git a/kmer_frequency_per_sequence.c b/kmer_frequency_per_sequence.c index 02a221b..62b98cb 100644 --- a/kmer_frequency_per_sequence.c +++ b/kmer_frequency_per_sequence.c @@ -27,7 +27,6 @@ int main(int argc, char **argv) { long kmer = 6; size_t len = 0; ssize_t read; - unsigned long width = 0; if(argc != 3) { printf("Please supply a filename and a kmer\n"); @@ -40,7 +39,7 @@ int main(int argc, char **argv) { exit(EXIT_FAILURE); } - width = (int)pow(4, kmer); + const unsigned long width = (unsigned long)1 << (kmer * 2); while ((read = getline(&line, &len, fh)) != -1) { if(line[0] != '>' && (read > kmer)) { diff --git a/kmer_total_count.c b/kmer_total_count.c index 113c979..345bceb 100644 --- a/kmer_total_count.c +++ b/kmer_total_count.c @@ -43,7 +43,7 @@ int main(int argc, char **argv) { const unsigned int kmer = atoi(argv[2]); // width is 4^kmer - const unsigned long width = pow(4, kmer); + const unsigned long width = (unsigned long)1 << (kmer * 2); // malloc our counts matrix unsigned long long * const counts = malloc((width+ 1) * sizeof(unsigned long long)); -- cgit v1.2.3