1 /* Copyright (c) 2014-2018. 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"
16 #include <boost/algorithm/string/classification.hpp>
17 #include <boost/algorithm/string/split.hpp>
19 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(surf_route_fat_tree, surf, "Routing for fat trees");
25 FatTreeZone::FatTreeZone(NetZone* father, std::string name) : ClusterZone(father, name)
27 XBT_DEBUG("Creating a new fat tree.");
30 FatTreeZone::~FatTreeZone()
32 for (unsigned int i = 0; i < this->nodes_.size(); i++) {
33 delete this->nodes_[i];
35 for (unsigned int i = 0; i < this->links_.size(); i++) {
36 delete this->links_[i];
40 bool FatTreeZone::isInSubTree(FatTreeNode* root, FatTreeNode* node)
42 XBT_DEBUG("Is %d(%u,%u) in the sub tree of %d(%u,%u) ?", node->id, node->level, node->position, root->id, root->level,
44 if (root->level <= node->level) {
47 for (unsigned int i = 0; i < node->level; i++) {
48 if (root->label[i] != node->label[i]) {
53 for (unsigned int i = root->level; i < this->levels_; i++) {
54 if (root->label[i] != node->label[i]) {
61 void FatTreeZone::get_local_route(NetPoint* src, NetPoint* dst, RouteCreationArgs* into, double* latency)
64 if (dst->is_router() || src->is_router())
67 /* Let's find the source and the destination in our internal structure */
68 auto searchedNode = this->compute_nodes_.find(src->id());
69 xbt_assert(searchedNode != this->compute_nodes_.end(), "Could not find the source %s [%u] in the fat tree",
70 src->get_cname(), src->id());
71 FatTreeNode* source = searchedNode->second;
73 searchedNode = this->compute_nodes_.find(dst->id());
74 xbt_assert(searchedNode != this->compute_nodes_.end(), "Could not find the destination %s [%u] in the fat tree",
75 dst->get_cname(), dst->id());
76 FatTreeNode* destination = searchedNode->second;
78 XBT_VERB("Get route and latency from '%s' [%u] to '%s' [%u] in a fat tree", src->get_cname(), src->id(),
79 dst->get_cname(), dst->id());
81 /* In case destination is the source, and there is a loopback, let's use it instead of going up to a switch */
82 if (source->id == destination->id && this->has_loopback_) {
83 into->link_list.push_back(source->loopback);
85 *latency += source->loopback->latency();
89 FatTreeNode* currentNode = source;
92 while (not isInSubTree(currentNode, destination)) {
93 int d = destination->position; // as in d-mod-k
95 for (unsigned int i = 0; i < currentNode->level; i++)
96 d /= this->num_parents_per_node_[i];
98 int k = this->num_parents_per_node_[currentNode->level];
100 into->link_list.push_back(currentNode->parents[d]->up_link_);
103 *latency += currentNode->parents[d]->up_link_->latency();
105 if (this->has_limiter_)
106 into->link_list.push_back(currentNode->limiter_link_);
107 currentNode = currentNode->parents[d]->up_node_;
110 XBT_DEBUG("%d(%u,%u) is in the sub tree of %d(%u,%u).", destination->id, destination->level, destination->position,
111 currentNode->id, currentNode->level, currentNode->position);
114 while (currentNode != destination) {
115 for (unsigned int i = 0; i < currentNode->children.size(); i++) {
116 if (i % this->num_children_per_node_[currentNode->level - 1] == destination->label[currentNode->level - 1]) {
117 into->link_list.push_back(currentNode->children[i]->down_link_);
119 *latency += currentNode->children[i]->down_link_->latency();
120 currentNode = currentNode->children[i]->down_node_;
121 if (this->has_limiter_)
122 into->link_list.push_back(currentNode->limiter_link_);
123 XBT_DEBUG("%d(%u,%u) is accessible through %d(%u,%u)", destination->id, destination->level,
124 destination->position, currentNode->id, currentNode->level, currentNode->position);
130 /* This function makes the assumption that parse_specific_arguments() and
131 * addNodes() have already been called
133 void FatTreeZone::seal()
135 if (this->levels_ == 0) {
138 this->generateSwitches();
140 if (XBT_LOG_ISENABLED(surf_route_fat_tree, xbt_log_priority_debug)) {
141 std::stringstream msgBuffer;
143 msgBuffer << "We are creating a fat tree of " << this->levels_ << " levels "
144 << "with " << this->nodes_by_level_[0] << " processing nodes";
145 for (unsigned int i = 1; i <= this->levels_; i++) {
146 msgBuffer << ", " << this->nodes_by_level_[i] << " switches at level " << i;
148 XBT_DEBUG("%s", msgBuffer.str().c_str());
150 msgBuffer << "Nodes are : ";
152 for (unsigned int i = 0; i < this->nodes_.size(); i++) {
153 msgBuffer << this->nodes_[i]->id << "(" << this->nodes_[i]->level << "," << this->nodes_[i]->position << ") ";
155 XBT_DEBUG("%s", msgBuffer.str().c_str());
158 this->generateLabels();
161 // Nodes are totally ordered, by level and then by position, in this->nodes
162 for (unsigned int i = 0; i < this->levels_; i++) {
163 for (unsigned int j = 0; j < this->nodes_by_level_[i]; j++) {
164 this->connectNodeToParents(this->nodes_[k]);
169 if (XBT_LOG_ISENABLED(surf_route_fat_tree, xbt_log_priority_debug)) {
170 std::stringstream msgBuffer;
171 msgBuffer << "Links are : ";
172 for (unsigned int i = 0; i < this->links_.size(); i++) {
173 msgBuffer << "(" << this->links_[i]->up_node_->id << "," << this->links_[i]->down_node_->id << ") ";
175 XBT_DEBUG("%s", msgBuffer.str().c_str());
179 int FatTreeZone::connectNodeToParents(FatTreeNode* node)
181 std::vector<FatTreeNode*>::iterator currentParentNode = this->nodes_.begin();
182 int connectionsNumber = 0;
183 const int level = node->level;
184 XBT_DEBUG("We are connecting node %d(%u,%u) to his parents.", node->id, node->level, node->position);
185 currentParentNode += this->getLevelPosition(level + 1);
186 for (unsigned int i = 0; i < this->nodes_by_level_[level + 1]; i++) {
187 if (this->areRelated(*currentParentNode, node)) {
188 XBT_DEBUG("%d(%u,%u) and %d(%u,%u) are related,"
189 " with %u links between them.",
190 node->id, node->level, node->position, (*currentParentNode)->id, (*currentParentNode)->level,
191 (*currentParentNode)->position, this->num_port_lower_level_[level]);
192 for (unsigned int j = 0; j < this->num_port_lower_level_[level]; j++) {
193 this->add_link(*currentParentNode, node->label[level] + j * this->num_children_per_node_[level], node,
194 (*currentParentNode)->label[level] + j * this->num_parents_per_node_[level]);
200 return connectionsNumber;
203 bool FatTreeZone::areRelated(FatTreeNode* parent, FatTreeNode* child)
205 std::stringstream msgBuffer;
207 if (XBT_LOG_ISENABLED(surf_route_fat_tree, xbt_log_priority_debug)) {
208 msgBuffer << "Are " << child->id << "(" << child->level << "," << child->position << ") <";
210 for (unsigned int i = 0; i < this->levels_; i++) {
211 msgBuffer << child->label[i] << ",";
215 msgBuffer << " and " << parent->id << "(" << parent->level << "," << parent->position << ") <";
216 for (unsigned int i = 0; i < this->levels_; i++) {
217 msgBuffer << parent->label[i] << ",";
220 msgBuffer << " related ? ";
221 XBT_DEBUG("%s", msgBuffer.str().c_str());
223 if (parent->level != child->level + 1) {
227 for (unsigned int i = 0; i < this->levels_; i++) {
228 if (parent->label[i] != child->label[i] && i + 1 != parent->level) {
235 void FatTreeZone::generateSwitches()
237 XBT_DEBUG("Generating switches.");
238 this->nodes_by_level_.resize(this->levels_ + 1, 0);
240 // Take care of the number of nodes by level
241 this->nodes_by_level_[0] = 1;
242 for (unsigned int i = 0; i < this->levels_; i++)
243 this->nodes_by_level_[0] *= this->num_children_per_node_[i];
245 if (this->nodes_by_level_[0] != this->nodes_.size()) {
246 surf_parse_error(std::string("The number of provided nodes does not fit with the wanted topology.") +
247 " Please check your platform description (We need " + std::to_string(this->nodes_by_level_[0]) +
248 "nodes, we got " + std::to_string(this->nodes_.size()));
252 for (unsigned int i = 0; i < this->levels_; i++) {
253 int nodesInThisLevel = 1;
255 for (unsigned int j = 0; j <= i; j++)
256 nodesInThisLevel *= this->num_parents_per_node_[j];
258 for (unsigned int j = i + 1; j < this->levels_; j++)
259 nodesInThisLevel *= this->num_children_per_node_[j];
261 this->nodes_by_level_[i + 1] = nodesInThisLevel;
264 // Create the switches
266 for (unsigned int i = 0; i < this->levels_; i++) {
267 for (unsigned int j = 0; j < this->nodes_by_level_[i + 1]; j++) {
268 FatTreeNode* newNode = new FatTreeNode(this->cluster_, --k, i + 1, j);
269 XBT_DEBUG("We create the switch %d(%u,%u)", newNode->id, newNode->level, newNode->position);
270 newNode->children.resize(this->num_children_per_node_[i] * this->num_port_lower_level_[i]);
271 if (i != this->levels_ - 1) {
272 newNode->parents.resize(this->num_parents_per_node_[i + 1] * this->num_port_lower_level_[i + 1]);
274 newNode->label.resize(this->levels_);
275 this->nodes_.push_back(newNode);
280 void FatTreeZone::generateLabels()
282 XBT_DEBUG("Generating labels.");
283 // TODO : check if nodesByLevel and nodes are filled
284 std::vector<int> maxLabel(this->levels_);
285 std::vector<int> currentLabel(this->levels_);
287 for (unsigned int i = 0; i <= this->levels_; i++) {
288 currentLabel.assign(this->levels_, 0);
289 for (unsigned int j = 0; j < this->levels_; j++) {
290 maxLabel[j] = j + 1 > i ? this->num_children_per_node_[j] : this->num_parents_per_node_[j];
293 for (unsigned int j = 0; j < this->nodes_by_level_[i]; j++) {
295 if (XBT_LOG_ISENABLED(surf_route_fat_tree, xbt_log_priority_debug)) {
296 std::stringstream msgBuffer;
298 msgBuffer << "Assigning label <";
299 for (unsigned int l = 0; l < this->levels_; l++) {
300 msgBuffer << currentLabel[l] << ",";
302 msgBuffer << "> to " << k << " (" << i << "," << j << ")";
304 XBT_DEBUG("%s", msgBuffer.str().c_str());
306 this->nodes_[k]->label.assign(currentLabel.begin(), currentLabel.end());
308 bool remainder = true;
309 unsigned int pos = 0;
310 while (remainder && pos < this->levels_) {
312 if (currentLabel[pos] >= maxLabel[pos]) {
313 currentLabel[pos] = 0;
326 int FatTreeZone::getLevelPosition(const unsigned int level)
328 xbt_assert(level <= this->levels_, "The impossible did happen. Yet again.");
329 int tempPosition = 0;
331 for (unsigned int i = 0; i < level; i++)
332 tempPosition += this->nodes_by_level_[i];
337 void FatTreeZone::add_processing_node(int id)
339 using std::make_pair;
340 static int position = 0;
341 FatTreeNode* newNode;
342 newNode = new FatTreeNode(this->cluster_, id, 0, position++);
343 newNode->parents.resize(this->num_parents_per_node_[0] * this->num_port_lower_level_[0]);
344 newNode->label.resize(this->levels_);
345 this->compute_nodes_.insert(make_pair(id, newNode));
346 this->nodes_.push_back(newNode);
349 void FatTreeZone::add_link(FatTreeNode* parent, unsigned int parentPort, FatTreeNode* child, unsigned int childPort)
351 FatTreeLink* newLink;
352 newLink = new FatTreeLink(this->cluster_, child, parent);
353 XBT_DEBUG("Creating a link between the parent (%u,%u,%u) and the child (%u,%u,%u)", parent->level, parent->position,
354 parentPort, child->level, child->position, childPort);
355 parent->children[parentPort] = newLink;
356 child->parents[childPort] = newLink;
358 this->links_.push_back(newLink);
361 void FatTreeZone::parse_specific_arguments(ClusterCreationArgs* cluster)
363 std::vector<std::string> parameters;
364 std::vector<std::string> tmp;
365 boost::split(parameters, cluster->topo_parameters, boost::is_any_of(";"));
366 const std::string error_msg {"Fat trees are defined by the levels number and 3 vectors, see the documentation for more information"};
368 // TODO : we have to check for zeros and negative numbers, or it might crash
369 if (parameters.size() != 4) {
370 surf_parse_error(error_msg);
373 // The first parts of topo_parameters should be the levels number
375 this->levels_ = std::stoi(parameters[0]);
376 } catch (std::invalid_argument& ia) {
377 throw std::invalid_argument(std::string("First parameter is not the amount of levels:") + parameters[0]);
380 // Then, a l-sized vector standing for the children number by level
381 boost::split(tmp, parameters[1], boost::is_any_of(","));
382 if (tmp.size() != this->levels_) {
383 surf_parse_error(error_msg);
385 for (size_t i = 0; i < tmp.size(); i++) {
387 this->num_children_per_node_.push_back(std::stoi(tmp[i]));
388 } catch (std::invalid_argument& ia) {
389 throw std::invalid_argument(std::string("Invalid lower level node number:") + tmp[i]);
393 // Then, a l-sized vector standing for the parents number by level
394 boost::split(tmp, parameters[2], boost::is_any_of(","));
395 if (tmp.size() != this->levels_) {
396 surf_parse_error(error_msg);
398 for (size_t i = 0; i < tmp.size(); i++) {
400 this->num_parents_per_node_.push_back(std::stoi(tmp[i]));
401 } catch (std::invalid_argument& ia) {
402 throw std::invalid_argument(std::string("Invalid upper level node number:") + tmp[i]);
406 // Finally, a l-sized vector standing for the ports number with the lower level
407 boost::split(tmp, parameters[3], boost::is_any_of(","));
408 if (tmp.size() != this->levels_) {
409 surf_parse_error(error_msg);
411 for (size_t i = 0; i < tmp.size(); i++) {
413 this->num_port_lower_level_.push_back(std::stoi(tmp[i]));
414 } catch (std::invalid_argument& ia) {
415 throw std::invalid_argument(std::string("Invalid lower level port number:") + tmp[i]);
418 this->cluster_ = cluster;
421 void FatTreeZone::generate_dot_file(const std::string& filename) const
424 file.open(filename, std::ios::out | std::ios::trunc);
425 xbt_assert(file.is_open(), "Unable to open file %s", filename.c_str());
427 file << "graph AsClusterFatTree {\n";
428 for (unsigned int i = 0; i < this->nodes_.size(); i++) {
429 file << this->nodes_[i]->id;
430 if (this->nodes_[i]->id < 0)
431 file << " [shape=circle];\n";
433 file << " [shape=hexagon];\n";
436 for (unsigned int i = 0; i < this->links_.size(); i++) {
437 file << this->links_[i]->down_node_->id << " -- " << this->links_[i]->up_node_->id << ";\n";
443 FatTreeNode::FatTreeNode(ClusterCreationArgs* cluster, int id, int level, int position)
444 : id(id), level(level), position(position)
446 LinkCreationArgs linkTemplate;
447 if (cluster->limiter_link) {
448 linkTemplate.bandwidth = cluster->limiter_link;
449 linkTemplate.latency = 0;
450 linkTemplate.policy = s4u::Link::SharingPolicy::SHARED;
451 linkTemplate.id = "limiter_"+std::to_string(id);
452 sg_platf_new_link(&linkTemplate);
453 this->limiter_link_ = resource::LinkImpl::byName(linkTemplate.id);
455 if (cluster->loopback_bw || cluster->loopback_lat) {
456 linkTemplate.bandwidth = cluster->loopback_bw;
457 linkTemplate.latency = cluster->loopback_lat;
458 linkTemplate.policy = s4u::Link::SharingPolicy::FATPIPE;
459 linkTemplate.id = "loopback_"+ std::to_string(id);
460 sg_platf_new_link(&linkTemplate);
461 this->loopback = resource::LinkImpl::byName(linkTemplate.id);
465 FatTreeLink::FatTreeLink(ClusterCreationArgs* cluster, FatTreeNode* downNode, FatTreeNode* upNode)
466 : up_node_(upNode), down_node_(downNode)
468 static int uniqueId = 0;
469 LinkCreationArgs linkTemplate;
470 linkTemplate.bandwidth = cluster->bw;
471 linkTemplate.latency = cluster->lat;
472 linkTemplate.policy = cluster->sharing_policy; // sthg to do with that ?
474 "link_from_" + std::to_string(downNode->id) + "_" + std::to_string(upNode->id) + "_" + std::to_string(uniqueId);
475 sg_platf_new_link(&linkTemplate);
477 if (cluster->sharing_policy == s4u::Link::SharingPolicy::SPLITDUPLEX) {
478 std::string tmpID = std::string(linkTemplate.id) + "_UP";
479 this->up_link_ = resource::LinkImpl::byName(tmpID); // check link?
480 tmpID = std::string(linkTemplate.id) + "_DOWN";
481 this->down_link_ = resource::LinkImpl::byName(tmpID); // check link ?
483 this->up_link_ = resource::LinkImpl::byName(linkTemplate.id);
484 this->down_link_ = this->up_link_;