aboutsummaryrefslogtreecommitdiff
path: root/src/c/quikr.c
diff options
context:
space:
mode:
authorCalvin Morrison <mutantturkey@gmail.com>2013-11-07 17:20:54 -0500
committerCalvin Morrison <mutantturkey@gmail.com>2013-11-07 17:20:54 -0500
commit5edd30aaa30ff843f42cc4f5a7d27253f85e5056 (patch)
treef4e8251aa721c0af056effa8427ceecad567864b /src/c/quikr.c
parentd7ed90dc545e54a09ef206634bb6042e5893c04b (diff)
remove debugging
Diffstat (limited to 'src/c/quikr.c')
-rw-r--r--src/c/quikr.c22
1 files changed, 2 insertions, 20 deletions
diff --git a/src/c/quikr.c b/src/c/quikr.c
index 1b8daa3..670cd3c 100644
--- a/src/c/quikr.c
+++ b/src/c/quikr.c
@@ -153,7 +153,6 @@ int main(int argc, char **argv) {
// load sensing matrix
struct matrix *sensing_matrix = load_sensing_matrix(sensing_matrix_filename, kmer);
-
// get our "rare" counts
while(1) {
rare_width = 0;
@@ -217,23 +216,6 @@ int main(int argc, char **argv) {
sensing_matrix_rare[x*rare_width] = 1.0;
}
- // DEBUG OUR MATRIX
- if(verbose) {
- printf("COUNT_MATRIX\n");
- for(x = 0; x < rare_width; x++)
- printf("%lf\t", count_matrix_rare[x]);
- printf("\n");
-
- printf("SENSING_MATRIX\n");
- for(y = 0; y < sensing_matrix->sequences - 1; y++) {
- for(x = 0; x < rare_width; x++) {
- printf("%lf\t", sensing_matrix_rare[x*rare_width + y]);
- }
- printf("%lf\n", sensing_matrix_rare[x*rare_width + (rare_width - 1)]);
- }
-
- }
-
// run NNLS
double *solution = nnls(sensing_matrix_rare, count_matrix_rare, sensing_matrix->sequences, rare_width);
@@ -246,9 +228,9 @@ int main(int argc, char **argv) {
fprintf(stderr, "Could not open %s for writing\n", output_filename);
exit(EXIT_FAILURE);
}
- for(x = 0; x < sensing_matrix->sequences; x++) {
+
+ for(x = 0; x < sensing_matrix->sequences; x++)
fprintf(output_fh, "%.10lf\n", solution[x]);
- }
fclose(output_fh);