# Patent

## Title
Estimating origin-destination flow entropy

## Abstract

The preferred embodiments of the present invention are directed to estimating entropy of origin-destination (OD) data flows in a network. To achieve this, first and second sketches are created corresponding to ingress (i.e. origin) and egress (i.e. destination) flows. The sketches allow estimating entropy associated with data streams as well as entropy associated with an intersection of two or more of the data streams, which provides a mechanism for estimating the entropy OD flows in a network.

## Background
BACKGROUND
1. Field of the Invention
The preferred embodiments of the present invention are directed to an approach for estimating entropies of origin-destination flows in a network.
2. Brief Description of the Related Art
Detection and diagnosis of network anomalies is important for network maintenance and troubleshooting. There are various automated processes for detecting and diagnosing network anomalies based on the status of or data stored by routers. Recently, the entropy of network traffic has been proposed, in many different contexts, as an effective and reliable metric for anomaly detection and diagnosis. Measuring this quantity exactly in real time, however, is typically not possible on high-speed links due to prohibitively high computational and memory requirements. For this reason, various data streaming algorithms have been proposed to approximate this quantity. Data streaming is concerned with processing a long stream of data items in one pass to answer a class of queries regarding the stream.
It is often important to know the entropies of origin-destination (OD) flows. An OD flow represents the traffic that enters an ingress point (called an origin) and exits at an egress point (called a destination). Knowing these quantities allows for more insight into the dynamics of traffic inside an Internet service provider (ISP) network. However, the conventional entropy estimation algorithms generally cannot be extended to solve the entropies of OD flows in a network.
SUMMARY OF THE INVENTION
Preferred embodiments are directed to a method, medium, and system, which can include a computing device, for estimating entropy in a network. A first sketch associated with at least one ingress flow in a network is generated and a second sketch associated with at least one egress flow in the network is generated. The first sketch can be generated by passing a packet to a bucket based on a flow label associated with the packet and incrementing at least one counter in the first sketch based on a packet identification associated with the packet. The at least one ingress flow and egress flow form a origin-destination flow pair. An origin-destination (OD) flow entropy based on the first and second sketches is determined.
To determine the OD flow entropy, a first Lp norm is extracted from the first sketch and a second Lp norm is extracted from the second sketch. An intersection between the first and second sketches is determined using the Intersection Measurable Property (IMP). An entropy norm can be determined using the intersection between the first and second sketches and an approximation for xlnx. A total number of packets in a stream can be determined using at least one of the first and second sketches, which can be used when determining the OD flow entropy. An estimate for the OD flow entropy can be computed using the entropy norm and the total number of packets.
Other objects and features of the present invention will become apparent from the following detailed description considered in conjunction with the accompanying drawings. It is to be understood, however, that the drawings are designed as an illustration only and not as a definition of the limits of the invention.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 shows a portion of an exemplary network;
FIG. 2 is a block diagram of a preferred embodiment of an entropy estimator;
FIG. 3 shows a comparison between an xlnx approximation using the xlnx estimator and an actual xlnx computation;
FIG. 4 is a flowchart showing a preferred implementation for estimating the entropy using the entropy estimator 200;
FIG. 5 is a flowchart showing a implementation for generating sketches and estimating the Lp norm ∥S∥p from the sketches;
FIG. 6 is flowchart showing a computation of the Lp norm using the bucketing technique in accordance with a preferred embodiment of the present invention
FIG. 7 shows an exemplary computing device for implementing entropy estimation of OD flows in a network;
FIG. 8 is a cumulative density plot showing an error of the entropy estimator 200 for a trace;
FIG. 9 shows the effect of increases a number of buckets used when creating a sketch;
FIG. 10 shows that the elephant detection mechanism allows good overall performance using the entropy estimator of the preferred embodiments of the present invention;
FIG. 11 compares the cumulative density plots of the error for between two traces having different measurement intervals; and
FIG. 12 shows an error plot for the entropy is comparable to that for the entropy norm.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
The preferred embodiments of the present invention are directed to estimating entropy of origin-destination (OD) data flows in a network. To achieve this, one or more data structures, referred to herein as sketches, are created for ingress (i.e. origin) and egress (i.e. destination) flows. The sketches allow estimating entropy associated with data streams as well as entropy associated with an intersection of two or more of the data streams, which provides a mechanism for estimating the entropy OD flows in a network.
A flow, as used herein, refers to data packets that have identical values in their headers, which may preferably include some combination of a source address, destination address, source port, destination port, and protocol. Entropy is a network metric that indicates changes in data stream flows in a network. A sketch, as used herein, is a lower dimensional vector used to represent higher dimensional data and preferably includes counters for tracking information associated with data stream flows in a network. The ability to compute the intersection of two or more data streams from their sketches is referred to herein as the intersection measurable property (IMP). Ingress nodes and egress nodes in the network can maintain a sketch of the traffic flows into and out of the network during the measurement interval.
The ability to estimate entropy for origin-destination (OD) flows in a network can be highly beneficial in detecting and diagnosing anomalies in a network. Anomalies, such as distributed denial of service (DDos) attacks, which inundate a network with an exorbitant amount of traffic thereby slowing the network performance down, network failures, flash crowds, and scheduled maintenance tasks, such as route Internetwork Operating System (IOS) updates and customer migration, can result in network performance degradation and service disruption. Many of these events occur in a distributed manner in terms of the signatures (i.e. patterns formed by the anomalies in the network) and impact. Detecting these events and evaluating their impact on network services often requires monitoring traffic from a number of locations across the network. Changes in traffic distribution as a result of anomalies and/or scheduled maintenance may not be visible in a traditional volume matrix. However, the preferred embodiments can capture these events in real time by examining OD flows in the network. The approaches described herein enable accurate estimations of OD flow entropies fast enough for processing traffic at very high speed links (e.g. 10 million packets per second).
For a stream S of packets that contains n (transport-layer) flows with sizes (number of packets) ai, a2, . . . , an, respectively, the empirical entropy H(S) is defined as follows:




  
    
      
        
          -
          
            
              ∑
              
                i
                =
                1
              
              n
            
            ⁢
            
              
                
                  a
                  i
                
                s
              
              ⁢
              
                
                  log
                  2
                
                ⁡
                
                  (
                  
                    
                      a
                      i
                    
                    s
                  
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        1
        )
      
    
  




where





  s
  =
  
    
      ∑
      
        i
        =
        1
      
      n
    
    ⁢
    
      a
      i
    
  




is the total number of packets in the stream S. To estimate the empirical entropy H(S) of a stream S, an estimate of a related quantity called entropy norm ∥S∥H can be used. The entropy norm ∥S∥H can be defined as Σiai ln(ai), since H(S) can be rewritten as follows:





  
    
      
        
          
            
              
                
                  H
                  ⁡
                  
                    (
                    S
                    )
                  
                
                =
                
                ⁢
                
                  -
                  
                    
                      ∑
                      i
                    
                    ⁢
                    
                      
                        
                          a
                          i
                        
                        s
                      
                      ⁢
                      
                        
                          log
                          2
                        
                        ⁡
                        
                          (
                          
                            
                              a
                              i
                            
                            s
                          
                          )
                        
                      
                    
                  
                
              
              ,
            
          
        
        
          
            
              
                =
                
                ⁢
                
                  
                    
                      log
                      2
                    
                    ⁡
                    
                      (
                      e
                      )
                    
                  
                  [
                  
                    
                      ln
                      ⁡
                      
                        (
                        s
                        )
                      
                    
                    -
                    
                      
                        1
                        s
                      
                      ⁢
                      
                        
                          ∑
                          i
                        
                        ⁢
                        
                          
                            a
                            i
                          
                          ⁢
                          
                            ln
                            ⁡
                            
                              (
                              
                                a
                                i
                              
                              )
                            
                          
                        
                      
                    
                  
                  ]
                
              
              ,
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  
                    log
                    2
                  
                  ⁡
                  
                    (
                    e
                    )
                  
                
                [
                
                  
                    ln
                    ⁡
                    
                      (
                      s
                      )
                    
                  
                  -
                  
                    
                      1
                      s
                    
                    ⁢
                    
                      
                        
                        S
                        
                      
                      H
                    
                  
                
                ]
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




To determine the entropy of OD flows, the entropy norm ∥S∥H and the total number of packets s in a stream S are determined. The entropy of an OD flow stream ODij between an ingress point i and an egress point j, is defined herein as the entropy of their intersection, i.e., H(ODij)≡H(Oi∩Dj). The total number of packets s in a stream S represents the volume of the OD flow and is typically an unknown quantity to be estimated separately. The preferred embodiments of the present invention include an approach to estimating the entropy norm ∥S∥H=∥Oi∩Dj∥H and the total number of packets s to facilitate a computation of the entropy of OD flows H(Oi∩Dj). As discussed below the entropy norm ∥S∥H can be approximated using an xln(x) approximation and an estimated Lp norm, which represents a distance measure of a vector, where p is preferably a number in the range of 0 to 2.
FIG. 1 shows an exemplary network 100 for which the entropy of OD flows can be determined. The network 100 can include provider edge (PE) routers 111-114, also referred to herein as edge routers, and provider routers (P-routers) 121-126, the implementation of each being generally know to those skilled in the art. Customers typically connect to the network 100 via the PE routers 111-114. For example, a customer router 141 can connect to the edge router 111 and a customer router 142 can connect to the edge router 114 to allow each customer router 141-142 to communicate with each other and with other customer routers connected to the network 100.
Streams S of packets can be transmitted over the network 100 from a PE router (i.e. an origin) to another PE router (i.e. a destination). Streams flowing into the network 100 via a PE router can represent an ingress flow and streams flowing out of the network via a PE router can represent an egress flow. These streams S of packets preferably include headers, where identical headers associate the packets with a flow. As the streams S of packets flow through the PE routers 111-114, the PE routers 111-114 preferably maintain sketches 150, including ingress sketches 152 for ingress flows and egress sketches 154 for egress flows. These sketches 150 can be used when determining the entropy between OD flows of the network 100.
In a preferred embodiment, a sketch is an array {right arrow over (Y)} (i.e. Y[1 . . . l]) of floating point counters initially set to 0.0 at the beginning of a measurement interval. The critical operator for this approach is a set of p-stable hash functions shi, i=1, . . ., , l, implemented by a p-stable hashing unit 160. The p-stable hash functions shi map a flow label packet identification (pkt.id) into a random value drawn from a stable distribution S(p) in such a way that the same flow label is preferably mapped to the same random value and that different flow labels are mapped to independent random values. In addition, these p-stable hash functions shi are independent of one another.
Stable distributions are natural generalizations of the Gaussian distribution and generally take a fairly complicated form with four parameters. However, the preferred embodiments of the present invention use standard (normalized) cases that have three fixed parameters and one free parameter p, resulting in a restricted/standardized family denoted S(p), p ε (0,2]. The S(p)'s (e.g., S(0), S(1), S(2)) are uniquely characterized as follows:

    
    
        Let X be a random variable that takes distribution S(p) with probability density function ƒs(p)(x). Then its characteristic function E[eitX] satisfies the following equation:
    
    






  
    
      
        
          
            E
            ⁡
            
              [
              
                e
                
                  i
                  ⁢
                  
                    
                  
                  ⁢
                  tX
                
              
              ]
            
          
          ≡
          
            
              ∫
              
                -
                ∞
              
              ∞
            
            ⁢
            
              
                
                  f
                  
                    s
                    ⁡
                    
                      (
                      p
                      )
                    
                  
                
                ⁡
                
                  (
                  x
                  )
                
              
              ⁢
              
                (
                
                  
                    cos
                    ⁡
                    
                      (
                      xt
                      )
                    
                  
                  +
                  
                    i
                    ·
                    
                      sin
                      ⁡
                      
                        (
                        xt
                        )
                      
                    
                  
                
                )
              
            
          
        
        =
        
          
            e
            
              -
              
                
                  
                  t
                  
                
                p
              
            
          
          .
        
      
    
    
      
        (
        3
        )
      
    
  




The existence and uniqueness of a probability density function satisfying the above equation is well known to those in the art. For any p, the probability density function ƒs(p)(x) of the distribution S(p) is continuous and infinitely differentiable on (−∞, +∞). However, fs(p)(x) takes a closed form for three p values (p=2, 1,0.5). S(p) with p=2 and 1 corresponds to two well-known examples of stable distributions to those skilled in the art. S(2) is the Gaussian distribution having a mean 0 and standard deviation 2, with probability density function




  
    f
    ⁡
    
      (
      x
      )
    
  
  =
  
    
      1
      
        2
        ⁢
        
          
            2
            ⁢
            π
          
        
      
    
    ⁢
    
      e
      
        
          -
          
            x
            2
          
        
        /
        8
      
    
  




and S(1) is the Cauchy distribution with the density function





  
    f
    ⁡
    
      (
      x
      )
    
  
  =
  
    
      1
      π
    
    ⁢
    
      
        1
        
          1
          +
          
            x
            2
          
        
      
      .
    
  




For p values other than 2, 1, and 0.5, the distribution S(p) is generated using a simulation formula. Some properties of the stable distribution include stability and symmetry.

To implement the stable hash functions shi, . . . , shl, a standard methodology for generating random variables with stable distributions S(p) is provided through the following simulation formula:




  
    
      
        
          X
          =
          
            
              [
              
                
                  
                    sin
                    ⁡
                    
                      (
                      
                        p
                        ⁢
                        
                          
                        
                        ⁢
                        θ
                      
                      )
                    
                  
                  
                    
                      cos
                      
                        1
                        /
                        p
                      
                    
                    ⁢
                    θ
                  
                
                ⁢
                
                  
                    (
                    
                      cos
                      ⁡
                      
                        (
                        
                          θ
                          ⁡
                          
                            (
                            
                              1
                              -
                              p
                            
                            )
                          
                        
                        )
                      
                    
                    )
                  
                  
                    
                      1
                      /
                      p
                    
                    -
                    1
                  
                
              
              ]
            
            [
            
              
                (
                
                  1
                  
                    
                      -
                      ln
                    
                    ⁢
                    
                      
                    
                    ⁢
                    r
                  
                
                )
              
              
                
                  1
                  /
                  p
                
                -
                1
              
            
            ]
          
        
        ,
      
    
    
      
        (
        4
        )
      
    
  




where θ is chosen uniformly in [−π/2, π/2] and r is chosen uniformly in [0, 1].

In one embodiment, the p-stable hash functions are implemented by fixing two uniform hash functions uhj1 and uhj2 that map a flow identifier pkt.id to a θ value uniformly distributed in [−π/2, π/2], and an r value uniformly distributed in [0, 1], respectively. These two values are inserted into equation (4). However, computing equation (4) typically requires thousands of CPU cycles, and it is not generally practical to perform 40 such computations for each incoming packet.
In a preferred embodiment, the rate at which these stable hash functions (i.e., sh′js) are computed is increased by performing memory lookups into precomputed tables 170. The term in the first bracket of equation (4) is a function of θ and the term in the second bracket is a function of r. For implementing each shj, two uniform hash functions uhj1 and uhj1 are fixed that map a flow identifier pkt.id to two index values uniformly distributed in [1 . . . N1] and [1 . . . N2] respectively. Two lookup tables T1 and T2 are allocated that contain N1 and N2 entries respectively, and table entries (for both T1 and T2) contain l=20 blocks of 4 bytes each. N1*20 i.i.d. random variables distributed as the term in the first bracket are precomputed and are filled into T1 and N2*20 i.i.d. random variables distributed as the term in the second bracket are computed and are filled into T2. For each incoming packet, l=20 random values T1 [uhj1 (pkt.id)][j]*T2[uhj2(pkt.id)][j], j=1, . . . , 20, are returned as the computation result for sh1(pkt.id), sh2(pkt.id), . . . , shl(pkt.id). Since each sub-sketch requires two tables, a total of four tables are used. In one implementation, both N1 and N2 are set to fairly large values, such as 1M. Simulations show that stable distribution values generated based on this approach are generally indistinguishable from real stable distribution values. However, the described preferred implementation is very fast: two memory reads (4 bytes each) and a floating point multiplication for computing each shj (pkt.id). Index values uhj1 (pkt.id) and uhj2(pkt.id) are generated for estimating the L1+anorm (p=1+a) and can be reused for the lookup operations performed in estimating the L1−anorm (p=1−a), since all entries in these four tables are mutually independent.
Packets are processed with the sketch such that for each incoming packet, the packet's flow label is hashed by shi and the result is added to the sketch Yi. The values of counters in the sketch are analyzed. For example, the value of the counter for Y1 can be analyzed as follows. Suppose that there are n flows with flow labels id1, . . . , idn, and of sizes a1, . . . , an. The p-stable hash functions are denoted as follows: sh1(id1) as X1, sh1(id2) as X2, . . . , sh1(idn) as Xn. By the property of sh1, we know that X1, . . . , Xn are independent and identically distributed (i.i.d.) random variables with a standard distribution S(p). After processing the stream, the counter value Y1 becomes a1*X1+a2*X2+. . . +an*Xn. Since each Xi has distribution S(p) and they are mutually independent, by the stability property, the counter value Y1 is distributed as (Σi|ai|p)1/pX, where (Σlai|p)1/p represents the Lp norm ∥S∥p and X is distributed as S(p). Therefore, counters Y1 , Y2, . . . , Yl are i.i.d. draws from distribution ∥S∥pX. The data streaming algorithm is able to “modulate” the signal to be estimated, (Σi|ai|p)1/p, into these counter values.
FIG. 2 shows a block diagram of a OD flow entropy estimator 200 (hereinafter “entropy estimator 200”) in accordance with the preferred embodiments of the present invention. The estimator 200 can include an entropy norm generator 210, a filtering unit 240, a total number of packets estimator 250 (hereinafter “packet estimator 250”), and an OD flow entropy generator 260 (hereinafter “entropy generator 260”). The entropy norm ∥S∥H=Oi∩Dj|H is estimated via the entropy norm generator 210 which in turn uses a Lp Norm generator 220 and an xln(x) approximator 230. An estimate for the entropy norm ∥S∥H can be calculated based on counter values in one or more sketches.
In a preferred embodiment, the entropy estimator 200 receives sketches for ingress and egress flows for the PE routers in the network 100. The theory of stable distributions is preferably used to estimate the distance measure Lp norm ∥S∥p of a stream S for values of p in (0, 2]. As discussed above, the Lp norm ∥S∥p for a stream S that contains n flows of sizes a1, . . . , an can be defined as (Σi|ai|p)1/p, which can be rewritten as follows:

∥S∥pp=Σi|ai|p.   (5)

Lp norms for p values slightly above or below 1 play a crucial role in facilitating an estimation of OD flow.

To extract the Lp norm ∥S∥p from the counter values in the sketches 150, the Lp norm generator 220 preferably implements a median estimator approach to determine the Lp norm ∥S∥p, which can be represented as follows:




  
    
      
        
          Λ
          (
          
            Y
            ->
          
          )
        
        ≡
        
          
            
              median
              (
              
                
                  
                  
                    Y
                    1
                  
                  
                
                ,
                ...
                ⁢
                
                  
                
                ,
                
                  
                  
                    Y
                    l
                  
                  
                
              
              )
            
            
              DMed
              p
            
          
          .
        
      
    
    
      
        (
        6
        )
      
    
  




Λ(.) is used to denote an operator that extracts Lp norm estimates from sketches. Here median(|Y1|, . . . , |Yl|) denotes the sample median of the absolute values of the counter values and DMedp denotes the distribution median of S+(p). S+(p)is used herein to refer to the probability distribution of a random variable |X|, where X has distribution S(p) because its probability distribution function (p.d.f) is exactly twice the positive half of the p.d.f. for S(p), due to the symmetry of S(p). In this case, DMedp is the unique x0 value such that Pr[|X|>x0]=0.5, where X has distribution S(p). Due to the symmetry of S(p), DMedp is exactly the three-quarter quantile of S(p). Although there is no closed form for DMedp for most of the p values, it can be numerically calculated by simulation or other programs known to those skilled in the art.
Since Y1/∥S∥p, . . . , Yl/∥S∥p are i.i.d. random variables with distribution S(p), taking absolute value gives us i.i.d. draws from S+(p). For large enough l, their median should be close to the distribution median of S+(p). Therefore, the median (|Y1|, . . . , |Yl|) is divided by the distribution median of S|(p) to get an estimator of ∥S∥p. The Lp norm is preferably generated by the Lp norm generator 220 based on the property of the median, however, those skilled in the art will recognize that it is possible to construct estimators based on other quantities.
Since the entropy norm ∥S∥H can be defined as Σiai ln(ai) and ai can be estimated by the Lp norm ∥S∥p, the xln(x) approximator 230 can be used to estimate x ln(x) using a linear combination of a small number of functions in the family {xp|p ε (0,2]} to compute the entropy norm ∥S∥H using the Lp norm ∥S∥p. The parameters c1, . . . , ck ε R, where R represents real numbers, and p1, . . . , pk ε (0, 2] are provided such that xln(x) is approximated as follows:




  
    
      
        
          x
          ⁢
          
            
          
          ⁢
          
            ln
            ⁡
            
              (
              x
              )
            
          
        
        ≈
        
          
            ∑
            
              j
              =
              1
            
            k
          
          ⁢
          
            
              c
              j
            
            *
            
              
                x
                pj
              
              .
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




Using this, the entropy norm ∥S∥H of an OD flow stream S can be estimated as follows:




  
    
      
        
          
            ∑
            
              j
              =
              1
            
            k
          
          ⁢
          
            
              c
              j
            
            *
            
              
                
                S
                
              
              pj
              pj
            
          
        
        ,
      
    
    
      
        (
        8
        )
      
    
  




where Lp norm ∥S∥pj, j=1, . . . , k can be estimated using the sketch for stream intersection described herein. Using a linear combination of only two functions in the family, in the form of





  
    
      f
      ⁡
      
        (
        x
        )
      
    
    =
    
      
        1
        
          2
          ⁢
          α
        
      
      ⁢
      
        (
        
          
            x
            
              1
              +
              α
            
          
          -
          
            x
            
              1
              -
              α
            
          
        
        )
      
    
  
  ,




the xln(x) approximator 230 can approximate xln(x) closely for all x values in a large interval (e.g., [1,1000], [1, 5000], etc.). Here a is a tunable parameter that takes small values. For example, FIG. 3 shows how closely xln(x) can be approximated when a is 0.05 using





  
    1
    
      2
      ⁢
      α
    
  
  ⁢
  
    (
    
      
        x
        
          1
          +
          α
        
      
      -
      
        x
        
          1
          -
          α
        
      
    
    )
  




within the interval [1, 1000]. In other words, if transport-layer flows in an OD flow stream S generally have less than 1000 packets, the entropy norm ∥S∥H of the OD flow can be estimated as 10(∥S∥1.051.05−∥S∥0.950.95). This “symmetry” of the exponents (1+α and 1−α) around 1 serves other important purposes described in more detail below. Using the Lp norm generator 220 and the xln(x) approximator 230, the entropy norm generator 210 can estimate the entropy norm ∥S∥H as discussed above.

Still referring to FIG. 2, very large transport-layer flows (e.g., 10,000, 100,000, etc. of packets) inside an OD flow stream S are identified and filtered out of the entropy estimation with a “sample and hold” algorithm implemented by a stream filtering unit 240 using methods known to those skilled in the art. The contributions to the OD flow entropy of these very large transport-layer flows are computed separately.
To compute the entropy H from the entropy norm ∥S∥H, the packet estimator 250 preferably estimates the total volume s of the OD flow stream S using the sketches for estimating Lp norms based on the formula




  
    
      1
      
        2
        ⁢
        α
      
    
    ⁢
    
      (
      
        
          x
          
            1
            +
            α
          
        
        -
        
          x
          
            1
            -
            α
          
        
      
      )
    
  
  ,




where α is a small number such as 0.05. For this case, the function x can be approximated by





  
    (
    
      
        x
        
          1
          +
          α
        
      
      +
      
        x
        
          1
          -
          α
        
      
    
    )
  
  2




and therefore the OD flow volume (i.e., the L1 norm) can be approximated by the average of ∥S∥1+α1+αand ∥S∥1−α1−α calculated from the OD flow stream S. Therefore, the sketch data structure described herein allows the estimated L1 norm for the total number of packets s and entropy norm ∥S∥H to be computed.

The entropy generator 260 can use the entropy norm ∥S∥H generated by the entropy norm generator 210 and the total number of packets s estimated using the packet estimator 250 to estimate the entropy of OD flows in the network 100. The entropy generator 260 can estimate the entropy of OD flows based on equation (2).
FIG. 4 shows a flowchart for estimating the entropy using the entropy estimator 200. A sketch is created that represents packets for a flow (step 400). Any flows that exceed a specified size can be filtered out by the filtering unit 240 to be analyzed separated. The entropy estimator 200 receives the sketch and extracts an estimate of the Lp norm ∥S∥p from the sketch using the Lp norm generator 220 (step 405). The xln(x) aproximator 230 generates the entropy norm ∥S∥H based on the computed Lp norms (step 410). A total number of packets s in a stream is determined using the packet estimator 250 (step 415). Subsequently, the entropy H(S) is calculated based on equation (2) using the entropy norm ∥S∥H and the total number of packets s (step 420).
FIG. 5 shows a flowchart for generating sketches 150 and estimating the Lp norm ∥S∥p from the sketches 150 in more detail. First a sketch is created as an array Y[1, . . . ,l] and the counters of the sketch are initialized to zero (step 500). P-stable hash functions sh1 through shl are fixed (step 505). As packets sequentially pass through a origin or destination PE router, the counter values of the sketch are incremented using the p-stable hash function based on the packet ID of the incoming packet (step 510). Once the sketch is computed for a given duration, the Lp norm ∥S∥p is extracted using equation (6) (step 515).
An analysis of the performance of the estimator for the Lp norm ∥S∥p can be based on a (ε,δ) bound for p=1.

    
    
        Theorem 1. Let {right arrow over (x)}=(x1, . . . x1) be i.i.d. samples form S(1),
    
    






  
    l
    =
    
      8
      ⁢
      
        
          (
          
            
              ln
              ⁢
              
                
              
              ⁢
              2
            
            +
            
              ln
              ⁢
              
                1
                δ
              
            
          
          )
        
        /
        
          ɛ
          2
        
      
    
  
  ,
  
    ɛ
    <
    0.2
  
  ,




the DMed1=1, and Pr[median(|X1|, . . . , |Xl|) ε [1−ε, 1+ε]]>1−δ. Thus (1) gives an (ε, δ) estimator for p=1.

For example, when p=1, δ=0.05, and ε=0.1, then l=2951. This is a very loose bound in the sense that we need a very large l. Thus an asymptotic normality of the median for some approximate analysis can be used as follows.




  
    
      
        
          Theorem
          ⁢
          
            
          
          ⁢
          2.
        
        ⁢
        
          
        
      
    
    
      
        
      
    
  
  
    
      
        
          
            Let
            ⁢
            
              
            
            ⁢
            f
          
          =
          
            f
            
              S
              +
              
                (
                p
                )
              
            
          
        
        ,
        
          m
          =
          
            DMed
            p
          
        
        ,
        
          l
          =
          
            
              (
              
                
                  
                    
                       
                      z
                    
                    ⁢
                    δ
                  
                  ⁢
                  
                    /
                  
                  ⁢
                  2
                
                
                  2
                  ⁢
                  
                    mf
                    ⁡
                    
                      (
                      m
                      )
                    
                  
                  ⁢
                  ɛ
                
              
              )
            
            2
          
        
        ,
        then
      
    
    
      
        (
        1
        )
      
    
  




    
    
        gives an estimator with asymptotic (ε, δ) bound. za is the number such that for standard normal distribution Z we have Pr[Z>za]=a.
    
    


The result for this is in the same order of O(1/ε2) as the above result, but the coefficient is much smaller as shown below. For example, when p=1, δ=0.05, and ε=0.1, then m=1, f(m)=1/π,zδ/2=2, and l=986 as compared with l=2951 from above.
For p=1.05, m=0.9938, f(m)=0.3324, l=916. For p=0.95, m=1.0078, f(m)=0.3030, and l=1072. Our simulations show that these are quite accurate bounds. We can see that mƒ(m) does not change much in a small neighborhood of p=1. Since we are only interested in p in a small neighborhood of 1, for rough arguments we may use mƒ(m) at p=1, which is 1/π.
As discussed above, the entropy of a stream of flows is computed based on approximating the entropy function xln(x) by a linear combination of expressions xp, p ε (0, 2]. This approximation can be performed up to an arbitrary relative error F. To make the formula simpler natural logarithm ln x instead of log2 x, noting that changing the base is simply a matter of multiplying by the appropriate constant, thus having no effect on relative error.

    
    
        Theorem 3. For any N>1, ε>0, there exists a ε (0, 1), c=1⁄2 a ε O(ln(N)/ε1/2, such that f(x)=c(x1+ax1−a) approximates the entropy function xlnx for x ε (1, N] within relative error bound ε, i.e.,
    
    






  
    
      
      
        
          
            f
            ⁡
            
              (
              x
              )
            
          
          -
          
            x
            ⁢
            
              
            
            ⁢
            ln
            ⁢
            
              
            
            ⁢
            x
          
        
        
          x
          ⁢
          
            
          
          ⁢
          ln
          ⁢
          
            
          
          ⁢
          x
        
      
      
    
    ≤
    ɛ
  
  ,




    
    
         which can be proved using the Taylor expansion.
    
    


An estimate for the entropy norm ∥S∥H can be computed using the xln(x) estimator 230 and the Lp norm estimator 220. For example, if α and c in Theorem 1 to get relative error bound ε0 on [1,N] and if l from Theorem 1 for p=1±α is chosen to achieve an asymptotic (ε, δ) error bound, then for p=1+α there are sketches {right arrow over (Y)} and for p=1−α there are sketches {right arrow over (Z)}. The entropy norm estimator 210 estimates the entropy norm ∥S∥H the entropy norm as follows:




  
    
      
        
          
            
              
                
                S
                
              
              ⁢
              n
            
            =
          
          ^
        
        ⁢
        
          1
          
            2
            ⁢
            α
          
        
        ⁢
        
          (
          
            
              
                Λ
                ⁡
                
                  (
                  
                    Y
                    →
                  
                  )
                
              
              
                1
                +
                α
              
            
            -
            
              
                Λ
                ⁡
                
                  (
                  
                    Z
                    →
                  
                  )
                
              
              
                1
                -
                α
              
            
          
          )
        
      
    
    
      
        (
        9
        )
      
    
  




The error of this estimation is based on an assumption that ai≦N. The estimate of the entropy norm ∥S∥H is within a relative error of roughly 2λcε+λ0ε0 with probability of roughly 1-2δ,




  
    c
    =
    
      1
      
        2
        ⁢
        a
      
    
  
  ,
  
    

  
  ⁢
  
    
      λ
      0
    
    =
    
      
        
          
          
            
              
                c
                ⁡
                
                  (
                  
                    
                      
                        
                        S
                        
                      
                      
                        1
                        +
                        a
                      
                      
                        1
                        +
                        a
                      
                    
                    -
                    
                      
                        
                        S
                        
                      
                      
                        1
                        +
                        a
                      
                      
                        1
                        +
                        a
                      
                    
                  
                  )
                
              
              ⁢
              
                /
              
              ⁢
              
                
                  
                  S
                  
                
                H
              
            
            -
            1
          
          
        
        ⁢
        
          /
        
        ⁢
        
          ɛ
          0
        
      
      ≤
      1
    
  
  ,
  
    

  
  ⁢
  
    λ
    -
    
      
        
          
          S
          
        
        
          1
          +
          a
        
        
          1
          +
          a
        
      
      ⁢
      
        /
      
      ⁢
      
        
          
          S
          
        
        H
      
    
  
  ,




and λ is typically less than 1.

For example, if N=1024, a=0.05, ε=0.001, δ=0.05, then ε0=0.023, l≈105. Assuming only that λ=λ0=1, then (2cε+ε0)≈0.04, i.e. entropy norm ∥S∥H can be approximated within 4% error with 90% probability using 105 samples. Assuming λ=0.5, then ε can be increased to 0.002, and l can be decreased to approximately 2.5×104; thereby achieving 4% error.
Using the same assumptions as above with a more aggressive error calculation, equation ( ) estimates the entropy norm ∥S∥H within a relative error of roughly √{square root over (2)}cλε+λ0ε0 with probability roughly 1-δ. As such, having l≈1.25×104 achieves the 4% error with probability 95%.
As discussed above, to compute the actual entropy




  
    H
    =
    
      
        
          log
          2
        
        ⁢
        s
      
      =
      
        
          1
          s
        
        ⁢
        
          
            ∑
            i
          
          ⁢
          
            
              a
              i
            
            ⁢
            
              log
              2
            
            ⁢
            
              a
              i
            
          
        
      
    
  
  ,




the volume of the traffic s or ∥S∥1 is first computed using the packet estimator 240. Since an objective of the preferred embodiments of the present invention is to calculate the entropy of OD pairs in a network, a traffic matrix for the network is computed. This is accomplished using the same sketches that are maintained. A single stream is discussed for this example, however, this is extended to a distributed case below for multiple streams.

    
    
        Theorem 4. Let a, N, and ε be as in Theorem 1. Then, the approximation (x1+a+x1−a)/2 approximates the function ƒ(x)=x with relative error at most 3ε. This approximation holds for all counts in the range [1,N] and an elephant sketch created by the filtering unit 240 to capture flows with size strictly greater than N.
    
    


The packet estimator 250 estimates the total number of packets s or ∥S∥1 as follows:




  
    
      
        
          
            (
            
              
                
                S
                
              
              ^
            
            )
          
          1
        
        ≡
        
          
            1
            2
          
          ⁢
          
            
              (
              
                
                  
                    Λ
                    (
                    
                      Y
                      ->
                    
                    )
                  
                  
                    1
                    +
                    α
                  
                
                +
                
                  
                    Λ
                    (
                    
                      Z
                      ->
                    
                    )
                  
                  
                    1
                    -
                    α
                  
                
              
              )
            
            .
          
        
      
    
    
      
        (
        10
        )
      
    
  




Equation (10) estimates ∥S∥1 roughly within relative error bound e+3λ′0e0 with probability 1-2δ, where λ′0=|(∥S∥1+α1+α+∥S∥1−α1−α)/∥S∥1−1|/3ε0≦1. An aggressive error estimate provides an error bound of roughly




  
    1
    
      2
    
  
  ∈
  
    
      +
      3
    
    ⁢
    
      λ
      θ
      ′
    
  
  ⁢
  
    ∈
    0
  




with probability 1-δ.

The above discussion can be extended to estimating entropy of OD flow pairs in the network based on the IMP property of the sketches using sketches at ingress nodes and egress nodes to estimate the Lp norm of the OD flows. For a given OD flow pair in a network, sketches for the ingress (origin) node and egress (destination) node of the pair are used to perform an estimation of the entropy for that OD pair. Subsequently, another OD pair can be selected and the appropriate sketches can be used to estimate the entropy of that OD pair. In the discussion that follows, the ingress stream is referred to as O and the egress stream is referred to as D. As discussed above, the entropy of OD pairs can be defined by the intersection between O and D (O∩D). Estimating the Lp norms for O∩D, p=1+α, using the Lp norm estimator 220 allows for estimating the entropy of OD pairs using the approach discussed above. The sketches at ingress nodes and egress nodes are preferably implemented with an identical number of counters l and an identical set of p-stable hash functions.
With a slight abuse of notation, {right arrow over (O)} denotes a sketch for an ingress node, and {right arrow over (D)} denotes a sketch for a corresponding egress node. {right arrow over (O)}+{right arrow over (D)} and {right arrow over (O)}−{right arrow over (D)} are the component-wise addition and difference of {right arrow over (O)} and {right arrow over (D)}, respectively, since the same values of 1 are used for the sketches.
The Lp norm estimator 220 for ∥O∩D∥p; ∥, can be defined as either




  
    
      
        
          
            Λ
            (
            
              
                O
                ->
              
              ,
              
                D
                ->
              
            
            )
          
          ≡
          
            
              (
              
                
                  
                    
                      Λ
                      (
                      
                        O
                        ->
                      
                      )
                    
                    p
                  
                  +
                  
                    
                      Λ
                      (
                      
                        D
                        ->
                      
                      )
                    
                    p
                  
                  -
                  
                    
                      Λ
                      (
                      
                        
                          O
                          ->
                        
                        -
                        
                          D
                          ->
                        
                      
                      )
                    
                    p
                  
                
                2
              
              )
            
            
              1
              /
              p
            
          
        
        ⁢
        
          

        
        ⁢
        or
      
    
    
      
        (
        11
        )
      
    
  
  
    
      
        
          
            
              Λ
              ′
            
            (
            
              
                O
                ->
              
              ,
              
                D
                ->
              
            
            )
          
          ≡
          
            
              (
              
                
                  
                    
                      Λ
                      (
                      
                        
                          O
                          ->
                        
                        +
                        
                          D
                          ->
                        
                      
                      )
                    
                    p
                  
                  -
                  
                    
                      Λ
                      (
                      
                        
                          O
                          ->
                        
                        -
                        
                          D
                          ->
                        
                      
                      )
                    
                    p
                  
                
                
                  2
                  p
                
              
              )
            
            
              s
              /
              p
            
          
        
        ,
      
    
    
      
        (
        12
        )
      
    
  




which are an extension of equation (6) to accommodate an OD pair.

The equations (11) and (12) described above provide good estimators for the Lp norm ∥O∩D∥p. The Lp norm estimator 220 preferably uses two copies of one of the above equations (11) or (12), one for p1=1−α and one for p2=1+α. Once the Lp norms are computed, the xln(x) estimator 230 is used by the entropy norm estimator 210 to compute a the entropy norm, which can be used by the entropy generator 260 to approximate the entropy between some or all of the OD pairs (i.e. ingress and egress node pairs) in the network.
The flows that enter through the ingress node or exit through the egress nodes can be portioned as follows:

    
    
        A=O−D=flows that enter at ingress but do not exit through the egress
        B=O∩D=flows that enter at ingress and exit through the egress
        C=D−O=flows that do not enter at ingress but exit through the egress
    
    


Based on equation (6), Λ({right arrow over (O)})is an estimator for entropy norm ∥O∥p and Λ({right arrow over (O)})p is an estimator for ∥O∥pp=∥A∪B∥pp=∥A∥pp+∥B∥pp. Similarly A({right arrow over (D)})p is an estimator for ∥B∥pp+∥C∥pp.
The sketch {right arrow over (O)}+{right arrow over (D)} holds the contributions of the flows in A, B and C, but with every packet from B contributing twice. B(2) is used to denote all the flows in B with packet counts doubled. Then Λ({right arrow over (O)}+{right arrow over (D)})p is an estimator for ∥A∪B(2)∪C∥pp=∥A∥pp+2p∥B∥pp+∥C∥pp, since the ingress and egress nodes are using the same sketch settings, as discussed above.
The sketch {right arrow over (O)}+{right arrow over (D)} exactly cancels out the contributions from all the flows in B, and leaves the contributions of flows from A and the negative of the contributions of flows from C. C(−1) is used to denote all the flows in C with packet counts multiplied by
−1. Then Λ({right arrow over (O)}−{right arrow over (D)})p is an estimator of ∥A∪C(−1)∥pp=∥A∥pp+∥C(−1)∥pp=∥A∥pp∥C∥p p. To sum up:

    
    
        Λ({right arrow over (O)})p estimates ∥A∥pp+∥B∥pp 
        Λ({right arrow over (D)})p estimates ∥B∥pp+∥C∥pp 
        Λ({right arrow over (O)}+{right arrow over (D)})p estimates ∥A∥pp+2p∥B∥pp+∥C∥pp 
        Λ({right arrow over (O)}−{right arrow over (D)})p estimates ∥A∥pp+∥C∥pp.
    
    


From the above, it can be determined that both equation (11) and equation (12) are reasonable estimators for ∥B∥p, i.e. ∥O∩D∥p.
Suppose we have chosen proper l such that equation (11) is roughly an (ε, δ) estimator. Suppose ∥O∩D∥pp=r1∥O∥pp, and ∥O∩D∥pp=r2∥D∥pp. Then equation (11) raised to the power p estimates ∥O∩D∥pp roughly within relative error bound




  
    (
    
      
        1
        
          r
          1
        
      
      +
      
        1
        
          r
          2
        
      
      -
      1
    
    )
  
  ∈




with probability at least 1-3δ. Similarly, equation (12) raised to the power p estimates ∥O∩D∥pp roughly within relative error bound





  
    
      2
      
        1
        -
        p
      
    
    ⁢
    
      (
      
        
          1
          
            r
            1
          
        
        +
        
          1
          
            r
            2
          
        
        +
        
          2
          
            p
            -
            r
          
        
        -
        2
      
      )
    
  
  ∈
  
    ≈
    
      (
      
        
          1
          
            r
            1
          
        
        +
        
          1
          
            r
            2
          
        
        -
        1
      
      )
    
  
  ∈




with probability at least 1-2δ. This gives a very loose rough upper bound on the relative error. The ratios r1 and r2 are related to, but not identical to, the ratio of OD flow traffic against the total traffic at the ingress and egress points.

To calculate OD flow entropy, Λ({right arrow over (Y)}) is replaced in equation (6) and (9) with Λ({right arrow over (O)}, {right arrow over (D)}) where {right arrow over (O)} and {right arrow over (D)} are L1+αsketches, and Λ({right arrow over (Z)}) is similarly replaced.
Incoming packets preferably trigger increments to all l counters for estimating one Lp norm, and the entropy estimator 200 uses two separate Lp norms computations. Such a large l is generally unacceptable for networking applications. For high-speed links, where each packet has tens of nanoseconds to process, it is generally not feasible to increment that many counters per packet, even if they are all in fast SRAM.
A methodology of bucketing is used to resolve feasibility concerns. With bucketing, a sketch becomes a two-dimensional array M[1 . . . k][1 . . . l]. For incoming packets, flow labels are hashed using a uniform hash, function uh, and the result uh(pkt.id) is the index of the bucket at which the packet should be processed. Then the counters are incremented M[uh(pkt.id)][1 . . . l] an a similar manner as discussed above with respect to FIG. 5. Finally, the Lp norm estimations are added up from all these buckets to obtain a final estimate.
Bucketing (i.e., with k buckets instead of 1) can reduce the standard deviation of the estimation of Lp norms by a factor slightly less than √{square root over (k)}, provided that the number of flows is much larger than the number of buckets k. A decrease in l can be compensated by an increase in k by a slightly larger factor. In a preferred implementation, the number of buckets k is typically on the order of 10,000. Such a large bucket size allows l to shrink to a small number such as 20 to achieve the same (or even better) estimation accuracy. Even on very high-speed links (e.g., 10M packets per second), a few tens of memory accesses per packet can be accommodated.
B is used to denote all the flows hashed to the ith bucket, and ∥B∥p, to denote the Lp norm of the flows in the ith bucket, similar to the definition of Lp norm ∥S∥p before (∥S∥pp=Σi∥Bi∥pp). Let Mi be the ith row vector of the sketch. Λ(Mi) as defined in equation (6) is an estimator of ∥Bi∥p, so naturally the estimator for ∥S∥p is:




  
    
      
        (
        
          
            
              
              S
              
            
            p
            p
          
          =
          
            
              [
              
                
                  ∑
                  
                    i
                    =
                    1
                  
                  k
                
                ⁢
                
                  
                    (
                    
                      Λ
                      ⁡
                      
                        (
                        
                          M
                          i
                        
                        )
                      
                    
                    )
                  
                  p
                
              
              ]
            
            
              1
              /
              p
            
          
        
        )
      
    
    
      
        (
        13
        )
      
    
  




In an ideal case of even distribution of flows into buckets, and all ∥Bi∥pp are the same, then ∥S∥pp=k∥Bi∥pp. Consider p=1, then the estimator Λ(Mi) is roughly Gaussian with mean ∥Bi∥p and standard deviation (1⁄2mƒ(m)√{square root over (l)})∥Bi∥1. By the Central Limit Theorem,




  
    Λ
    ⁡
    
      (
      M
      )
    
  
  =
  
    
      ∑
      
        i
        =
        1
      
      k
    
    ⁢
    
      Λ
      ⁡
      
        (
        
          M
          i
        
        )
      
    
  




is asymptotically Gaussian, and its standard deviation is roughly √{square root over (k)}(1⁄2mƒ(m)√{square root over (l)})∥B1∥1=(1⁄2mƒ(m)√{square root over (lk)})∥S∥1. If equation (6) was used without bucketing, the standard deviation would be (1⁄2mƒ(m)√{square root over (l)})∥S∥1. So lk is performing the same role as l did when bucketing is not used, or in other words, k buckets reduces standard error roughly by a factor of √{square root over (k)}.

Generally there is not even distribution of flows into various buckets. However, when the number of flows is far larger than the number of buckets, which will be the case with parameter settings and intended workload of the preferred embodiments of the present invention, the factor of error reduction is only slightly less than √{square root over (k)}. When k is increased to be on the same order of the number of flows, however, the factor of error reduction typically no longer scale as √{square root over (k)} since there can be many empty buckets that do not contribute to the reduction of estimation error, and since there can be a more uneven distribution of flows into buckets. Therefore, when l is fixed, the estimation error cannot generally be brought down arbitrarily close to 0 by increasing k arbitrarily.
Another issue is the bias of median estimator in equation (6), that is, the expected value of the sample median of l samples is typically not equal to the distribution median DMedp. When buckets are not implemented, the asymptotic normality implies that the bias is much smaller than the standard error, so we could ignore the issue. However, when buckets are implement (e.g., k buckets) to reduce the standard error by a factor of √{square root over (k)}, the bias becomes significant. Let EMedp,l denote the expected value of the median of l samples from distribution S+(p). So the estimator for ∥Biμp is redefined as follows




  
    
      
        
          
            Λ
            →
          
          ⁡
          
            (
            
              M
              i
            
            )
          
        
        ≡
        
          
            
              median
              ⁡
              
                (
                
                  
                    
                    
                      
                        M
                        ⁡
                        
                          [
                          i
                          ]
                        
                      
                      ⁡
                      
                        [
                        1
                        ]
                      
                    
                    
                  
                  ,
                  ...
                  ⁢
                  
                    
                  
                  ,
                  
                    
                    
                      
                        M
                        ⁡
                        
                          [
                          i
                          ]
                        
                      
                      ⁡
                      
                        [
                        l
                        ]
                      
                    
                    
                  
                
                )
              
            
            
              EMed
              
                p
                ,
                l
              
            
          
          .
        
      
    
    
      
        (
        14
        )
      
    
  




Equation (14) replaces Λ(Mi) in (13) and gives Lp norm estimation used be the Lp norm estimator 220 when buckets are used. Note that equation (14) is an unbiased estimator, but equation (13) still may be biased.
It can be assumed that M and N are Lp sketches with buckets at ingress node O and egress node D respectively, and the two sketches use the same settings. {right arrow over (O)} and {right arrow over (D)} in equations (11) and (12) can be replaced with M and N, which as discussed above provide estimators of OD flow Lp norms.
EMedp,l can be numerically calculated using the p.d.f. formula for order statistics when ƒs(p) has closed form. Or it can be derived via simulation. VEMedp,l denotes the variance of the sample median of l samples. For example, when p=1, l=20, then EMed1,20=1.069, VMed1,20=0.149, so the standard deviation is 0.386. The distribution median is DMed1=1, so we can see the bias 0.069 is much smaller than the standard deviation 0.386. Also the asymptotic standard deviation is 0.351, which is close to the actual value of 0.385.
FIG. 6 is a flowchart showing computation of the Lp norm using the bucketing technique in accordance with a preferred embodiment of the present invention. A sketch M[1 . . . k][1 . . . l] is initialized to zero (step 600). The p-stable hash functions sh1 through shl are fixed (step 605). The hash function uh maps flow labels to buckets { 1, . . . , k} (step 610) and the expected sample median EMedp,l is calculated (step 615). For incoming packets (pkt), increment counters in the sketch for j from 1 to l(M[uh(pkt.id)][j]+=shj(pkt.id)) (step 620). Finally, the Lp norm is generated based on equation 14 (step 625).
FIG. 7 depicts an exemplary computing device 700 for implementing entropy estimation of OD flows in a network in accordance with the preferred embodiments of the present invention. The computing device 700 can be a mainframe, personal computer (PC), laptop computer, workstation, handheld device, such as a PDA, or the like. In the illustrated embodiment, the computing device 700 includes a central processing unit (CPU) 702 and preferably a display device 704. The display device 704 enables the computing device 700 to communicate directly with a user through a visual display. The computing device 700 can further include data entry device(s) 706, such as a keyboard, touch screen, and/or mouse.
The computing device 700 can include storage 708 for storing data, such as sketches, Lp norms ∥S∥p, total number of packets s, entropy norms ∥S∥H, entropy, etc., and instructions, such as instruction for computing the Lp norms ∥S∥p, total number of packets s, entropy norms ∥S∥H, entropy, etc. The storage 708 can include such technologies as a floppy drive, hard drive, tape drive, Flash drive, optical drive, read only memory (ROM), random access memory (RAM), and the like.
Applications 710, such as the entropy estimator 200 or components thereof for performing the process described above, can be resident in the storage 708 or may be distributed among multiple computing devices. The storage 708 can be local or remote to the computing device 700 and can include a database 720 to which the computing device 700 can connect via the network 100 or another network. The computing device 700 includes a network interface 712 for communicating over the network 100 to query and receive trace information from the nodes of the network 100 and/or to retrieve or store information in the database 720. The CPU 702 operates to run the application in storage 708 by performing instructions therein and storing data resulting from the performed instructions, which may be depicted via the display 704 or by other mechanisms known to those skilled in the art, such a print out from a printer.
The preferred embodiments can work with high link speeds of up to 10 million packets per second using commodity CPU/memory at a reasonable cost. As discussed above, two sub-sketches for estimating the L1+a and the L1−a norms respectively, each of which preferably consists of k*l real-valued counters. The number of counters per bucket l is set to 20 in a preferred implementation. Since the sketches can implemented using inexpensive high-throughput DRAM, the number of buckets k can be very large (say up to millions).
Each incoming packet, l=20 counters are incremented per sketch (two sub-sketches). Single-precision real number counters (4 bytes each) are used to minimize memory I/O (space not an issue), as 7 decimal digits of precision is accurate enough for the preferred embodiments. This involves 320 bytes of memory reads or writes, since each counter increment preferably involves a memory read and a memory write. Generating realizations of p-stable distributions from two pre-computed tables in order to compute sh preferably involves another 320 bytes of memory reads. In total, each incoming packet preferably triggers 640 bytes of memory reads/writes. Using commodity RDRAM DIMM 6400, the sketches of the preferred embodiments can deliver a throughput of 10 million packets per second.
The sketches can be implemented using RDRAM DIMM 6400 (named after its 6400 MB/s sustained throughput for burst accesses), except that the elephant detection module is preferably implemented using a small amount of SRAM. RDRAM can deliver a very high throughput for read/write in burst mode (a series of accesses to consecutive memory locations). Since 640 bytes of memory accesses triggered by each incoming packet consist of 4 large contiguous blocks of 160 bytes each, the 6400 MB/s through-put provided by RDRAM DIMM 6400 can be used efficiently. Memory space/consumption is generally not a concern, since the sketches are implemented using DRAM. Even if we need millions of buckets, only hundreds of MBs are consumed.
A packet monitor was deployed on a 1 Gbit/second ingress link from a data center into a large tier-1 ISP backbone network. The data center hosts tens of thousands of Web sites as well as servers providing a wide range of services such as multimedia, mail, etc. The link monitored is one of multiple links connecting this data center to the Internet. All the traffic carried by this link enters the ISP backbone network at the same ingress router. For each observed packet, the monitor captured its IP header fields as well as UDP/TCP and ICMP information for defining flows.
Five-minute traces and a one-hour trace were collected. A routing table dumped at the ingress router was used to determine the egress router for each packet, thus determining the OD flows to each possible egress router. Because the packet monitoring capability implemented to gather data was not available at egress routers, synthetic traffic traces are generated at egress routers so that they contain corresponding OD flows observed at the ingress router and the flow size distribution at egress routers are dictated so that they preferably match the flow size distribution of the ingress trace.
The following two traces are used to illustrate results for the entropy estimation approach described herein.

    
    
        Trace 1: A one hour trace collected that contains 0.4 billion packets which belong to 1.8 million flows. We chose one egress router so that the traffic between the origin and destination comprised of 5% of the total traffic arriving at the ingress router.
        Trace 2: A five minute trace that is similar to Trace 1, where traffic between the origin and our chosen egress router comprised of 5% of the total traffic arriving at the ingress router. The traffic in this trace is purposely chosen as being a subset of the traffic for Trace 1 so that we can directly compare the performance of entropy estimator 200 for five minutes and one hour intervals.
    
    


For each trace, the experiment is repeated 1000 times with independently generated sketch data structures for which the cumulative density function of the relative error are computed. Unless stated otherwise, the parameters we used for each experiment were: number of buckets k=50, 000, number of registers in each bucket l=20, sample and hold sampling rate P=0.001, and one million entries in each lookup table.
An elephant flow (for which the contribution to the entropy are computed separately) is defined to be any flow with at least N=1000 packets in this illustrative implementation. α=0.05 is used, which satisfies the constraint α<1/ln N. Hence, at every ingress and egress point we have a pair of sketches computing the L1.05 and L0.95 norms of the traffic.

    
    
        Results from equations (11) and (12) are compared using the sketches created for the above traces. The cumulative density plots for the error of the entropy estimator 200 using equations (11) and (12) for Trace 1 are given in FIG. 8. Implementations using both equations have reasonably small and comparable error values. This also holds on all five minute traces and hence we fix and use equation (11) for the remaining evaluation.
    
    


The effect of varying the number of buckets on the performance of the entropy estimator is studied. Keeping all other parameters fixed at reasonable values, the number of buckets is varied between k=5000 and k=80000, with increasing factors of two. FIG. 9 shows the results of Trace 2. It is shown that increasing the number of buckets increases the accuracy (as expected), but with diminishing returns.
As discussed above, the entropy computation large (elephant) flows are separated by means of a sample and hold technique. By varying the sampling rate the probability with which flows that are above our elephant threshold (i.e., flows of size greater than 1000) are sampled can increase or decrease. It was shown that the sampling rate did not affect the performance of our algorithm significantly. FIG. 10 shows that, even with a small sampling rate (e.g., 1 in 1000), the elephant detection mechanism allows good overall performance for our algorithm.
FIG. 11 compares the cumulative density plots of the error for the five minute, Trace 2, and the one hour, Trace 1, which have the same origin and destination nodes and similar traffic distributions. Even though there is an order of magnitude difference in the size of these traces, not only does the error remain comparable, but also does the distribution of the error. Different trace sizes show that the entropy estimator 200 is robust to changes in the size of traces as long as the fraction of cross-traffic is held constant.
The performance of the entropy estimator 200 in computing the actual entropy (as opposed to the entropy norm) of the OD flows is evaluated. The computation of entropy can have additional error since we need to make use of sketches to estimate the total volume of traffic between a source and destination. In FIG. 12 the error plot for the entropy is comparable to that for the entropy norm. Hence, this confirms that the entropy estimator 200 is a robust estimator of the entropy of OD flows.
While preferred embodiments of the present invention have been described herein, it is expressly noted that the present invention is not limited to these embodiments, but rather the intention is that additions and modifications to what is expressly described herein are also included within the scope of the invention. Moreover, it is to be understood that the features of the various embodiments described herein are not mutually exclusive and can exist in various combinations and permutations, even if such combinations or permutations are not made explicit herein, without departing from the spirit and scope of the invention.

## Claims

1. A method of estimating entropy in a network comprising:
generating a first sketch associated with at least one ingress flow in a network by a first router;
generating a second sketch associated with at least one egress flow in the network, the at least one ingress flow and the at least one egress flow forming an origin-destination (OD) flow pair by a second router; and
determining an OD flow entropy based on the first and second sketches using an entropy estimator implemented by a computing device,
wherein determining an OD flow entropy comprises:
extracting a first Lp norm from the first sketch and a second Lp norm from the second sketch; and
determining an intersection between the first and second sketches using the Intersection Measurable Property (IMP).

2. The method of claim 1, wherein determining an OD flow entropy comprises determining an entropy norm using the intersection between the first and second sketches and an approximation for x*ln(x) for all x from about 1 to about 1000.

3. The method of claim 1, wherein determining an OD flow entropy comprises determining a total number of packets in a stream using at least one of the first and second sketches.

4. The method of claim 3, wherein determining an origin-destination (OD) flow entropy comprises computing an estimate for the OD flow entropy using the entropy norm and the total number of packets.

5. The method of claim 1, wherein generating a first sketch comprises:
passing a packet to a bucket based on a flow label associated with the packet; and
incrementing at least one counter in the first sketch based on a packet identification associated with the packet.

6. A non-transitory computer-readable medium comprising instructions executable by a computing device for estimating entropy in a network by:
generating a first sketch associated with at least one ingress flow in a network;
generating a second sketch associated with at least one egress flow in the network, the at least one ingress flow and the at least one egress flow forming an origin-destination (OD) flow pair; and
determining an OD flow entropy based on the first and second sketches,
wherein determining an OD flow entropy comprises:
extracting a first Lp norm from the first sketch and a second Lp norm from the second sketch; and
determining an intersection between the first and second sketches using the Intersection Measurable Property (IMP).

7. The medium of claim 6, wherein determining an OD flow entropy comprises determining an entropy norm using the intersection between the first and second sketches and an approximation for x*ln(x) for all x from about 1 to about 1000.

8. The medium of claim 6, wherein determining an OD flow entropy comprises determining a total number of packets in a stream using at least one of the first and second sketches.

9. The medium of claim 8, wherein determining an origin-destination (OD) flow entropy comprises computing an estimate for the OD flow entropy using the entropy norm and the total number of packets.

10. The medium of claim 6, wherein generating a first sketch comprises:
passing a packet to a bucket based on a flow label associated with the packet; and
incrementing at least one counter in the first sketch based on a packet identification associated with the packet.

11. A system for estimating entropy in a network comprising:
a computing device configured to generate a first sketch associated with at least one ingress flow in a network, generate a second sketch associated with at least one egress flow in the network, the at least one ingress flow and the at least one egress flow forming an origin-destination (OD) flow pair, and to determine an OD flow entropy based on the first and second sketches,
wherein the computing device is configured to extract a first Lp norm from the first sketch and a second Lp norm from the second sketch and determine an intersection between the first and second sketches the Intersection Measurable Property (IMP).

12. The system of claim 11, wherein the computing device is configured to determine an entropy norm using the intersection between the first and second sketches and an approximation for x*ln(x) for all x from about 1 to about 1000.

13. The system of claim 11, wherein the computing device is configured to determine a total number of packets in a stream using at least one of the first and second sketches.

14. The system of claim 13, wherein the computing device is configured to compute an estimate for the OD flow entropy using the entropy norm and the total number of packets.

15. The system of claim 11, wherein the computing device is configured to pass a packet to a bucket based on a flow label associated with the packet and increment at least one counter in the first sketch based on a packet identification associated with the packet.