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