From 8d711d36102ae88c429096c740b4d54f09a7c934 Mon Sep 17 00:00:00 2001 From: Prashant Gaurav Date: Mon, 27 Jun 2011 18:03:36 -0400 Subject: [PATCH] fixed issue 17 --- gtfold-mfe/src/main.cc | 13 +++++++------ gtfold-mfe/src/traceback.c | 2 +- 2 files changed, 8 insertions(+), 7 deletions(-) diff --git a/gtfold-mfe/src/main.cc b/gtfold-mfe/src/main.cc index 5a7f707..280ed9e 100644 --- a/gtfold-mfe/src/main.cc +++ b/gtfold-mfe/src/main.cc @@ -243,10 +243,13 @@ int main(int argc, char** argv) { printf("Done.\n\n"); printf("Results:\n"); - printf("- Minimum Free Energy: %12.4f kcal/mol\n", energy/100.00); + if (energy >= MAXENG) + printf("- Minimum Free Energy: %12.4f kcal/mol\n", 0.00); + else + printf("- Minimum Free Energy: %12.4f kcal/mol\n", energy/100.00); printf("- MFE runtime: %9.6f seconds\n", t1); - - + + if (SUBOPT_ENABLED) { t1 = get_seconds(); ss_map_t subopt_data = subopt_traceback(seq.length(), 100.0*suboptDelta); @@ -263,17 +266,15 @@ int main(int argc, char** argv) { t1 = get_seconds(); trace(seq.length(), VERBOSE, UNAMODE, T_MISMATCH); t1 = get_seconds() - t1; - + printf("\n"); print_sequence(seq.length()); print_structure(seq.length()); if (CONS_ENABLED) print_constraints(seq.length()); - if (SHAPE_ENABLED && VERBOSE) print_shapeArray(seq.length()); - save_ct_file(outputFile, seq, energy); printf("\nMFE structure saved in .ct format to %s\n", outputFile.c_str()); diff --git a/gtfold-mfe/src/traceback.c b/gtfold-mfe/src/traceback.c index 1a9511f..6cea549 100644 --- a/gtfold-mfe/src/traceback.c +++ b/gtfold-mfe/src/traceback.c @@ -45,7 +45,7 @@ void trace(int len, int vv, int mode, int mismatch) { structure[i] = 0; if (W[len] >= MAXENG) { - printf("\n No Structure "); + printf("- No Structure \n"); return; }