Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Let's trust the mathematics here: sqrt + fabs + fabs IS positive or nul :)
[simgrid.git] / src / surf / surf_routing_vivaldi.c
1 /* Copyright (c) 2009, 2010, 2011. The SimGrid Team.
2  * All rights reserved.                                                     */
3
4 /* This program is free software; you can redistribute it and/or modify it
5  * under the terms of the license (GNU LGPL) which comes with this package. */
6 #include "surf_routing_private.h"
7
8 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(surf_route_vivaldi, surf, "Routing part of surf");
9
10 /* Business methods */
11 static route_t vivaldi_get_route(AS_t rc, const char *src, const char *dst)
12 {
13           xbt_assert(rc && src
14                       && dst,
15                       "Invalid params for \"get_route\" function at AS \"%s\"",
16                       rc->name);
17
18           route_t new_e_route = xbt_new0(s_route_t, 1);
19           new_e_route->src_gateway = ROUTER_PEER(src);
20           new_e_route->dst_gateway = ROUTER_PEER(dst);
21           new_e_route->link_list =  xbt_dynar_new(0, NULL);
22           return new_e_route;
23 }
24
25 static double euclidean_dist_comp(int index, xbt_dynar_t src, xbt_dynar_t dst)
26 {
27         double src_coord, dst_coord;
28
29         // FIXME converting from string to float each time we need a coordinate is ... suboptimal
30         src_coord = atof(xbt_dynar_get_as(src, index, char *));
31         dst_coord = atof(xbt_dynar_get_as(dst, index, char *));
32
33         return (src_coord-dst_coord)*(src_coord-dst_coord);
34
35 }
36
37 static double base_vivaldi_get_latency (const char *src, const char *dst)
38 {
39   double euclidean_dist;
40   xbt_dynar_t src_ctn, dst_ctn;
41   src_ctn = xbt_lib_get_or_null(host_lib, src, COORD_HOST_LEVEL);
42   if(!src_ctn) src_ctn = xbt_lib_get_or_null(as_router_lib, src, COORD_ASR_LEVEL);
43   dst_ctn = xbt_lib_get_or_null(host_lib, dst, COORD_HOST_LEVEL);
44   if(!dst_ctn) dst_ctn = xbt_lib_get_or_null(as_router_lib, dst, COORD_ASR_LEVEL);
45
46   if(dst_ctn == NULL || src_ctn == NULL)
47   xbt_die("Coord src '%s' :%p   dst '%s' :%p",src,src_ctn,dst,dst_ctn);
48
49   euclidean_dist = sqrt (euclidean_dist_comp(0,src_ctn,dst_ctn)+euclidean_dist_comp(1,src_ctn,dst_ctn))
50                                  + fabs(atof(xbt_dynar_get_as(src_ctn, 2, char *)))+fabs(atof(xbt_dynar_get_as(dst_ctn, 2, char *)));
51
52   //From .ms to .s
53   return euclidean_dist / 1000;
54 }
55
56 static double vivaldi_get_latency (AS_t rc,const char *src, const char *dst, route_t e_route)
57 {
58   if(routing_get_network_element_type(src) == SURF_NETWORK_ELEMENT_AS) {
59           int need_to_clean = e_route?0:1;
60           double latency;
61           e_route = e_route ? e_route : rc->get_route(rc, src, dst);
62           latency = base_vivaldi_get_latency(e_route->src_gateway,e_route->dst_gateway);
63           if(need_to_clean) generic_free_route(e_route);
64           return latency;
65   } else {
66           return base_vivaldi_get_latency(src,dst);
67   }
68 }
69
70 /* Creation routing model functions */
71 AS_t model_vivaldi_create(void)
72 {
73           AS_t new_component = model_rulebased_create();
74           new_component->get_route = vivaldi_get_route;
75           new_component->get_latency = vivaldi_get_latency;
76           return new_component;
77 }