TransitionInfo.cc revision 10447
110447Snilay@cs.wisc.edu#include "model/TransitionInfo.h" 210447Snilay@cs.wisc.edu 310447Snilay@cs.wisc.edunamespace DSENT 410447Snilay@cs.wisc.edu{ 510447Snilay@cs.wisc.edu TransitionInfo::TransitionInfo() 610447Snilay@cs.wisc.edu : m_number_transitions_00_(0.25), m_number_transitions_01_(0.25), 710447Snilay@cs.wisc.edu m_number_transitions_11_(0.25), m_frequency_multiplier_(1.0) 810447Snilay@cs.wisc.edu { 910447Snilay@cs.wisc.edu update(); 1010447Snilay@cs.wisc.edu } 1110447Snilay@cs.wisc.edu 1210447Snilay@cs.wisc.edu TransitionInfo::TransitionInfo(double number_transitions_00_, double number_transitions_01_, 1310447Snilay@cs.wisc.edu double number_transitions_11_) 1410447Snilay@cs.wisc.edu : m_number_transitions_00_(number_transitions_00_), m_number_transitions_01_(number_transitions_01_), 1510447Snilay@cs.wisc.edu m_number_transitions_11_(number_transitions_11_) 1610447Snilay@cs.wisc.edu { 1710447Snilay@cs.wisc.edu m_frequency_multiplier_ = m_number_transitions_00_ + 2.0 * m_number_transitions_01_ + m_number_transitions_11_; 1810447Snilay@cs.wisc.edu update(); 1910447Snilay@cs.wisc.edu } 2010447Snilay@cs.wisc.edu 2110447Snilay@cs.wisc.edu TransitionInfo::~TransitionInfo() 2210447Snilay@cs.wisc.edu {} 2310447Snilay@cs.wisc.edu 2410447Snilay@cs.wisc.edu void TransitionInfo::update() 2510447Snilay@cs.wisc.edu { 2610447Snilay@cs.wisc.edu ASSERT(m_number_transitions_00_ >= 0.0, "[Error] Number of 0->0 transitions (" + 2710447Snilay@cs.wisc.edu (String)m_number_transitions_00_ + ") must be >= 0.0!"); 2810447Snilay@cs.wisc.edu ASSERT(m_number_transitions_01_ >= 0.0, "[Error] Number of 0->1 transitions (" + 2910447Snilay@cs.wisc.edu (String)m_number_transitions_01_ + ") must be >= 0.0!"); 3010447Snilay@cs.wisc.edu ASSERT(m_number_transitions_11_ >= 0.0, "[Error] Number of 1->1 transitions (" + (String)m_number_transitions_11_ + ") must be >= 0.0!"); 3110447Snilay@cs.wisc.edu 3210447Snilay@cs.wisc.edu m_probability_1_ = (m_number_transitions_01_ + m_number_transitions_11_) / m_frequency_multiplier_; 3310447Snilay@cs.wisc.edu m_probability_0_ = 1.0 - m_probability_1_; 3410447Snilay@cs.wisc.edu return; 3510447Snilay@cs.wisc.edu } 3610447Snilay@cs.wisc.edu 3710447Snilay@cs.wisc.edu TransitionInfo TransitionInfo::scaleFrequencyMultiplier(double frequency_multiplier_) const 3810447Snilay@cs.wisc.edu { 3910447Snilay@cs.wisc.edu // A short cut if frequency_multiplier_ == m_frequency_multiplier_ to avoid excess calculations 4010447Snilay@cs.wisc.edu if(frequency_multiplier_ == m_frequency_multiplier_) 4110447Snilay@cs.wisc.edu { 4210447Snilay@cs.wisc.edu return TransitionInfo(m_number_transitions_00_, m_number_transitions_01_, m_number_transitions_11_); 4310447Snilay@cs.wisc.edu } 4410447Snilay@cs.wisc.edu else if (frequency_multiplier_ > m_frequency_multiplier_) 4510447Snilay@cs.wisc.edu { 4610447Snilay@cs.wisc.edu double freq_ratio = frequency_multiplier_ / m_frequency_multiplier_; 4710447Snilay@cs.wisc.edu double diff_num_trans_01 = m_number_transitions_01_ * (freq_ratio - 1.0); 4810447Snilay@cs.wisc.edu double norm_num_trans_00 = m_number_transitions_00_ * freq_ratio + diff_num_trans_01; 4910447Snilay@cs.wisc.edu double norm_num_trans_01 = m_number_transitions_01_; 5010447Snilay@cs.wisc.edu double norm_num_trans_11 = m_number_transitions_11_ * freq_ratio + diff_num_trans_01; 5110447Snilay@cs.wisc.edu 5210447Snilay@cs.wisc.edu return TransitionInfo(norm_num_trans_00, norm_num_trans_01, norm_num_trans_11); 5310447Snilay@cs.wisc.edu } 5410447Snilay@cs.wisc.edu else 5510447Snilay@cs.wisc.edu { 5610447Snilay@cs.wisc.edu ASSERT(false, "[Error] Cannot scale to frequency multiplier (" + (String) frequency_multiplier_ + 5710447Snilay@cs.wisc.edu ") since current frequency multiplier (" + (String) m_frequency_multiplier_ + ") is bigger!"); 5810447Snilay@cs.wisc.edu } 5910447Snilay@cs.wisc.edu } 6010447Snilay@cs.wisc.edu 6110447Snilay@cs.wisc.edu void TransitionInfo::print(std::ostream& ost_) const 6210447Snilay@cs.wisc.edu { 6310447Snilay@cs.wisc.edu ost_ << "[" << m_number_transitions_00_ << ", " << m_number_transitions_01_ << ", " << m_number_transitions_11_ << "]"; 6410447Snilay@cs.wisc.edu ost_ << " [" << m_number_transitions_00_ / m_frequency_multiplier_; 6510447Snilay@cs.wisc.edu ost_ << ", " << m_number_transitions_01_*2.0 / m_frequency_multiplier_; 6610447Snilay@cs.wisc.edu ost_ << ", " << m_number_transitions_11_ / m_frequency_multiplier_ << "]" << endl; 6710447Snilay@cs.wisc.edu return; 6810447Snilay@cs.wisc.edu } 6910447Snilay@cs.wisc.edu} // namespace DSENT 7010447Snilay@cs.wisc.edu 71