Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Merge commit '045db1657e870c721be490b411868f4181a12ced' into surf++
[simgrid.git] / src / surf / cpu_cas01.cpp
1 /* Copyright (c) 2009-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
7 #include "cpu_cas01.hpp"
8 #include "cpu_ti.hpp"
9 #include "surf.hpp"
10 #include "maxmin_private.h"
11 #include "simgrid/sg_config.h"
12
13 extern "C" {
14 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(surf_cpu_cas, surf_cpu,
15                                 "Logging specific to the SURF CPU IMPROVED module");
16 }
17
18 static xbt_swag_t
19     cpu_running_action_set_that_does_not_need_being_checked = NULL;
20
21 /*************
22  * CallBacks *
23  *************/
24
25 static void parse_cpu_init(sg_platf_host_cbarg_t host){
26   ((CpuCas01ModelPtr)surf_cpu_model)->parseInit(host);
27 }
28
29 static void cpu_add_traces_cpu(){
30   surf_cpu_model->addTraces();
31 }
32
33 static void cpu_define_callbacks()
34 {
35   sg_platf_host_add_cb(parse_cpu_init);
36   sg_platf_postparse_add_cb(cpu_add_traces_cpu);
37 }
38
39 /*********
40  * Model *
41  *********/
42 void surf_cpu_model_init_Cas01()
43 {
44   char *optim = xbt_cfg_get_string(_sg_cfg_set, "cpu/optim");
45
46   if (surf_cpu_model)
47     return;
48
49   if (!strcmp(optim, "TI")) {
50     surf_cpu_model_init_ti();
51     return;
52   }
53
54   surf_cpu_model = new CpuCas01Model();
55   cpu_define_callbacks();
56   ModelPtr model = static_cast<ModelPtr>(surf_cpu_model);
57   xbt_dynar_push(model_list, &model);
58 }
59
60 CpuCas01Model::CpuCas01Model() : CpuModel("cpu")
61 {
62   ActionPtr action;
63   ActionLmmPtr actionlmm;
64
65   char *optim = xbt_cfg_get_string(_sg_cfg_set, "cpu/optim");
66   int select = xbt_cfg_get_boolean(_sg_cfg_set, "cpu/maxmin_selective_update");
67
68   if (!strcmp(optim, "Full")) {
69     p_updateMechanism = UM_FULL;
70     m_selectiveUpdate = select;
71   } else if (!strcmp(optim, "Lazy")) {
72     p_updateMechanism = UM_LAZY;
73     m_selectiveUpdate = 1;
74     xbt_assert((select == 1)
75                ||
76                (xbt_cfg_is_default_value
77                 (_sg_cfg_set, "cpu/maxmin_selective_update")),
78                "Disabling selective update while using the lazy update mechanism is dumb!");
79   } else {
80     xbt_die("Unsupported optimization (%s) for this model", optim);
81   }
82
83   cpu_running_action_set_that_does_not_need_being_checked =
84       xbt_swag_new(xbt_swag_offset(*action, p_stateHookup));
85
86   if (p_updateMechanism == UM_LAZY) {
87         shareResources = &CpuCas01Model::shareResourcesLazy;
88         updateActionsState = &CpuCas01Model::updateActionsStateLazy;
89
90   } else if (p_updateMechanism == UM_FULL) {
91         shareResources = &CpuCas01Model::shareResourcesFull;
92         updateActionsState = &CpuCas01Model::updateActionsStateFull;
93   } else
94     xbt_die("Invalid cpu update mechanism!");
95
96   if (!p_maxminSystem) {
97     p_maxminSystem = lmm_system_new(m_selectiveUpdate);
98   }
99
100   if (p_updateMechanism == UM_LAZY) {
101     p_actionHeap = xbt_heap_new(8, NULL);
102     xbt_heap_set_update_callback(p_actionHeap,  surf_action_lmm_update_index_heap);
103     p_modifiedSet = xbt_swag_new(xbt_swag_offset(*actionlmm, p_actionListHookup));
104     p_maxminSystem->keep_track = p_modifiedSet;
105   }
106 }
107
108 CpuCas01Model::~CpuCas01Model()
109 {
110   lmm_system_free(p_maxminSystem);
111   p_maxminSystem = NULL;
112
113   if (p_actionHeap)
114     xbt_heap_free(p_actionHeap);
115   xbt_swag_free(p_modifiedSet);
116
117   surf_cpu_model = NULL;
118
119   xbt_swag_free(cpu_running_action_set_that_does_not_need_being_checked);
120   cpu_running_action_set_that_does_not_need_being_checked = NULL;
121 }
122
123 void CpuCas01Model::parseInit(sg_platf_host_cbarg_t host)
124 {
125   createResource(host->id,
126         host->power_peak,
127         host->pstate,
128         host->power_scale,
129         host->power_trace,
130         host->core_amount,
131         host->initial_state,
132         host->state_trace,
133         host->properties);
134 }
135
136 CpuCas01LmmPtr CpuCas01Model::createResource(const char *name, xbt_dynar_t power_peak,
137                                   int pstate, double power_scale,
138                           tmgr_trace_t power_trace, int core,
139                           e_surf_resource_state_t state_initial,
140                           tmgr_trace_t state_trace,
141                           xbt_dict_t cpu_properties)
142 {
143   CpuPtr cpu = NULL;
144   xbt_assert(!surf_cpu_resource_priv(surf_cpu_resource_by_name(name)),
145              "Host '%s' declared several times in the platform file",
146              name);
147   xbt_assert(power_peak > 0, "Power has to be >0");
148   xbt_assert(core > 0, "Invalid number of cores %d", core);
149
150   cpu = new CpuCas01Lmm(this, name, power_peak, pstate, power_scale, power_trace, core, state_initial, state_trace, cpu_properties);
151   xbt_lib_set(host_lib, name, SURF_CPU_LEVEL, static_cast<ResourcePtr>(cpu));
152
153   return (CpuCas01LmmPtr) xbt_lib_get_elm_or_null(host_lib, name);
154 }
155
156 double CpuCas01Model::shareResourcesFull(double now)
157 {
158   return Model::shareResourcesMaxMin(p_runningActionSet,
159                              p_maxminSystem, lmm_solve);
160 }
161
162 void CpuCas01Model::addTraces()
163 {
164   xbt_dict_cursor_t cursor = NULL;
165   char *trace_name, *elm;
166   static int called = 0;
167   if (called)
168     return;
169   called = 1;
170
171   /* connect all traces relative to hosts */
172   xbt_dict_foreach(trace_connect_list_host_avail, cursor, trace_name, elm) {
173     tmgr_trace_t trace = (tmgr_trace_t) xbt_dict_get_or_null(traces_set_list, trace_name);
174     CpuCas01LmmPtr host = static_cast<CpuCas01LmmPtr>(surf_cpu_resource_priv(surf_cpu_resource_by_name(elm)));
175
176     xbt_assert(host, "Host %s undefined", elm);
177     xbt_assert(trace, "Trace %s undefined", trace_name);
178
179     host->p_stateEvent =
180         tmgr_history_add_trace(history, trace, 0.0, 0, static_cast<ResourcePtr>(host));
181   }
182
183   xbt_dict_foreach(trace_connect_list_power, cursor, trace_name, elm) {
184     tmgr_trace_t trace = (tmgr_trace_t) xbt_dict_get_or_null(traces_set_list, trace_name);
185     CpuCas01LmmPtr host = dynamic_cast<CpuCas01LmmPtr>(static_cast<ResourcePtr>(surf_cpu_resource_priv(surf_cpu_resource_by_name(elm))));
186
187     xbt_assert(host, "Host %s undefined", elm);
188     xbt_assert(trace, "Trace %s undefined", trace_name);
189
190     host->p_powerEvent =
191         tmgr_history_add_trace(history, trace, 0.0, 0, static_cast<ResourcePtr>(host));
192   }
193 }
194
195 /************
196  * Resource *
197  ************/
198 CpuCas01Lmm::CpuCas01Lmm(CpuCas01ModelPtr model, const char *name, xbt_dynar_t powerPeak,
199                 int pstate, double powerScale, tmgr_trace_t powerTrace, int core,
200         e_surf_resource_state_t stateInitial, tmgr_trace_t stateTrace,
201         xbt_dict_t properties) :
202         CpuLmm(model, name, properties), Resource(model, name, properties) {
203   m_powerPeak = xbt_dynar_get_as(powerPeak, pstate, double);
204   p_powerPeakList = powerPeak;
205   m_pstate = pstate;
206
207   p_energy = xbt_new(s_energy_cpu_cas01_t, 1);
208   p_energy->total_energy = 0;
209   p_energy->power_range_watts_list = getWattsRangeList();
210   p_energy->last_updated = surf_get_clock();
211
212   XBT_DEBUG("CPU create: peak=%f, pstate=%d", m_powerPeak, m_pstate);
213
214   m_powerScale = powerScale;
215   m_core = core;
216   p_stateCurrent = stateInitial;
217   if (powerTrace)
218     p_powerEvent = tmgr_history_add_trace(history, powerTrace, 0.0, 0, static_cast<ResourcePtr>(this));
219
220   if (stateTrace)
221     p_stateEvent = tmgr_history_add_trace(history, stateTrace, 0.0, 0, static_cast<ResourcePtr>(this));
222
223   p_constraint = lmm_constraint_new(p_model->p_maxminSystem, this, m_core * m_powerScale * m_powerPeak);
224 }
225
226 CpuCas01Lmm::~CpuCas01Lmm(){
227   unsigned int iter;
228   xbt_dynar_t power_tuple = NULL;
229   xbt_dynar_foreach(p_energy->power_range_watts_list, iter, power_tuple)
230     xbt_dynar_free(&power_tuple);
231   xbt_dynar_free(&p_energy->power_range_watts_list);
232   xbt_dynar_free(&p_powerPeakList);
233   xbt_free(p_energy);
234   return;
235 }
236
237 bool CpuCas01Lmm::isUsed()
238 {
239   return lmm_constraint_used(p_model->p_maxminSystem, p_constraint);
240 }
241
242 void CpuCas01Lmm::updateState(tmgr_trace_event_t event_type, double value, double date)
243 {
244   lmm_variable_t var = NULL;
245   lmm_element_t elem = NULL;
246
247   if (event_type == p_powerEvent) {
248     m_powerScale = value;
249     lmm_update_constraint_bound(surf_cpu_model->p_maxminSystem, p_constraint,
250                                 m_core * m_powerScale *
251                                 m_powerPeak);
252 #ifdef HAVE_TRACING
253     TRACE_surf_host_set_power(date, m_name,
254                               m_core * m_powerScale *
255                               m_powerPeak);
256 #endif
257     while ((var = lmm_get_var_from_cnst
258             (surf_cpu_model->p_maxminSystem, p_constraint, &elem))) {
259       CpuCas01ActionLmmPtr action = static_cast<CpuCas01ActionLmmPtr>(static_cast<ActionLmmPtr>(lmm_variable_id(var)));
260
261       lmm_update_variable_bound(surf_cpu_model->p_maxminSystem,
262                                 action->p_variable,
263                                 m_powerScale * m_powerPeak);
264     }
265     if (tmgr_trace_event_free(event_type))
266       p_powerEvent = NULL;
267   } else if (event_type == p_stateEvent) {
268     if (value > 0) {
269       if(p_stateCurrent == SURF_RESOURCE_OFF)
270         xbt_dynar_push_as(host_that_restart, char*, (char *)m_name);
271       p_stateCurrent = SURF_RESOURCE_ON;
272     } else {
273       lmm_constraint_t cnst = p_constraint;
274
275       p_stateCurrent = SURF_RESOURCE_OFF;
276
277       while ((var = lmm_get_var_from_cnst(surf_cpu_model->p_maxminSystem, cnst, &elem))) {
278         ActionLmmPtr action = static_cast<ActionLmmPtr>(lmm_variable_id(var));
279
280         if (action->getState() == SURF_ACTION_RUNNING ||
281             action->getState() == SURF_ACTION_READY ||
282             action->getState() == SURF_ACTION_NOT_IN_THE_SYSTEM) {
283           action->m_finish = date;
284           action->setState(SURF_ACTION_FAILED);
285         }
286       }
287     }
288     if (tmgr_trace_event_free(event_type))
289       p_stateEvent = NULL;
290   } else {
291     XBT_CRITICAL("Unknown event ! \n");
292     xbt_abort();
293   }
294
295   return;
296 }
297
298 ActionPtr CpuCas01Lmm::execute(double size)
299 {
300
301   XBT_IN("(%s,%g)", m_name, size);
302   CpuCas01ActionLmmPtr action = new CpuCas01ActionLmm(surf_cpu_model, size, p_stateCurrent != SURF_RESOURCE_ON);
303
304   action->m_suspended = 0;     /* Should be useless because of the
305                                                    calloc but it seems to help valgrind... */
306
307   action->p_variable =
308       lmm_variable_new(surf_cpu_model->p_maxminSystem, static_cast<ActionLmmPtr>(action),
309                        action->m_priority,
310                        m_powerScale * m_powerPeak, 1);
311   if (surf_cpu_model->p_updateMechanism == UM_LAZY) {
312     action->m_indexHeap = -1;
313     action->m_lastUpdate = surf_get_clock();
314     action->m_lastValue = 0.0;
315   }
316   lmm_expand(surf_cpu_model->p_maxminSystem, p_constraint,
317              action->p_variable, 1.0);
318   XBT_OUT();
319   return action;
320 }
321
322 ActionPtr CpuCas01Lmm::sleep(double duration)
323 {
324   if (duration > 0)
325     duration = MAX(duration, MAXMIN_PRECISION);
326
327   XBT_IN("(%s,%g)", m_name, duration);
328   CpuCas01ActionLmmPtr action = dynamic_cast<CpuCas01ActionLmmPtr>(execute(1.0));
329
330   // FIXME: sleep variables should not consume 1.0 in lmm_expand
331   action->m_maxDuration = duration;
332   action->m_suspended = 2;
333   if (duration == NO_MAX_DURATION) {
334     /* Move to the *end* of the corresponding action set. This convention
335        is used to speed up update_resource_state  */
336     xbt_swag_remove(static_cast<ActionPtr>(action), action->p_stateSet);
337     action->p_stateSet = cpu_running_action_set_that_does_not_need_being_checked;
338     xbt_swag_insert(static_cast<ActionPtr>(action), action->p_stateSet);
339   }
340
341   lmm_update_variable_weight(surf_cpu_model->p_maxminSystem,
342                              action->p_variable, 0.0);
343   if (surf_cpu_model->p_updateMechanism == UM_LAZY) {     // remove action from the heap
344     action->heapRemove(surf_cpu_model->p_actionHeap);
345     // this is necessary for a variable with weight 0 since such
346     // variables are ignored in lmm and we need to set its max_duration
347     // correctly at the next call to share_resources
348     xbt_swag_insert_at_head(static_cast<ActionLmmPtr>(action), surf_cpu_model->p_modifiedSet);
349   }
350
351   XBT_OUT();
352   return action;
353 }
354
355 xbt_dynar_t CpuCas01Lmm::getWattsRangeList()
356 {
357         xbt_dynar_t power_range_list;
358         xbt_dynar_t power_tuple;
359         int i = 0, pstate_nb=0;
360         xbt_dynar_t current_power_values;
361         double min_power, max_power;
362
363         if (m_properties == NULL)
364                 return NULL;
365
366         char* all_power_values_str = (char*)xbt_dict_get_or_null(m_properties, "power_per_state");
367
368         if (all_power_values_str == NULL)
369                 return NULL;
370
371
372         power_range_list = xbt_dynar_new(sizeof(xbt_dynar_t), NULL);
373         xbt_dynar_t all_power_values = xbt_str_split(all_power_values_str, ",");
374
375         pstate_nb = xbt_dynar_length(all_power_values);
376         for (i=0; i< pstate_nb; i++)
377         {
378                 /* retrieve the power values associated with the current pstate */
379                 current_power_values = xbt_str_split(xbt_dynar_get_as(all_power_values, i, char*), ":");
380                 xbt_assert(xbt_dynar_length(current_power_values) > 1,
381                                 "Power properties incorrectly defined - could not retrieve min and max power values for host %s",
382                                 m_name);
383
384                 /* min_power corresponds to the idle power (cpu load = 0) */
385                 /* max_power is the power consumed at 100% cpu load       */
386                 min_power = atof(xbt_dynar_get_as(current_power_values, 0, char*));
387                 max_power = atof(xbt_dynar_get_as(current_power_values, 1, char*));
388
389                 power_tuple = xbt_dynar_new(sizeof(double), NULL);
390                 xbt_dynar_push_as(power_tuple, double, min_power);
391                 xbt_dynar_push_as(power_tuple, double, max_power);
392
393                 xbt_dynar_push_as(power_range_list, xbt_dynar_t, power_tuple);
394                 xbt_dynar_free(&current_power_values);
395         }
396         xbt_dynar_free(&all_power_values);
397         return power_range_list;
398 }
399
400 /**
401  * Computes the power consumed by the host according to the current pstate and processor load
402  *
403  */
404 double CpuCas01Lmm::getCurrentWattsValue(double cpu_load)
405 {
406         xbt_dynar_t power_range_list = p_energy->power_range_watts_list;
407
408         if (power_range_list == NULL)
409         {
410                 XBT_DEBUG("No power range properties specified for host %s", m_name);
411                 return 0;
412         }
413         xbt_assert(xbt_dynar_length(power_range_list) == xbt_dynar_length(p_powerPeakList),
414                                                 "The number of power ranges in the properties does not match the number of pstates for host %s",
415                                                 m_name);
416
417     /* retrieve the power values associated with the current pstate */
418     xbt_dynar_t current_power_values = xbt_dynar_get_as(power_range_list, m_pstate, xbt_dynar_t);
419
420     /* min_power corresponds to the idle power (cpu load = 0) */
421     /* max_power is the power consumed at 100% cpu load       */
422     double min_power = xbt_dynar_get_as(current_power_values, 0, double);
423     double max_power = xbt_dynar_get_as(current_power_values, 1, double);
424     double power_slope = max_power - min_power;
425
426     double current_power = min_power + cpu_load * power_slope;
427
428         XBT_DEBUG("[get_current_watts] min_power=%f, max_power=%f, slope=%f", min_power, max_power, power_slope);
429     XBT_DEBUG("[get_current_watts] Current power (watts) = %f, load = %f", current_power, cpu_load);
430
431         return current_power;
432 }
433
434 /**
435  * Updates the total energy consumed as the sum of the current energy and
436  *                                               the energy consumed by the current action
437  */
438 void CpuCas01Lmm::updateEnergy(double cpu_load)
439 {
440   double start_time = p_energy->last_updated;
441   double finish_time = surf_get_clock();
442
443   XBT_DEBUG("[cpu_update_energy] action time interval=(%f-%f), current power peak=%f, current pstate=%d",
444                   start_time, finish_time, m_powerPeak, m_pstate);
445   double current_energy = p_energy->total_energy;
446   double action_energy = getCurrentWattsValue(cpu_load)*(finish_time-start_time);
447
448   p_energy->total_energy = current_energy + action_energy;
449   p_energy->last_updated = finish_time;
450
451   XBT_DEBUG("[cpu_update_energy] old_energy_value=%f, action_energy_value=%f", current_energy, action_energy);
452 }
453
454 double CpuCas01Lmm::getCurrentPowerPeak()
455 {
456   return m_powerPeak;
457 }
458
459 double CpuCas01Lmm::getPowerPeakAt(int pstate_index)
460 {
461   xbt_dynar_t plist = p_powerPeakList;
462   xbt_assert((pstate_index <= xbt_dynar_length(plist)), "Invalid parameters (pstate index out of bounds)");
463
464   return xbt_dynar_get_as(plist, pstate_index, double);
465 }
466
467 int CpuCas01Lmm::getNbPstates()
468 {
469   return xbt_dynar_length(p_powerPeakList);
470 }
471
472 void CpuCas01Lmm::setPowerPeakAt(int pstate_index)
473 {
474   xbt_dynar_t plist = p_powerPeakList;
475   xbt_assert((pstate_index <= xbt_dynar_length(plist)), "Invalid parameters (pstate index out of bounds)");
476
477   double new_power_peak = xbt_dynar_get_as(plist, pstate_index, double);
478   m_pstate = pstate_index;
479   m_powerPeak = new_power_peak;
480 }
481
482 double CpuCas01Lmm::getConsumedEnergy()
483 {
484   return p_energy->total_energy;
485 }
486
487 /**********
488  * Action *
489  **********/
490
491 /**
492  * Update the CPU total energy for a finished action
493  *
494  */
495 void CpuCas01ActionLmm::updateEnergy()
496 {
497   CpuCas01LmmPtr cpu  = static_cast<CpuCas01LmmPtr>(lmm_constraint_id(lmm_get_cnst_from_var
498                                                                                   (p_model->p_maxminSystem,
499                                                                                                   p_variable, 0)));
500
501   if(cpu->p_energy->last_updated < surf_get_clock()) {
502         double load = lmm_constraint_get_usage(cpu->p_constraint) / cpu->m_powerPeak;
503     cpu->updateEnergy(load);
504   }
505 }