A HYBRID GENETIC ALGORITHM FOR THE WEIGHT SETTING PROBLEM IN OSPF/IS-IS ROUTING
|
|
- Conrad Stevens
- 6 years ago
- Views:
Transcription
1 A HYBRID GENETIC ALGORITHM FOR THE WEIGHT SETTING PROBLEM IN OSPF/IS-IS ROUTING L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP Abstract. Intra-domain traffic engineering aims to make more efficient use of network resources within an autonomous system. Interior Gateway Protocols such as OSPF (Open Shortest Path First) and IS-IS (Intermediate System- Intermediate System) are commonly used to select the paths along which traffic is routed within an autonomous system. These routing protocols direct traffic based on link weights assigned by the network operator. Each router in the autonomous system computes shortest paths and creates destination tables used to direct each packet to the next router on the path to its final destination. Given a set of traffic demands between origin-destination pairs, the OSPF weight setting problem consists in determining weights to be assigned to the links so as to optimize a cost function, typically associated with a network congestion measure. In this paper, we propose a genetic algorithm with a local improvement procedure for the OSPF weight setting problem. The local improvement procedure makes use of an efficient dynamic shortest path algorithm to recompute shortest paths after the modification of link weights. We test the algorithm on a set of real and synthetic test problems and show that it produces near-optimal solutions. We compare the hybrid algorithm with other algorithms for this problem illustrating its efficiency and robustness. 1. Introduction The Internet is divided into many routing domains, called autonomous systems (ASes). These ASes interact to control and deliver IP traffic. They typically fall under the administration of a single institution, such as a company, a university, or a service provider. Neighboring ASes use the Border Gateway Protocol (BGP) to route traffic [21]. The goal of intra-domain traffic engineering [9] consists in improving user performance and making more efficient use of network resources within an AS. Interior Gateway Protocols (IGPs) such as OSPF (Open Shortest Path First) and IS-IS (Intermediate System-Intermediate System) are commonly used to select the paths along which traffic is routed within an AS. These routing protocols direct traffic based on link weights assigned by the network operator. Each router in the AS computes shortest paths and creates destination tables used to direct each IP packet to the next router on the path to its final destination. OSPF calculates routes as follows. To each link is assigned an integer weight ranging from 1 to (= ). The weight of a path is the Date: June 24, Key words and phrases. OSPF routing, IS-IS routing, Internet, metaheuristics, genetic algorithm, optimized crossover, local search. This research was done while the first author was a visiting scholar at the Internet and Network Systems Research Center at AT&T Labs Research. AT&T Labs Research Technical Report TD-5NTN5G. 1
2 2 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP sum of the link weights on the path. OSPF mandates that each router computes a graph of shortest paths with itself as the root [13]. This graph gives the least weight routes (including multiple routes in case of ties) to all destinations in the AS. In the case of multiple shortest paths originating at a router, OSPF is usually implemented so that it will accomplish load balancing by splitting the traffic flow over all shortest paths leaving from each router [16]. In this paper, we consider that traffic is split evenly between all outgoing links on the shortest paths to the destination IP address. OSPF requires routers to exchange routing information with all the other routers in the AS. Complete network topology knowledge is required for the computation of the shortest paths. Given a set of traffic demands between origin-destination pairs [8], the OSPF weight setting problem consists in determining weights to be assigned to the links so as to optimize a cost function, typically associated with a network congestion measure. The NP-hardness of the OSPF weight setting problem was established in [10]. Previous work on optimizing OSPF weights have either chosen weights so as to avoid multiple shortest paths from source to destination or applied a protocol for breaking ties, thus selecting a unique shortest path for each source-destination pair [19, 14, 3]. Fortz and Thorup [10] were the first to consider even traffic splitting in OSPF weight setting. They proposed a local search heuristic and tested it on a realistic AT&T backbone network and on synthetic networks. Ericsson, Resende, and Pardalos [7] proposed a genetic algorithm and used the set of test problems considered in [10]. Sridharan, Guérin, and Diot [20] developed another heuristic for a slightly different version of the problem, in which flow is split among a subset of the outgoing links on the shortest paths to the destination IP address. In this paper, we propose a hybrid genetic algorithm incorporating a local improvement procedure to the crossover operator of the genetic algorithm proposed in [7]. The local improvement procedure makes use of an efficient dynamic shortest path algorithm to recompute shortest paths after the modification of link weights. We compare the hybrid algorithm with the genetic algorithm as well as with the local search procedure in [10]. In the next section, we give the mathematical formulation of the OSPF weight setting problem. The hybrid genetic algorithm is described in Section 3 and the local improvement procedure in Section 4. Section 5 describes efficient algorithms for solution update used in the local improvement procedure. Computational results are reported in Section 6. Concluding remarks are made in the last section. 2. Problem formulation In a data communication network, nodes and arcs represent routers and transmission links, respectively. Let N and A denote, respectively, the sets of nodes and arcs. Data packets are routed along links, which have fixed capacities. Consider a directed network graph G = (N, A) with a capacity c a for each a A, and a demand matrix D that, for each pair (s, t) N N, gives the demand d st in traffic flow from node s to node t. Then, the OSPF weight setting problem consists in assigning positive integer weights w a [1, w max ] to each arc a A, such that a measure of routing cost is optimized when the demands are routed according to the rules of the OSPF protocol. The OSPF protocol allows for w max
3 HYBRID GA FOR OSPF/IS-IS ROUTING cost per unit of capacity trunk utilization rate Figure 1. Piecewise linear function Φ a (l a ). For each pair (s, t) and each arc a, let f a (st) indicate how much of the traffic flow from s to t goes over arc a. Let l a be the total load on arc a, i.e. the sum of the flows going over a, and let the trunk utilization rate u a = l a /c a. The routing cost in each arc a A is taken as the piecewise linear function Φ a (l a ), proposed by Fortz and Thorup [10] and depicted in Figure 1, which increasingly penalizes flows approaching or violating the capacity limits: (1) Φ a (l a ) = u a, u a [0, 1/3) 3 u a 2/3, u a [1/3, 2/3), 10 u a 16/3, u a [2/3, 9/10), 70 u a 178/3, u a [9/10, 1), 500 u a 1468/3, u a [1, 11/10), 5000 u a 16318/3, u a [11/10, ). OSP F (w) Given a weight assignment w and the loads la associated with each arc a A corresponding to the routes obtained with OSPF, we denote its routing cost by Φ OSP F (w) = a A Φ OSP F (w) a(la ). The OSPF weight setting problem is then equivalent to finding arc weights w [1, w max ] such that Φ OSP F (w) is minimized. The general routing problem can be formulated as the following linear programming problem with a piecewise linear objective function: (2) Φ OP T = min Φ = a A Φ a (l a )
4 4 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP subject to (3) (4) (5) (6) (7) (8) (9) (10) (11) u:(u,v) A f (st) (u,v) u:(v,u) A f (st) (v,u) = d st if v = s, d st if v = t, 0 otherwise, l a = (s,t) N N v, s, t N, f (st) a, a A, Φ a (l a ) l a, a A, Φ a (l a ) 3l a 2/3c a, a A, Φ a (l a ) 10l a 16/3c a, a A, Φ a (l a ) 70l a 178/3c a, a A, Φ a (l a ) 500l a 1468/3c a, a A, Φ a (l a ) 5000l a 16318/3c a, a A, f (st) a 0, a A; s, t N. Constraints (3) are flow conservation constraints that ensure routing of the desired traffic. Constraints (4) define the load on each arc a and constraints (5 10) define the cost on each arc a according to the cost function Φ a (l a ). The above is a relaxation of OSPF routing, as it allows for arbitrary routing of traffic. Then, Φ OP T is a lower bound on the optimal OSPF routing cost Φ OSP F (w ). Also, if Φ OSP F (1) denotes the optimal OSPF routing cost when unit weights are used, then Φ OSP F (w ) Φ OSP F (1). Fortz and Thorup [10] proposed a normalizing scaling factor for the routing cost which makes possible comparisons across different network sizes and topologies: Φ UNCAP = d st h st, (s,t) N N where h st is the minimum hop count between nodes s and t. For any routing cost Φ, the scaled routing cost is defined as Φ = Φ/Φ UNCAP. Using this notation, the following results hold: The optimal routing costs satisfy 1 Φ OP T Φ OSP F (w ) Φ OSP F (1) Given any solution to (2-11) with normalized routing cost Φ, then Φ = 1 if and only if all arc loads are below 1/3 of their capacities and all demands are routed on minimum hop routes. Given any solution to (2-11) where all arcs are at their maximum capacity, then the normalized routing cost Φ = We say that a routing congests a network if Φ Hybrid genetic algorithm for OSPF weight setting In this section, we summarize the detailed description of the genetic algorithm given in [7] and propose a hybrid genetic algorithm by adding a local improvement procedure after the crossover. A genetic algorithm is a population-based metaheuristic for combinatorial optimization. In this context, a population is simply a set of feasible solutions. Solutions
5 HYBRID GA FOR OSPF/IS-IS ROUTING 5 in a population are combined (through crossover) and perturbed (by mutation) to produce a new generation of solutions. When solutions are combined, attributes of higher-quality solutions have a greater probability to be passed down to the next generation. This process is repeated over many generations as long as the quality of the solutions in the new population improves over time. We next show how this idea can be explored for weight setting in OSPF routing. Each solution is represented by an array of integer weights, where each component corresponds to the weight of an arc of the network. Each individual weight belongs to the interval [1, w max ]. Each solution w is associated with a fitness value defined by the OSPF routing cost Φ OSP F (w). The initial population is randomly generated, with arc weights selected from a uniform distribution in the interval [1, w max /3]. The population is partitioned into three sets A, B, and C. The best solutions are kept in A, while the worst ones are in C. All solutions in A are promoted to the next generation. Solutions in B are replaced by crossover of one parent from A with another from B C using the random keys crossover scheme of Bean [2]. All solutions in C are replaced by new randomly generated solutions with arc weights selected in the interval [1, w max ]. In the random keys scheme, crossover is carried out on a selected pair of parent solutions to produce an offspring solution. Each selected pair consists of an elite parent and a non-elite parent. The elite parent is selected, at random, uniformly from solutions in set A, while the non-elite parent is selected, at random, uniformly from solutions in set B C. Each weight in the offspring solution is either inherited from one of its parents or is reset by mutation. With mutation probability p m, the weight is reset to a value selected at random in the interval [1, w max ]. If mutation does not occur, then the child inherits the weight from its elite parent with a given probability p A > 1/2. Otherwise, it inherits the weight from its non-elite parent. The hybrid genetic algorithm proposed in this paper, applies a local improvement procedure to each offspring solution obtained by crossover. This local improvement procedure is described in the next section. 4. Local Improvement Procedure In this section, we describe the local improvement procedure. Starting from an initial solution, the local improvement procedure analyzes solutions in the neighborhood of a current solution w in the search for a solution having a smaller routing cost. If such a solution exists, then it replaces the current solution. Otherwise, the current solution is returned as a local minimum. The local improvement procedure is incorporated in the genetic algorithm, described in Section 3, to enhance its ability to find better-quality solutions with less computational effort. Local improvement is applied to each solution generated by the crossover operator. Besides being computationally demanding, the use of large neighborhoods in a hybrid genetic algorithm can lead to loss of population diversity, and consequently premature convergence to low-quality local minima. We next describe the local improvement procedure using a reduced neighborhood. As before, let l a denote the total load on arc a A in the solution defined by the current weight settings w. We recall that Φ a (l a ) denotes the routing cost on this arc. The local improvement procedure examines the effect of increasing the weights of a subset of the arcs. These candidate arcs are selected among those with the highest routing costs and whose weight is smaller than w max. To reduce the routing
6 6 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP procedure LocalImprovement(q, w) 1 dontlook a 0, a A; 2 i 1; 3 while i q do 4 Renumber the arc indices such that Φ a (l a ) Φ a+1 (l a+1 ), a = 1,..., A 1; 5 a 0; 6 for a = 1,..., A while a = 0 do 7 if dontlook a = 1 then dontlook a 0; 8 else if w a < w max then a a; 9 end for; 10 if a = 0 then return; 11 dontlook a 1; 12 for ŵ = w a + 1,..., w a + (w max w a )/4 do 13 w a w a, a A, a a ; 14 w a ŵ; 15 if Φ OSP F (w ) < Φ OSP F (w) then 16 w w ; 17 dontlook a 0; 18 i 0; 19 end if 20 end for 21 i i + 1; 22 end while end LocalImprovement. Figure 2. Pseudo-code of procedure LocalImprovement. cost of a candidate arc, the procedure attempts to increase its weight to induce a reduction on its load. If this leads to a reduction in the overall routing cost, the change is accepted and the procedure is restarted. The procedure stops at a local minimum when no improvement results from changing the weights of the candidate arcs. The pseudo-code in Figure 2 describes the local improvement procedure in detail. The procedure LocalImprovement takes as input parameters the current solution defined by the weights w and a parameter q which specifies the maximum number of candidate arcs to be examined at each local improvement iteration. To speed up the search, we disallow the weight increase of arcs for which no weight increase leads to an improvement in the routing cost in the previous iteration. To implement this strategy, we make use of a don t look bit for each arc. The don t look bits are set unmarked in line 1 and the counter of candidate arcs is initialized in line 2. The loop in lines 3 to 22 investigates at most q selected candidate arcs for weight increase in the current solution. The arc indices are renumbered in line 4 such that the arcs are considered in non-increasing order of routing cost. The loop in lines 6 to 9 searches for an unmarked arc with weight less than w max. Marked arcs which cannot be selected at the current iteration are unmarked in line 7 for future investigation. Arc a is selected in line 8. If no arc
7 HYBRID GA FOR OSPF/IS-IS ROUTING 7 satisfying these conditions is found, the procedure stops in line 10 returning the current weights w as the local minimum. In line 11, arc a is temporarily marked to disallow its investigation in the next iteration, unless a weight change in w a results in a better solution. The loop in lines 12 to 20 examines all possible weight changes for arc a in the range [w a + 1, w a + (w max w a )/4 ]. A neighbor solution w, keeping all arc weights unchanged except for arc a, is built in lines 13 and 14. If the new solution w has a smaller routing cost than the current solution (test in line 15), then the current solution is updated in line 16, arc a is unmarked in line 17, and the arc counter i is reset in line 18. In line 21, we increment the candidate arc counter i. The routing cost Φ OSP F (w ) associated with the neighbor solution w must be evaluated in line 15. Instead of computing it from scratch, we use fast update procedures for recomputing the shortest path graphs as well as the arc loads. These procedures are considered in the next section of the paper. Once the new arc loads are known, the total routing cost is computed as the sum of the individual arc routing costs. 5. Fast updates of arc loads and routing costs In this section, we describe the procedures used for fast update of the cost (line 15 of procedure LocalImprovement) and arc loads. We are in the situation where l are the loads associated with the current weight settings w and the weight of a unique arc a is increased by exactly a unit. Let T be the set of destination nodes and denote by g t = (N, A t ) the shortest path graph associated with each destination node t T. One or more of these shortest path graphs will be affected by the change of the weight of arc a from w a to w a + 1. Consequently, the loads of some of the arcs in each graph will change. The procedures described in this section use a set of data structures that work like a memory for the solution. With a weight change, the shortest path graph and the loads can change, and the memories are updated instead of recomputed from scratch. Each shortest path graph with node t as destination has an A -vector A t indicating the arcs in g t. If arc a is in the shortest path graph, then A t a = 1. Otherwise, A t a = 0. Another A -vector, l t, associated with the arcs, stores the partial loads flowing to t traversing each arc a A. The total load from each arc is represented in the A -vector l a which stores the total load traversing each arc a A. The N -vectors π t and δ t are associated with the nodes. The distance from each node to the destination t is stored in π t, while δ t keeps the number of arcs outgoing from each node in g t. All these structures are populated in the beginning and set free at the end of procedure LocalImprovement. For simplicity, they were omitted from this procedure and from the parameter list of the procedures described in this section. The pseudo-code in Figure 3 summarizes the main steps of the update procedure. The new load l a on each arc a A is set to zero in line 1. The loop in lines 2 to 6 considers each destination node t T. For each one of them, the shortest path graph g t is updated in line 3 and the partial arc loads are updated in line 4. The arc loads are updated in line 5. Lines 1 and 5 are removed from this procedure in case the arc loads be updated inside procedure UpdateLoads. Finally, the cost Φ OSP F (w ) of the new solution is computed in line 7. In the remainder of this
8 8 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP procedure UpdateCost(a, d, l) 1 forall a A do l a 0; 2 forall t T do 3 UpdateShortestPaths(a, w, t); 4 UpdateLoads(d, t); 5 forall a A t do l a l a + l t a; 6 end forall 7 Φ OSP F (w ) a A Φ a(l a ); end UpdateCost. Figure 3. Pseudo-code of procedure UpdateCost. section we describe the procedures UpdateShortestPaths and UpdateLoads used in lines 3 and Dynamic reverse shortest path algorithm. We denote by g t = (N, A t ) the shortest paths graph associated with each destination node t T. Since the weight of a unique arc a was changed, the graph g t does not have to be recomputed from scratch. Instead, we update the part of it which is affected by the weight change. Ramalingam and Reps [17] and Frigioni et al. [12] proposed efficient algorithms for these dynamic computations in Dijkstra s algorithm. These two algorithms are compared experimentally in [11]. Although the algorithm of Frigioni et al. is theoretically better, the algorithm of Ramalingam and Reps usually runs faster in practice. Due to the nature of the OSPF weight setting problem, we use the reversed version of Dijkstra s shortest path algorithm. The pseudo-code of the specialized dynamic shortest path algorithm for unit weight increases is given in Figure 4. Buriol et al. [4] showed empirically that this algorithm is faster than the general dynamic reverse shortest path algorithm of Ramalingam and Reps [17]. First, it identifies the set Q of nodes whose distance labels change due to the increased weight. Next, the shortest paths graph is updated by deleting and adding arcs for which at least one of its extremities belongs to Q. Algorithm UpdateShortestPaths takes as parameters the arc a = ( u, v) whose weight changed, the current setting of weights w, and a destination node t T. The algorithm checks in line 1 if arc a does not belong to the shortest paths graph g t, in which case the weight change does not affect the latter and the procedure returns. Arc a is eliminated from the shortest paths graph g t in line 2. In line 3 the tail node of arc a is inserted in a heap containing all nodes for which the load of its outgoing links might change. This heap will be used later by the load update procedure. The distances to the destination node t are used as priority keys and the root contains the node with maximum distance. The outdegree δ u + of the tail node of arc a is updated in line 4. In line 5, we check if there is an alternative path to the destination starting from u. If this is the case, the procedure returns in line 5, since no further change is needed. The set Q of nodes affected by the weight change in arc a is initialized with node u in line 6. The loop in lines 7 to 15 builds the set Q. For each node identified in this set (line 7), its distance π t (v) to the destination node is increased by 1 in line 8. The loop in lines 16 to 24 updates the shortest paths. Each node u in set Q is investigated one-by-one in line 16 and each outgoing arc a is scanned in line 17. We check in line 18 if arc a belongs to the
9 HYBRID GA FOR OSPF/IS-IS ROUTING 9 procedure UpdateShortestPaths(a = ( u, v), w, t) 1 if a / A t return; 2 A t A t \ {a }; 3 HeapInsertMax(H, u, π t (u)); 4 δ u + δ u + 1; 5 if δ u + > 0 then return 6 Q = {u}; 7 forall v Q do 8 π t (v) π t (v) + 1; 9 forall a = (u, v) IN(v) A t do 10 A t A t \ {a}; 11 HeapInsertMax(H, u, π t (u)); 12 δ + (u) δ + (u) 1; 13 if δ + (u) = 0 then Q Q {u}; 14 end forall 15 end forall 16 forall u Q do 17 forall a = (u, v) OUT(u) do 18 if π t (u) = w a + π t (v) then 19 A t A t {a}; 20 HeapInsertMax(H, u, π t (u)); 21 δ + (u) δ + (u) + 1; 22 end if 23 end forall 24 end forall end UpdateShortestPaths. Figure 4. Pseudo-code of procedure UpdateShortestPaths. new shortest path to the destination. If so, arc a is inserted in the shortest paths graph in line 19, its tail u is inserted in the heap in line 20, and its outdegree δ + (u) is updated in line Dynamic load update. We first recall that procedure UpdateShortestPaths built a heap H containing all nodes for which the set of outgoing arcs was modified in the shortest paths graph. The pseudo-code in Figure 5 summarizes the main steps of the load update procedure. We denote by l t a the load on arc a A associated with the destination node t T. Procedure UpdateLoads takes as parameters the demands d and a destination node t. The loop in lines 1 to 9 removes nodes from the heap until the heap becomes empty. The node u with maximum distance to the destination node is removed in line 2. The total load flowing through node u is equal d ut + a=(v,u) A t l a. The load in each arc leaving node u is computed in line 3. The loop in lines 4 to 7 scans all arcs leaving node u in the current shortest paths graph for which the partial load l t a has to be updated. The new partial load l t a is set in line 5 and the head v of arc a is inserted in the heap H in line 6.
10 10 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP procedure UpdateLoads(H, d, t) 1 while HeapSize(H) > 0 do 2 u HeapExtractMax(H); 3 load (d ut + a=(v,u) A t lt a)/δ + (u); 4 forall a = (u, v) A t : l t a load do 5 l t a load; 6 HeapInsertMax(H, v, π(v)); 7 end forall 8 end while end UpdateLoads. Figure 5. Pseudo-code of procedure UpdateLoads. 6. Computational Results In this section, we describe the experimental results using the hybrid genetic algorithm introduced in this paper. We describe the computer environment, list the values of the algorithm parameters, present the test problems, and outline the experimental setup. In the experiments, we compare the hybrid genetic algorithm with the lower bound associated with the linear program (2 11) and other heuristics The setup. The experiments were done on an SGI Challenge computer ( MHz MIPS R10000 processors) with 7.6 Gb of memory. Each run used a single processor. The algorithms were implemented in C and compiled with the MIPSpro cc compiler, version 7.30, using flag -O3. Running times were measured with the getrusage function. Random numbers were generated in the hybrid genetic algorithm as well as in the pure genetic algorithm using Matsumoto and Nishimura s Mersenne Twister [15]. The following parameters were set in both the pure and the hybrid genetic algorithms: Population size: 50. Weight range: [1, w max = 20]. Population partitioning placed the top 25% of the solutions (rounded up to 13) in set A, the bottom 5% of the solutions (rounded up to 3) in set C, and the remaining solutions in set B. Probability that mutation occurs: p m = Probability that an offspring inherits the weight from the elite parent during crossover: p A = 0.7. The number of generations varied according to the type of experiment. In addition, the maximum number of candidate arcs is set to q = 5 in the local improvement procedure in the hybrid genetic algorithm. The experiments were done on 13 networks from four classes proposed by Fortz and Thorup [10] and also used in [7]. The networks are summarized in Table 1. The AT&T Worldnet backbone is a proposed real-world network of 90 routers and 274 links. The 2-level hierarchical networks are generated using the GT-ITM generator [24], based on a model of Calvert et al. [5] and Zegura et al. [25]. This model
11 HYBRID GA FOR OSPF/IS-IS ROUTING 11 Table 1. Network characteristics: Class name, instance name, number of nodes ( N ), number of arcs ( A ), number of destination nodes ( T ), number of demand pairs (o-d pairs), total demand ( d uv ), and demand scaling factor (ρ). Class Name N A T o-d pairs duv ρ AT&T backbone att level hierarchical hier hier100a hier50a hier50b Random rand rand100b rand rand50a Waxman wax wax100a wax wax50a uses two types of arcs: local access arcs have capacities equal to 200, while long distance arcs have capacities equal to For the class of random networks, the probability of having an arc between two nodes is given by a parameter that controls the density of the network. All arc capacities are set to In Waxman networks, the nodes are points uniformly distributed in the unit square. The probability of having an arc between nodes u and v is ηe (u,v)/(2θ), where η is a parameter used to control the density of the network, (u, v) is the Euclidean distance between u and v, and θ is the maximum distance between any two nodes in the network [23]. All arc capacities are set to Fortz and Thorup generated the demands to force some nodes to be more active senders or receivers than others, thus modeling hot spots on the network. Their generation assigns higher demands to closely located nodes pairs. For each instance, 12 distinct demand matrices D 1, D 2,..., D 12 are generated. Starting from demand matrix D 1, the other demand matrices are generated by repeatedly multiplying D 1 by a scaling factor: D k = ρ k 1 D 1, k = 1,..., 12. Table 1 summarizes the test problems. For each problem, the table lists its class, name, number of nodes, number of arcs, number of destination nodes, number of origin-destination pairs, the total demand of D 1, and the scaling factor Fixed time comparison. In this subsection, we compare the hybrid genetic algorithm with three heuristics: InvCap: weights are set proportional to the inverse of the link capacity, i.e. w a = c max /c a, where c max is the maximum link capacity; GA: the basic genetic algorithm without the local search used by the hybrid genetic algorithm; LS: the local search algorithm of Fortz and Thorup [10]; as well as with LPLB, the linear programming lower bound Φ OP T. InvCap is used in Cisco IOS 10.3 and later by default [6, 22]. GA is derived from the genetic algorithm in [7], which will be also compared with the hybrid genetic algorithm later in Section 6.4. LS is the implementation used in [10].
12 12 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP Table 2. Routing costs for att with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total att 20 InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP-LB 15 cost demand Figure 6. InvCap, GA, HGA, LS, and LP lower bound on att. For each of the 13 networks, 12 increasingly loaded traffic demand matrices are considered. InvCap and LPLB were run a single time. Ten one-hour runs were done with GA, HGA, and LS on each instance and average routing costs computed. Tables 2 to 14 and Figures 6 to 18 summarize these results. For each demand level, the tables list the normalized costs for InvCap and LPLB as well as the average
13 HYBRID GA FOR OSPF/IS-IS ROUTING 13 Table 3. Routing costs for hier50a with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total hier50a 20 InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP-LB 15 cost demand Figure 7. InvCap, GA, HGA, LS, and LP lower bound on hier50a. normalized costs over ten one-hour runs for GA, HGA, and LS. The last row in each table lists the sum of the normalized average costs for each algorithm. Normalized cost values less than (i.e., recall that when the routing cost exceeds we say that the routing congests the network; see Section 2) are separated from those for which the network is congested by a line segment in the tables. The distribution of the costs can be seen in the figures, where all ten cost values for each algorithm and each demand point are plotted together with the average costs.
14 14 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP Table 4. Routing costs for hier50b with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP-LB hier50b 10 cost demand Figure 8. InvCap, GA, HGA, LS, and LP lower bound on hier50b.
15 HYBRID GA FOR OSPF/IS-IS ROUTING 15 Table 5. Routing costs for hier100 with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total hier cost 10 InvCap GA 5 GA (mean) LS LS (mean) HGA HGA (mean) LP-LB demand Figure 9. InvCap, GA, HGA, LS, and LP lower bound on hier100.
16 16 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP Table 6. Routing costs for hier100a with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP-LB hier100a cost demand Figure 10. InvCap, GA, HGA, LS, and LP lower bound on hier100a.
17 HYBRID GA FOR OSPF/IS-IS ROUTING 17 Table 7. Routing costs for rand50 with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total rand cost demand InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP-LB Figure 11. InvCap, GA, HGA, LS, and LP lower bound on rand50.
18 18 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP Table 8. Routing costs for rand50a with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP-LB rand50a cost demand Figure 12. InvCap, GA, HGA, LS, and LP lower bound on rand50a.
19 HYBRID GA FOR OSPF/IS-IS ROUTING 19 Table 9. Routing costs for rand100 with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP-LB rand cost demand Figure 13. InvCap, GA, HGA, LS, and LP lower bound on rand100.
20 20 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP Table 10. Routing costs for rand100b with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP-LB rand100b 120 cost demand Figure 14. InvCap, GA, HGA, LS, and LP lower bound on rand100b.
21 HYBRID GA FOR OSPF/IS-IS ROUTING 21 Table 11. Routing costs for wax50 with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total wax InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP-LB 8 cost demand Figure 15. InvCap, GA, HGA, LS, and LP lower bound on wax50.
22 22 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP Table 12. Routing costs for wax50a with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP-LB wax50a 8 cost demand Figure 16. InvCap, GA, HGA, LS, and LP lower bound on wax50a.
23 HYBRID GA FOR OSPF/IS-IS ROUTING 23 Table 13. Routing costs for wax100 with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP-LB wax100 cost demand Figure 17. InvCap, GA, HGA, LS, and LP lower bound on wax100.
24 24 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP Table 14. Routing costs for wax100a with scaled projected demands. Solutions are averaged over ten one-hour runs. Demand InvCap GA HGA LS LPLB Total InvCap GA GA (mean) LS LS (mean) HGA HGA (mean) LP wax100a 20 cost demand Figure 18. InvCap, GA, HGA, LS, and LP lower bound on wax100a. We make the following remarks about the computational results. The pure genetic algorithm (GA) consistently found better solutions than InvCap. Solution differences increased with traffic intensity. The relative differences Φ InvCap Φ GA /Φ InvCap between the solution values obtained by GA and InvCap varied from 13.6% on rand100b to 94.8% on hier50b. HGA found solutions at least as good as GA on all networks and for all demand levels. Solution differences increased with traffic intensity. The relative differences
25 HYBRID GA FOR OSPF/IS-IS ROUTING att/demand= cost 100 HGA GA LP lower bound time (secs) Figure 19. Cost as a function of time on 1-hour run: HGA versus GA on att with demand Φ GA Φ HGA /Φ GA varied from as little as 1.9% on hier100 to as large as 88.0% on rand100b. HGA not only found better-quality solutions, but did so in less CPU time (see Figures 19 and 20, which compare one run of HGA and GA each on the att network with demand D = ). Figure 19 shows the value of the best-quality solution in the population as a function of CPU time, while Figure 20 shows the value of the best-quality solution in the population as a function of the generation of the algorithm. These figures illustrate how close to the LP lower bound the HGA comes and how much faster HGA is to converge compared to GA. Of the 13 network classes, HGA found the best average solutions in seven classes, while LS found the best in the remaining six. On classes rand100 and rand100b, solution values found by LS were 37.7% and 32.2% smaller with respect to those found by HGA. On the other hand, on classes hier50a and hier100a, solution values found by HGA were 6.5% and 8.1% smaller with respect to those found by LS. On the remaining nine network instances, the relative difference between average solution values found by HGA and LS was small, varying from 0.09% to 3.5%. HGA and LS on average found solution values that varied from 0.6% to 105.5% of the linear programming lower bound. HGA found solution values with a relative gap of less than 5% for five of the 13 network classes, while LS did so for three classes. Both HGA and LS found solution values with less than 10% relative gap on six classes Distribution of time-to-target-solution-value. To study and compare the computation times, we used the methodology proposed by Aiex et al.[1] and further explored by Resende and Ribeiro e.g. in [18].
26 26 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP 1000 att/demand= cost 100 HGA GA LP lower bound generations Figure 20. Cost as a function of generations on 1-hour run: HGA versus GA on att with demand Without loss of generality, we considered networks att with the demand equal to , hier50a with the demand equal to , and rand50 with the demand equal to to illustrate the general behavior observed for most instances. For each of them, we performed one hundred independent runs with different seeds of each algorithm GA, HGA, and LS, considering a given parameter value look4. Each execution was terminated when a solution of value less than or equal to the target value look4 was found or when the time limit of one hour was reached. Three different values (corresponding to easy, medium, and difficult cases) of look4 were investigated for each network: 2.89, 2.77 and 2.64 for network att, 2.32, 2.21 and 2.11 for network hier50a, and 2.93, 2.81 and 2.68 for network rand50. Different target values were used since the networks and demands were different. Empirical probability distributions for the time-to-target-solution-value are plotted in Figures 21 to 24. Runs which failed to find a solution of value less than or equal to the target value look4 within the one-hour time limit were discarded in these plots. To plot the empirical distribution for each algorithm and each instance, we associate with the i-th smallest running time t i a probability p i = (i 1 2 )/100, and plot the points z i = (t i, p i ), for i = 1,..., n r, where n r 100 is the number of runs which found a solution of value less than or equal to the target value look4 within the one-hour time limit. HGA and LS found solutions with value less than or equal to the target in all runs associated with network att (Figures 21 and 22). GA failed to find solutions at least as good as the target value on eight runs with the easiest target, 19 runs with the medium target, and 59 runs with the hardest target. HGA is not only much faster than GA, but also the computation times of the former are more predictable than those of the latter. As we can see from Figure 21, in many runs of GA
27 HYBRID GA FOR OSPF/IS-IS ROUTING cumulative probability HGA: look4 = 2.89 HGA: look4 = 2.77 HGA: look4 = GA: look4 = 2.89 GA: look4 = 2.77 GA: look4 = time to target value (seconds) Figure 21. Time to target solution value: HGA versus GA on network att with demand cumulative probability HGA: look4 = 2.89 HGA: look4 = 2.77 HGA: look4 = LS: look4 = 2.89 LS: look4 = 2.77 LS: look4 = time to target value (seconds) Figure 22. Time to target solution value: HGA versus LS on network att with demand
28 28 L.S. BURIOL, M.G.C. RESENDE, C.C. RIBEIRO, AND M. THORUP HGA: look4 = 2.32 HGA: look4 = 2.21 HGA: look4 = 2.11 GA: look4 = 2.32 GA: look4 = cumulative probability time to target value (seconds) Figure 23. Time to target solution value: HGA versus GA on network rand50 with demand cumulative probability HGA: look4 = 2.32 HGA: look4 = 2.21 HGA: look4 = LS: look4 = 2.32 LS: look4 = 2.21 LS: look4 = time to target value (seconds) Figure 24. Time to target solution value: HGA versus LS on network hier50a with demand
29 HYBRID GA FOR OSPF/IS-IS ROUTING 29 the computation times are several orders of magnitude larger than those of HGA. Considering Figure 22, we notice that the computation times of HGA are more predictable than those of LS. The latter are several times larger than the former in many runs. In more than 30% of the runs LS encounters difficulties to converge and requires very long computation times. The figure seems to suggest that this is related to the fact that LS frequently gets stuck at a local minimum and the escape mechanism often needs to be applied repeatedly to succeed. In the case of network hier50a (Figures 23 and 24), HGA performed much better than GA, as we can see from Figure 23. HGA found solutions at least as good as the target within the one-hour time limit in all runs for the easy and medium target values, and in 94 out of the 100 runs for the more difficult target. Contrarily, GA rarely found solutions at least as good as the target: 44 times for the easy target, 14 for the medium target, and only once in the hardest case (which cannot be plotted in the figure). HGA also did better than LS for this network, as shown in Figure 24. LS found solutions at least as good as the target in all runs for the easier target value. However, it failed in two runs with the medium target. For the harder target value, it missed the target in 48 runs. The instances associated with random graphs are among the hardest ones for the genetic algorithms. Network rand50 is particularly difficult for GA, which failed to find a solution at least as good as the target within the one-hour time limit in all 300 runs. Contrarily, HGA found solutions at least as good as the target in all runs. Although HGA performs better than LS for many classes of instances (as illustrated in Figures 22 and 24), LS did better in this case. Both algorithms are robust and found solutions at least as good as the target in all runs, but LS was usually faster Comparison with the original GA. We conclude the experimental results section by comparing the GA used in the experiments described above with the GA presented in [7] to show that both implementations produce comparable results. Consequently, the conclusions we make regarding the comparison of HGA and GA should carry over to the GA of [7], which we refer to as GA 0. Although GA and GA 0 are based on the same C code, they differ in some aspects: GA uses a population of size 50, while GA 0 uses one of size 200. Note that the size 200 population was used in [7] for the instances in this experiment. For larger instances, the experiments in [7] used a population of size 100. GA randomly generates an initial population with weights in the interval [1, 7], while GA 0 randomly generates all but one element of the initial population with weights in the interval [1, 20] and completes the remaining element of the population with the weights corresponding to the heuristic InvCap. Define the population parameters α = A / A B C and β = C / A B C. GA uses population parameters (α, β) = (.25,.05) while GA 0 uses population parameters (α, β) = (.20,.10), We ran both genetic algorithm implementations for ten independent one-hour trials on three networks, att, hier50a, and rand50, with two demand matrices each. We used the matrices with the two largest demands in the experiments described earlier in this section. For network att, we use demands and
SURVIVABLE IP NETWORK DESIGN WITH OSPF ROUTING
SURVIVABLE IP NETWORK DESIGN WITH OSPF ROUTING L.S. BURIOL, M.G.C. RESENDE, AND M. THORUP Abstract. Internet protocol (IP) traffic follows rules established by routing protocols. Shortest path based protocols,
More informationA GRASP HEURISTIC FOR THE COOPERATIVE COMMUNICATION PROBLEM IN AD HOC NETWORKS
A GRASP HEURISTIC FOR THE COOPERATIVE COMMUNICATION PROBLEM IN AD HOC NETWORKS C. COMMANDER, C.A.S. OLIVEIRA, P.M. PARDALOS, AND M.G.C. RESENDE ABSTRACT. Ad hoc networks are composed of a set of wireless
More informationA GRASP heuristic for the Cooperative Communication Problem in Ad Hoc Networks
MIC2005: The Sixth Metaheuristics International Conference??-1 A GRASP heuristic for the Cooperative Communication Problem in Ad Hoc Networks Clayton Commander Carlos A.S. Oliveira Panos M. Pardalos Mauricio
More informationMetaheuristics in network design
Metaheuristics in network design Talk given at the Tippie School of Management University of Iowa, Iowa City, IA Mauricio G. C. Resende AT&T Labs Research Florham Park, New Jersey mgcr@att.com Summary
More informationGateways Placement in Backbone Wireless Mesh Networks
I. J. Communications, Network and System Sciences, 2009, 1, 1-89 Published Online February 2009 in SciRes (http://www.scirp.org/journal/ijcns/). Gateways Placement in Backbone Wireless Mesh Networks Abstract
More informationOSPF Fundamentals. Agenda. OSPF Principles. L41 - OSPF Fundamentals. Open Shortest Path First Routing Protocol Internet s Second IGP
OSPF Fundamentals Open Shortest Path First Routing Protocol Internet s Second IGP Agenda OSPF Principles Introduction The Dijkstra Algorithm Communication Procedures LSA Broadcast Handling Splitted Area
More informationOSPF - Open Shortest Path First. OSPF Fundamentals. Agenda. OSPF Topology Database
OSPF - Open Shortest Path First OSPF Fundamentals Open Shortest Path First Routing Protocol Internet s Second IGP distance vector protocols like RIP have several dramatic disadvantages: slow adaptation
More informationLink State Routing. Brad Karp UCL Computer Science. CS 3035/GZ01 3 rd December 2013
Link State Routing Brad Karp UCL Computer Science CS 33/GZ 3 rd December 3 Outline Link State Approach to Routing Finding Links: Hello Protocol Building a Map: Flooding Protocol Healing after Partitions:
More informationShort course: Some applications of combinatorial optimization in telecommunications
Advanced School and Workshop on Mathematical Techniques and Problems in Telecommunications Tomar, Portugal September 8-12, 2003 Short course: Some applications of combinatorial optimization in telecommunications
More informationLink-state protocols and Open Shortest Path First (OSPF)
Fixed Internetworking Protocols and Networks Link-state protocols and Open Shortest Path First (OSPF) Rune Hylsberg Jacobsen Aarhus School of Engineering rhj@iha.dk 0 ITIFN Objectives Describe the basic
More informationLink State Routing. Stefano Vissicchio UCL Computer Science CS 3035/GZ01
Link State Routing Stefano Vissicchio UCL Computer Science CS 335/GZ Reminder: Intra-domain Routing Problem Shortest paths problem: What path between two vertices offers minimal sum of edge weights? Classic
More informationConfiguring OSPF. Information About OSPF CHAPTER
CHAPTER 22 This chapter describes how to configure the ASASM to route data, perform authentication, and redistribute routing information using the Open Shortest Path First (OSPF) routing protocol. The
More informationCS 457 Lecture 16 Routing Continued. Spring 2010
CS 457 Lecture 16 Routing Continued Spring 2010 Scaling Link-State Routing Overhead of link-state routing Flooding link-state packets throughout the network Running Dijkstra s shortest-path algorithm Introducing
More informationDesign of Parallel Algorithms. Communication Algorithms
+ Design of Parallel Algorithms Communication Algorithms + Topic Overview n One-to-All Broadcast and All-to-One Reduction n All-to-All Broadcast and Reduction n All-Reduce and Prefix-Sum Operations n Scatter
More informationA Study of Dynamic Routing and Wavelength Assignment with Imprecise Network State Information
A Study of Dynamic Routing and Wavelength Assignment with Imprecise Network State Information Jun Zhou Department of Computer Science Florida State University Tallahassee, FL 326 zhou@cs.fsu.edu Xin Yuan
More informationConfiguring the maximum number of external LSAs in LSDB 27 Configuring OSPF exit overflow interval 28 Enabling compatibility with RFC Logging
Contents Configuring OSPF 1 Overview 1 OSPF packets 1 LSA types 1 OSPF areas 2 Router types 4 Route types 5 Route calculation 6 OSPF network types 6 DR and BDR 6 Protocols and standards 8 OSPF configuration
More informationOSPF Domain / OSPF Area. OSPF Advanced Topics. OSPF Domain / OSPF Area. Agenda
OSPF Domain / OSPF Area OSPF Advanced Topics Areas,, Backbone, Summary-LSA, ASBR, Stub Area, Route Summarization, Virtual Links, Header Details OSPF domain can be divided in multiple OSPF areas to improve
More informationIntroduction to OSPF. ISP Workshops. Last updated 11 November 2013
Introduction to OSPF ISP Workshops Last updated 11 November 2013 1 OSPF p Open Shortest Path First p Open: n Meaning an Open Standard n Developed by IETF (OSPF Working Group) for IP RFC1247 n Current standard
More informationTransportation Timetabling
Outline DM87 SCHEDULING, TIMETABLING AND ROUTING 1. Sports Timetabling Lecture 16 Transportation Timetabling Marco Chiarandini 2. Transportation Timetabling Tanker Scheduling Air Transport Train Timetabling
More informationDistance-Vector Routing
Distance-Vector Routing Antonio Carzaniga Faculty of Informatics University of Lugano June 8, 2007 c 2005 2007 Antonio Carzaniga 1 Recap on link-state routing Distance-vector routing Bellman-Ford equation
More informationChapter 12. Cross-Layer Optimization for Multi- Hop Cognitive Radio Networks
Chapter 12 Cross-Layer Optimization for Multi- Hop Cognitive Radio Networks 1 Outline CR network (CRN) properties Mathematical models at multiple layers Case study 2 Traditional Radio vs CR Traditional
More informationA Factorial Representation of Permutations and Its Application to Flow-Shop Scheduling
Systems and Computers in Japan, Vol. 38, No. 1, 2007 Translated from Denshi Joho Tsushin Gakkai Ronbunshi, Vol. J85-D-I, No. 5, May 2002, pp. 411 423 A Factorial Representation of Permutations and Its
More informationA Comparative Study of Quality of Service Routing Schemes That Tolerate Imprecise State Information
A Comparative Study of Quality of Service Routing Schemes That Tolerate Imprecise State Information Xin Yuan Wei Zheng Department of Computer Science, Florida State University, Tallahassee, FL 330 {xyuan,zheng}@cs.fsu.edu
More informationTravel time uncertainty and network models
Travel time uncertainty and network models CE 392C TRAVEL TIME UNCERTAINTY One major assumption throughout the semester is that travel times can be predicted exactly and are the same every day. C = 25.87321
More informationEmpirical Probability Based QoS Routing
Empirical Probability Based QoS Routing Xin Yuan Guang Yang Department of Computer Science, Florida State University, Tallahassee, FL 3230 {xyuan,guanyang}@cs.fsu.edu Abstract We study Quality-of-Service
More informationPart VII: VRP - advanced topics
Part VII: VRP - advanced topics c R.F. Hartl, S.N. Parragh 1/32 Overview Dealing with TW and duration constraints Solving VRP to optimality c R.F. Hartl, S.N. Parragh 2/32 Dealing with TW and duration
More informationOn Optimal Routing with Multiple Traffic Matrices
On Optimal Routing with Multiple Traffic Matrices Chun Zhang a, Yong Liu a, Weibo Gong b, Jim Kurose a, Robert Moll a and Don Towsley a a Dept. of Computer Science, b Dept. of Electrical & Computer Engineering,
More informationShuffled Complex Evolution
Shuffled Complex Evolution Shuffled Complex Evolution An Evolutionary algorithm That performs local and global search A solution evolves locally through a memetic evolution (Local search) This local search
More informationTRAFFIC SIGNAL CONTROL WITH ANT COLONY OPTIMIZATION. A Thesis presented to the Faculty of California Polytechnic State University, San Luis Obispo
TRAFFIC SIGNAL CONTROL WITH ANT COLONY OPTIMIZATION A Thesis presented to the Faculty of California Polytechnic State University, San Luis Obispo In Partial Fulfillment of the Requirements for the Degree
More informationIntroduction to Local and Wide Area Networks
Introduction to Local and Wide Area Networks Lecturers Amnach Khawne Jirasak Sittigorn Chapter 1 1 Routing Protocols and Concepts Chapter 10 : Link-State Routing Protocols Chapter 11 : OSPF Chapter 1 2
More informationCSE/EE 461. Link State Routing. Last Time. This Lecture. Routing Algorithms Introduction Distance Vector routing (RIP)
CSE/EE 46 Link State Routing Last Time Routing Algorithms Introduction Distance Vector routing (RIP) Application Presentation Session Transport Network Data Link Physical This Lecture Routing Algorithms
More informationCreative Commons: Attribution 3.0 Hong Kong License
Title A simultaneous bus route design and frequency setting problem for Tin Shui Wai, Hong Kong Author(s) Szeto, WY; Wu, Y Citation European Journal Of Operational Research, 2011, v. 209 n. 2, p. 141-155
More informationFast Placement Optimization of Power Supply Pads
Fast Placement Optimization of Power Supply Pads Yu Zhong Martin D. F. Wong Dept. of Electrical and Computer Engineering Dept. of Electrical and Computer Engineering Univ. of Illinois at Urbana-Champaign
More informationEnergy Saving Routing Strategies in IP Networks
Energy Saving Routing Strategies in IP Networks M. Polverini; M. Listanti DIET Department - University of Roma Sapienza, Via Eudossiana 8, 84 Roma, Italy 2 june 24 [scale=.8]figure/logo.eps M. Polverini
More informationFOUR TOTAL TRANSFER CAPABILITY. 4.1 Total transfer capability CHAPTER
CHAPTER FOUR TOTAL TRANSFER CAPABILITY R structuring of power system aims at involving the private power producers in the system to supply power. The restructured electric power industry is characterized
More informationComputer Networks II
ipartimento di Informatica e Sistemistica omputer Networks II Routing protocols Overview Luca Becchetti Luca.Becchetti@dis.uniroma.it.. 2009/200 Goals escribe approaches and give overview of mechanisms
More informationA new mixed integer linear programming formulation for one problem of exploration of online social networks
manuscript No. (will be inserted by the editor) A new mixed integer linear programming formulation for one problem of exploration of online social networks Aleksandra Petrović Received: date / Accepted:
More informationUtilization-Aware Adaptive Back-Pressure Traffic Signal Control
Utilization-Aware Adaptive Back-Pressure Traffic Signal Control Wanli Chang, Samarjit Chakraborty and Anuradha Annaswamy Abstract Back-pressure control of traffic signal, which computes the control phase
More informationOptimizing Client Association in 60 GHz Wireless Access Networks
Optimizing Client Association in 60 GHz Wireless Access Networks G Athanasiou, C Weeraddana, C Fischione, and L Tassiulas KTH Royal Institute of Technology, Stockholm, Sweden University of Thessaly, Volos,
More informationDistributed Collaborative Path Planning in Sensor Networks with Multiple Mobile Sensor Nodes
7th Mediterranean Conference on Control & Automation Makedonia Palace, Thessaloniki, Greece June 4-6, 009 Distributed Collaborative Path Planning in Sensor Networks with Multiple Mobile Sensor Nodes Theofanis
More informationVehicle routing problems with road-network information
50 Dominique Feillet Mines Saint-Etienne and LIMOS, CMP Georges Charpak, F-13541 Gardanne, France Vehicle routing problems with road-network information ORBEL - Liège, February 1, 2018 Vehicle Routing
More informationHow Much Can Sub-band Virtual Concatenation (VCAT) Help Static Routing and Spectrum Assignment in Elastic Optical Networks?
How Much Can Sub-band Virtual Concatenation (VCAT) Help Static Routing and Spectrum Assignment in Elastic Optical Networks? (Invited) Xin Yuan, Gangxiang Shen School of Electronic and Information Engineering
More informationNetwork Layer (Routing)
Network Layer (Routing) Where we are in the ourse Moving on up to the Network Layer! Application Transport Network Link Physical SE 61 University of Washington Topics Network service models Datagrams (packets),
More informationThe Capability of Error Correction for Burst-noise Channels Using Error Estimating Code
The Capability of Error Correction for Burst-noise Channels Using Error Estimating Code Yaoyu Wang Nanjing University yaoyu.wang.nju@gmail.com June 10, 2016 Yaoyu Wang (NJU) Error correction with EEC June
More informationOptical Networks with Limited Wavelength Conversion.
Practical Routing and Wavelength Assignment algorithms for All Optical Networks with Limited Wavelength Conversion M.D. Swaminathan*, Indian Institute of Science, Bangalore, India. Abstract We present
More informationTable of Contents. OSPF Configuration 1
Table of Contents OSPF Configuration 1 Introduction to OSPF 1 Basic Concepts 2 Area Based OSPF Network Partition 3 Router Types 7 Classification of OSPF Networks 9 DR and BDR 9 OSPF Packet Formats 11 Supported
More informationComputational aspects of two-player zero-sum games Course notes for Computational Game Theory Section 3 Fall 2010
Computational aspects of two-player zero-sum games Course notes for Computational Game Theory Section 3 Fall 21 Peter Bro Miltersen November 1, 21 Version 1.3 3 Extensive form games (Game Trees, Kuhn Trees)
More informationEconomic Design of Control Chart Using Differential Evolution
Economic Design of Control Chart Using Differential Evolution Rukmini V. Kasarapu 1, Vijaya Babu Vommi 2 1 Assistant Professor, Department of Mechanical Engineering, Anil Neerukonda Institute of Technology
More informationLocalization (Position Estimation) Problem in WSN
Localization (Position Estimation) Problem in WSN [1] Convex Position Estimation in Wireless Sensor Networks by L. Doherty, K.S.J. Pister, and L.E. Ghaoui [2] Semidefinite Programming for Ad Hoc Wireless
More informationA Location-Aware Routing Metric (ALARM) for Multi-Hop, Multi-Channel Wireless Mesh Networks
A Location-Aware Routing Metric (ALARM) for Multi-Hop, Multi-Channel Wireless Mesh Networks Eiman Alotaibi, Sumit Roy Dept. of Electrical Engineering U. Washington Box 352500 Seattle, WA 98195 eman76,roy@ee.washington.edu
More informationTABLE OF CONTENTS CHAPTER NO. TITLE PAGE NO. LIST OF TABLES LIST OF FIGURES LIST OF SYMBOLS AND ABBREVIATIONS
vi TABLE OF CONTENTS CHAPTER TITLE PAGE ABSTRACT LIST OF TABLES LIST OF FIGURES LIST OF SYMBOLS AND ABBREVIATIONS iii viii x xiv 1 INTRODUCTION 1 1.1 DISK SCHEDULING 1 1.2 WINDOW-CONSTRAINED SCHEDULING
More informationLink State Routing. In particular OSPF. dr. C. P. J. Koymans. Informatics Institute University of Amsterdam. March 4, 2008
Link State Routing In particular OSPF dr. C. P. J. Koymans Informatics Institute University of Amsterdam March 4, 2008 dr. C. P. J. Koymans (UvA) Link State Routing March 4, 2008 1 / 70 1 Link State Protocols
More information1 This work was partially supported by NSF Grant No. CCR , and by the URI International Engineering Program.
Combined Error Correcting and Compressing Codes Extended Summary Thomas Wenisch Peter F. Swaszek Augustus K. Uht 1 University of Rhode Island, Kingston RI Submitted to International Symposium on Information
More informationLecture 8 Link-State Routing
6998-02: Internet Routing Lecture 8 Link-State Routing John Ioannidis AT&T Labs Research ji+ir@cs.columbia.edu Copyright 2002 by John Ioannidis. All Rights Reserved. Announcements Lectures 1-5, 7-8 are
More informationSkip Lists S 3 S 2 S 1. 2/6/2016 7:04 AM Skip Lists 1
Skip Lists S 3 15 15 23 10 15 23 36 2/6/2016 7:04 AM Skip Lists 1 Outline and Reading What is a skip list Operations Search Insertion Deletion Implementation Analysis Space usage Search and update times
More informationEnergy-Efficient MANET Routing: Ideal vs. Realistic Performance
Energy-Efficient MANET Routing: Ideal vs. Realistic Performance Paper by: Thomas Knuz IEEE IWCMC Conference Aug. 2008 Presented by: Farzana Yasmeen For : CSE 6590 2013.11.12 Contents Introduction Review:
More informationLecture 20: Combinatorial Search (1997) Steven Skiena. skiena
Lecture 20: Combinatorial Search (1997) Steven Skiena Department of Computer Science State University of New York Stony Brook, NY 11794 4400 http://www.cs.sunysb.edu/ skiena Give an O(n lg k)-time algorithm
More informationGeneralized Game Trees
Generalized Game Trees Richard E. Korf Computer Science Department University of California, Los Angeles Los Angeles, Ca. 90024 Abstract We consider two generalizations of the standard two-player game
More informationAn Optimization Approach for Real Time Evacuation Reroute. Planning
An Optimization Approach for Real Time Evacuation Reroute Planning Gino J. Lim and M. Reza Baharnemati and Seon Jin Kim November 16, 2015 Abstract This paper addresses evacuation route management in the
More informationJunos Intermediate Routing
Junos Intermediate Routing Chapter 4: Open Shortest Path First 2012 Juniper Networks, Inc. All rights reserved. www.juniper.net Worldwide Education Services Chapter Objectives After successfully completing
More informationFrequency Hopping Pattern Recognition Algorithms for Wireless Sensor Networks
Frequency Hopping Pattern Recognition Algorithms for Wireless Sensor Networks Min Song, Trent Allison Department of Electrical and Computer Engineering Old Dominion University Norfolk, VA 23529, USA Abstract
More informationLink and Link Impedance 2018/02/13. VECTOR DATA ANALYSIS Network Analysis TYPES OF OPERATIONS
VECTOR DATA ANALYSIS Network Analysis A network is a system of linear features that has the appropriate attributes for the flow of objects. A network is typically topology-based: lines (arcs) meet at intersections
More informationGraphs and Network Flows IE411. Lecture 14. Dr. Ted Ralphs
Graphs and Network Flows IE411 Lecture 14 Dr. Ted Ralphs IE411 Lecture 14 1 Review: Labeling Algorithm Pros Guaranteed to solve any max flow problem with integral arc capacities Provides constructive tool
More informationSearch then involves moving from state-to-state in the problem space to find a goal (or to terminate without finding a goal).
Search Can often solve a problem using search. Two requirements to use search: Goal Formulation. Need goals to limit search and allow termination. Problem formulation. Compact representation of problem
More informationA Hybrid Evolutionary Approach for Multi Robot Path Exploration Problem
A Hybrid Evolutionary Approach for Multi Robot Path Exploration Problem K.. enthilkumar and K. K. Bharadwaj Abstract - Robot Path Exploration problem or Robot Motion planning problem is one of the famous
More informationMulti-Dimensional Conflict Graph Based Computing for Optimal Capacity in MR-MC Wireless Networks
Multi-Dimensional Conflict Graph Based Computing for Optimal Capacity in MR-MC Wireless Networks Hongkun Li, Yu Cheng, Chi Zhou Department of Electrical and Computer Engineering Illinois Institute of Technology
More informationC Commands. Send comments to
This chapter describes the Cisco NX-OS Open Shortest Path First (OSPF) commands that begin with C. UCR-583 clear ip ospf neighbor clear ip ospf neighbor To clear neighbor statistics and reset adjacencies
More informationScheduling. Radek Mařík. April 28, 2015 FEE CTU, K Radek Mařík Scheduling April 28, / 48
Scheduling Radek Mařík FEE CTU, K13132 April 28, 2015 Radek Mařík (marikr@fel.cvut.cz) Scheduling April 28, 2015 1 / 48 Outline 1 Introduction to Scheduling Methodology Overview 2 Classification of Scheduling
More informationISudoku. Jonathon Makepeace Matthew Harris Jamie Sparrow Julian Hillebrand
Jonathon Makepeace Matthew Harris Jamie Sparrow Julian Hillebrand ISudoku Abstract In this paper, we will analyze and discuss the Sudoku puzzle and implement different algorithms to solve the puzzle. After
More informationChapter 5 OPTIMIZATION OF BOW TIE ANTENNA USING GENETIC ALGORITHM
Chapter 5 OPTIMIZATION OF BOW TIE ANTENNA USING GENETIC ALGORITHM 5.1 Introduction This chapter focuses on the use of an optimization technique known as genetic algorithm to optimize the dimensions of
More informationControl of the Contract of a Public Transport Service
Control of the Contract of a Public Transport Service Andrea Lodi, Enrico Malaguti, Nicolás E. Stier-Moses Tommaso Bonino DEIS, University of Bologna Graduate School of Business, Columbia University SRM
More informationHeuristic Search with Pre-Computed Databases
Heuristic Search with Pre-Computed Databases Tsan-sheng Hsu tshsu@iis.sinica.edu.tw http://www.iis.sinica.edu.tw/~tshsu 1 Abstract Use pre-computed partial results to improve the efficiency of heuristic
More informationProblem 1 (15 points: Graded by Shahin) Recall the network structure of our in-class trading experiment shown in Figure 1
Solutions for Homework 2 Networked Life, Fall 204 Prof Michael Kearns Due as hardcopy at the start of class, Tuesday December 9 Problem (5 points: Graded by Shahin) Recall the network structure of our
More informationA GRAPH THEORETICAL APPROACH TO SOLVING SCRAMBLE SQUARES PUZZLES. 1. Introduction
GRPH THEORETICL PPROCH TO SOLVING SCRMLE SQURES PUZZLES SRH MSON ND MLI ZHNG bstract. Scramble Squares puzzle is made up of nine square pieces such that each edge of each piece contains half of an image.
More informationAd Hoc and Neighborhood Search Methods for Placement of Mesh Routers in Wireless Mesh Networks
29 29th IEEE International Conference on Distributed Computing Systems Workshops Ad Hoc and Neighborhood Search Methods for Placement of Mesh Routers in Wireless Mesh Networks Fatos Xhafa Department of
More informationCCO Commun. Comb. Optim.
Communications in Combinatorics and Optimization Vol. 2 No. 2, 2017 pp.149-159 DOI: 10.22049/CCO.2017.25918.1055 CCO Commun. Comb. Optim. Graceful labelings of the generalized Petersen graphs Zehui Shao
More informationData Gathering. Chapter 4. Ad Hoc and Sensor Networks Roger Wattenhofer 4/1
Data Gathering Chapter 4 Ad Hoc and Sensor Networks Roger Wattenhofer 4/1 Environmental Monitoring (PermaSense) Understand global warming in alpine environment Harsh environmental conditions Swiss made
More informationWireless Internet Routing. IEEE s
Wireless Internet Routing IEEE 802.11s 1 Acknowledgments Cigdem Sengul, Deutsche Telekom Laboratories 2 Outline Introduction Interworking Topology discovery Routing 3 IEEE 802.11a/b/g /n /s IEEE 802.11s:
More informationFree Cell Solver. Copyright 2001 Kevin Atkinson Shari Holstege December 11, 2001
Free Cell Solver Copyright 2001 Kevin Atkinson Shari Holstege December 11, 2001 Abstract We created an agent that plays the Free Cell version of Solitaire by searching through the space of possible sequences
More informationAccess point selection algorithms for maximizing throughputs in wireless LAN environment
Access point selection algorithms for maximizing throughputs in wireless LAN environment Akihiro Fujiwara Yasuhiro Sagara Masahiko Nakamura Department of Computer Science and Electronics Kyushu Institute
More informationTraffic Control for a Swarm of Robots: Avoiding Group Conflicts
Traffic Control for a Swarm of Robots: Avoiding Group Conflicts Leandro Soriano Marcolino and Luiz Chaimowicz Abstract A very common problem in the navigation of robotic swarms is when groups of robots
More informationQuestion No: 2 In an OSPF Hello packet, which of the following fields must match for all neighbor routers on the segment? Choose three answers.
Volume: 335 Questions Question No: 1 What is the default preference value for a static route in the Alcatel-Lucent 7750 SR? A. 0 B. 5 C. 10 D. 15 Answer: B Question No: 2 In an OSPF Hello packet, which
More informationRating and Generating Sudoku Puzzles Based On Constraint Satisfaction Problems
Rating and Generating Sudoku Puzzles Based On Constraint Satisfaction Problems Bahare Fatemi, Seyed Mehran Kazemi, Nazanin Mehrasa International Science Index, Computer and Information Engineering waset.org/publication/9999524
More informationLink State Routing. In particular OSPF. Karst Koymans. Informatics Institute University of Amsterdam. (version 16.3, 2017/03/09 11:25:31)
Link State Routing In particular OSPF Karst Koymans Informatics Institute University of Amsterdam (version 16.3, 2017/03/09 11:25:31) Tuesday, March 7, 2017 Karst Koymans (UvA) Link State Routing Tuesday,
More informationisudoku Computing Solutions to Sudoku Puzzles w/ 3 Algorithms by: Gavin Hillebrand Jamie Sparrow Jonathon Makepeace Matthew Harris
isudoku Computing Solutions to Sudoku Puzzles w/ 3 Algorithms by: Gavin Hillebrand Jamie Sparrow Jonathon Makepeace Matthew Harris What is Sudoku? A logic-based puzzle game Heavily based in combinatorics
More informationGuide to OSPF Application on the CSS 11000
Guide to OSPF Application on the CSS 11000 Document ID: 12638 Contents Introduction Before You Begin Conventions Prerequisites Components Used Description OSPF Configuration Task List Configuration Global
More informationIntroduction. Introduction ROBUST SENSOR POSITIONING IN WIRELESS AD HOC SENSOR NETWORKS. Smart Wireless Sensor Systems 1
ROBUST SENSOR POSITIONING IN WIRELESS AD HOC SENSOR NETWORKS Xiang Ji and Hongyuan Zha Material taken from Sensor Network Operations by Shashi Phoa, Thomas La Porta and Christopher Griffin, John Wiley,
More informationJoint QoS Multicast Routing and Channel Assignment in Multiradio Multichannel Wireless Mesh Networks using Intelligent Computational Methods
Joint QoS Multicast Routing and Channel Assignment in Multiradio Multichannel Wireless Mesh Networks using Intelligent Computational Methods Hui Cheng,a, Shengxiang Yang b a Department of Computer Science,
More information4.4 Shortest Paths in a Graph Revisited
4.4 Shortest Paths in a Graph Revisited shortest path from computer science department to Einstein's house Algorithm Design by Éva Tardos and Jon Kleinberg Slides by Kevin Wayne Copyright 2004 Addison
More informationLecture 13 Register Allocation: Coalescing
Lecture 13 Register llocation: Coalescing I. Motivation II. Coalescing Overview III. lgorithms: Simple & Safe lgorithm riggs lgorithm George s lgorithm Phillip. Gibbons 15-745: Register Coalescing 1 Review:
More informationStanford University CS261: Optimization Handout 9 Luca Trevisan February 1, 2011
Stanford University CS261: Optimization Handout 9 Luca Trevisan February 1, 2011 Lecture 9 In which we introduce the maximum flow problem. 1 Flows in Networks Today we start talking about the Maximum Flow
More informationA Practical Approach to Bitrate Control in Wireless Mesh Networks using Wireless Network Utility Maximization
A Practical Approach to Bitrate Control in Wireless Mesh Networks using Wireless Network Utility Maximization EE359 Course Project Mayank Jain Department of Electrical Engineering Stanford University Introduction
More informationNovel Placement Mesh Router Approach for Wireless Mesh Network
Novel Placement Mesh Router Approach for Wireless Mesh Network Mohsen Rezaei 1, Mehdi Agha Sarram 2,Vali Derhami 3,and Hossein Mahboob Sarvestani 4 Electrical and Computer Engineering Department, Yazd
More informationM U LT I C A S T C O M M U N I C AT I O N S. Tarik Cicic
M U LT I C A S T C O M M U N I C AT I O N S Tarik Cicic 9..08 O V E R V I E W One-to-many communication, why and how Algorithmic approach: Steiner trees Practical algorithms Multicast tree types Basic
More informationA Performance Comparison of Multi-Hop Wireless Ad Hoc Network Routing Protocols
A Performance Comparison of Multi-Hop Wireless Ad Hoc Network Routing Protocols Josh Broch, David Maltz, David Johnson, Yih-Chun Hu and Jorjeta Jetcheva Computer Science Department Carnegie Mellon University
More informationSimulative Comparison of MPLS Protection Switching vs. OSPF Re-routing
C O R P O R A T E T E C H N O L O Y Simulative Comparison of MPLS Protection Switching vs. OSPF Re-routing nformation & Sandrine PASQUALINI Antoine FROT Andreas Iselt Andreas Kirstädter C O R P O R A T
More informationCHANNEL ASSIGNMENT AND LOAD DISTRIBUTION IN A POWER- MANAGED WLAN
CHANNEL ASSIGNMENT AND LOAD DISTRIBUTION IN A POWER- MANAGED WLAN Mohamad Haidar Robert Akl Hussain Al-Rizzo Yupo Chan University of Arkansas at University of Arkansas at University of Arkansas at University
More informationInterference-Aware Joint Routing and TDMA Link Scheduling for Static Wireless Networks
Interference-Aware Joint Routing and TDMA Link Scheduling for Static Wireless Networks Yu Wang Weizhao Wang Xiang-Yang Li Wen-Zhan Song Abstract We study efficient interference-aware joint routing and
More informationOptimisation and Operations Research
Optimisation and Operations Research Lecture : Graph Problems and Dijkstra s algorithm Matthew Roughan http://www.maths.adelaide.edu.au/matthew.roughan/ Lecture_notes/OORII/
More informationCitation for published version (APA): Nutma, T. A. (2010). Kac-Moody Symmetries and Gauged Supergravity Groningen: s.n.
University of Groningen Kac-Moody Symmetries and Gauged Supergravity Nutma, Teake IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please
More information