1 /* Copyright (c) 2013-2017. The SimGrid Team. All rights reserved. */
3 /* This program is free software; you can redistribute it and/or modify it
4 * under the terms of the license (GNU LGPL) which comes with this package. */
6 #include <boost/algorithm/string.hpp>
8 #include "simgrid/s4u/Engine.hpp"
9 #include "simgrid/s4u/Host.hpp"
11 #include "src/kernel/routing/NetPoint.hpp"
12 #include "src/kernel/routing/VivaldiZone.hpp"
13 #include "src/surf/network_interface.hpp"
15 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(surf_route_vivaldi, surf, "Routing part of surf");
21 simgrid::xbt::Extension<NetPoint, Coords> Coords::EXTENSION_ID;
23 Coords::Coords(NetPoint* netpoint, std::string coordStr)
25 if (not Coords::EXTENSION_ID.valid())
26 Coords::EXTENSION_ID = NetPoint::extension_create<Coords>();
28 std::vector<std::string> string_values;
29 boost::split(string_values, coordStr, boost::is_any_of(" "));
30 xbt_assert(string_values.size() == 3, "Coordinates of %s must have 3 dimensions", netpoint->cname());
32 for (auto str : string_values)
34 coords.push_back(std::stod(str));
35 } catch (std::invalid_argument const& ia) {
36 throw std::invalid_argument(std::string("Invalid coordinate: ") + ia.what());
38 coords.shrink_to_fit();
40 netpoint->extension_set<Coords>(this);
41 XBT_DEBUG("Coords of %s %p: %s", netpoint->cname(), netpoint, coordStr.c_str());
43 }; // namespace vivaldi
45 static inline double euclidean_dist_comp(int index, std::vector<double>* src, std::vector<double>* dst)
47 double src_coord = src->at(index);
48 double dst_coord = dst->at(index);
50 return (src_coord - dst_coord) * (src_coord - dst_coord);
53 static std::vector<double>* getCoordsFromNetpoint(NetPoint* np)
55 simgrid::kernel::routing::vivaldi::Coords* coords = np->extension<simgrid::kernel::routing::vivaldi::Coords>();
56 xbt_assert(coords, "Please specify the Vivaldi coordinates of %s %s (%p)",
57 (np->isNetZone() ? "Netzone" : (np->isHost() ? "Host" : "Router")), np->cname(), np);
58 return &coords->coords;
60 VivaldiZone::VivaldiZone(NetZone* father, const char* name) : ClusterZone(father, name)
64 void VivaldiZone::setPeerLink(NetPoint* netpoint, double bw_in, double bw_out, std::string coord)
66 xbt_assert(netpoint->netzone() == this, "Cannot add a peer link to a netpoint that is not in this netzone");
68 new simgrid::kernel::routing::vivaldi::Coords(netpoint, coord);
70 std::string link_up = "link_" + netpoint->name() + "_UP";
71 std::string link_down = "link_" + netpoint->name() + "_DOWN";
72 surf::LinkImpl* linkUp = surf_network_model->createLink(link_up.c_str(), bw_out, 0, SURF_LINK_SHARED);
73 surf::LinkImpl* linkDown = surf_network_model->createLink(link_down.c_str(), bw_in, 0, SURF_LINK_SHARED);
74 privateLinks_.insert({netpoint->id(), {linkUp, linkDown}});
77 void VivaldiZone::getLocalRoute(NetPoint* src, NetPoint* dst, sg_platf_route_cbarg_t route, double* lat)
79 XBT_DEBUG("vivaldi getLocalRoute from '%s'[%u] '%s'[%u]", src->cname(), src->id(), dst->cname(), dst->id());
81 if (src->isNetZone()) {
82 std::string srcName = "router_" + src->name();
83 std::string dstName = "router_" + dst->name();
84 route->gw_src = simgrid::s4u::Engine::getInstance()->getNetpointByNameOrNull(srcName.c_str());
85 route->gw_dst = simgrid::s4u::Engine::getInstance()->getNetpointByNameOrNull(dstName.c_str());
88 /* Retrieve the private links */
90 std::pair<surf::LinkImpl*, surf::LinkImpl*> info = privateLinks_.at(src->id());
92 route->link_list->push_back(info.first);
94 *lat += info.first->latency();
96 } catch (std::out_of_range& unfound) {
97 XBT_DEBUG("Source of private link (%u) doesn't exist", src->id());
101 std::pair<surf::LinkImpl*, surf::LinkImpl*> info = privateLinks_.at(dst->id());
103 route->link_list->push_back(info.second);
105 *lat += info.second->latency();
107 } catch (std::out_of_range& unfound) {
108 XBT_DEBUG("Destination of private link (%u) doesn't exist", dst->id());
111 /* Compute the extra latency due to the euclidean distance if needed */
113 std::vector<double>* srcCoords = getCoordsFromNetpoint(src);
114 std::vector<double>* dstCoords = getCoordsFromNetpoint(dst);
116 double euclidean_dist =
117 sqrt(euclidean_dist_comp(0, srcCoords, dstCoords) + euclidean_dist_comp(1, srcCoords, dstCoords)) +
118 fabs(srcCoords->at(2)) + fabs(dstCoords->at(2));
120 XBT_DEBUG("Updating latency %f += %f", *lat, euclidean_dist);
121 *lat += euclidean_dist / 1000.0; // From .ms to .s