AccessibilityDisabled.h

Go to the documentation of this file.

/*
 * AccessibilityDisabled.h
 *
 *  Created on: 25.06.2014
 *      Author: Mmann
 */
 
#ifndef INTARNA_ACCESSIBILITYDISABLED_H_
#define INTARNA_ACCESSIBILITYDISABLED_H_
 
#include "IntaRNA/Accessibility.h"
 
 
namespace IntaRNA {
 
class AccessibilityDisabled: public Accessibility {
 
public:
 
    AccessibilityDisabled( const RnaSequence& sequence
                            , const size_t maxLength
                            , const AccessibilityConstraint * const accConstr
                        );
 
    virtual ~AccessibilityDisabled();
 
    virtual
    E_type
    getED( const size_t from, const size_t to ) const;
 
};
 
 
 
 
inline
AccessibilityDisabled::AccessibilityDisabled(const RnaSequence& seq
                , const size_t maxLength
                , const AccessibilityConstraint * const accConstr)
 :
    Accessibility(seq, maxLength, accConstr)
{
#if INTARNA_MULITHREADING
    #pragma omp critical(intarna_omp_logOutput)
#endif
    { VLOG(2) <<"no accessibility requested..."; }
}
 
 
inline
AccessibilityDisabled::~AccessibilityDisabled()
{
}
 
 
inline
E_type
AccessibilityDisabled::
getED( const size_t from, const size_t to ) const
{
    // input check
    checkIndices(from,to);
 
    if ((to-from+1) <= getMaxLength()+1) { // "getMaxLength()+1" to allow for dangling-end probability computation
        // check for constrained end positions
        if (!getAccConstraint().isAccessible(from) || !getAccConstraint().isAccessible(to)) {
            // end position blocked --> omit accessibility
            return ED_UPPER_BOUND;
        }
        // else: no accessibility computation done --> always zero
        return (E_type)0;
    } else {
        // region length exceeds maximally allowed length -> no value
        return ED_UPPER_BOUND;
    }
}
 
 
} // namespace
 
#endif /* ACCESSIBILITYDISABLED_H_ */