TransitionInfo.cc revision 10448
110448Snilay@cs.wisc.edu/* Copyright (c) 2012 Massachusetts Institute of Technology 210448Snilay@cs.wisc.edu * 310448Snilay@cs.wisc.edu * Permission is hereby granted, free of charge, to any person obtaining a copy 410448Snilay@cs.wisc.edu * of this software and associated documentation files (the "Software"), to deal 510448Snilay@cs.wisc.edu * in the Software without restriction, including without limitation the rights 610448Snilay@cs.wisc.edu * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell 710448Snilay@cs.wisc.edu * copies of the Software, and to permit persons to whom the Software is 810448Snilay@cs.wisc.edu * furnished to do so, subject to the following conditions: 910448Snilay@cs.wisc.edu * 1010448Snilay@cs.wisc.edu * The above copyright notice and this permission notice shall be included in 1110448Snilay@cs.wisc.edu * all copies or substantial portions of the Software. 1210448Snilay@cs.wisc.edu * 1310448Snilay@cs.wisc.edu * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 1410448Snilay@cs.wisc.edu * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 1510448Snilay@cs.wisc.edu * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE 1610448Snilay@cs.wisc.edu * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER 1710448Snilay@cs.wisc.edu * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, 1810448Snilay@cs.wisc.edu * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN 1910448Snilay@cs.wisc.edu * THE SOFTWARE. 2010448Snilay@cs.wisc.edu */ 2110448Snilay@cs.wisc.edu 2210447Snilay@cs.wisc.edu#include "model/TransitionInfo.h" 2310447Snilay@cs.wisc.edu 2410447Snilay@cs.wisc.edunamespace DSENT 2510447Snilay@cs.wisc.edu{ 2610447Snilay@cs.wisc.edu TransitionInfo::TransitionInfo() 2710447Snilay@cs.wisc.edu : m_number_transitions_00_(0.25), m_number_transitions_01_(0.25), 2810447Snilay@cs.wisc.edu m_number_transitions_11_(0.25), m_frequency_multiplier_(1.0) 2910447Snilay@cs.wisc.edu { 3010447Snilay@cs.wisc.edu update(); 3110447Snilay@cs.wisc.edu } 3210447Snilay@cs.wisc.edu 3310447Snilay@cs.wisc.edu TransitionInfo::TransitionInfo(double number_transitions_00_, double number_transitions_01_, 3410447Snilay@cs.wisc.edu double number_transitions_11_) 3510447Snilay@cs.wisc.edu : m_number_transitions_00_(number_transitions_00_), m_number_transitions_01_(number_transitions_01_), 3610447Snilay@cs.wisc.edu m_number_transitions_11_(number_transitions_11_) 3710447Snilay@cs.wisc.edu { 3810447Snilay@cs.wisc.edu m_frequency_multiplier_ = m_number_transitions_00_ + 2.0 * m_number_transitions_01_ + m_number_transitions_11_; 3910447Snilay@cs.wisc.edu update(); 4010447Snilay@cs.wisc.edu } 4110447Snilay@cs.wisc.edu 4210447Snilay@cs.wisc.edu TransitionInfo::~TransitionInfo() 4310447Snilay@cs.wisc.edu {} 4410447Snilay@cs.wisc.edu 4510447Snilay@cs.wisc.edu void TransitionInfo::update() 4610447Snilay@cs.wisc.edu { 4710447Snilay@cs.wisc.edu ASSERT(m_number_transitions_00_ >= 0.0, "[Error] Number of 0->0 transitions (" + 4810447Snilay@cs.wisc.edu (String)m_number_transitions_00_ + ") must be >= 0.0!"); 4910447Snilay@cs.wisc.edu ASSERT(m_number_transitions_01_ >= 0.0, "[Error] Number of 0->1 transitions (" + 5010447Snilay@cs.wisc.edu (String)m_number_transitions_01_ + ") must be >= 0.0!"); 5110447Snilay@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!"); 5210447Snilay@cs.wisc.edu 5310447Snilay@cs.wisc.edu m_probability_1_ = (m_number_transitions_01_ + m_number_transitions_11_) / m_frequency_multiplier_; 5410447Snilay@cs.wisc.edu m_probability_0_ = 1.0 - m_probability_1_; 5510447Snilay@cs.wisc.edu return; 5610447Snilay@cs.wisc.edu } 5710447Snilay@cs.wisc.edu 5810447Snilay@cs.wisc.edu TransitionInfo TransitionInfo::scaleFrequencyMultiplier(double frequency_multiplier_) const 5910447Snilay@cs.wisc.edu { 6010447Snilay@cs.wisc.edu // A short cut if frequency_multiplier_ == m_frequency_multiplier_ to avoid excess calculations 6110447Snilay@cs.wisc.edu if(frequency_multiplier_ == m_frequency_multiplier_) 6210447Snilay@cs.wisc.edu { 6310447Snilay@cs.wisc.edu return TransitionInfo(m_number_transitions_00_, m_number_transitions_01_, m_number_transitions_11_); 6410447Snilay@cs.wisc.edu } 6510447Snilay@cs.wisc.edu else if (frequency_multiplier_ > m_frequency_multiplier_) 6610447Snilay@cs.wisc.edu { 6710447Snilay@cs.wisc.edu double freq_ratio = frequency_multiplier_ / m_frequency_multiplier_; 6810447Snilay@cs.wisc.edu double diff_num_trans_01 = m_number_transitions_01_ * (freq_ratio - 1.0); 6910447Snilay@cs.wisc.edu double norm_num_trans_00 = m_number_transitions_00_ * freq_ratio + diff_num_trans_01; 7010447Snilay@cs.wisc.edu double norm_num_trans_01 = m_number_transitions_01_; 7110447Snilay@cs.wisc.edu double norm_num_trans_11 = m_number_transitions_11_ * freq_ratio + diff_num_trans_01; 7210447Snilay@cs.wisc.edu 7310447Snilay@cs.wisc.edu return TransitionInfo(norm_num_trans_00, norm_num_trans_01, norm_num_trans_11); 7410447Snilay@cs.wisc.edu } 7510447Snilay@cs.wisc.edu else 7610447Snilay@cs.wisc.edu { 7710447Snilay@cs.wisc.edu ASSERT(false, "[Error] Cannot scale to frequency multiplier (" + (String) frequency_multiplier_ + 7810447Snilay@cs.wisc.edu ") since current frequency multiplier (" + (String) m_frequency_multiplier_ + ") is bigger!"); 7910447Snilay@cs.wisc.edu } 8010447Snilay@cs.wisc.edu } 8110447Snilay@cs.wisc.edu 8210447Snilay@cs.wisc.edu void TransitionInfo::print(std::ostream& ost_) const 8310447Snilay@cs.wisc.edu { 8410447Snilay@cs.wisc.edu ost_ << "[" << m_number_transitions_00_ << ", " << m_number_transitions_01_ << ", " << m_number_transitions_11_ << "]"; 8510447Snilay@cs.wisc.edu ost_ << " [" << m_number_transitions_00_ / m_frequency_multiplier_; 8610447Snilay@cs.wisc.edu ost_ << ", " << m_number_transitions_01_*2.0 / m_frequency_multiplier_; 8710447Snilay@cs.wisc.edu ost_ << ", " << m_number_transitions_11_ / m_frequency_multiplier_ << "]" << endl; 8810447Snilay@cs.wisc.edu return; 8910447Snilay@cs.wisc.edu } 9010447Snilay@cs.wisc.edu} // namespace DSENT 9110447Snilay@cs.wisc.edu 92