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