diff --git a/ChangeLog b/ChangeLog index 74222b7..f2164a9 100644 --- a/ChangeLog +++ b/ChangeLog @@ -16,6 +16,11 @@ ################################################################################ ################################################################################ +240430 Martin Raden + * IntaRNA/PredictorMfeEns : + * updateZ() + * BUGFIX : no handling of partial Z including ED values (so far not used) + 240429 Martin Raden * IntaRNA/Interaction * Seed::operator"<" @@ -28,9 +33,6 @@ * IntaRNA/SeedHandler : * addSeeds() * BUGFIX : if internal seed base pairs were not found seed was still added - * IntaRNA/PredictorMfeEns : - * updateZ() - * BUGFIX : no handling of partial Z including ED values (so far not used) ################################################################################ ### version 3.4.0 diff --git a/src/IntaRNA/PredictorMfeEns2dSeedExtension.cpp b/src/IntaRNA/PredictorMfeEns2dSeedExtension.cpp index 4b19565..250339a 100644 --- a/src/IntaRNA/PredictorMfeEns2dSeedExtension.cpp +++ b/src/IntaRNA/PredictorMfeEns2dSeedExtension.cpp @@ -102,13 +102,13 @@ predict( const IndexRange & r1, const IndexRange & r2 ) for (size_t l2 = 0; l2< hybridZ_left.size2(); l2++) { // check complementarity of boundary if ( Z_equal(hybridZ_left(l1,l2), 0.0) ) continue; - // ensure max interaction length in seq 1 + // iterate extension right of seed in seq 1 for (size_t r1 = 0; r1 < hybridZ_right.size1() ; r1++) { - // check interaction length + // ensure max interaction length in seq 1 if (sj1+r1-si1+l1 >= energy.getAccessibility1().getMaxLength()) break; - // ensure max interaction length in seq 2 + // iterate extension right of seed in seq 2 for (size_t r2 = 0; r2 < hybridZ_right.size2() ; r2++) { - // check interaction length + // ensure max interaction length in seq 2 if (sj2+r2-si2+l2 >= energy.getAccessibility2().getMaxLength()) break; // check complementarity of boundary if (Z_equal(hybridZ_right(r1,r2),0.0)) continue;