+ cleaning up code base
[qpalma.git] / dyn_prog / result_align.cpp
index d0c6944..2386360 100644 (file)
@@ -15,17 +15,19 @@ void increaseFeatureCount(penalty_struct* qparam, int dnaChar, int estChar, doub
 
    //printf("Current index %d est/dna: %d %d score %f\n",currentPos,estChar,dnaChar,estprb);
 
-   //printf("before\n");
-   //int p_idx;
-   //for(p_idx=0;p_idx<currentStruct.len;p_idx++) {
-   //   printf("%f ",currentStruct.limits[p_idx]);
+   //if( estChar == 1 and dnaChar == 1) {
+   //   printf("before\n");
+   //   int p_idx;
+   //   for(p_idx=0;p_idx<currentStruct.len;p_idx++) {
+   //      printf("%f ",currentStruct.limits[p_idx]);
+   //   }
+   //   printf("\n");
+
+   //   for(p_idx=0;p_idx<currentStruct.len;p_idx++) {
+   //      printf("%f ",currentStruct.penalties[p_idx]);
+   //   }
+   //   printf("\n");
    //}
-   //printf("\n");
-
-   //for(p_idx=0;p_idx<currentStruct.len;p_idx++) {
-   //   printf("%f ",currentStruct.penalties[p_idx]);
-   //}
-   //printf("\n");
 
    double value = estprb;
    int Lower = 0;
@@ -56,15 +58,18 @@ void increaseFeatureCount(penalty_struct* qparam, int dnaChar, int estChar, doub
    currentStruct.penalties[Upper] += weightup;
    currentStruct.penalties[Lower] += weightlow; 
 
-   //printf("estprb/Lower/Upper %f %d %d\n",estprb,Lower,Upper);
-   //printf("after\n");
-   //for(p_idx=0;p_idx<currentStruct.len;p_idx++) {
-   //   printf("%f ",currentStruct.limits[p_idx]);
+   //if( estChar == 1 and dnaChar == 1) {
+   //   printf("estprb/Lower/Upper %f %d %d\n",estprb,Lower,Upper);
+   //   int p_idx;
+   //   printf("after\n");
+   //   for(p_idx=0;p_idx<currentStruct.len;p_idx++) {
+   //      printf("%f ",currentStruct.limits[p_idx]);
+   //   }
+   //   printf("\n");
+   //   for(p_idx=0;p_idx<currentStruct.len;p_idx++)
+   //      printf("%f ",currentStruct.penalties[p_idx]);
+   //   printf("\n");
    //}
-   //printf("\n");
-   //for(p_idx=0;p_idx<currentStruct.len;p_idx++)
-   //   printf("%f ",currentStruct.penalties[p_idx]);
-   //printf("\n");
    qparam[currentPos] = currentStruct;
 }