1 /* Copyright (c) 2014-2021. 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. */
10 #include "simgrid/kernel/routing/FatTreeZone.hpp"
11 #include "simgrid/kernel/routing/NetPoint.hpp"
12 #include "src/surf/network_interface.hpp"
13 #include "src/surf/xml/platf_private.hpp"
15 #include <boost/algorithm/string/classification.hpp>
16 #include <boost/algorithm/string/split.hpp>
18 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(surf_route_fat_tree, surf, "Routing for fat trees");
24 FatTreeZone::~FatTreeZone()
26 for (FatTreeNode const* node : this->nodes_)
28 for (FatTreeLink const* link : this->links_)
32 bool FatTreeZone::is_in_sub_tree(FatTreeNode* root, FatTreeNode* node) const
34 XBT_DEBUG("Is %d(%u,%u) in the sub tree of %d(%u,%u) ?", node->id, node->level, node->position, root->id, root->level,
36 if (root->level <= node->level) {
39 for (unsigned int i = 0; i < node->level; i++) {
40 if (root->label[i] != node->label[i]) {
45 for (unsigned int i = root->level; i < this->levels_; i++) {
46 if (root->label[i] != node->label[i]) {
53 void FatTreeZone::get_local_route(NetPoint* src, NetPoint* dst, RouteCreationArgs* into, double* latency)
55 if (dst->is_router() || src->is_router())
58 /* Let's find the source and the destination in our internal structure */
59 auto searchedNode = this->compute_nodes_.find(src->id());
60 xbt_assert(searchedNode != this->compute_nodes_.end(), "Could not find the source %s [%u] in the fat tree",
61 src->get_cname(), src->id());
62 FatTreeNode* source = searchedNode->second;
64 searchedNode = this->compute_nodes_.find(dst->id());
65 xbt_assert(searchedNode != this->compute_nodes_.end(), "Could not find the destination %s [%u] in the fat tree",
66 dst->get_cname(), dst->id());
67 FatTreeNode* destination = searchedNode->second;
69 XBT_VERB("Get route and latency from '%s' [%u] to '%s' [%u] in a fat tree", src->get_cname(), src->id(),
70 dst->get_cname(), dst->id());
72 /* In case destination is the source, and there is a loopback, let's use it instead of going up to a switch */
73 if (source->id == destination->id && has_loopback()) {
74 into->link_list.push_back(source->loopback_);
76 *latency += source->loopback_->get_latency();
80 FatTreeNode* currentNode = source;
83 while (not is_in_sub_tree(currentNode, destination)) {
84 int d = destination->position; // as in d-mod-k
86 for (unsigned int i = 0; i < currentNode->level; i++)
87 d /= this->num_parents_per_node_[i];
89 int k = this->num_parents_per_node_[currentNode->level];
91 into->link_list.push_back(currentNode->parents[d]->up_link_);
94 *latency += currentNode->parents[d]->up_link_->get_latency();
97 into->link_list.push_back(currentNode->limiter_link_);
98 currentNode = currentNode->parents[d]->up_node_;
101 XBT_DEBUG("%d(%u,%u) is in the sub tree of %d(%u,%u).", destination->id, destination->level, destination->position,
102 currentNode->id, currentNode->level, currentNode->position);
105 while (currentNode != destination) {
106 for (unsigned int i = 0; i < currentNode->children.size(); i++) {
107 if (i % this->num_children_per_node_[currentNode->level - 1] == destination->label[currentNode->level - 1]) {
108 into->link_list.push_back(currentNode->children[i]->down_link_);
110 *latency += currentNode->children[i]->down_link_->get_latency();
111 currentNode = currentNode->children[i]->down_node_;
113 into->link_list.push_back(currentNode->limiter_link_);
114 XBT_DEBUG("%d(%u,%u) is accessible through %d(%u,%u)", destination->id, destination->level,
115 destination->position, currentNode->id, currentNode->level, currentNode->position);
121 /* This function makes the assumption that parse_specific_arguments() and
122 * addNodes() have already been called
124 void FatTreeZone::do_seal()
126 if (this->levels_ == 0) {
129 this->generate_switches();
131 if (XBT_LOG_ISENABLED(surf_route_fat_tree, xbt_log_priority_debug)) {
132 std::stringstream msgBuffer;
134 msgBuffer << "We are creating a fat tree of " << this->levels_ << " levels "
135 << "with " << this->nodes_by_level_[0] << " processing nodes";
136 for (unsigned int i = 1; i <= this->levels_; i++) {
137 msgBuffer << ", " << this->nodes_by_level_[i] << " switches at level " << i;
139 XBT_DEBUG("%s", msgBuffer.str().c_str());
141 msgBuffer << "Nodes are : ";
143 for (FatTreeNode const* node : this->nodes_) {
144 msgBuffer << node->id << "(" << node->level << "," << node->position << ") ";
146 XBT_DEBUG("%s", msgBuffer.str().c_str());
149 this->generate_labels();
152 // Nodes are totally ordered, by level and then by position, in this->nodes
153 for (unsigned int i = 0; i < this->levels_; i++) {
154 for (unsigned int j = 0; j < this->nodes_by_level_[i]; j++) {
155 this->connect_node_to_parents(this->nodes_[k]);
160 if (XBT_LOG_ISENABLED(surf_route_fat_tree, xbt_log_priority_debug)) {
161 std::stringstream msgBuffer;
162 msgBuffer << "Links are : ";
163 for (FatTreeLink const* link : this->links_) {
164 msgBuffer << "(" << link->up_node_->id << "," << link->down_node_->id << ") ";
166 XBT_DEBUG("%s", msgBuffer.str().c_str());
170 int FatTreeZone::connect_node_to_parents(FatTreeNode* node)
172 auto currentParentNode = this->nodes_.begin();
173 int connectionsNumber = 0;
174 const int level = node->level;
175 XBT_DEBUG("We are connecting node %d(%u,%u) to his parents.", node->id, node->level, node->position);
176 currentParentNode += this->get_level_position(level + 1);
177 for (unsigned int i = 0; i < this->nodes_by_level_[level + 1]; i++) {
178 if (this->are_related(*currentParentNode, node)) {
179 XBT_DEBUG("%d(%u,%u) and %d(%u,%u) are related,"
180 " with %u links between them.",
181 node->id, node->level, node->position, (*currentParentNode)->id, (*currentParentNode)->level,
182 (*currentParentNode)->position, this->num_port_lower_level_[level]);
183 for (unsigned int j = 0; j < this->num_port_lower_level_[level]; j++) {
184 this->add_link(*currentParentNode, node->label[level] + j * this->num_children_per_node_[level], node,
185 (*currentParentNode)->label[level] + j * this->num_parents_per_node_[level]);
191 return connectionsNumber;
194 bool FatTreeZone::are_related(FatTreeNode* parent, FatTreeNode* child) const
196 std::stringstream msgBuffer;
198 if (XBT_LOG_ISENABLED(surf_route_fat_tree, xbt_log_priority_debug)) {
199 msgBuffer << "Are " << child->id << "(" << child->level << "," << child->position << ") <";
201 for (unsigned int i = 0; i < this->levels_; i++) {
202 msgBuffer << child->label[i] << ",";
206 msgBuffer << " and " << parent->id << "(" << parent->level << "," << parent->position << ") <";
207 for (unsigned int i = 0; i < this->levels_; i++) {
208 msgBuffer << parent->label[i] << ",";
211 msgBuffer << " related ? ";
212 XBT_DEBUG("%s", msgBuffer.str().c_str());
214 if (parent->level != child->level + 1) {
218 for (unsigned int i = 0; i < this->levels_; i++) {
219 if (parent->label[i] != child->label[i] && i + 1 != parent->level) {
226 void FatTreeZone::generate_switches()
228 XBT_DEBUG("Generating switches.");
229 this->nodes_by_level_.resize(this->levels_ + 1, 0);
231 // Take care of the number of nodes by level
232 this->nodes_by_level_[0] = 1;
233 for (unsigned int i = 0; i < this->levels_; i++)
234 this->nodes_by_level_[0] *= this->num_children_per_node_[i];
236 if (this->nodes_by_level_[0] != this->nodes_.size()) {
237 surf_parse_error(std::string("The number of provided nodes does not fit with the wanted topology.") +
238 " Please check your platform description (We need " + std::to_string(this->nodes_by_level_[0]) +
239 "nodes, we got " + std::to_string(this->nodes_.size()));
242 for (unsigned int i = 0; i < this->levels_; i++) {
243 int nodesInThisLevel = 1;
245 for (unsigned int j = 0; j <= i; j++)
246 nodesInThisLevel *= this->num_parents_per_node_[j];
248 for (unsigned int j = i + 1; j < this->levels_; j++)
249 nodesInThisLevel *= this->num_children_per_node_[j];
251 this->nodes_by_level_[i + 1] = nodesInThisLevel;
254 // Create the switches
256 for (unsigned int i = 0; i < this->levels_; i++) {
257 for (unsigned int j = 0; j < this->nodes_by_level_[i + 1]; j++) {
258 auto* newNode = new FatTreeNode(--k, i + 1, j, nullptr, nullptr);
259 XBT_DEBUG("We create the switch %d(%u,%u)", newNode->id, newNode->level, newNode->position);
260 newNode->children.resize(this->num_children_per_node_[i] * this->num_port_lower_level_[i]);
261 if (i != this->levels_ - 1) {
262 newNode->parents.resize(this->num_parents_per_node_[i + 1] * this->num_port_lower_level_[i + 1]);
264 newNode->label.resize(this->levels_);
265 this->nodes_.push_back(newNode);
270 void FatTreeZone::generate_labels()
272 XBT_DEBUG("Generating labels.");
273 // TODO : check if nodesByLevel and nodes are filled
274 std::vector<int> maxLabel(this->levels_);
275 std::vector<int> currentLabel(this->levels_);
277 for (unsigned int i = 0; i <= this->levels_; i++) {
278 currentLabel.assign(this->levels_, 0);
279 for (unsigned int j = 0; j < this->levels_; j++) {
280 maxLabel[j] = j + 1 > i ? this->num_children_per_node_[j] : this->num_parents_per_node_[j];
283 for (unsigned int j = 0; j < this->nodes_by_level_[i]; j++) {
284 if (XBT_LOG_ISENABLED(surf_route_fat_tree, xbt_log_priority_debug)) {
285 std::stringstream msgBuffer;
287 msgBuffer << "Assigning label <";
288 for (unsigned int l = 0; l < this->levels_; l++) {
289 msgBuffer << currentLabel[l] << ",";
291 msgBuffer << "> to " << k << " (" << i << "," << j << ")";
293 XBT_DEBUG("%s", msgBuffer.str().c_str());
295 this->nodes_[k]->label.assign(currentLabel.begin(), currentLabel.end());
297 bool remainder = true;
298 unsigned int pos = 0;
299 while (remainder && pos < this->levels_) {
301 if (currentLabel[pos] >= maxLabel[pos]) {
302 currentLabel[pos] = 0;
315 int FatTreeZone::get_level_position(const unsigned int level)
317 xbt_assert(level <= this->levels_, "The impossible did happen. Yet again.");
318 int tempPosition = 0;
320 for (unsigned int i = 0; i < level; i++)
321 tempPosition += this->nodes_by_level_[i];
326 void FatTreeZone::add_processing_node(int id, resource::LinkImpl* limiter, resource::LinkImpl* loopback)
328 using std::make_pair;
329 static int position = 0;
330 auto* newNode = new FatTreeNode(id, 0, position++, limiter, loopback);
331 newNode->parents.resize(this->num_parents_per_node_[0] * this->num_port_lower_level_[0]);
332 newNode->label.resize(this->levels_);
333 this->compute_nodes_.insert(make_pair(id, newNode));
334 this->nodes_.push_back(newNode);
337 void FatTreeZone::add_link(FatTreeNode* parent, unsigned int parentPort, FatTreeNode* child, unsigned int childPort)
339 static int uniqueId = 0;
340 const s4u::Link* linkup;
341 const s4u::Link* linkdown;
343 "link_from_" + std::to_string(child->id) + "_" + std::to_string(parent->id) + "_" + std::to_string(uniqueId);
345 if (cluster_->sharing_policy == s4u::Link::SharingPolicy::SPLITDUPLEX) {
346 linkup = create_link(id + "_UP", std::vector<double>{cluster_->bw})->set_latency(cluster_->lat)->seal();
347 linkdown = create_link(id + "_DOWN", std::vector<double>{cluster_->bw})->set_latency(cluster_->lat)->seal();
349 linkup = create_link(id, std::vector<double>{cluster_->bw})->set_latency(cluster_->lat)->seal();
354 auto* newLink = new FatTreeLink(child, parent, linkup->get_impl(), linkdown->get_impl());
355 XBT_DEBUG("Creating a link between the parent (%u,%u,%u) and the child (%u,%u,%u)", parent->level, parent->position,
356 parentPort, child->level, child->position, childPort);
357 parent->children[parentPort] = newLink;
358 child->parents[childPort] = newLink;
360 this->links_.push_back(newLink);
363 void FatTreeZone::parse_specific_arguments(ClusterCreationArgs* cluster)
365 std::vector<std::string> parameters;
366 std::vector<std::string> tmp;
367 boost::split(parameters, cluster->topo_parameters, boost::is_any_of(";"));
369 // TODO : we have to check for zeros and negative numbers, or it might crash
371 parameters.size() == 4,
372 "Fat trees are defined by the levels number and 3 vectors, see the documentation for more information.");
374 // The first parts of topo_parameters should be the levels number
376 this->levels_ = std::stoi(parameters[0]);
377 } catch (const std::invalid_argument&) {
378 surf_parse_error(std::string("First parameter is not the amount of levels: ") + parameters[0]);
381 // Then, a l-sized vector standing for the children number by level
382 boost::split(tmp, parameters[1], boost::is_any_of(","));
383 surf_parse_assert(tmp.size() == this->levels_, std::string("You specified ") + std::to_string(this->levels_) +
384 " levels but the child count vector (the first one) contains " +
385 std::to_string(tmp.size()) + " levels.");
387 for (std::string const& level : tmp) {
389 this->num_children_per_node_.push_back(std::stoi(level));
390 } catch (const std::invalid_argument&) {
391 surf_parse_error(std::string("Invalid child count: ") + level);
395 // Then, a l-sized vector standing for the parents number by level
396 boost::split(tmp, parameters[2], boost::is_any_of(","));
397 surf_parse_assert(tmp.size() == this->levels_, std::string("You specified ") + std::to_string(this->levels_) +
398 " levels but the parent count vector (the second one) contains " +
399 std::to_string(tmp.size()) + " levels.");
400 for (std::string const& parent : tmp) {
402 this->num_parents_per_node_.push_back(std::stoi(parent));
403 } catch (const std::invalid_argument&) {
404 surf_parse_error(std::string("Invalid parent count: ") + parent);
408 // Finally, a l-sized vector standing for the ports number with the lower level
409 boost::split(tmp, parameters[3], boost::is_any_of(","));
410 surf_parse_assert(tmp.size() == this->levels_, std::string("You specified ") + std::to_string(this->levels_) +
411 " levels but the port count vector (the third one) contains " +
412 std::to_string(tmp.size()) + " levels.");
413 for (std::string const& port : tmp) {
415 this->num_port_lower_level_.push_back(std::stoi(port));
416 } catch (const std::invalid_argument&) {
417 throw std::invalid_argument(std::string("Invalid lower level port number:") + port);
420 this->cluster_ = cluster;
423 void FatTreeZone::generate_dot_file(const std::string& filename) const
426 file.open(filename, std::ios::out | std::ios::trunc);
427 xbt_assert(file.is_open(), "Unable to open file %s", filename.c_str());
429 file << "graph AsClusterFatTree {\n";
430 for (FatTreeNode const* node : this->nodes_) {
433 file << " [shape=circle];\n";
435 file << " [shape=hexagon];\n";
438 for (FatTreeLink const* link : this->links_) {
439 file << link->down_node_->id << " -- " << link->up_node_->id << ";\n";
444 } // namespace routing
445 } // namespace kernel
448 NetZone* create_fatTree_zone(const std::string& name)
450 return (new kernel::routing::FatTreeZone(name))->get_iface();
454 } // namespace simgrid