# Patent

## Title
Calculations for admission control

## Abstract

The present invention relates to on-line admission-control decisions. Specifically, the invention concerns general delay bounds for both deterministic and statistical cases for Differentiated Services (DiffServ) networks. More specifically, a detailed method of calculation in each case is followed by simpler methods of calculation that are more appropriate for on-line admission-control decisions. Relatively involved occupancy bound calculations for various service classes take place only at the time of network configuration or reconfiguration. At the time of admission control only simple occupancy threshold compliance calculations need to be performed. Concrete illustrations are provided for deriving bounds for the EF and AF classes provided by DiffServ. These results are applicable to both layer-3 networks that support DiffServ and layer-2 networks that support the more restricted class of service functions.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application claims the benefit of the filing date of U.S. Provisional Patent Application No. 60/650,272, filed Feb. 4, 2005, the disclosure of which is hereby incorporated herein by reference.


GOVERNMENT LICENSE RIGHTS
This invention was made with Government support under contract NBCHC030132 awarded by DARPA. The Government may have certain rights in this invention.


FIELD OF THE INVENTION
The present invention relates to on-line admission-control decisions. Specifically, the invention concerns general delay bounds for both deterministic and statistical cases for Differentiated Services (DiffServ) networks. More specifically, a detailed method of calculation in each case is followed by simpler methods of calculation that are more appropriate for on-line admission-control decisions. Relatively involved occupancy bound calculations for various service classes take place only at the time of network configuration or reconfiguration. At the time of admission control only simple occupancy threshold compliance calculations need to be performed.
BACKGROUND OF THE INVENTION
The seminal work of Parekh and Gallager, entitled “A Generalized Processor Sharing Approach to Flow Control in Integrated Services Networks, The multiple node case”, EEEE/ACM Trans. Networking, vol. 2-2, pages 137-150, April 1994, has shown that it is possible to provide Quality of service (QoS) guarantees in an Internet Protocol (IP) network. Their techniques require certain fundamental additions to the networking infrastructure, including (1) admission control, (2) marking of admitted traffic so that the network routers and switches can associate it with an admitted flow, (3) policing of admitted traffic to make sure that it does not exceed its promised demand on network resources, and (4) managing and queuing traffic from admitted flows in certain ways in the core of the network to ensure that the service promised to admitted flows is not compromised by congestion in other portions of the network.
The initial QoS solution for IP networks, Integrated Services (IntServ), was directly inspired by the work of Parekh and Gallager, supra. That approach requires that network elements provide separate queuing for each admitted flow and entails the maintenance of per-flow state in each network element. IntServ approaches also tend to make admission decisions locally, rather than based on a network-wide policy.
Differentiated Services (DiffServ) is an approach to providing QoS to aggregate flows. It may provide satisfactory QoS without needing per-flow state in each network element. Additionally, in conjunction with a Bandwidth Broker for admission control, the DiffServ-based approach can support admission control that is driven by network-wide policies. Our delay bound techniques are readily applicable to layer-2 networks that support class of service (CoS), as CoS and DiffServ provide similar differentiated aggregate flow traffic treatment.
DiffServ uses aggregate flows in the sense that multiple flows share the same queue in the network core. Additional analysis is needed to bound the impact of aggregation of flows. Providing-end-to-end QoS requires the presence of individual network elements that have sufficiently strong QoS behavior, together with an analysis of network-wide accumulation of traffic bursts.
The analysis for aggregated flows relies on approaches where it is known that in the absence of bounds on the total number of network elements and on the maximum path length, it is possible to construct pathological scenarios in which no delay bound can be met. Fortunately, there are many situations where bounded delay can be achieved. These include networks where there is a known bound on the maximum number of hops taken by any flow. Many important networks, such as enterprise networks or networks confined to a ship or other vehicle are of bounded size and hence have a determinable upper bound on the number of hops taken by any flow.
Within a bounded network, there may be a small fraction of traffic that is so time-fragile that deterministic bounds on delay are important, a much larger fraction of less fragile, but still time-sensitive traffic for which approximate, probabilistic QoS is sufficient, and an even larger fraction for which best-effort delivery is completely satisfactory. To support this scenario of differentiated service requirements, without resorting to separate queues for individual flows, the network core would treat packets according to an aggregate classification. To respect the requirements of a more sensitive aggregate, the treatment could be one based on a strict priority, or it could be based on weighted sharing discipline. In any case, admissions need to be controlled appropriately while accounting correctly for the treatment of the aggregates.
Motivated by this scenario, the admission calculations in the present invention address this need within a general framework of any number of priority classes and, within each priority class, any number of weighted-share sub-classes.
SUMMARY OF THE INVENTION
Each session is imagined as beginning with an admission request that provides enough information to determine the session's service requirements and some features of the session's traffic, such as policing or shaping parameters (σ, ρ) bounding bursts and rates, with the meaning that the amount of that session's traffic entering the network during any interval of duration t will never exceed σ+ρt. In general, admission decisions are combined with routing decisions, but the following description will be phrased as if an admission decision is to be determined for a given session over a given route.
For the most efficient use of the network's resources that serve the traffic directly, the admission decision should be based on predictions of service quality that are as accurate as possible. Yet, because of the need not to delay these decisions too greatly, accuracy must be sacrificed, and the corresponding calculations are basically split between relatively complicated off-line determinations of suitable thresholds and, during the admission decisions themselves, relatively simple on-line comparisons of these thresholds with elementary combinations of traffic descriptions. Naturally, the off-line calculations must ensure that the thresholds are chosen so that service quality is assured as long as the admission decisions respect the thresholds. Of course, in this case, the loss of accuracy, that is created by the introduction of these thresholds, takes the form of loss of network efficiency: more admission requests are denied.
Different sorts of calculations arise for dealing with deterministic worst-case service requirements and for dealing with more statistical requirements. Naturally, the deterministic version concerns a guarantee about actual delays (GAD). To avoid impractical complexities, however, the statistical version is neither a guarantee about actual statistics (GAS), which would just be a more complicated worst-case guarantee, nor a guarantee about actual statistical expectations (GASE), which would require complex calculations of the implications of multitudes of cases beyond just the worst case. Instead, the statistical version described below is merely a guarantee about certain estimates of statistical expectations (GESE), where the estimates do include some pessimistic features in the conservative direction to take serious account of the traffic fluctuations that are possible.
The primary advance made by the present invention both for GAD and GESE is that the novel formulation is set in a more general context supporting any number of priority classes and within each priority class any number of weighted-fair queuing service classes. The scheme explicitly takes into account all the individual traffic streams that are presently in progress from all classes together with the stream whose admission is being considered. The main technique used in an earlier work to increase utilization bounds is one of dividing a class into several subclasses, each subclass supporting traffic with different deadline values in place of a single more stringent deadline for the entire class. However, implementation of multiple priorities in a single class is not readily realizable in commercial routers and switches.
Calculations according to the present invention for GESEs may be compared with earlier work. In comparison with those earlier works that come closest to ours, our calculations take more seriously the fluctuations that can arise while conforming to traffic descriptors. Thus, to deal with the possible fluctuations in the arrivals over a particular interval from independent streams of an aggregate, the present analysis relies on exponential moments for certain on-off distributions for the individual streams, where the on-off distribution for a given stream is the one that maximizes the exponential moments, among all distributions that respect the given bounds on the peak and the mean for the given stream in the given interval. In contrast, the prior art analysis, that comes closest to ours, relies on Gaussian approximations, where, for each stream, the corresponding variance is taken from the distribution that maximizes the entropy, among the distributions matching the given peak and mean. First, of course, the use of Gaussian approximations is more hopeful than conservative. Yet, a more disquieting concern is that these maximum-entropy variances are all smaller than the variances of other distributions that are consistent with the given peaks and means. In more detail, the maximum-entropy distribution is continuous, so it is not an on-off distribution, and the on-off distribution is the one that maximizes the variance. In this sense, the prior analysis is less serious about the fluctuations that can possibly arise.



FIG. 1 is a flow diagram illustrating operations performed at a device in accordance with a representative embodiment.
FIG. 2 is a flow diagram illustrating operations performed at a device in accordance with a representative embodiment.



DETAILED DESCRIPTION
In the following description, a network is imagined as serving sessions that can be partitioned into a few priority classes, which, in turn, can be partitioned into a few weighted fair queuing (WFQ) sub-classes. For instance, the priorities might be high, medium, and low, with high priority intended for real-time applications that could require a deterministic guarantee on delay, with low priority being basically best-effort service that requires no guarantee at all, and with medium priority intended for non-real-time applications that still require some kind of guarantee, such as a statistical guarantee. In this instance, the high-priority packets would have absolute, though non-pre-emptive, priority over the packets of the lower-priority classes, while the medium-priority packets would have absolute priority over low-priority packets. Moreover, continuing with this instance, the medium-priority applications are imagined as being partitioned into WFQ sub-classes or types, with each type getting a guaranteed fraction of the medium-priority effort of any congested resource, typically a link. More generally and abstractly, let L be the set of links in a network, where, for this description, the links are imagined as being the only network resources that could ever become congested and create queuing delays. For each link lεL, let C(l) be the capacity or service rate of the link. In addition, for the link l, it is convenient to introduce the notion of an effective propagation delay for the various packets traversing the link l, that represents the time between a packet ending its transmission onto the source end of link l and the same packet arriving at the queue for the next link on its route. So this time can vary from packet to packet, as the different packets deal with quirks in operating-system latencies at the next node. The purpose of including these latencies into the effective propagation delays is because the variations in these latencies are practically independent of the traffic, and they turn out to have the same effect on the analysis. Accordingly, let Δ(l) be the minimum possible value for the effective propagation delay of link l, and let J(l), the “jitter”, be the difference between the maximum and minimum possible values for the effective propagation delay of link l. Incidentally, the importance of distinguishing the portion J(l) of the bound Δ(l)+J(l) on the effective propagation time of link l is that J(l), unlike Δ(l), can contribute to increasing the sizes of bursts in streams traversing l. In addition, let M be the maximum size of a packet.
Let K be the set of priority classes. This set of priority classes is totally ordered, with each packet of one class having absolute service-order priority over packets of lower classes. For each class kεK, let Tk be the set of service types within class k. Reflecting the fact that a given service type belongs to just one priority class, the sets Tk are disjoint. Within class k, the packets of the different types in Tk are served according to a WFQ allocation of the class-k effort at any link. For each link lεL, for each type τεTk, let w(τ, l) be the weight that type τ gets at link l, so that the same type can get different weights at different links (which makes sense when different links could be handling different mixtures of the same traffic types). Here, the weights are taken to be normalized, so that




  
    
      ∑
      
        τ
        ∈
        
          T
          k
        
      
    
    ⁢
    
      w
      ⁡
      
        (
        
          τ
          ,
          l
        
        )
      
    
  
  =
  1




for any kεK and any lεL. For each type τεTk, let Sτ be the set of traffic streams of type τ that are imagined as being in progress for the purpose of an admission-control decision (so these streams include the session requesting admission, as well as the sessions that are really already in progress). Reflecting the fact that a given stream belongs to just one service type, the sets Sτ are disjoint. For each stream sεSτ, let H(s) be the number of hops in the route of s, and for any j with 0≦j<H(s), let λ(j,s) be the j-th link on the route of s. Also, let σs and ρs be the burst and rate shaping parameters that bound the traffic of stream s, at its entry into the network, i.e. as it arrives to queue at link λ(0,s).

Given this information about the network's configuration and traffic, the task for the calculations is to determine whether each session's service requirement is met.
The following are two calculation schemes that could be used, in principle, by an admission-control procedure to support a GAD (guarantee about actual delay). The first scheme is closely related to earlier work, though the earlier work deals with an abstract rate function whose identification for a service type at a low priority class implicitly requires some preliminary calculations for the performance of higher priority classes. Here, however, the relevant calculations for class-dependent performance are given explicitly, except that no reference is ever actually made to any abstract rate function. In any case, one of the calculation schemes being presented here takes explicit account of all the individual traffic streams that are presently in progress, together with the stream whose admission is being considered, though the complexity of the scheme renders it impractical for on-line admission control. A second scheme is obtained from the first by the introduction of thresholds, such as occupancy thresholds, where the idea is that most of the calculations are to be done off-line to identify the thresholds, so that only simple comparisons with thresholds need to be done on-line during admission control.
The next step is to try to identify several kinds of suitable bounds. For any stream s, for any j<H(s), one of the desired quantities is a value σs,j that can be used as a bound on the burstiness of stream s as it arrives at the queue for link λ(j,s). Naturally, for j=0, σs,0=σs is a suitable choice, but for larger values of j, some account must be taken for the possibility of stream s having acquired additional burstiness in traversing the hops leading up to hop j from the point at which stream s entered the network. Next, for any type τ of stream, and for any link lεL, the other desired quantity is a value Dτ(l) that can be used as a bound on the queuing delay that a type-τ packet could encounter at link l. Intuitively, one can understand how potential delays along the early hops taken by a stream can increase the burstiness of the stream as it arrives at later hops, and one can understand how the burstiness of the streams arriving at a link can increase the queuing delays at that link. Next is formalization of these intuitions into relations for governing suitable choices for the quantities σs,j and Dτ(l).
The easier intuition to formalize is the one relating to a stream's burstiness increasing as it traverses more hops. Thus, consider a stream s of type τ, and consider some hop j<H(s). Let l=λ(j,s), the link being traversed by stream s at hop j. Suppose that σs,j bounds the burstiness of stream s as it enters the queue for link l, and suppose that Dτ(l) is a bound on the queuing delay that a type-τ packet can encounter at link l. To see how bursty stream s can be as it leaves link l, suppose that one packet of stream s is delayed in queue by the maximal amount Dτ(l) and also encounters the maximal jitter J(l), while a later packet encounters neither queuing nor jitter delay at l. Comparing the interval separating the packets upon entering the queue with the interval after exiting the influence of the link, one finds that the latter interval is shorter than the former by Dτ(l)+J(l), which translates into a potential increase in burstiness of ρs(Dτ(l)+J(l)).
Thus,

σs,j+1=σs,j+ρs(Dτ(λ(j,s))+J(λ(j,s)))   (1)

is a suitable formalization of the intuition of how variable delays can increase burstiness.

The intuition of how bursts contribute to delays is much more difficult to formalize in the case of multiple types and classes of traffic, because of various complications. One complication is that a packet at one priority class can be delayed by later arrivals at higher priorities. The delay can even be from later arrivals at the same priority, if these arrivals belong to competing WFQ types. Finally, the fact, that traffic comes in discrete packets, rather than in continuous fluid, adds two more complications, one complication being in the merely approximate success of WFQ in isolating delays to one service type from the effects of traffic bursts in other service types, and the other complication being that a high-priority packet can be delayed by a low-priority packet that begins its transmission just before the high-priority packet arrives. Eventually, a formula, that expresses how bursts contribute to delays, must account for all of these complications.
To begin to express these complications quantitatively, it helps first to identify aggregate rate and burst bounds, based on bounds for the individual streams. Thus, for any service type τ and any link l, let



  
    
      
        
          ρ
          τ
        
        ⁡
        
          (
          l
          )
        
      
      =
      
        
          ∑
          
            s
            ∈
            
              S
              τ
            
          
        
        ⁢
        
          
            ∑
            
              
                
                  0
                  ≤
                  j
                  <
                  
                    H
                    ⁡
                    
                      (
                      s
                      )
                    
                  
                
                :
                
                  λ
                  ⁡
                  
                    (
                    
                      j
                      ,
                      s
                    
                    )
                  
                
              
              =
              l
            
          
          ⁢
          
            ρ
            s
          
        
      
    
    ,
    
      

    
    ⁢
    and
  




  
    
      
        σ
        τ
      
      ⁡
      
        (
        l
        )
      
    
    =
    
      
        ∑
        
          s
          ∈
          
            S
            τ
          
        
      
      ⁢
      
        
          ∑
          
            
              
                0
                ≤
                j
                <
                
                  H
                  ⁡
                  
                    (
                    s
                    )
                  
                
              
              :
              
                λ
                ⁡
                
                  (
                  
                    j
                    ,
                    s
                  
                  )
                
              
            
            =
            l
          
        
        ⁢
        
          
            σ
            
              s
              ,
              j
            
          
          .
        
      
    
  



Then, for any priority class k, let



  
    
      
        
          ρ
          
            >
            k
          
        
        ⁡
        
          (
          l
          )
        
      
      =
      
        
          ∑
          
            
              k
              ′
            
            ∈
            
              K
              :
              
                
                  k
                  ′
                
                >
                k
              
            
          
        
        ⁢
        
          
            ∑
            
              
                τ
                ′
              
              ∈
              
                T
                
                  κ
                  ′
                
              
            
          
          ⁢
          
            
              ρ
              
                τ
                ′
              
            
            ⁡
            
              (
              l
              )
            
          
        
      
    
    ,
    
      

    
    ⁢
    and
  




  
    
      
        σ
        
          >
          k
        
      
      ⁡
      
        (
        l
        )
      
    
    =
    
      
        ∑
        
          
            k
            ′
          
          ∈
          
            K
            :
            
              
                k
                ′
              
              >
              k
            
          
        
      
      ⁢
      
        
          ∑
          
            
              τ
              ′
            
            ∈
            
              T
              
                k
                ′
              
            
          
        
        ⁢
        
          
            
              σ
              
                τ
                ′
              
            
            ⁡
            
              (
              l
              )
            
          
          .
        
      
    
  



Then the service rate available to priority class k is Ck(l)=C(l)−ρ>k(l), which is strictly greater than zero, if admission control is doing its job. Of course, only the amount w(τ,l)Ck(l) is reliably available to a particular service type τεTk.

For the subtler complications, it helps to identify corresponding indicators. First, for any priority class k, let l′k=I (k is not the lowest priority), where, for any clause c, I(c) is the indicator of clause c, which is 1 if c is true and 0 if c is false. This indicator is related to the subtle complication that even high-priority packets can experience up to a packet's worth of delay from the influence of a low-priority packet at the start of a high-priority busy period, if the higher-priority packets cannot preempt the lower-priority packet, provided that there is a lower priority class. Next, for any priority class k and any service type τεTk, let {tilde over (l)}τ=I({τ}≠Tk). The corresponding subtle complication here is that WFQ, as a discrete approximation to an idealized partition of continuous service, can miss its target allocation by up to a packet's worth, provided that there is more than one service type in the priority class. Combining all the considerations above, one finds that the intuition of how suitable choices for delay bounds can be derived from suitable burstiness bounds can be expressed as




  
    
      
        
          
            D
            τ
          
          ⁡
          
            (
            l
            )
          
        
        =
        
          
            
              
                M
                ⁢
                
                  
                
                ⁢
                
                  ι
                  k
                  ′
                
              
              +
              
                
                  σ
                  
                    >
                    k
                  
                
                ⁡
                
                  (
                  l
                  )
                
              
            
            
              
                C
                k
              
              ⁡
              
                (
                l
                )
              
            
          
          +
          
            
              
                
                  M
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      ι
                      ~
                    
                    τ
                  
                
                +
                
                  
                    σ
                    τ
                  
                  ⁡
                  
                    (
                    l
                    )
                  
                
              
              
                
                  w
                  ⁡
                  
                    (
                    
                      τ
                      ,
                      l
                    
                    )
                  
                
                ⁢
                
                  
                    C
                    k
                  
                  ⁡
                  
                    (
                    l
                    )
                  
                
              
            
            .
          
        
      
    
    
      
        (
        2
        )
      
    
  




The intuitive relations between bursts and delays have now been formalized in the relations (1) and (2). The problem with these relations, however, is the circularity that they typically generate for the dependencies among the quantities of interest. In other words, there is typically no simple order in which one can first identify the values for some “earlier” quantities and then use the relations successively to identify values for the “later” quantities one step at a time. Instead, many of the identifications have to be made simultaneously.

The following description uses the relations identified previously to construct a scheme of calculations in terms of individual streams. Unfortunately, the scheme seems too complicated for on-line admission control, which is why the following description adds a modification in terms of thresholds, so that on-line calculations can be simple. The task is to identify choices for the quantities σs,j and Dτ(l) that solve simultaneously the relations of the previous paragraphs, primarily (1) and (2).
To organize this kind of simultaneous solution, it helps first to use the relations (1) to express all the burst bounds σs,j explicitly in terms of the delay bounds Dτ(l). Thus, recalling that σs,0=σs, one finds that, for any service type τ, for any stream sεSτ, for any j<H(s),




  
    
      
        
          σ
          
            s
            ,
            j
          
        
        =
        
          
            σ
            s
          
          +
          
            
              ρ
              s
            
            ⁢
            
              
                ∑
                
                  0
                  ≤
                  
                    j
                    ′
                  
                  <
                  j
                
              
              ⁢
              
                
                  [
                  
                    
                      
                        D
                        τ
                      
                      ⁡
                      
                        (
                        
                          λ
                          ⁡
                          
                            (
                            
                              
                                j
                                ′
                              
                              ,
                              s
                            
                            )
                          
                        
                        )
                      
                    
                    +
                    
                      J
                      ⁡
                      
                        (
                        
                          λ
                          ⁡
                          
                            (
                            
                              
                                j
                                ′
                              
                              ,
                              s
                            
                            )
                          
                        
                        )
                      
                    
                  
                  ]
                
                .
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




The idea now is to use (3) to eliminate all the burst bounds σs,j in favor of the delay bounds Dτ(l). As a result, (2) is modified to express Dτ(l) in terms of other delay bounds Dτ(l′) for the same service type τ at other links l′, other delay bounds for higher priority classes, and still other influences that are independent of any delay bounds. So the delay bounds can indeed be identified in a manner that is at least partially ordered: the delay bounds for higher priorities can be identified before any delay bounds for lower priorities are considered. Moreover, within a single priority class, the delay bounds for the different service types can be identified independently.
In accordance with this procedure, the resulting modification of (2) takes the form




  
    
      
        
          
            
              D
              τ
            
            ⁡
            
              (
              l
              )
            
          
          =
          
            
              
                b
                τ
              
              ⁡
              
                (
                l
                )
              
            
            +
            
              
                ∑
                
                  
                    l
                    '
                  
                  ∈
                  L
                
              
              ⁢
              
                
              
              ⁢
              
                
                  
                    Γ
                    τ
                  
                  ⁡
                  
                    (
                    
                      l
                      ,
                      
                        l
                        '
                      
                    
                    )
                  
                
                ⁢
                
                  
                    D
                    τ
                  
                  ⁡
                  
                    (
                    
                      l
                      '
                    
                    )
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        4
        )
      
    
  




where the quantities bτ(l) capture the influence of intrinsic bursts and jitters, as well as delay bounds for higher priorities, while Γτ(l,l′) captures the influence of Dτ(l′) on Dτ(l). Explicitly, these quantities are given by




  
    
      
        
          b
          τ
        
        ⁡
        
          (
          l
          )
        
      
      =
      
        
          
            
              M
              ⁢
              
                
              
              ⁢
              
                ι
                k
                ′
              
            
            +
            
              
                σ
                
                  >
                  k
                
              
              ⁡
              
                (
                l
                )
              
            
          
          
            
              C
              k
            
            ⁡
            
              (
              l
              )
            
          
        
        +
        
          
            
              M
              ⁢
              
                
              
              ⁢
              
                
                  ι
                  ~
                
                τ
              
            
            +
            
              
                σ
                τ
                J
              
              ⁡
              
                (
                l
                )
              
            
          
          
            
              w
              ⁡
              
                (
                
                  τ
                  ,
                  l
                
                )
              
            
            ⁢
            
              
                C
                k
              
              ⁡
              
                (
                l
                )
              
            
          
        
      
    
    ,
    
      

    
    ⁢
    where
  




  
    
      
        σ
        τ
        J
      
      =
      
        
          
            ∑
            
              s
              ∈
              
                S
                τ
              
            
          
          ⁢
          
            
              ∑
              
                
                  
                    0
                    ≤
                    j
                    <
                    
                      H
                      ⁡
                      
                        (
                        s
                        )
                      
                    
                  
                  :
                  
                    λ
                    ⁡
                    
                      (
                      
                        j
                        ,
                        s
                      
                      )
                    
                  
                
                =
                l
              
            
            ⁢
            
              σ
              s
            
          
        
        +
        
          
            ∑
            
              s
              ∈
              
                S
                τ
              
            
          
          ⁢
          
            
              ∑
              
                
                  
                    0
                    ≤
                    
                      j
                      ′
                    
                    <
                    j
                    <
                    
                      H
                      ⁡
                      
                        (
                        s
                        )
                      
                    
                  
                  :
                  
                    λ
                    ⁡
                    
                      (
                      
                        j
                        ,
                        s
                      
                      )
                    
                  
                
                =
                l
              
            
            ⁢
            
              
                ρ
                s
              
              ⁢
              
                J
                ⁡
                
                  (
                  
                    λ
                    ⁡
                    
                      (
                      
                        
                          j
                          ′
                        
                        ,
                        s
                      
                      )
                    
                  
                  )
                
              
            
          
        
      
    
    ,
    
      

    
    ⁢
    and
  




  
    
      
        Γ
        τ
      
      ⁡
      
        (
        
          l
          ,
          
            l
            ′
          
        
        )
      
    
    =
    
      
        
          
            ∑
            
              s
              ∈
              
                S
                τ
              
            
          
          ⁢
          
            
              ∑
              
                0
                ≤
                
                  j
                  ′
                
                <
                j
                <
                
                  H
                  ⁡
                  
                    (
                    s
                    )
                  
                
              
            
            ⁢
            
              
                ρ
                s
              
              ⁢
              
                I
                ⁡
                
                  (
                  
                    
                      λ
                      ⁡
                      
                        (
                        
                          
                            j
                            ′
                          
                          ,
                          s
                        
                        )
                      
                    
                    =
                    
                      
                        
                          
                            l
                            ′
                          
                          &
                        
                        ⁢
                        
                          
                        
                        ⁢
                        
                          λ
                          ⁡
                          
                            (
                            
                              j
                              ,
                              s
                            
                            )
                          
                        
                      
                      =
                      l
                    
                  
                  )
                
              
            
          
        
        
          
            w
            ⁡
            
              (
              
                τ
                ,
                l
              
              )
            
          
          ⁢
          
            
              C
              k
            
            ⁡
            
              (
              l
              )
            
          
        
      
      .
    
  



The formulation (4) can be regarded as the matrix equation

Dτ=bτ+ΓτDτ,  (5)

if Dτ and bτ are regarded as column vectors with index set L, and if Γτ is regarded as a square matrix with index set L×L. Hence, if the matrix IεΓτ is invertible, where I here is the L×L identity matrix, then the solution to (5) is Dτ=(I−Γτ)−1bτ. Yet, invertibility alone is insufficient for the solution to be meaningful, because, with all the entries of Γτ being non-negative, the equation (5) can be used to establish finite, positive bounds on delays only if all the entries of the matrix (I−Γτ)−1 are also non-negative, which turns out to be equivalent to the condition that the power series





  
    ∑
    
      n
      =
      0
    
    ∞
  
  ⁢
  
    Γ
    τ
    n
  




converge to (I−Γτ)−1. That is, meaningful finite positive bounds cannot be established in this way unless the relevant occupancies entering the terms of Γτ are sufficiently small (though the precise version of the condition is, as was just noted, the convergence of the power series).

Assuming that the desired bounds Dτ(l) on individual link delays are produced, they are easily combined into bounds on end-to-end delays, at least in a conceptual sense. Explicitly, for any individual session S of type τ within priority class k, one can guarantee that the maximum end-to-end delay will not exceed




  
    
      D
      s
    
    =
    
      
        ∑
        
          0
          ≤
          j
          ≤
          
            H
            ⁡
            
              (
              s
              )
            
          
        
      
      ⁢
      
        [
        
          
            
              D
              τ
            
            ⁡
            
              (
              
                λ
                ⁡
                
                  (
                  
                    j
                    ,
                    s
                  
                  )
                
              
              )
            
          
          +
          
            Δ
            ⁡
            
              (
              
                λ
                ⁡
                
                  (
                  
                    j
                    ,
                    s
                  
                  )
                
              
              )
            
          
          +
          
            J
            ⁡
            
              (
              
                λ
                ⁡
                
                  (
                  
                    j
                    ,
                    s
                  
                  )
                
              
              )
            
          
        
        ]
      
    
  
  ,




and the maximum end-to-end delay variation will not exceed





  
    D
    s
    ′
  
  =
  
    
      ∑
      
        0
        ≤
        j
        <
        
          H
          ⁡
          
            (
            s
            )
          
        
      
    
    ⁢
    
      
        [
        
          
            
              D
              τ
            
            ⁡
            
              (
              
                λ
                ⁡
                
                  (
                  
                    j
                    ,
                    s
                  
                  )
                
              
              )
            
          
          +
          
            J
            ⁡
            
              (
              
                λ
                ⁡
                
                  (
                  
                    j
                    ,
                    s
                  
                  )
                
              
              )
            
          
        
        ]
      
      .
    
  




Presumably, the criteria for satisfactory service to a stream s involve tolerances δs and δ′s for maximal delays and delay variations. Accordingly, as long as the inequalities Ds≦δs and D′s≦δ′s are all satisfied, the calculations justify admitting all the streams.
As an alternative second scheme to the first scheme that was just described, thresholds are introduced to allow for simple on-line comparisons in making admission decisions, where the idea is that the complicated off-line calculations identifying the thresholds would be done ahead of time.
Primarily, the thresholds considered here are occupancy thresholds, but the scheme here also relies on a bound h on all the hop counts H(s) of all the streams s that will ever be admitted, as well as bounds Tτ on the burst times for the streams of type τ, with the meaning that, for any stream s of type τ that is ever admitted, σs≦ρsTτ. In a strict sense, the weights w(τ,l) are also thresholds that should be tuned to optimize the admission of anticipated traffic, but the emphasis here is on the new thresholds that do not appear in the calculations of the first scheme that was described above. Accordingly, the emphasis is on the occupancy thresholds αk where αk can be regarded as the occupancy allocated to priority class k, but its more precise meaning is that, during an on-line admission decision for a stream s of type τ of priority class k, the stream s would be rejected if, for any link l on the route of s, the admission of s would cause ρτ(l) to exceed αk w(τ,l)C(l).
An implication of basing admission decisions on thresholds for aggregates rather than on the individual streams is that, once one has to deal with a situation for a given link in which the streams routed through that link could have been routed almost arbitrarily through other links in the network, it ceases to make much sense to try to identify distinct delay bounds Dτ(l) for different links l. Instead, for a given service type τ, it is appropriate to identify just a single delay bound Dτ that applies to all the links. An implication of this change is that the modified version of the matrix equation (5) is a scalar equation.
The specific quantities, that enter the calculation of the delay bound (Dτ) for a service type τ at priority class k, are the following. First, for any priority class k′ higher than k, let Dk′=maxτ′εTRDτ′, and Tk′=maxτ′εTRTτ′.
Next, let




  
    α
    
      >
      k
    
  
  =
  
    
      ∑
      
        
          k
          ′
        
        ∈
        K
      
    
    ⁢
    
      
        α
        
          k
          ′
        
      
      .
    
  




Next, introduce the “averages”




  
    
      D
      
        >
        k
      
    
    =
    
      
        α
        
          >
          k
        
        
          -
          1
        
      
      ⁢
      
        
          ∑
          
            
              k
              ′
            
            ∈
            
              
            
            ⁢
            
              K
              :
              
                
                  k
                  ′
                
                >
                k
              
            
          
        
        ⁢
        
          
            α
            
              k
              ′
            
          
          ⁢
          
            D
            
              k
              ′
            
          
        
      
    
  




  and




  
    
      T
      
        >
        k
      
    
    =
    
      
        α
        
          >
          k
        
        
          -
          1
        
      
      ⁢
      
        
          ∑
          
            
              k
              ′
            
            ∈
            
              
            
            ⁢
            
              K
              :
              
                
                  k
                  ′
                
                >
                k
              
            
          
        
        ⁢
        
          
            α
            
              k
              ′
            
          
          ⁢
          
            
              T
              
                k
                ′
              
            
            .
          
        
      
    
  



Also, let J=maxlεLJ(l). Finally, for each link lεL, let




  
    
      
        b
        ~
      
      τ
    
    ⁡
    
      (
      l
      )
    
  
  =
  
  ⁢
  
    
      
        M
        
          
            C
            ⁡
            
              (
              l
              )
            
          
          ⁢
          
            (
            
              1
              -
              
                α
                
                  >
                  k
                
              
            
            )
          
        
      
      ⁡
      
        [
        
          
            ι
            k
            ′
          
          +
          
            
              
                ι
                ~
              
              τ
            
            
              w
              ⁡
              
                (
                
                  τ
                  ,
                  l
                
                )
              
            
          
        
        ]
      
    
    +
    
      
        
          α
          
            >
            k
          
        
        
          1
          -
          
            α
            
              >
              k
            
          
        
      
      [
      
        
          
            T
            
              >
              k
            
          
          +
          
            
              (
              
                h
                -
                1
              
              )
            
            ⁢
            
              (
              
                
                  D
                  
                    
                      
                    
                    ⁢
                    
                      >
                      
                        
                      
                      ⁢
                      k
                    
                  
                
                +
                
                  J
                  ~
                
              
              )
            
            ]
          
          +
          
            
              
                
                  
                
                ⁢
                
                  α
                  
                    
                      
                    
                    ⁢
                    k
                  
                
              
              
                
                  
                
                ⁢
                
                  1
                  ⁢
                  
                    
                  
                  -
                  
                    
                  
                  ⁢
                  
                    α
                    
                      
                        
                      
                      ⁢
                      
                        >
                        
                          
                        
                        ⁢
                        k
                      
                    
                  
                
              
            
            ⁡
            
              [
              
                
                  T
                  
                    
                      
                    
                    ⁢
                    τ
                  
                
                +
                
                  
                    (
                    
                      h
                      -
                      1
                    
                    )
                  
                  ⁢
                  
                    J
                    ~
                  
                
              
              ]
            
          
        
        ,
      
    
  




and set





  
    
      b
      ^
    
    τ
  
  =
  
    
      max
      
        l
        ∈
        L
      
    
    ⁢
    
      
    
    ⁢
    
      
        
          
            b
            ~
          
          τ
        
        ⁡
        
          (
          l
          )
        
      
      .
    
  




These quantities now allow for an easy determination of the delay bound Dτ, because the modified version of the matrix equation (5) is the scalar equation




  
    D
    τ
  
  =
  
    
      
        b
        ^
      
      τ
    
    +
    
      
        
          
            α
            k
          
          ⁡
          
            (
            
              h
              -
              1
            
            )
          
        
        
          1
          -
          
            α
            
              >
              k
            
          
        
      
      ⁢
      
        
          D
          τ
        
        .
      
    
  




Thus, as long as αk(h−1)/(1−α>k)<1, then the individual link delays for type-τ streams are all bounded by




  
    D
    τ
  
  =
  
    
      
        [
        
          1
          -
          
            
              
                α
                k
              
              ⁡
              
                (
                
                  h
                  -
                  1
                
                )
              
            
            
              1
              -
              
                α
                
                  >
                  k
                
              
            
          
        
        ]
      
      
        -
        1
      
    
    ⁢
    
      
        
          b
          ^
        
        τ
      
      .
    
  




Bounds on the end-to-end delays are then identified rather routinely. First, let




  
    
      Δ
      ^
    
    =
    
      
        max
        s
      
      ⁢
      
        
          ∑
          
            0
            ≤
            j
            ≤
            
              H
              ⁡
              
                (
                s
                )
              
            
          
        
        ⁢
        
          Δ
          ⁡
          
            (
            
              λ
              ⁡
              
                (
                
                  j
                  ,
                  s
                
                )
              
            
            )
          
        
      
    
  
  ,




where an assumption being made here is that the longest route is the same for all service types. Then a bound on the longest end-to-end delay for all type-τ streams is given by {circumflex over (D)}τ={circumflex over (Δ)}+h(Dτ+ J), while a bound on the longest delay variation is {circumflex over (D)}′τ=h(Dτ+ J). Accordingly, if there are tolerances δτ and τ′τ on the maximal end-to-end delay and delay variation for type-τ streams, then the point of the off-line calculations is to select choices for the thresholds like αk to ensure that {circumflex over (D)}τ≦δτ and {circumflex over (D)}′τ≦δ′τ for all service types τ.

The following illustrations deal with 4 service types in total, distributed among 3 priority classes: a high-priority Expedited Forwarding (EF) type, two medium-priority Assured Forwarding (AF) types, AF0 and AF1, and a low-priority best-effort type that really has no delay target. The service rate of each link will be taken to be C(l)=1 Gb/s. The maximum packet size will be taken to be M=12,000 bits, corresponding to 1500 bytes. The bound on the number of hops in a route will be taken to be h=10, so h−1=9. The bound on operating-system jitters at each hop will be taken to be J=1 ms. The bound on the propagation delay along the longest route in the network will be taken to be {circumflex over (Δ)}=4 μs, corresponding to an image of a 300-meter ship networked with fiber having an index of refraction of 2, but with sloppy enough wiring and routing to double the propagations that might otherwise be obtained. The burst times of the various service types will be taken to be TEF=2 ms, TAF0=10 ms, and TAF1=50 ms. As for the WFQ weights distinguishing the AF types, they will be taken to be w(AF0,l)=0.6 and w(AF1,l)=0.4 at all links l. Finally, the target delays will be taken to be δEF=125 ms, δAF0=500 ms, and δAF1=1 second.
The computations for the high-priority EF type are rather direct. First, in terms of milliseconds, {circumflex over (b)}EF=11αEF+0.024. Hence,



  
    
      
        D
        EF
      
      =
      
        
          
            
              b
              ^
            
            EF
          
          
            1
            -
            
              9
              ⁢
              
                
              
              ⁢
              
                α
                EF
              
            
          
        
        =
        
          
            
              11
              ⁢
              
                
              
              ⁢
              
                α
                EF
              
            
            +
            0.024
          
          
            1
            -
            
              9
              ⁢
              
                
              
              ⁢
              
                α
                EF
              
            
          
        
      
    
    ,
    
      

    
    ⁢
    so
  




  
    
      
        D
        ^
      
      EF
    
    =
    
      
        
          Δ
          ^
        
        +
        
          h
          ⁡
          
            [
            
              
                D
                EF
              
              +
              
                J
                _
              
            
            ]
          
        
      
      =
      
        10.004
        +
        
          
            
              
                110
                ⁢
                
                  α
                  EF
                
              
              +
              0.24
            
            
              1
              -
              
                9
                ⁢
                
                  
                
                ⁢
                
                  α
                  EF
                
              
            
          
          .
        
      
    
  



With αEF=0.1, the result is {circumflex over (D)}EF=122.404 ms, which does not exceed the target δEF=125 ms.
The computations for the AF types are similar, except for dealing with the extra complication of dealing with the effects of the higher-priority EF type. In particular, one must account for




  
    
      1
      -
      
        α
        
          >
          AF
        
      
    
    =
    0.9
  
  ,
  
    
      
        a
        
          >
          AF
        
      
      
        1
        -
        
          α
          
            >
            AF
          
        
      
    
    =
    
      1
      9
    
  
  ,




and DEF=11.24 ms. First, for each of the two AF types τ,





  
    
      b
      ^
    
    τ
  
  =
  
    
      
        0.12
        9
      
      ⁢
      
        (
        
          1
          +
          
            1
            
              w
              ⁡
              
                (
                τ
                )
              
            
          
        
        )
      
    
    +
    12.24
    +
    
      2
      9
    
    +
    
      10
      ⁢
      
        
      
      ⁢
      
        
          
            α
            AF
          
          ⁡
          
            (
            
              
                
                  T
                  τ
                
                9
              
              +
              1
            
            )
          
        
        .
      
    
  




Next, Dτ={circumflex over (b)}τ/(1−10αAF). Finally, {circumflex over (D)}τ=10.004+10{circumflex over (b)}τ/(1−10αAF). With αAF=0.06, {circumflex over (D)}AF0≈354 ms, while {circumflex over (D)}AF1≈421 ms, which meet the service targets (and type AF1 even meets the target for AF0). Incidentally, notice that the choice of αAF=0.07 would violate the target for AF0.
Analogous with the previous discussion of the calculations for GADs there will now be a similar discussion of calculations for GESEs (guarantees about estimates of statistical expectations). The aim here is not to provide a guarantee on actual statistical expectations; instead, for the sake of practicality, the aim is merely to provide a guarantee for rough approximations to such statistical expectations. Before identifying the specific approximations entering the guarantee, however, the first step is to specify the service targets that will be considered. In particular, the calculations here will be developed on the conception that for a stream s of service type τ, most of the end-to-end delays are to be bounded by some target δτ, with some tolerance ετ for the fraction of packets with end-to-end delays above δτ, where ετ might be on the order of 10−3. Now, for each stream s, the maximum value, that can be attained by the effective end-to-end propagation delay, is




  
    
      Δ
      _
    
    s
  
  =
  
    
      ∑
      
        0
        ≤
        j
        <
        
          H
          ⁡
          
            (
            s
            )
          
        
      
    
    ⁢
    
      
        [
        
          
            Δ
            ⁡
            
              (
              
                λ
                ⁡
                
                  (
                  
                    j
                    ,
                    s
                  
                  )
                
              
              )
            
          
          +
          
            J
            ⁡
            
              (
              
                λ
                ⁡
                
                  (
                  
                    j
                    ,
                    s
                  
                  )
                
              
              )
            
          
        
        ]
      
      .
    
  




Thus, for a stream s of type τ, it makes sense to admit s only if Δs<δτ, and then only if one expects that the end-to-end queuing delay encountered by a packet of stream s to exceed {hacek over (δ)}s=δτ− Δs with a probability no greater than ετ. This formulation, of comparing ετ with an estimate of the probability for a queuing delay to exceed {hacek over (δ)}s, is basically the version of the service target to be supported by the calculations to be developed here.
The next step is to list the simplifying approximations to be made. As one example of such a simplifying approximation, the analysis here will assume that the queue sizes at different links are independent. Another simplification is to assume that the different streams are independent, which is at least more reasonable, though perhaps not entirely accurate, if different streams turn out to support different aspects of a single application e.g. voice and video channels of a single conversation. In addition, another simplifying approximation is to treat the traffic being generated by the different streams as forming, not just independent stochastic processes, but stationary stochastic processes. Yet another aspect of the simplifying approach to be used in this analysis, is to employ the philosophy of large-deviations analysis in paying attention only to the most likely ways of having bad things happen in the service of a stream, or, more formally, of estimating the probability of a union as being roughly equal to the maximum of the individual probabilities.
Another simplifying approximation, that will be made here, is to assume that the burstiness of a stream's traffic as it arrives at one link on its route is the same as the burstiness at other links along the route. Indeed, it is the tracing of the potential increase in the burstiness of a stream along the successive links of the route that creates most of the complexity in the worst-case analysis. So the analysis here will basically ignore this complexity.
The simplifications listed above are all in the optimistic direction (i.e. they tend to underestimate the service-degrading possibilities of the real network), so the next simplifications will be in the conservative direction. One of these simplifications is to replace probabilities by certain exponential moments, where the simplification occurs if the probability has the form P(S>s), where s=ΣiXi is the sum of independent random variables Xi whose distributions are relatively easy to manipulate. Formally, the relevant bound is P(S≧s)≦Eep(S−s) for any p≧0, which is valid because, whenever S≧s, 1≦ep(S−s), and whenever S<s, 0<ep(S−s). Here, the simplification is that P(S≧s) is typically very complicated to compute, whereas EepS is simply the product of the EepXi.
In the calculations for the GESE, Xi will often be the random amount of traffic being generated by a particular stream s during a particular time interval of duration t, so the knowledge available about Xi will include certain bounds on Xi, but will not include all the details of the distribution of Xi. Basically, if the burst-rate pair (σ, ρ) is known for the stream s, then the values of the random variable Xi are known to be confined to the interval [0,σ+ρt]. Of course, if Xi is taken from a stationary process that is governed by the pair (σ,ρ), then EXi≦ρt. Now, regardless of the value of p>0, of all the distributions of random variables X that are confined to some interval [0,B] with B>0, and that have a mean bounded by some prescribed bound mε[0,B], the distribution, that maximizes the exponential moment EepX, is the on-off distribution that assigns probability




  
    m
    B
  
  ⁢
  
    
  
  ⁢
  to
  ⁢
  
    
  
  ⁢
  B




and





  1
  -
  
    
      m
      B
    
    ⁢
    
      
    
    ⁢
    t
    ⁢
    
      
    
    ⁢
    o
    ⁢
    
      
    
    ⁢
    0.
  




Accordingly, in the calculations to be presented later, the exponential moment EepXi will be estimated by the quantity





  
    
      
        
          μ
          ⁡
          
            (
            
              
                p
                ;
                t
              
              ,
              σ
              ,
              ρ
            
            )
          
        
        =
        
          
            
              
                ρ
                ⁢
                
                  
                
                ⁢
                t
              
              
                σ
                +
                
                  ρ
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
              
            
            ⁢
            
              e
              
                p
                ⁡
                
                  (
                  
                    σ
                    +
                    
                      
                    
                    ⁢
                    
                      p
                      ⁢
                      
                        
                      
                      ⁢
                      t
                    
                  
                  )
                
              
            
          
          +
          
            σ
            
              σ
              +
              
                ρ
                ⁢
                
                  
                
                ⁢
                t
              
            
          
        
      
    
  
  
    
      
        =
        
          1
          +
          
            
              
                ρ
                ⁢
                
                  
                
                ⁢
                t
              
              
                σ
                +
                
                  
                
                ⁢
                
                  ρ
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
              
            
            ⁢
            
              
                (
                
                  
                    e
                    
                      p
                      ⁡
                      
                        (
                        
                          σ
                          +
                          
                            
                          
                          ⁢
                          
                            ρ
                            ⁢
                            
                              
                            
                            ⁢
                            t
                          
                        
                        )
                      
                    
                  
                  -
                  1
                
                )
              
              .
            
          
        
      
    
  




As for the specific considerations motivating specific calculations, recall first that the service target for a stream s of type τ is that the end-to-end queuing delay not exceed {hacek over (δ)}s with a probability greater than ετ. Evidently, the main task is to identify the probabilities for large queuing delays. Of course, the end-to-end queuing delay along a route is simply the sum of the individual queuing delays at each of the links along the route. Formally, let Dj be the random queuing delay at link λ(j,s). Then ΣjDj is the end-to-end queuing delay. Hence, if the end-to-end queuing delay exceeds {hacek over (δ)}s, then for some sequence of values djε[0,∞) with 0≦j<H(s), where Σjdj={hacek over (δ)}s, Dj must be at least as large as dj for each j. In view of the simplifying approximation procedure of treating the performance at different links as being independent, the corresponding probability will be approximated as ΠjP(Dj≧dj). Moreover, in view of the large-deviations procedure of approximating the probability of a union as the maximum of the individual probabilities, the procedure here will be to make the approximation




  
    
      P
      ⁡
      
        (
        
          
            
              ∑
              j
            
            ⁢
            
              D
              j
            
          
          >
          
            
              δ
              ⋓
            
            s
          
        
        )
      
    
    ≈
    
      
        sup
        
          
            
              d
              ∈
              
                
                  [
                  
                    0
                    ,
                    ∞
                  
                  )
                
                
                  H
                  ⁡
                  
                    (
                    s
                    )
                  
                
              
            
            :
            
              ∑
              
                
              
              ⁢
              d
            
          
          =
          
            
              δ
              ⋓
            
            s
          
        
      
      ⁢
      
        
          ∏
          
            0
            ≤
            j
            <
            
              H
              ⁡
              
                (
                s
                )
              
            
          
        
        ⁢
        
          
        
        ⁢
        
          P
          ⁡
          
            (
            
              
                D
                j
              
              ≥
              
                d
                j
              
            
            )
          
        
      
    
  
  ,




i.e. to pick the sequence d={dj} to maximize the expression ΠjP(Dj≧dj). In any case, most of the rest of the development of this approximation will be concerned with estimating the individual factors P(Dj≧dj).

In considering the queuing delay to a “test” packet of stream s of type T in priority class k at a link l=λ(j,s) on its route, let time 0 be the time at which the packet arrives at the queue for link l. Naturally, this test packet is delayed by all the type-τ packets in queue ahead of the test packet, as well as all the higher-priority packets that either are in queue ahead of the test packet or arrive later, but still before the test packet is served. As for the other WFQ types τ′ of class k that are competing with type τ, these arrivals can also delay the test packet even by arriving after the test packet, but only to the extent of the bandwidth allocated to these other types. Let tj be the age of the busy period, relative to classes k and higher, that is in progress at time 0. Let A>k be the random amount of higher-priority traffic that arrives at link l during the interval (−tj, dj]. If A>k≧C(l)(tj+dj), then the link serves higher-priority traffic exclusively during this interval, unless the amount A>k arrives primarily during the latter part of the interval (thereby allowing the link to serve some lower-priority traffic before the bulk of the amount A>k arrives). Within this relatively rare event of A>k≧C(l) (tj+dj), however, the most likely pattern of arrival is close to the pattern of a nearly constant arrival rate of A>k/(tj+dj), which does not allow any lower-priority traffic to be served. Thus, with this kind of large-deviations excuse, it is reasonable to presume that, whenever A>k≧C(l) (tj+dj), no low-priority traffic is served, and, in particular, the test packet is delayed more than dj.
A more likely way of delaying the test packet more than dj, however, is to share the burden of occupying the link. The effective service Sk=C(l)(tj+dj)−A>k available to priority class k during the interval (−tj, dj], specifically, would be positive, but this service would be occupied by other class-k competitors of the test packet. Let Aτ be the amount of arrivals of type τ during the interval (−tj,0]. Similarly, let Ak−τ be the amount of traffic of priority class k, but outside type τ, that arrives during the interval (−tj,dj]. So Ak−τ can be regarded as the “external” traffic of class k that competes with type τ and can contribute to the delay of the test packet, unless the WFQ weights mitigate the influence of this external traffic. Now, if Ak−τ≦(1−w(τ,l))Sk, then there is at least a decent chance that all of Ak−τ can contribute to the delay of the test packet (unless individual classes τ′ generate more than their share of the traffic). In a strict sense, one should consider all possible subsets of Tk−{τ}, the set of WFQ types competing with type τ, as the set of possible types whose every arrival contributes to the delay of the test packet (so the types in the complement have even greater arrivals, but these types in the complement contribute to the delay of the test packet only their fair share in proportion to w(τ′,l)). As a rougher approximation, however, one can imagine the event Dj≧dj as arising mostly in just two ways: either

Aτ+Ak−τ≧Sk&Ak−τ≦(1−w(τ,l))Sk,  (6)

or

Aτ≧w(τ,l)Sk&Ak−τ>(1−w(τ,l)Sk.  (7)

A complication arises in trying to estimate the joint probability of the combination (6) with the techniques that have been described so far. The problem is that the techniques estimate probabilities of the form of random variables being far above their means, but the combination (6) includes a condition of one random variable, Ak−τ, being not too far above its mean. Recall now that the basic technique for estimating the probability P(X≧x) is with an exponential moment of the form Eep(X−x) with p≧0. Naturally, the tightest bound of this form, that can be obtained on the probability, is to take p to minimize the expression Eep(X−x). Letting ΛX(p)=log(EepX), notice that the minimizing value of p is the one that makes ΛX′(p)=x. It turns out that, if X is a sum of independent random variables Xi and if each of these random variables has a continuous distribution (though this is not the situation of the analysis here), then the most likely values of the individual random variables Xi, given that the sum X≧x, are given by xi=ΛXi′(p), where p is the minimizing value that makes ΛX′(p)=x. The idea is to use this minimizing value of p to identify the most likely values of the components of a sum that is barely attaining a particular rare value.
Pursuing this idea, it is useful to introduce the following quantities. First, let μ′(p;t,σ,ρ)=ρtep(σ+ρt), the derivative of μ(p;t,σ,ρ) with respect top. Then, let x(p;t,σ,ρ)=μ′(p;t,σ,ρ)/μ(p;t,σ,ρ). Next, for any stream s′, let As′(p;t)=x(p;t,σs′,ρs′). Also, for any stream s′, let H(s′) be the set of links traversed by s′. Then, let



  
    
      
        
          A
          τ
        
        ⁡
        
          (
          
            
              p
              ;
              l
            
            ,
            t
          
          )
        
      
      =
      
        
          ∑
          
            
              s
              ′
            
            ∈
            
              
                S
                τ
              
              :
              
                l
                ∈
                
                  H
                  ⁡
                  
                    (
                    
                      s
                      ′
                    
                    )
                  
                
              
            
          
        
        ⁢
        
          
            A
            
              s
              ′
            
          
          ⁡
          
            (
            
              p
              ;
              t
            
            )
          
        
      
    
    ,
    
      

    
    ⁢
    
      
        
          A
          
            k
            -
            τ
          
        
        ⁡
        
          (
          
            
              p
              ;
              l
            
            ,
            t
          
          )
        
      
      =
      
        
          ∑
          
            
              τ
              ′
            
            ∈
            
              
                T
                k
              
              -
              
                {
                τ
                }
              
            
          
        
        ⁢
        
          
            ∑
            
              
                s
                ′
              
              ∈
              
                
                  S
                  τ
                
                :
                
                  l
                  ∈
                  
                    H
                    ⁡
                    
                      (
                      
                        s
                        ′
                      
                      )
                    
                  
                
              
            
          
          ⁢
          
            
              A
              
                s
                ′
              
            
            ⁡
            
              (
              
                p
                ;
                t
              
              )
            
          
        
      
    
    ,
    
      

    
    ⁢
    and
  




  
    
      
        A
        
          >
          k
        
      
      ⁡
      
        (
        
          
            p
            ;
            l
          
          ,
          t
        
        )
      
    
    =
    
      
        ∑
        
          
            k
            ′
          
          ∈
          
            K
            :
            
              
                k
                ′
              
              >
              k
            
          
        
      
      ⁢
      
        
          ∑
          
            
              τ
              ′
            
            ∈
            
              T
              
                k
                ′
              
            
          
        
        ⁢
        
          
            ∑
            
              
                s
                ′
              
              ∈
              
                
                  S
                  
                    τ
                    ′
                  
                
                :
                
                  l
                  ∈
                  
                    H
                    ⁡
                    
                      (
                      
                        s
                        ′
                      
                      )
                    
                  
                
              
            
          
          ⁢
          
            
              A
              
                s
                ′
              
            
            ⁢
            
              
                (
                
                  p
                  ;
                  t
                
                )
              
              .
            
          
        
      
    
  



With the aid of these quantities, consider the event Aτ+Ak−τ≧Sk, a basic event that occurs in the definition of the combination (6). Notice that this basic event is equivalent to Aτ+Ak−τ+A>k≧C(l)(tj+dj). Let p be the corresponding minimizing value, the value that minimizes Eep[Aτ+Ak−τA>k−C(l)(tj+dj)], or, more precisely, the version of this expectation that is based on the on-off distributions for the individual streams making up the arrival quantities Aτ, Ak−τ, and A>k. For this minimizing value p, make note of the corresponding values Ak−τ(p;l,tj+dj) and A>k(p;l,tj+dj), and check whether

Ak−τ(p;l,tj+dj)≦(1−w(τ,l))[C(l)(tj+dj)−A>k(p;l,tj+dj)].  (8)

If so, then the combination (6) makes the dominant contribution to P(Dj≧dj); otherwise, the other combination (7) is the dominant contributor. In this way, one of these two contributors can be identified as the appropriate estimate of P(Dj≧dj), at least if the quantities dj and tj are specified.
Now, if the second contributor turns out to be the dominant one, then that contributor is the probability of the intersection of two rare events, but these two events are dependent. So the probability of the intersection is not just the product of the probabilities. To handle the dependence, one must simply consider all relevant values a for the amount A>k that enters both of the conditions of (7). Of course, values a, that are below the mean A>k(0;tj+dj,l) of A>k are not relevant toward contributing the most likely situations for delaying the test packet more than dj, so attention can be restricted to a≧A>k(0;tj+dj,l). Now, given A>k=a≧A>k(0;tj+dj,l), the random variables Aτ and Ak−τ are independent. Accordingly, for a given value a of A>k, which corresponds to Sk(a)=C(l)(tj+dj)−a, the usual techniques can be used to estimate

P(Aτ>w(τ,l)Sk(a)&Ak−τ>(1−w(τ,l))Sk(a)),

which is then the estimate for the conditional probability of Dj>dj, given A>k=a. Hence, multiplying by P(A>k=a), which is estimated by P(A>k≧a) for a≧A>k (0;tj+dj) in the large-deviations procedure, and then maximizing over a, one obtains the final estimate for P(Dj≧dj), at least when based on the idea of the busy period having age tj.

Of course, the age tj of the busy period could be anything, at least in principle. Continuing with the large-deviations idea of approximating the probability of a union with the maximum of the individual probabilities, the actual estimate of P(Dj≧dj) should be the maximum, over all choices of tjε[0,∞) of the values obtained for each of these choices. Naturally, as was stated when the Dj were introduced, the final estimate of P(ΣjDj>{hacek over (δ)}s) is obtained by maximizing over all choices of the decomposition of {hacek over (δ)}s into the components dj.
It is convenient to introduce some auxiliary functions that represent various exponential moments for certain sums of on-off random variables that are taken to represent certain random arrivals. Explicitly, for any service type τ, and any link lεL, for any p,t≧0, let




  
    
      μ
      τ
    
    ⁡
    
      (
      
        
          p
          ;
          t
        
        ,
        l
      
      )
    
  
  =
  
    
      ∏
      
        s
        ∈
        
          
            S
            τ
          
          :
          
            l
            ∈
            
              H
              ⁡
              
                (
                s
                )
              
            
          
        
      
    
    ⁢
    
      
    
    ⁢
    
      
        μ
        ⁡
        
          (
          
            
              p
              ;
              t
            
            ,
            
              σ
              s
            
            ,
            
              ρ
              s
            
          
          )
        
      
      .
    
  




    
    
        Similarly, for any priority class kεK, and any link lεL, for any p, t≧0, let
    
    






  
    
      μ
      k
    
    ⁡
    
      (
      
        
          p
          ;
          t
        
        ,
        l
      
      )
    
  
  =
  
    
      ∏
      
        τ
        ∈
        
          T
          k
        
      
    
    ⁢
    
      
    
    ⁢
    
      
        
          μ
          τ
        
        ⁡
        
          (
          
            
              p
              ;
              t
            
            ,
            l
          
          )
        
      
      .
    
  




    
    
        Next, for any kεK, for any lεL, for any p, t≧0, let
    
    






  
    
      μ
      
        >
        k
      
    
    ⁡
    
      (
      
        
          p
          ;
          t
        
        ,
        l
      
      )
    
  
  ⁢
  
    
      ∏
      
        
          k
          ′
        
        ∈
        
          K
          :
          
            
              k
              ′
            
            >
            k
          
        
      
    
    ⁢
    
      
    
    ⁢
    
      
        
          μ
          
            k
            ′
          
        
        ⁡
        
          (
          
            
              p
              ;
              t
            
            ,
            l
          
          )
        
      
      .
    
  




    
    
        Finally, for any kεK, for any τεTk for any lεL, for any p, t≧0, let
    
    






  
    
      μ
      
        k
        -
        τ
      
    
    ⁡
    
      (
      
        
          p
          ;
          t
        
        ,
        l
      
      )
    
  
  =
  
    
      ∏
      
        
          τ
          ′
        
        ∈
        
          
            T
            k
          
          -
          
            {
            τ
            }
          
        
      
    
    ⁢
    
      
    
    ⁢
    
      
        
          μ
          
            τ
            ′
          
        
        ⁡
        
          (
          
            
              p
              ;
              t
            
            ,
            l
          
          )
        
      
      .
    
  




Basically, μτ(p;tj,l) will be used to estimate EepAτ, μk−τ(p;tj+dj,l) will be used to estimate EepAk−τ, and μ>k(p;tj+dj,l) will be used to estimate EepA>k.
The next step is to introduce a certain combination of these quantities, where this combination would serve as the estimate of P(Dj>dj) if the contribution of (6) would turn out to be dominant. More precisely, for any kεK, for any τεTk, for any lεL, for any p,t,d≧0, let

{tilde over (μ)}τ(p;t,d,l)=μτ(p;t;l)·μk−τ(p;t+d,l)·μ>k(p;t+d,l)e−pC(l)(t+d).

Next, for any kεK, for any τεTk, for any lεL, for any t,d≧0, let




  
    
      
        p
        τ
      
      ⁡
      
        (
        
          t
          ,
          d
          ,
          l
        
        )
      
    
    =
    
      arg
      ⁢
      
        
      
      ⁢
      
        
          min
          p
        
        ⁢
        
          
        
        ⁢
        
          
            
              μ
              ~
            
            τ
          
          ⁡
          
            (
            
              
                p
                ;
                t
              
              ,
              d
              ,
              l
            
            )
          
        
      
    
  
  ,




the value of p that minimizes {tilde over (μ)}τ(p;t,d,l). Naturally, the minimum value itself is then {tilde over (μ)}τ(t,d,l)={tilde over (μ)}τ(pτ(t,d,l);t,d,l). Now, the quantity, that summarizes whether the condition (8) holds, is





  
    
      
        
          
            θ
            τ
          
          ⁡
          
            (
            
              t
              ,
              d
              ,
              l
            
            )
          
        
        =
        
        ⁢
        
          
            (
            
              1
              -
              
                w
                ⁡
                
                  (
                  
                    τ
                    ,
                    l
                  
                  )
                
              
            
            )
          
          [
          
            
              
                C
                ⁡
                
                  (
                  l
                  )
                
              
              ⁢
              
                (
                
                  t
                  +
                  d
                
                )
              
            
            -
          
        
      
    
  
  
    
      
        
          
          ⁢
          
            
              A
              
                >
                k
              
            
            ⁡
            
              (
              
                
                  
                    
                      p
                      τ
                    
                    ⁡
                    
                      (
                      
                        t
                        ,
                        d
                        ,
                        l
                      
                      )
                    
                  
                  ;
                  l
                
                ,
                
                  t
                  +
                  d
                
              
              )
            
          
          ]
        
        -
      
    
  
  
    
      
        
        ⁢
        
          
            
              A
              
                k
                -
                τ
              
            
            ⁡
            
              (
              
                
                  
                    
                      p
                      τ
                    
                    ⁡
                    
                      (
                      
                        t
                        ,
                        d
                        ,
                        l
                      
                      )
                    
                  
                  ;
                  l
                
                ,
                
                  t
                  +
                  d
                
              
              )
            
          
          .
        
      
    
  




Specifically, if Θτ(t,d,l)≧0, then {hacek over (μ)}τ(t,d,l) is to be used to estimate P(Dj≧dj), at least when t is regarded as the age of the busy period.
To deal with the other possibility, however, some additional quantities must be introduced. First, corresponding to the idea of estimating P(A>k>a), define, for any kεK, for any lεL, for any p,t≧0, for any a≧A>k(0;t,l), π>k(a;p,t,l)=μ>k(p;t,l)epa. Then let p>k(a;t,l)=arg minπ>k(a;p,t,l), the minimizing value of p, so the minimum value itself is π>k(a;t,l)=π>k(a;t,l),t,l).
The next quantities are for estimating certain conditional probabilities. Thus, for any service type τ, for any link lεL, for any p,t≧0, and for any y, let πτ(p;t,y,l)=μτ(p;t,l)e−py. Then, for any service type τ, for any link lεL, for any t≧0, and for any y, let πτ(t,y,l)=minpπ(p;t,y,l). Similarly, for any kεK, for any τεTk, for any lεL, for any p,t≧0, for any y, let πk−τ(p;t,y,l)=μk−τ(p;t,l)e−py, and let πk−τ(t,y,l)=minpπk−τ(p;t,y,l).
The next step is to multiply these estimates of conditional and unconditional probabilities. Formally, for any kεK, for any τεTk, for any lεL, for any t,d≧0, for any a≧A>k(0,t+d,l), let



  
    
      
        
          
            
              
                π
                ~
              
              τ
            
            ⁡
            
              (
              
                a
                ,
                t
                ,
                d
                ,
                l
              
              )
            
          
          =
          
          ⁢
          
            
              
                
                  π
                  _
                
                
                  k
                  -
                  τ
                
              
              ⁡
              
                (
                
                  
                    t
                    +
                    d
                  
                  ,
                  
                    
                      (
                      
                        1
                        -
                        
                          w
                          ⁡
                          
                            (
                            
                              τ
                              ,
                              l
                            
                            )
                          
                        
                      
                      )
                    
                    ⁡
                    
                      [
                      
                        
                          
                            C
                            ⁡
                            
                              (
                              l
                              )
                            
                          
                          ⁢
                          
                            (
                            
                              t
                              +
                              d
                            
                            )
                          
                        
                        -
                        a
                      
                      ]
                    
                  
                  ,
                  l
                
                )
              
            
            ·
          
        
      
    
    
      
        
          
          ⁢
          
            
              
                
                  π
                  _
                
                τ
              
              ⁡
              
                (
                
                  t
                  ,
                  
                    
                      w
                      ⁡
                      
                        (
                        
                          τ
                          ,
                          l
                        
                        )
                      
                    
                    ⁡
                    
                      [
                      
                        
                          
                            C
                            ⁡
                            
                              (
                              l
                              )
                            
                          
                          ⁢
                          
                            (
                            
                              t
                              +
                              d
                            
                            )
                          
                        
                        -
                        a
                      
                      ]
                    
                  
                  ,
                  l
                
                )
              
            
            ·
          
        
      
    
    
      
        
          
          ⁢
          
            
              
                
                  π
                  _
                
                
                  >
                  k
                
              
              ⁡
              
                (
                
                  
                    a
                    ;
                    
                      t
                      +
                      d
                    
                  
                  ,
                  l
                
                )
              
            
            .
          
        
      
    
  




  
    Then
    ⁢
    
      
    
    ⁢
    let
  




  
    
      
        
          π
          ^
        
        τ
      
      ⁡
      
        (
        
          t
          ,
          d
          ,
          l
        
        )
      
    
    =
    
      
        sup
        
          a
          ≥
          
            
              A
              
                >
                k
              
            
            ⁡
            
              (
              
                
                  0
                  ;
                  l
                
                ,
                
                  t
                  +
                  d
                
              
              )
            
          
        
      
      ⁢
      
        
          
            
              π
              ~
            
            τ
          
          ⁡
          
            (
            
              a
              ,
              t
              ,
              d
              ,
              l
            
            )
          
        
        .
      
    
  



The significance of this quantity is that, if Θτ(t,d,l)<0, then πτ(t,d,l) is to be used to estimate P(Dj≧dj), at least when t is regarded as the age of the busy period.
The next step is to specify explicitly the version of the quantity that depends on Θτ(t,d,l). Accordingly, for any service type τ, for any link lεL, for any t,d≧0, let




  
    
      q
      τ
    
    ⁡
    
      (
      
        t
        ,
        d
        ,
        l
      
      )
    
  
  =
  
    {
    
      
        
          
            
              
                μ
                ⋓
              
              τ
            
            ⁡
            
              (
              
                t
                ,
                d
                ,
                l
              
              )
            
          
        
        
          
            
              if
              ⁢
              
                
              
              ⁢
              
                
                  θ
                  τ
                
                ⁡
                
                  (
                  
                    t
                    ,
                    d
                    ,
                    l
                  
                  )
                
              
            
            ≥
            0
          
        
      
      
        
          
            
              
                π
                ^
              
              τ
            
            ⁢
            
              (
              
                t
                ,
                d
                ,
                l
              
              )
            
          
        
        
          
            
              if
              ⁢
              
                
              
              ⁢
              
                
                  θ
                  τ
                
                ⁡
                
                  (
                  
                    t
                    ,
                    d
                    ,
                    l
                  
                  )
                
              
            
            <
            0.
          
        
      
    
  




The next step is to maximize over t. Thus, let qτ(d,l)=supr≧0qτ(t,d,l), which is the estimate for P(Dj>dj) without any disclaimers about the age of the busy period.
Finally, for any service type τ, for any stream sεSτ, the estimate for the probability of the end-to-end queuing delay of stream s exceeding any delay δ>0 is given by




  
    
      
        q
        ~
      
      s
    
    ⁡
    
      (
      δ
      )
    
  
  =
  
    
      sup
      
        
          
            
              d
              ∈
              
                
                  [
                  
                    0
                    ,
                    ∞
                  
                  )
                
                
                  H
                  ⁡
                  
                    (
                    s
                    )
                  
                
              
            
            :
            
              ∑
              
                
              
              ⁢
              d
            
          
          =
          δ
        
        ⁢
        
          
        
      
    
    ⁢
    
      
        ∏
        
          0
          ≤
          j
          <
          
            H
            ⁡
            
              (
              s
              )
            
          
        
      
      ⁢
      
        
          
            
              
                
              
              ⁢
              q
            
            _
          
          τ
        
        ⁢
        
          
            (
            
              
                d
                j
              
              ,
              
                λ
                ⁡
                
                  (
                  
                    j
                    ,
                    s
                  
                  )
                
              
            
            )
          
          .
        
      
    
  




In particular, the estimate for the probability of the end-to-end queuing delay of stream s exceeding {hacek over (δ)}s is {tilde over (q)}s({hacek over (δ)}s).
Taking into consideration that the on-line calculations in an admission-control procedure should be reasonably fast, the following is a simpler version of the calculation scheme. As with the GAD calculations, the idea is to introduce some isolating thresholds, with the purpose that the on-line computations during an admission decision should only be checking whether the admission of the new session would continue to respect the isolating thresholds, where the off-line computations, that determine the values for these thresholds, would be choosing these values in such a way that, as long as the thresholds are respected, all the admitted sessions would be receiving service that meets their standards. For definiteness, the isolating thresholds here will be the same as for the GAD case: for each priority class k, a bound αk on the occupancy due to priority class k, and for each service type τ, a bound Tτ on the burst times σs/ρs of all the streams s admitted as type τ. The on-line computations for admission control are then very quick, because they involve only two kinds of simple checks: first checking that a stream s requesting admission satisfies the burst-time bound Tτ, and then checking that, for each hop on the stream's route, the relevant occupancy bound αk continues to be respected at that hop. The next part of the discussion identifies the off-line computations that need to be done for dealing with the isolating thresholds αk and Tτ.
Mainly, the off-line computations here follow the same pattern as in the discussion above, except for a few modifications. Indeed, the principal change is to revise the definition of the function μτ(p;t,l), that summarizes the type-d traffic arriving at a link lεL, to reflect the fact that, in the off-line situation, one does not know the details of the individual streams. Instead, one only has the aggregate bounds. Accordingly, in this situation, the appropriate definition of the quantity is μτ(p;t,l)=μ(p,t,αk, w(τ,l)C(l)Tτ,αkw(τ,l)C(l)), where kεK is the priority class of service type τ. Similarly, the analogous changes should be made to the definitions of Aτ(p;l,t), Ak−τ(p;l,t), and A>k(p;l,t). Explicitly, the revised definitions are the following.



  
    
      
        
          A
          τ
        
        ⁡
        
          (
          
            
              p
              ;
              l
            
            ,
            t
          
          )
        
      
      =
      
        
          ∂
          
            ∂
            p
          
        
        ⁢
        
          log
          ⁡
          
            (
            
              
                μ
                τ
              
              ⁡
              
                (
                
                  
                    p
                    ;
                    t
                  
                  ,
                  l
                
                )
              
            
            )
          
        
      
    
    ,
    
      

    
    ⁢
    
      
        
          A
          
            k
            -
            τ
          
        
        ⁡
        
          (
          
            
              p
              ;
              l
            
            ,
            t
          
          )
        
      
      =
      
        
          ∑
          
            
              τ
              ′
            
            ∈
            
              
            
            ⁢
            
              
                T
                k
              
              -
              
                {
                τ
                }
              
            
          
        
        ⁢
        
          
            A
            
              τ
              ′
            
          
          ⁡
          
            (
            
              
                p
                ;
                l
              
              ,
              t
            
            )
          
        
      
    
    ,
    and
  




  
    
      
        A
        
          >
          k
        
      
      ⁡
      
        (
        
          
            p
            ;
            l
          
          ,
          t
        
        )
      
    
    =
    
      
        ∑
        
          
            k
            ′
          
          ∈
          
            
          
          ⁢
          K
        
      
      ⁢
      
        
          ∑
          
            
              τ
              ′
            
            ∈
            
              T
              
                k
                ′
              
            
          
        
        ⁢
        
          
            
              A
              
                τ
                ′
              
            
            ⁡
            
              (
              
                
                  p
                  ;
                  l
                
                ,
                t
              
              )
            
          
          .
        
      
    
  



After these preliminary changes the rest of the calculations of the earlier discussion remain unchanged until after the quantities qτ(d,l) are identified, which basically characterize queuing delays at the individual links. At this point, the remaining task is to combine the link delays into end-to-end delays, which, in the earlier discussion, was represented by the calculation of the quantities {tilde over (q)}s({hacek over (δ)}s) and the comparison with the targets ετ. Here, however, the individual streams s, that will be present at any one instant, are not known. Accordingly, the off-line computation must perform these comparisons for all possible routes. Overall, the full task for the off-line computation is to tune the isolating thresholds to be as “permissive” as possible, in the sense of tending to lead to the admission of more of the sessions that are likely to request admission, while missing none of the service targets for delay.
Some illustrative concrete calculations will now be presented, based on the statistical calculations that have just been described, where the values to be obtained for the isolating parameters αk will turn out to differ from the values obtained in the earlier illustrative concrete calculations aiming at worst-case targets (and the differences will be in the direction of admitting more sessions). Yet, the meaning of these parameters as bounds on certain sums of the rate parameters of the streams being admitted to the network is the same. The calculations are for the same network as above, with one difference being the switch from worst-case targets for the service types to statistical targets for these same service types. Accordingly, the basic parameter values, that define the network environment, such as the link capacities C(l), the delay targets δτ, and the bound h on the hop count of any route, will basically be the same with the one exception being that now, since some of the service targets are statistical, there will be some tolerances ετ for the probabilities of missing the delay targets δτ.
Now, in a real network, the expectation is that deterministic targets will be employed for higher priorities, while statistical targets will be employed for lower priorities (except that the lowest priority would typically have no target at all). To illustrate the implications of the statistical calculations, the following discussion will be phrased as if all the higher-priority service types, EF, AF0, and AF1, have statistical targets. For definiteness, the corresponding tolerances will be taken to be εEF=10−6, εAF0=10−3, and εAF1=10−2, though the precise values will turn out not to matter in the calculations to be presented here. If the above calculations were followed in complete detail, then the resulting values for the occupancy thresholds αk would depend on the tolerances ετ, and, as long as the sum of the occupancy thresholds remained below the value 1, the thresholds would rise as the tolerances rose. Yet, the complexity of the calculations would tend to hide the main point to be made here that, by adopting the statistical approach, one can “justify” the admission of more sessions. Indeed, this point can be made with a simplification of the detailed calculations, so it is this simplified calculation that will be followed here, even though the simplified calculation ignores the values of the tolerances ετ.
Basically, the simplified calculation is just the application of the more detailed calculation to the special case of zero values for all the tolerances. So occupancy thresholds αk, that are adequate according to the simplified calculation, will certainly be adequate for satisfying strictly positive tolerances. The simplification of this special case of zero tolerances is that, instead of dealing with all the minima, maxima, and exponential moments, one can simply deal with worst cases. In contrast to the earlier worst-case calculations for GADs, however, the simplified “worst-case” calculations here inherit the assumption of no growth in a stream's burstiness as it traverses its route. Thus, it is this assumption that allows the statistical calculations to “justify” the admission of more streams for the special case of zero tolerances, where these calculations are equivalent to the version of worst-case calculations that assume that there is no growth in burstiness along a stream's route.
Now, for the network being considered, the point to be made by these calculations is that the statistical calculations can justify raising the isolating thresholds from αEF=0.1 and αAF=0.06 to, say, αEF=0.5 and αAF=0.38. Ultimately, of course, one has the possibility of raising one at the cost of dropping the other, so the values just identified can be changed to adapt to different expectations for the traffic mix to be served. As for establishing the adequacy of these values, under the assumption of burstiness parameters that are impervious to queuing variations, consider first the implications for the high-priority EF traffic. Under the assumptions being imposed, the EF traffic, that can arrive at any given link l, can have a burst of size no greater than αEFC(l)TEF. Thus, as long as αEF≦1, so that high-priority queues cannot grow any greater than would be produced by a maximal burst, the maximum queuing delay that can be suffered by an EF packet at a link l is no greater than αEFC(l)TEF/C(l)=αEFTEF, where use is made here of the further assumption in the statistical calculations that no leftover low-priority packet can contribute any delay to higher-priority traffic. Evidently, with the choices of the parameter values for the network, the maximal end-to-end delay for an EF packet is bounded, in milliseconds, by {circumflex over (Δ)}+h(αEFTEF+ J)=10.004+20αEF which is always less than δEF=125 ms, as long as αEF≦1. In particular, the EF service target is met with the value αEF=0.5 advertised above.
The analysis for the lower-priority AF service types is more complicated. Suppose a “test” packet of some AF type τ arrives at a link l at time 0 as part of a burst of type-τ traffic, along with a burst of EF traffic, suppose that this burst begins a busy period for the link in which all the traffic streams are generating traffic at the rates given by their rate parameters, and suppose that the test packet loses all the “tie-breaking” procedures to all the competing packets in the initial burst. This is the situation that creates the greatest delay for the test packet, assuming that the link is not “oversubscribed”, i.e. as long as αEF+αAF≦1. Of course, the link spends the initial phase of the busy period serving EF traffic alone. The size of the initial burst of EF traffic could be as large as αEFC(l)TEF, and with further EF traffic arriving at rate αEFC(l), this initial phase could last as long as




  
    
      D
      ^
    
    
      >
      AF
    
  
  =
  
    
      
        
          α
          EF
        
        ⁢
        
          C
          ⁡
          
            (
            l
            )
          
        
        ⁢
        
          T
          EF
        
      
      
        
          C
          ⁡
          
            (
            l
            )
          
        
        -
        
          
            α
            EF
          
          ⁢
          
            C
            ⁡
            
              (
              l
              )
            
          
        
      
    
    =
    
      
        
          
            α
            EF
          
          ⁢
          
            T
            EF
          
        
        
          1
          -
          
            α
            EF
          
        
      
      .
    
  




After this phase, the link can serve the lower-priority type-τ traffic, but, if the EF traffic continues to arrive at its maximal allowed rate, and if the other AF service type occupies its share of the link's remaining capacity, then the link will be delivering service to the type-τ traffic at rate w(τ,l)C(l)(1−αEF). In particular, the initial burst of type-τ traffic, of size no greater than αAFC(l)Tτ, will have been served within a time




  
    
      
        α
        AF
      
      ⁢
      
        C
        ⁡
        
          (
          l
          )
        
      
      ⁢
      
        T
        τ
      
    
    
      
        w
        ⁡
        
          (
          
            τ
            ,
            l
          
          )
        
      
      ⁢
      
        C
        ⁡
        
          (
          l
          )
        
      
      ⁢
      
        (
        
          1
          -
          
            α
            EF
          
        
        )
      
    
  
  =
  
    
      
        α
        AF
      
      ⁢
      
        T
        τ
      
    
    
      
        w
        ⁡
        
          (
          
            τ
            ,
            l
          
          )
        
      
      ⁢
      
        (
        
          1
          -
          
            α
            EF
          
        
        )
      
    
  




after {circumflex over (D)}>AF. Since the test packet is part of this burst, the maximum delay, that a type-τ packet could experience on a single link, is bounded by





  
    
      
        
          D
          ^
        
        
          >
          AF
        
      
      +
      
        
          
            α
            AF
          
          ⁢
          
            T
            τ
          
        
        
          
            w
            ⁡
            
              (
              
                τ
                ,
                l
              
              )
            
          
          ⁢
          
            (
            
              1
              -
              
                α
                EF
              
            
            )
          
        
      
    
    =
    
      
        
          
            α
            EF
          
          ⁢
          
            T
            EF
          
        
        
          1
          -
          
            α
            EF
          
        
      
      +
      
        
          
            α
            AF
          
          ⁢
          
            T
            τ
          
        
        
          
            w
            ⁡
            
              (
              
                τ
                ,
                l
              
              )
            
          
          ⁢
          
            (
            
              1
              -
              
                α
                EF
              
            
            )
          
        
      
    
  
  ;




provided that αEF+αAF≦1.

The next step is to use this information to identify bounds on end-to-end delays for AF traffic. First, for a particular type τ of AF traffic, the maximal end-to-end delay in milliseconds is




  
    
      
        
          
            D
            ^
          
          τ
        
        =
        
          
            Δ
            ^
          
          +
          
            h
            ⁡
            
              [
              
                
                  J
                  _
                
                +
                
                  
                    
                      α
                      EF
                    
                    ⁢
                    
                      T
                      EF
                    
                  
                  
                    1
                    -
                    
                      α
                      EF
                    
                  
                
                +
                
                  
                    
                      α
                      AF
                    
                    ⁢
                    
                      T
                      τ
                    
                  
                  
                    
                      w
                      ⁡
                      
                        (
                        
                          τ
                          ,
                          l
                        
                        )
                      
                    
                    ⁢
                    
                      (
                      
                        1
                        -
                        
                          α
                          EF
                        
                      
                      )
                    
                  
                
              
              ]
            
          
        
      
    
  
  
    
      
        
          =
          
            30.004
            +
            
              
                20
                ⁢
                
                  α
                  AF
                
                ⁢
                
                  T
                  τ
                
              
              
                w
                ⁡
                
                  (
                  
                    τ
                    ,
                    l
                  
                  )
                
              
            
          
        
        ,
      
    
  




where use has been made of the known parameter values, including the choice αEF=0.5. In particular,




  
    
      
        
          D
          ^
        
        
          AF
          ⁢
          
            
          
          ⁢
          0
        
      
      =
      
        30.004
        +
        
          
            200
            ⁢
            
              α
              AF
            
          
          0.6
        
      
    
    ,
    while
  




  
    
      
        D
        ^
      
      
        AF
        ⁢
        
          
        
        ⁢
        1
      
    
    =
    
      
        30.004
        +
        
          
            1000
            ⁢
            
              α
              AF
            
          
          0.4
        
      
      =
      
        30.004
        +
        
          2500
          ⁢
          
            
              α
              AF
            
            .
          
        
      
    
  



On the basis of these expressions, one can choose αAF appropriately. First, notice that {circumflex over (D)}AF0≦500=δAF0 as long as αAF≦1. On the other hand, {circumflex over (D)}AF1 needs αAF to be rather smaller to avoid exceeding δAF1=1000 ms. Indeed, αAF has to be strictly smaller than 0.4. In particular, the value αAF=0.38 advertised above is adequate. One point to notice, however, is that, at least when αEF is as high as 0.5, the value of αAF must be restricted to be so small as not fully to occupy the link, in the sense that αAF+αEF must be strictly smaller than 1. The more important point, though, is that based on the statistical calculations, one can “justify” raising the occupancy bounds significantly (to αEF=0.5 and αAF=0.38 in this example from the values αEF=0.1 and αAF=0.06 that could be justified by the worst-case calculations).
FIG. 1 is a flow diagram illustrating operations performed at a device in accordance with a representative embodiment. Additional, fewer, or different operations may be performed depending on the particular implementation. In addition, the operations may be performed in an alternative order depending on the particular embodiment. In an operation 100, a request for a new session in a communications network is received at the device. In an embodiment, the request includes a priority class, a service type, and session characteristics. In an operation 110, an occupancy threshold value associated with the priority class is identified. A burst time threshold value associated with the priority class and the service type is identified in an operation 120. In a representative embodiment, both the occupancy threshold value and the burst time threshold value are predetermined. In an operation 130, a target delay threshold value associated with the priority class and the service type is identified. In an operation 140, an end-to-end delay value is calculated based on the identified occupancy threshold value, the burst time threshold value, and the session characteristics. In an operation 150, the calculated end-to-end value is compared with the identified target delay threshold value to determine if a delay threshold is satisfied. In an operation 160, the new session in the communications network is admitted if the delay threshold is satisfied.
FIG. 2 is a flow diagram illustrating operations performed at a device in accordance with a representative embodiment. Additional, fewer, or different operations may be performed depending on the particular implementation. In addition, the operations may be performed in an alternative order depending on the particular embodiment. In an operation 200, a request for a new session in a communications network is received at the device. In a representative embodiment, the request includes a priority class, a service type, and session characteristics. A burst time threshold value associated with the priority class and the service type is identified in an operation 210. In an operation 220, a burst time value calculated for the new session based on the session characteristics is compared with the identified burst time threshold value to determine if a burst time threshold is satisfied. An occupancy threshold value associated with the priority class is identified in an operation 230. In a representative embodiment, both the identified occupancy threshold value and the identified burst time threshold value are predetermined. In an operation 240, it is determined if the identified occupancy threshold value is satisfied at each hop on a route identified for the new session based on the session characteristics. In an operation 250, the new session in the communications network is admitted if the burst time threshold is satisfied and the identified occupancy threshold value is satisfied at each hop.
While there has been described a method of calculating admission control in an internet protocol network, it will be apparent to those skilled in the art that modifications and variations are possible without deviating from the teachings and broad scope of the present invention which shall be limited solely by the scope of the claims appended hereto.

## Claims

1. A method for providing admission control in a communications network, the method comprising:
at a first device,
receiving a request for a new session in a communications network, the request comprising a priority class, a service type, and session characteristics;
identifying an occupancy threshold value associated with the priority class, wherein the occupancy threshold value is predetermined;
identifying a burst time threshold value associated with the priority class and the service type, wherein the burst time threshold value is predetermined;
identifying a target delay threshold value associated with the priority class and the service type;
calculating an end-to-end delay value based on the identified occupancy threshold value, the burst time threshold value, and the session characteristics;
comparing the calculated end-to-end delay value with the identified target delay threshold value to determine if a delay threshold is satisfied; and
admitting the new session in the communications network if the delay threshold is satisfied.

2. The method of claim 1, further comprising rejecting the new session if the delay threshold is not satisfied.

3. The method of claim 1, further comprising:
identifying a target delay variation threshold value associated with the priority class and the service type;
calculating an end-to-end delay variation value based on the identified occupancy threshold value, the burst time threshold value, and the session characteristics;
comparing the calculated end-to-end delay value variation with the identified target delay variation threshold value to determine if a delay variation threshold is satisfied; and
admitting the new session in the communications network if the delay variation threshold is satisfied.

4. The method of claim 3, further comprising rejecting the new session if the delay variation threshold is not satisfied.

5. The method of claim 3, wherein the occupancy threshold value is predetermined to insure that a worst case end-to-end delay variation value satisfies the delay variation threshold for all service types.

6. The method of claim 3, wherein the occupancy threshold value is predetermined to insure that a statistical end-to-end delay variation value satisfies the delay variation threshold for all service types.

7. The method of claim 3, wherein the occupancy threshold value for a first priority class is predetermined to insure that a worst case end-to-end delay variation value satisfies the delay variation threshold for all service types, and further wherein the occupancy threshold value for a second priority class is predetermined to insure that a statistical end-to-end delay variation value satisfies the delay variation threshold for all service types.

8. The method of claim 1, wherein the calculated end-to-end delay value includes a delay value resulting from a higher priority class of the priority class.

9. The method of claim 8, wherein the calculated end-to-end delay value includes a second delay value resulting from a higher service type of the service type.

10. The method of claim 1, wherein the occupancy threshold value is predetermined to insure that a worst case end-to-end delay value satisfies the delay threshold for all service types.

11. The method of claim 1, wherein the occupancy threshold value is predetermined to insure that a statistical end-to-end delay value satisfies the delay threshold for all service types.

12. The method of claim 1, wherein the occupancy threshold value for a first priority class is predetermined to insure that a worst case end-to-end delay value satisfies the delay threshold for all service types, and further wherein the occupancy threshold value for a second priority class is predetermined to insure that a statistical end-to-end delay value satisfies the delay threshold for all service types.

13. A method for providing admission control in a communications network, the method comprising:
at a first device,
receiving a request for a new session in a communications network, the request comprising a priority class, a service type, and session characteristics;
identifying a burst time threshold value associated with the priority class and the service type, wherein the burst time threshold value is predetermined;
comparing a burst time value calculated for the new session based on the session characteristics with the identified burst time threshold value to determine if a burst time threshold is satisfied;
identifying an occupancy threshold value associated with the priority class, wherein the occupancy threshold value is predetermined;
determining if the identified occupancy threshold value is satisfied at each hop on a route identified for the new session based on the session characteristics; and
admitting the new session in the communications network if the burst time threshold is satisfied and the identified occupancy threshold value is satisfied at each hop.

14. The method of claim 13, further comprising calculating the occupancy threshold value.

15. The method of claim 14, wherein calculating the occupancy threshold value for a first priority class comprises:
identifying a target delay threshold value associated with the first priority class and the service type;
calculating a worst case end-to-end delay variation value; and
selecting a value for the occupancy threshold value for the first priority class such that the calculated worst case end-to-end delay variation value satisfies the identified target delay threshold value.

16. The method of claim 15, wherein calculating the occupancy threshold value for a second priority class comprises:
identifying a second target delay threshold value associated with the second priority class and the service type;
calculating a statistical end-to-end delay variation value; and
selecting a second value for the occupancy threshold value for the second priority class such that the calculated statistical end-to-end delay variation value satisfies the identified second target delay threshold value.

17. The method of claim 14, wherein calculating the occupancy threshold value comprises:
identifying a target delay threshold value associated with the priority class and the service type;
calculating a statistical end-to-end delay variation value; and
selecting a value for the occupancy threshold value such that the calculated statistical end-to-end delay variation value satisfies the identified target delay threshold value.