TransitionInfo.cc revision 10447
1#include "model/TransitionInfo.h" 2 3namespace DSENT 4{ 5 TransitionInfo::TransitionInfo() 6 : m_number_transitions_00_(0.25), m_number_transitions_01_(0.25), 7 m_number_transitions_11_(0.25), m_frequency_multiplier_(1.0) 8 { 9 update(); 10 } 11 12 TransitionInfo::TransitionInfo(double number_transitions_00_, double number_transitions_01_, 13 double number_transitions_11_) 14 : m_number_transitions_00_(number_transitions_00_), m_number_transitions_01_(number_transitions_01_), 15 m_number_transitions_11_(number_transitions_11_) 16 { 17 m_frequency_multiplier_ = m_number_transitions_00_ + 2.0 * m_number_transitions_01_ + m_number_transitions_11_; 18 update(); 19 } 20 21 TransitionInfo::~TransitionInfo() 22 {} 23 24 void TransitionInfo::update() 25 { 26 ASSERT(m_number_transitions_00_ >= 0.0, "[Error] Number of 0->0 transitions (" + 27 (String)m_number_transitions_00_ + ") must be >= 0.0!"); 28 ASSERT(m_number_transitions_01_ >= 0.0, "[Error] Number of 0->1 transitions (" + 29 (String)m_number_transitions_01_ + ") must be >= 0.0!"); 30 ASSERT(m_number_transitions_11_ >= 0.0, "[Error] Number of 1->1 transitions (" + (String)m_number_transitions_11_ + ") must be >= 0.0!"); 31 32 m_probability_1_ = (m_number_transitions_01_ + m_number_transitions_11_) / m_frequency_multiplier_; 33 m_probability_0_ = 1.0 - m_probability_1_; 34 return; 35 } 36 37 TransitionInfo TransitionInfo::scaleFrequencyMultiplier(double frequency_multiplier_) const 38 { 39 // A short cut if frequency_multiplier_ == m_frequency_multiplier_ to avoid excess calculations 40 if(frequency_multiplier_ == m_frequency_multiplier_) 41 { 42 return TransitionInfo(m_number_transitions_00_, m_number_transitions_01_, m_number_transitions_11_); 43 } 44 else if (frequency_multiplier_ > m_frequency_multiplier_) 45 { 46 double freq_ratio = frequency_multiplier_ / m_frequency_multiplier_; 47 double diff_num_trans_01 = m_number_transitions_01_ * (freq_ratio - 1.0); 48 double norm_num_trans_00 = m_number_transitions_00_ * freq_ratio + diff_num_trans_01; 49 double norm_num_trans_01 = m_number_transitions_01_; 50 double norm_num_trans_11 = m_number_transitions_11_ * freq_ratio + diff_num_trans_01; 51 52 return TransitionInfo(norm_num_trans_00, norm_num_trans_01, norm_num_trans_11); 53 } 54 else 55 { 56 ASSERT(false, "[Error] Cannot scale to frequency multiplier (" + (String) frequency_multiplier_ + 57 ") since current frequency multiplier (" + (String) m_frequency_multiplier_ + ") is bigger!"); 58 } 59 } 60 61 void TransitionInfo::print(std::ostream& ost_) const 62 { 63 ost_ << "[" << m_number_transitions_00_ << ", " << m_number_transitions_01_ << ", " << m_number_transitions_11_ << "]"; 64 ost_ << " [" << m_number_transitions_00_ / m_frequency_multiplier_; 65 ost_ << ", " << m_number_transitions_01_*2.0 / m_frequency_multiplier_; 66 ost_ << ", " << m_number_transitions_11_ / m_frequency_multiplier_ << "]" << endl; 67 return; 68 } 69} // namespace DSENT 70 71