OpticalWavelength.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 2310447Snilay@cs.wisc.edu#include "model/optical_graph/OpticalWavelength.h" 2410447Snilay@cs.wisc.edu#include "model/optical_graph/OpticalNode.h" 2510447Snilay@cs.wisc.edu#include "model/optical_graph/OpticalLaser.h" 2610447Snilay@cs.wisc.edu#include "model/optical_graph/OpticalModulator.h" 2710447Snilay@cs.wisc.edu#include "model/optical_graph/OpticalFilter.h" 2810447Snilay@cs.wisc.edu#include "model/optical_graph/OpticalDetector.h" 2910447Snilay@cs.wisc.edu#include "model/optical_graph/OpticalWavelength.h" 3010447Snilay@cs.wisc.edu#include <list> 3110447Snilay@cs.wisc.edu#include <cmath> 3210447Snilay@cs.wisc.edu 3310447Snilay@cs.wisc.edunamespace DSENT 3410447Snilay@cs.wisc.edu{ 3510447Snilay@cs.wisc.edu using std::list; 3610447Snilay@cs.wisc.edu using std::min; 3710447Snilay@cs.wisc.edu 3810447Snilay@cs.wisc.edu OpticalWavelength::OpticalWavelength(const String& instance_name_, const WavelengthGroup& wavelengths_) 3910447Snilay@cs.wisc.edu : m_instance_name_(instance_name_), m_wavelengths_(wavelengths_) 4010447Snilay@cs.wisc.edu { 4110447Snilay@cs.wisc.edu m_data_paths_ = new vector<OpticalDataPath>; 4210447Snilay@cs.wisc.edu } 4310447Snilay@cs.wisc.edu 4410447Snilay@cs.wisc.edu OpticalWavelength::~OpticalWavelength() 4510447Snilay@cs.wisc.edu { 4610447Snilay@cs.wisc.edu delete m_data_paths_; 4710447Snilay@cs.wisc.edu } 4810447Snilay@cs.wisc.edu 4910447Snilay@cs.wisc.edu const String& OpticalWavelength::getInstanceName() const 5010447Snilay@cs.wisc.edu { 5110447Snilay@cs.wisc.edu return m_instance_name_; 5210447Snilay@cs.wisc.edu } 5310447Snilay@cs.wisc.edu 5410447Snilay@cs.wisc.edu void OpticalWavelength::addDataPath(OpticalLaser* laser_, OpticalModulator* modulator_, OpticalDetector* detector_, double loss_) 5510447Snilay@cs.wisc.edu { 5610447Snilay@cs.wisc.edu // Expected wavelengths check 5710447Snilay@cs.wisc.edu ASSERT(laser_->isExpected(getWavelengths()), "[Error] " + getInstanceName() + 5810447Snilay@cs.wisc.edu " -> " + laser_->getInstanceName() + " is not expecting the set wavelengths!"); 5910447Snilay@cs.wisc.edu ASSERT(modulator_->isExpected(getWavelengths()), "[Error] " + getInstanceName() + 6010447Snilay@cs.wisc.edu " -> " + modulator_->getInstanceName() + " is not expecting the set wavelengths!"); 6110447Snilay@cs.wisc.edu ASSERT(detector_->isExpected(getWavelengths()), "[Error] " + getInstanceName() + 6210447Snilay@cs.wisc.edu " -> " + detector_->getInstanceName() + " is not expecting the set wavelengths!"); 6310447Snilay@cs.wisc.edu 6410447Snilay@cs.wisc.edu // Check to see if the modulator and laser already have a data path entry 6510447Snilay@cs.wisc.edu bool entry_exists = false; 6610447Snilay@cs.wisc.edu for (unsigned int i = 0; i < m_data_paths_->size(); ++i) 6710447Snilay@cs.wisc.edu { 6810447Snilay@cs.wisc.edu OpticalDataPath& current = m_data_paths_->at(i); 6910447Snilay@cs.wisc.edu bool current_laser = current.laser == laser_; 7010447Snilay@cs.wisc.edu bool current_modulator = current.modulator == modulator_; 7110447Snilay@cs.wisc.edu 7210447Snilay@cs.wisc.edu ASSERT((current_modulator && current_laser) || !current_modulator, "[Error] " + 7310447Snilay@cs.wisc.edu getInstanceName() + " -> Modulator is the same, but laser is different?"); 7410447Snilay@cs.wisc.edu 7510447Snilay@cs.wisc.edu // If it is already in the table 7610447Snilay@cs.wisc.edu if (current_modulator) 7710447Snilay@cs.wisc.edu { 7810447Snilay@cs.wisc.edu entry_exists = true; 7910447Snilay@cs.wisc.edu current.detectors.push_back(detector_); 8010447Snilay@cs.wisc.edu current.losses.push_back(loss_); 8110447Snilay@cs.wisc.edu } 8210447Snilay@cs.wisc.edu } 8310447Snilay@cs.wisc.edu 8410447Snilay@cs.wisc.edu // If it wasn't found, add the entry 8510447Snilay@cs.wisc.edu if (!entry_exists) 8610447Snilay@cs.wisc.edu m_data_paths_->push_back(OpticalDataPath(laser_, modulator_, detector_, loss_)); 8710447Snilay@cs.wisc.edu return; 8810447Snilay@cs.wisc.edu } 8910447Snilay@cs.wisc.edu 9010447Snilay@cs.wisc.edu const vector<OpticalDataPath>* OpticalWavelength::getDataPaths() const 9110447Snilay@cs.wisc.edu { 9210447Snilay@cs.wisc.edu return (const vector<OpticalDataPath>*) m_data_paths_; 9310447Snilay@cs.wisc.edu } 9410447Snilay@cs.wisc.edu 9510447Snilay@cs.wisc.edu WavelengthGroup OpticalWavelength::getWavelengths() const 9610447Snilay@cs.wisc.edu { 9710447Snilay@cs.wisc.edu return m_wavelengths_; 9810447Snilay@cs.wisc.edu } 9910447Snilay@cs.wisc.edu 10010447Snilay@cs.wisc.edu double OpticalWavelength::getLaserPower(unsigned int number_detectors_) const 10110447Snilay@cs.wisc.edu { 10210447Snilay@cs.wisc.edu ASSERT(number_detectors_ > 0, "[Error] " + getInstanceName() + 10310447Snilay@cs.wisc.edu " -> Number of detectors must be non-zero!"); 10410447Snilay@cs.wisc.edu // Find the number of actual wavelengths 10510447Snilay@cs.wisc.edu int number_wavelengths = getWavelengths().second - getWavelengths().first + 1; 10610447Snilay@cs.wisc.edu // Laser power sum 10710447Snilay@cs.wisc.edu double laser_power_sum = 0; 10810447Snilay@cs.wisc.edu // Loop through all data paths 10910447Snilay@cs.wisc.edu for (unsigned int i = 0; i < getDataPaths()->size(); ++i) 11010447Snilay@cs.wisc.edu { 11110447Snilay@cs.wisc.edu // Get the current data_path 11210447Snilay@cs.wisc.edu const OpticalDataPath& current_path = getDataPaths()->at(i); 11310447Snilay@cs.wisc.edu // Create data structure holding the worstcase detectors 11410447Snilay@cs.wisc.edu list<double>* detectors = new list<double>(); 11510447Snilay@cs.wisc.edu // Get the extinction ratio of the modulator 11610447Snilay@cs.wisc.edu double ER_dB = current_path.modulator->getExtinctionRatio(); 11710447Snilay@cs.wisc.edu // Get the insertion loss of the modulator 11810447Snilay@cs.wisc.edu double IR_dB = current_path.modulator->getInsertionLoss(); 11910447Snilay@cs.wisc.edu // Walk through all detectors in a data path 12010447Snilay@cs.wisc.edu for (unsigned int j = 0; j < current_path.detectors.size(); ++j) 12110447Snilay@cs.wisc.edu { 12210447Snilay@cs.wisc.edu // Convert sensitivity, extinction ratio, and path loss to a required laser power 12310447Snilay@cs.wisc.edu double current_laser_power = current_path.detectors[j]->getSensitivity(ER_dB) * 12410447Snilay@cs.wisc.edu std::pow(10.0, (current_path.losses[j] + IR_dB) / 10.0) * 12510447Snilay@cs.wisc.edu 1.0 / (1.0 - pow(10, -ER_dB / 10)); 12610447Snilay@cs.wisc.edu 12710447Snilay@cs.wisc.edu // Add the laser power 12810447Snilay@cs.wisc.edu detectors->push_back(current_laser_power); 12910447Snilay@cs.wisc.edu } 13010447Snilay@cs.wisc.edu // Cap the number of detectors 13110447Snilay@cs.wisc.edu number_detectors_ = std::min(number_detectors_, (unsigned int) current_path.detectors.size()); 13210447Snilay@cs.wisc.edu // Sort the detectors list in ascending order, only necessary if the number 13310447Snilay@cs.wisc.edu // of detectors is < total number of detectors 13410447Snilay@cs.wisc.edu if (number_detectors_ < detectors->size()) 13510447Snilay@cs.wisc.edu detectors->sort(); 13610447Snilay@cs.wisc.edu // Sum up the laser power from the worst-case detectors 13710447Snilay@cs.wisc.edu list<double>::reverse_iterator iter = detectors->rbegin(); 13810447Snilay@cs.wisc.edu for (unsigned int j = 0; j < number_detectors_; ++j) 13910447Snilay@cs.wisc.edu { 14010447Snilay@cs.wisc.edu laser_power_sum += (*iter) / current_path.laser->getEfficiency(); 14110447Snilay@cs.wisc.edu ++iter; 14210447Snilay@cs.wisc.edu } 14310447Snilay@cs.wisc.edu delete detectors; 14410447Snilay@cs.wisc.edu } 14510447Snilay@cs.wisc.edu return number_wavelengths * laser_power_sum; 14610447Snilay@cs.wisc.edu } 14710447Snilay@cs.wisc.edu 14810447Snilay@cs.wisc.edu} // namespace DSENT 14910447Snilay@cs.wisc.edu 15010447Snilay@cs.wisc.edu 151