1 #include "network_smpi.hpp"
2 #include "simgrid/sg_config.h"
4 XBT_LOG_EXTERNAL_DEFAULT_CATEGORY(surf_network);
6 xbt_dynar_t smpi_bw_factor = NULL;
7 xbt_dynar_t smpi_lat_factor = NULL;
9 typedef struct s_smpi_factor *smpi_factor_t;
10 typedef struct s_smpi_factor {
15 xbt_dict_t gap_lookup = NULL;
17 static int factor_cmp(const void *pa, const void *pb)
19 return (((s_smpi_factor_t*)pa)->factor > ((s_smpi_factor_t*)pb)->factor) ? 1 :
20 (((s_smpi_factor_t*)pa)->factor < ((s_smpi_factor_t*)pb)->factor) ? -1 : 0;
24 static xbt_dynar_t parse_factor(const char *smpi_coef_string)
27 unsigned int iter = 0;
29 xbt_dynar_t smpi_factor, radical_elements, radical_elements2 = NULL;
31 smpi_factor = xbt_dynar_new(sizeof(s_smpi_factor_t), NULL);
32 radical_elements = xbt_str_split(smpi_coef_string, ";");
33 xbt_dynar_foreach(radical_elements, iter, value) {
35 radical_elements2 = xbt_str_split(value, ":");
36 if (xbt_dynar_length(radical_elements2) != 2)
37 xbt_die("Malformed radical for smpi factor!");
38 fact.factor = atol(xbt_dynar_get_as(radical_elements2, 0, char *));
39 fact.value = atof(xbt_dynar_get_as(radical_elements2, 1, char *));
40 xbt_dynar_push_as(smpi_factor, s_smpi_factor_t, fact);
41 XBT_DEBUG("smpi_factor:\t%ld : %f", fact.factor, fact.value);
42 xbt_dynar_free(&radical_elements2);
44 xbt_dynar_free(&radical_elements);
46 xbt_dynar_sort(smpi_factor, &factor_cmp);
47 xbt_dynar_foreach(smpi_factor, iter, fact) {
48 XBT_DEBUG("ordered smpi_factor:\t%ld : %f", fact.factor, fact.value);
58 /************************************************************************/
59 /* New model based on LV08 and experimental results of MPI ping-pongs */
60 /************************************************************************/
61 /* @Inproceedings{smpi_ipdps, */
62 /* author={Pierre-Nicolas Clauss and Mark Stillwell and Stéphane Genaud and Frédéric Suter and Henri Casanova and Martin Quinson}, */
63 /* title={Single Node On-Line Simulation of {MPI} Applications with SMPI}, */
64 /* booktitle={25th IEEE International Parallel and Distributed Processing Symposium (IPDPS'11)}, */
65 /* address={Anchorage (Alaska) USA}, */
69 void surf_network_model_init_SMPI(void)
72 if (surf_network_model)
74 surf_network_model = new NetworkSmpiModel();
75 net_define_callbacks();
76 xbt_dynar_push(model_list, &surf_network_model);
77 //network_solve = lmm_solve;
79 xbt_cfg_setdefault_double(_sg_cfg_set, "network/sender_gap", 10e-6);
80 xbt_cfg_setdefault_double(_sg_cfg_set, "network/weight_S", 8775);
83 NetworkSmpiModel::NetworkSmpiModel()
84 : NetworkCm02Model() {
88 NetworkSmpiModel::~NetworkSmpiModel(){
90 xbt_dict_free(&gap_lookup);
93 xbt_dynar_free(&smpi_bw_factor);
94 smpi_bw_factor = NULL;
96 if (smpi_lat_factor) {
97 xbt_dynar_free(&smpi_lat_factor);
98 smpi_lat_factor = NULL;
102 void NetworkSmpiModel::gapAppend(double size, const NetworkLinkPtr link, NetworkActionPtr action)
104 const char *src = link->getName();
106 //surf_action_network_CM02_t last_action;
109 if (sg_sender_gap > 0.0) {
111 gap_lookup = xbt_dict_new_homogeneous(NULL);
113 fifo = (xbt_fifo_t) xbt_dict_get_or_null(gap_lookup, src);
114 action->m_senderGap = 0.0;
115 if (fifo && xbt_fifo_size(fifo) > 0) {
116 /* Compute gap from last send */
118 (surf_action_network_CM02_t)
119 xbt_fifo_get_item_content(xbt_fifo_get_last_item(fifo));*/
120 // bw = net_get_link_bandwidth(link);
121 action->m_senderGap = sg_sender_gap;
122 /* max(sg_sender_gap,last_action->sender.size / bw);*/
123 action->m_latency += action->m_senderGap;
125 /* Append action as last send */
126 /*action->sender.link_name = link->lmm_resource.generic_resource.name;
128 (xbt_fifo_t) xbt_dict_get_or_null(gap_lookup,
129 action->sender.link_name);
131 fifo = xbt_fifo_new();
132 xbt_dict_set(gap_lookup, action->sender.link_name, fifo, NULL);
134 action->sender.fifo_item = xbt_fifo_push(fifo, action);*/
135 action->m_senderSize = size;
139 void NetworkSmpiModel::gapRemove(ActionPtr lmm_action)
143 NetworkCm02ActionPtr action = static_cast<NetworkCm02ActionPtr>(lmm_action);
145 if (sg_sender_gap > 0.0 && action->p_senderLinkName
146 && action->p_senderFifoItem) {
148 (xbt_fifo_t) xbt_dict_get_or_null(gap_lookup,
149 action->p_senderLinkName);
150 xbt_fifo_remove_item(fifo, action->p_senderFifoItem);
151 size = xbt_fifo_size(fifo);
154 xbt_dict_remove(gap_lookup, action->p_senderLinkName);
155 size = xbt_dict_length(gap_lookup);
157 xbt_dict_free(&gap_lookup);
163 double NetworkSmpiModel::bandwidthFactor(double size)
167 parse_factor(sg_cfg_get_string("smpi/bw_factor"));
169 unsigned int iter = 0;
170 s_smpi_factor_t fact;
172 xbt_dynar_foreach(smpi_bw_factor, iter, fact) {
173 if (size <= fact.factor) {
174 XBT_DEBUG("%f <= %ld return %f", size, fact.factor, current);
179 XBT_DEBUG("%f > %ld return %f", size, fact.factor, current);
184 double NetworkSmpiModel::latencyFactor(double size)
186 if (!smpi_lat_factor)
188 parse_factor(sg_cfg_get_string("smpi/lat_factor"));
190 unsigned int iter = 0;
191 s_smpi_factor_t fact;
193 xbt_dynar_foreach(smpi_lat_factor, iter, fact) {
194 if (size <= fact.factor) {
195 XBT_DEBUG("%f <= %ld return %f", size, fact.factor, current);
200 XBT_DEBUG("%f > %ld return %f", size, fact.factor, current);
205 double NetworkSmpiModel::bandwidthConstraint(double rate, double bound, double size)
207 return rate < 0 ? bound : min(bound, rate * bandwidthFactor(size));