pareto.cc
上传用户:rrhhcc
上传日期:2015-12-11
资源大小:54129k
文件大小:6k
源码类别:

通讯编程

开发平台:

Visual C++

  1. /* -*- Mode:C++; c-basic-offset:8; tab-width:8; indent-tabs-mode:t -*- */
  2. /*
  3.  * Copyright (c) Xerox Corporation 1997. All rights reserved.
  4.  *  
  5.  * This program is free software; you can redistribute it and/or modify it
  6.  * under the terms of the GNU General Public License as published by the
  7.  * Free Software Foundation; either version 2 of the License, or (at your
  8.  * option) any later version.
  9.  *
  10.  * This program is distributed in the hope that it will be useful, but
  11.  * WITHOUT ANY WARRANTY; without even the implied warranty of
  12.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13.  * General Public License for more details.
  14.  *
  15.  * You should have received a copy of the GNU General Public License along
  16.  * with this program; if not, write to the Free Software Foundation, Inc.,
  17.  * 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
  18.  *
  19.  * Linking this file statically or dynamically with other modules is making
  20.  * a combined work based on this file.  Thus, the terms and conditions of
  21.  * the GNU General Public License cover the whole combination.
  22.  *
  23.  * In addition, as a special exception, the copyright holders of this file
  24.  * give you permission to combine this file with free software programs or
  25.  * libraries that are released under the GNU LGPL and with code included in
  26.  * the standard release of ns-2 under the Apache 2.0 license or under
  27.  * otherwise-compatible licenses with advertising requirements (or modified
  28.  * versions of such code, with unchanged license).  You may copy and
  29.  * distribute such a system following the terms of the GNU GPL for this
  30.  * file and the licenses of the other code concerned, provided that you
  31.  * include the source code of that other code when and as the GNU GPL
  32.  * requires distribution of source code.
  33.  *
  34.  * Note that people who make modified versions of this file are not
  35.  * obligated to grant this special exception for their modified versions;
  36.  * it is their choice whether to do so.  The GNU General Public License
  37.  * gives permission to release a modified version without this exception;
  38.  * this exception also makes it possible to release a modified version
  39.  * which carries forward this exception.
  40.  */
  41. #ifndef lint
  42. static const char rcsid[] =
  43.     "@(#) $Header: /cvsroot/nsnam/ns-2/tools/pareto.cc,v 1.9 2005/08/26 05:05:31 tomh Exp $ (Xerox)";
  44. #endif
  45.  
  46. #include "random.h"
  47. #include "trafgen.h"
  48. /* implement an on/off source with average on and off times taken
  49.  * from a pareto distribution.  (enough of these sources multiplexed
  50.  * produces aggregate traffic that is LRD).  It is parameterized
  51.  * by the average burst time, average idle time, burst rate, and
  52.  * pareto shape parameter and packet size.
  53.  */
  54. class POO_Traffic : public TrafficGenerator {
  55.  public:
  56. POO_Traffic();
  57. virtual double next_interval(int&);
  58. int on()  { return on_ ; }
  59.   // Added by Debojyoti Dutta October 12th 2000
  60. int command(int argc, const char*const* argv);
  61. protected:
  62. void init();
  63. double ontime_;  /* average length of burst (sec) */
  64. double offtime_; /* average idle period (sec) */
  65. double rate_;    /* send rate during burst (bps) */
  66. double interval_; /* inter-packet time at burst rate */
  67. double burstlen_; /* average # packets/burst */
  68. double shape_;    /* pareto shape parameter */
  69. unsigned int rem_; /* number of packets remaining in current burst */
  70. double p1_;       /* parameter for pareto distribution to compute
  71.    * number of packets in burst.
  72.            */
  73. double p2_;       /* parameter for pareto distribution to compute
  74.          * length of idle period.
  75.          */
  76. int on_;          /* denotes whether in the on or off state */
  77. // Added by Debojyoti Dutta 13th October 2000
  78. RNG * rng_; /* If the user wants to specify his own RNG object */
  79. };
  80. static class POOTrafficClass : public TclClass {
  81.  public:
  82. POOTrafficClass() : TclClass("Application/Traffic/Pareto") {}
  83.   TclObject* create(int, const char*const*) {
  84. return (new POO_Traffic());
  85. }
  86. } class_poo_traffic;
  87. // Added by Debojyoti Dutta October 12th 2000
  88. // This is a new command that allows us to use 
  89. // our own RNG object for random number generation
  90. // when generating application traffic
  91. int POO_Traffic::command(int argc, const char*const* argv){
  92.         
  93. Tcl& tcl = Tcl::instance();
  94.         if(argc==3){
  95.                 if (strcmp(argv[1], "use-rng") == 0) {
  96. rng_ = (RNG*)TclObject::lookup(argv[2]);
  97. if (rng_ == 0) {
  98. tcl.resultf("no such RNG %s", argv[2]);
  99. return(TCL_ERROR);
  100. }                        
  101. return (TCL_OK);
  102.                 }
  103.         }
  104.         return Application::command(argc,argv);
  105. }
  106. POO_Traffic::POO_Traffic() : rng_(NULL)
  107. {
  108. bind_time("burst_time_", &ontime_);
  109. bind_time("idle_time_", &offtime_);
  110. bind_bw("rate_", &rate_);
  111. bind("shape_", &shape_);
  112. bind("packetSize_", &size_);
  113. }
  114. void POO_Traffic::init()
  115. {
  116. interval_ = (double)(size_ << 3)/(double)rate_;
  117. burstlen_ = ontime_/interval_;
  118. rem_ = 0;
  119. on_ = 0;
  120. p1_ = burstlen_ * (shape_ - 1.0)/shape_;
  121. p2_ = offtime_ * (shape_ - 1.0)/shape_;
  122. if (agent_)
  123. agent_->set_pkttype(PT_PARETO);
  124. }
  125. double POO_Traffic::next_interval(int& size)
  126. {
  127. double t = interval_;
  128. on_ = 1;
  129. if (rem_ == 0) {
  130. /* compute number of packets in next burst */
  131. if(rng_ == 0){
  132. rem_ = int(Random::pareto(p1_, shape_) + .5);
  133. }
  134. else{
  135. // Added by Debojyoti Dutta 13th October 2000
  136. rem_ = int(rng_->pareto(p1_, shape_) + .5);
  137. }
  138. /* make sure we got at least 1 */
  139. if (rem_ == 0)
  140. rem_ = 1;
  141. /* start of an idle period, compute idle time */
  142. if(rng_ == 0){
  143. t += Random::pareto(p2_, shape_);
  144. }
  145. else{
  146. // Added by Debojyoti Dutta 13th October 2000
  147. t += rng_->pareto(p2_, shape_);
  148. }
  149. on_ = 0;
  150. }
  151. rem_--;
  152. size = size_;
  153. return(t);
  154. }