src/surf/cpu_ti.c

   1
   2 /* Copyright (c) 2009, 2010. The SimGrid Team.
   3  * All rights reserved.                                                     */
   4
   5 /* This program is free software; you can redistribute it and/or modify it
   6  * under the terms of the license (GNU LGPL) which comes with this package. */
   7
   8 /*
   9         commit: e2d6799c4182f00443b3013aadb1c2412372460f
  10         This commit retrieves the old implementation of CPU_TI with multi-levels.
  11 */
  12
  13 #include "surf_private.h"
  14 #include "trace_mgr_private.h"
  15 #include "cpu_ti_private.h"
  16 #include "xbt/heap.h"
  17
  18
  19 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(surf_cpu_ti, surf,
  20                                 "Logging specific to the SURF CPU TRACE INTEGRATION module");
  21
  22
  23 static xbt_swag_t
  24     cpu_ti_running_action_set_that_does_not_need_being_checked = NULL;
  25 static xbt_swag_t cpu_ti_modified_cpu = NULL;
  26 static xbt_heap_t cpu_ti_action_heap;
  27
  28 /* prototypes of new trace functions */
  29 static double surf_cpu_ti_integrate_trace(surf_cpu_ti_tgmr_t trace,
  30                                           double a, double b);
  31
  32
  33 static double surf_cpu_ti_solve_trace(surf_cpu_ti_tgmr_t trace, double a,
  34                                       double amount);
  35 static double surf_cpu_ti_solve_trace_somewhat_simple(surf_cpu_ti_tgmr_t
  36                                                       trace, double a,
  37                                                       double amount);
  38
  39 static void surf_cpu_ti_free_tmgr(surf_cpu_ti_tgmr_t trace);
  40
  41 static double surf_cpu_ti_integrate_trace_simple(surf_cpu_ti_trace_t trace,
  42                                                  double a, double b);
  43 static double surf_cpu_ti_integrate_trace_simple_point(surf_cpu_ti_trace_t
  44                                                        trace, double a);
  45 static double surf_cpu_ti_solve_trace_simple(surf_cpu_ti_trace_t trace,
  46                                              double a, double amount);
  47 static int surf_cpu_ti_binary_search(double *array, double a, int low,
  48                                      int high);
  49 /* end prototypes */
  50
  51 static void surf_cpu_ti_free_trace(surf_cpu_ti_trace_t trace)
  52 {
  53   if (trace->time_points)
  54     xbt_free(trace->time_points);
  55   if (trace->integral)
  56     xbt_free(trace->integral);
  57   xbt_free(trace);
  58 }
  59
  60 static void surf_cpu_ti_free_tmgr(surf_cpu_ti_tgmr_t trace)
  61 {
  62   if (trace->trace)
  63     surf_cpu_ti_free_trace(trace->trace);
  64   xbt_free(trace);
  65 }
  66
  67 static surf_cpu_ti_trace_t surf_cpu_ti_trace_new(tmgr_trace_t power_trace)
  68 {
  69   surf_cpu_ti_trace_t trace;
  70   s_tmgr_event_t val;
  71   unsigned int cpt;
  72   double integral = 0;
  73   double time = 0;
  74   int i = 0;
  75   trace = xbt_new0(s_surf_cpu_ti_trace_t, 1);
  76   trace->time_points =
  77       xbt_malloc0(sizeof(double) *
  78                   (xbt_dynar_length(power_trace->event_list) + 1));
  79   trace->integral =
  80       xbt_malloc0(sizeof(double) *
  81                   (xbt_dynar_length(power_trace->event_list) + 1));
  82   trace->nb_points = xbt_dynar_length(power_trace->event_list);
  83   xbt_dynar_foreach(power_trace->event_list, cpt, val) {
  84     trace->time_points[i] = time;
  85     trace->integral[i] = integral;
  86     integral += val.delta * val.value;
  87     time += val.delta;
  88     i++;
  89   }
  90   trace->time_points[i] = time;
  91   trace->integral[i] = integral;
  92   return trace;
  93 }
  94
  95 /**
  96 * \brief Creates a new integration trace from a tmgr_trace_t
  97 *
  98 * \param        power_trace             CPU availability trace
  99 * \param        value                                   Percentage of CPU power available (useful to fixed tracing)
 100 * \param        spacing                         Initial spacing
 101 * \return       Integration trace structure
 102 */
 103 static surf_cpu_ti_tgmr_t cpu_ti_parse_trace(tmgr_trace_t power_trace,
 104                                              double value)
 105 {
 106   surf_cpu_ti_tgmr_t trace;
 107   double total_time = 0.0;
 108   s_tmgr_event_t val;
 109   unsigned int cpt;
 110   trace = xbt_new0(s_surf_cpu_ti_tgmr_t, 1);
 111
 112 /* no availability file, fixed trace */
 113   if (!power_trace) {
 114     trace->type = TRACE_FIXED;
 115     trace->value = value;
 116     XBT_DEBUG("No availabily trace. Constant value = %lf", value);
 117     return trace;
 118   }
 119
 120   /* only one point available, fixed trace */
 121   if (xbt_dynar_length(power_trace->event_list) == 1) {
 122     xbt_dynar_get_cpy(power_trace->event_list, 0, &val);
 123     trace->type = TRACE_FIXED;
 124     trace->value = val.value;
 125     return trace;
 126   }
 127
 128   trace->type = TRACE_DYNAMIC;
 129   trace->power_trace = power_trace;
 130
 131   /* count the total time of trace file */
 132   xbt_dynar_foreach(power_trace->event_list, cpt, val) {
 133     total_time += val.delta;
 134   }
 135   trace->trace = surf_cpu_ti_trace_new(power_trace);
 136   trace->last_time = total_time;
 137   trace->total =
 138       surf_cpu_ti_integrate_trace_simple(trace->trace, 0, total_time);
 139
 140   XBT_DEBUG("Total integral %lf, last_time %lf ",
 141          trace->total, trace->last_time);
 142
 143   return trace;
 144 }
 145
 146
 147 static void* cpu_ti_create_resource(char *name, double power_peak,
 148                            double power_scale,
 149                            tmgr_trace_t power_trace,
 150                            int core,
 151                            e_surf_resource_state_t state_initial,
 152                            tmgr_trace_t state_trace,
 153                            xbt_dict_t cpu_properties)
 154 {
 155   tmgr_trace_t empty_trace;
 156   s_tmgr_event_t val;
 157   cpu_ti_t cpu = xbt_new0(s_cpu_ti_t, 1);
 158   s_surf_action_cpu_ti_t ti_action;
 159   xbt_assert(core==1,"Multi-core not handled with this model yet");
 160   xbt_assert(!surf_cpu_resource_by_name(name),
 161               "Host '%s' declared several times in the platform file",
 162               name);
 163   xbt_assert(core==1,"Multi-core not handled with this model yet");
 164   cpu->action_set =
 165       xbt_swag_new(xbt_swag_offset(ti_action, cpu_list_hookup));
 166   cpu->generic_resource.model = surf_cpu_model;
 167   cpu->generic_resource.name = name;
 168   cpu->generic_resource.properties = cpu_properties;
 169   cpu->power_peak = power_peak;
 170   xbt_assert(cpu->power_peak > 0, "Power has to be >0");
 171   XBT_DEBUG("power scale %lf", power_scale);
 172   cpu->power_scale = power_scale;
 173   cpu->avail_trace = cpu_ti_parse_trace(power_trace, power_scale);
 174   cpu->state_current = state_initial;
 175   if (state_trace)
 176     cpu->state_event =
 177         tmgr_history_add_trace(history, state_trace, 0.0, 0, cpu);
 178   if (power_trace && xbt_dynar_length(power_trace->event_list) > 1) {
 179     /* add a fake trace event if periodicity == 0 */
 180     xbt_dynar_get_cpy(power_trace->event_list,
 181                       xbt_dynar_length(power_trace->event_list) - 1, &val);
 182     if (val.delta == 0) {
 183       empty_trace = tmgr_empty_trace_new();
 184       cpu->power_event =
 185           tmgr_history_add_trace(history, empty_trace,
 186                                  cpu->avail_trace->last_time, 0, cpu);
 187     }
 188   }
 189   xbt_lib_set(host_lib, name, SURF_CPU_LEVEL, cpu);
 190
 191   return cpu;
 192 }
 193
 194
 195 static void parse_cpu_ti_init(void)
 196 {
 197   double power_peak = 0.0;
 198   double power_scale = 0.0;
 199   int core = 0;
 200   tmgr_trace_t power_trace = NULL;
 201   e_surf_resource_state_t state_initial = SURF_RESOURCE_OFF;
 202   tmgr_trace_t state_trace = NULL;
 203
 204   power_peak = get_cpu_power(A_surfxml_host_power);
 205   surf_parse_get_double(&power_scale, A_surfxml_host_availability);
 206   power_trace = tmgr_trace_new(A_surfxml_host_availability_file);
 207   surf_parse_get_int(&core, A_surfxml_host_core);
 208
 209   xbt_assert((A_surfxml_host_state == A_surfxml_host_state_ON) ||
 210               (A_surfxml_host_state == A_surfxml_host_state_OFF),
 211               "Invalid state");
 212   if (A_surfxml_host_state == A_surfxml_host_state_ON)
 213     state_initial = SURF_RESOURCE_ON;
 214   if (A_surfxml_host_state == A_surfxml_host_state_OFF)
 215     state_initial = SURF_RESOURCE_OFF;
 216   state_trace = tmgr_trace_new(A_surfxml_host_state_file);
 217
 218   cpu_ti_create_resource(xbt_strdup(A_surfxml_host_id), power_peak, power_scale,
 219              power_trace, core, state_initial, state_trace,
 220              current_property_set);
 221   current_property_set = NULL;
 222
 223 }
 224
 225 static void add_traces_cpu_ti(void)
 226 {
 227   xbt_dict_cursor_t cursor = NULL;
 228   char *trace_name, *elm;
 229
 230   static int called = 0;
 231
 232   if (called)
 233     return;
 234   called = 1;
 235
 236 /* connect all traces relative to hosts */
 237   xbt_dict_foreach(trace_connect_list_host_avail, cursor, trace_name, elm) {
 238     tmgr_trace_t trace = xbt_dict_get_or_null(traces_set_list, trace_name);
 239     cpu_ti_t cpu = surf_cpu_resource_by_name(elm);
 240
 241     xbt_assert(cpu, "Host %s undefined", elm);
 242     xbt_assert(trace, "Trace %s undefined", trace_name);
 243
 244     if (cpu->state_event) {
 245       XBT_DEBUG("Trace already configured for this CPU(%s), ignoring it",
 246              elm);
 247       continue;
 248     }
 249     XBT_DEBUG("Add state trace: %s to CPU(%s)", trace_name, elm);
 250     cpu->state_event = tmgr_history_add_trace(history, trace, 0.0, 0, cpu);
 251   }
 252
 253   xbt_dict_foreach(trace_connect_list_power, cursor, trace_name, elm) {
 254     tmgr_trace_t trace = xbt_dict_get_or_null(traces_set_list, trace_name);
 255     cpu_ti_t cpu = surf_cpu_resource_by_name(elm);
 256
 257     xbt_assert(cpu, "Host %s undefined", elm);
 258     xbt_assert(trace, "Trace %s undefined", trace_name);
 259
 260     XBT_DEBUG("Add power trace: %s to CPU(%s)", trace_name, elm);
 261     if (cpu->avail_trace)
 262       surf_cpu_ti_free_tmgr(cpu->avail_trace);
 263
 264     cpu->avail_trace = cpu_ti_parse_trace(trace, cpu->power_scale);
 265
 266     /* add a fake trace event if periodicity == 0 */
 267     if (trace && xbt_dynar_length(trace->event_list) > 1) {
 268       s_tmgr_event_t val;
 269       xbt_dynar_get_cpy(trace->event_list,
 270                         xbt_dynar_length(trace->event_list) - 1, &val);
 271       if (val.delta == 0) {
 272         tmgr_trace_t empty_trace;
 273         empty_trace = tmgr_empty_trace_new();
 274         cpu->power_event =
 275             tmgr_history_add_trace(history, empty_trace,
 276                                    cpu->avail_trace->last_time, 0, cpu);
 277       }
 278     }
 279   }
 280 }
 281
 282 static void cpu_ti_define_callbacks(const char *file)
 283 {
 284   surfxml_add_callback(ETag_surfxml_host_cb_list, parse_cpu_ti_init);
 285   surfxml_add_callback(ETag_surfxml_platform_cb_list, &add_traces_cpu_ti);
 286 }
 287
 288 static int cpu_ti_resource_used(void *resource_id)
 289 {
 290   cpu_ti_t cpu = resource_id;
 291   return xbt_swag_size(cpu->action_set);
 292 }
 293
 294 static int cpu_ti_action_unref(surf_action_t action)
 295 {
 296   action->refcount--;
 297   if (!action->refcount) {
 298     xbt_swag_remove(action, action->state_set);
 299     /* remove from action_set */
 300     xbt_swag_remove(action, ACTION_GET_CPU(action)->action_set);
 301     /* remove from heap */
 302     xbt_heap_remove(cpu_ti_action_heap,
 303                     ((surf_action_cpu_ti_t) action)->index_heap);
 304     xbt_swag_insert(ACTION_GET_CPU(action), cpu_ti_modified_cpu);
 305     surf_action_free(&action);
 306     return 1;
 307   }
 308   return 0;
 309 }
 310
 311 static void cpu_ti_action_cancel(surf_action_t action)
 312 {
 313   surf_action_state_set(action, SURF_ACTION_FAILED);
 314   xbt_heap_remove(cpu_ti_action_heap,
 315                   ((surf_action_cpu_ti_t) action)->index_heap);
 316   xbt_swag_insert(ACTION_GET_CPU(action), cpu_ti_modified_cpu);
 317   return;
 318 }
 319
 320 static void cpu_ti_action_state_set(surf_action_t action,
 321                                     e_surf_action_state_t state)
 322 {
 323   surf_action_state_set(action, state);
 324   xbt_swag_insert(ACTION_GET_CPU(action), cpu_ti_modified_cpu);
 325   return;
 326 }
 327
 328 /**
 329 * \brief Update the remaining amount of actions
 330 *
 331 * \param        cpu             Cpu on which the actions are running
 332 * \param        now             Current time
 333 */
 334 static void cpu_ti_update_remaining_amount(cpu_ti_t cpu, double now)
 335 {
 336 #define GENERIC_ACTION(action) action->generic_action
 337   double area_total;
 338   surf_action_cpu_ti_t action;
 339
 340 /* already updated */
 341   if (cpu->last_update >= now)
 342     return;
 343
 344 /* calcule the surface */
 345   area_total =
 346       surf_cpu_ti_integrate_trace(cpu->avail_trace, cpu->last_update,
 347                                   now) * cpu->power_peak;
 348   XBT_DEBUG("Flops total: %lf, Last update %lf", area_total,
 349          cpu->last_update);
 350
 351   xbt_swag_foreach(action, cpu->action_set) {
 352     /* action not running, skip it */
 353     if (GENERIC_ACTION(action).state_set !=
 354         surf_cpu_model->states.running_action_set)
 355       continue;
 356
 357     /* bogus priority, skip it */
 358     if (GENERIC_ACTION(action).priority <= 0)
 359       continue;
 360
 361     /* action suspended, skip it */
 362     if (action->suspended != 0)
 363       continue;
 364
 365     /* action don't need update */
 366     if (GENERIC_ACTION(action).start >= now)
 367       continue;
 368
 369     /* skip action that are finishing now */
 370     if (GENERIC_ACTION(action).finish >= 0
 371         && GENERIC_ACTION(action).finish <= now)
 372       continue;
 373
 374     /* update remaining */
 375     double_update(&(GENERIC_ACTION(action).remains),
 376                   area_total / (cpu->sum_priority *
 377                                 GENERIC_ACTION(action).priority));
 378     XBT_DEBUG("Update remaining action(%p) remaining %lf", action,
 379            GENERIC_ACTION(action).remains);
 380   }
 381   cpu->last_update = now;
 382 #undef GENERIC_ACTION
 383 }
 384
 385 /**
 386 * \brief Update the finish date of action if necessary
 387 *
 388 * \param        cpu             Cpu on which the actions are running
 389 * \param        now             Current time
 390 */
 391 static void cpu_ti_update_action_finish_date(cpu_ti_t cpu, double now)
 392 {
 393 #define GENERIC_ACTION(action) action->generic_action
 394   surf_action_cpu_ti_t action;
 395   double sum_priority = 0.0, total_area, min_finish = -1;
 396
 397 /* update remaning amount of actions */
 398   cpu_ti_update_remaining_amount(cpu, now);
 399
 400   xbt_swag_foreach(action, cpu->action_set) {
 401     /* action not running, skip it */
 402     if (GENERIC_ACTION(action).state_set !=
 403         surf_cpu_model->states.running_action_set)
 404       continue;
 405
 406     /* bogus priority, skip it */
 407     if (GENERIC_ACTION(action).priority <= 0)
 408       continue;
 409
 410     /* action suspended, skip it */
 411     if (action->suspended != 0)
 412       continue;
 413
 414     sum_priority += 1.0 / GENERIC_ACTION(action).priority;
 415   }
 416   cpu->sum_priority = sum_priority;
 417
 418   xbt_swag_foreach(action, cpu->action_set) {
 419     min_finish = -1;
 420     /* action not running, skip it */
 421     if (GENERIC_ACTION(action).state_set !=
 422         surf_cpu_model->states.running_action_set)
 423       continue;
 424
 425     /* verify if the action is really running on cpu */
 426     if (action->suspended == 0 && GENERIC_ACTION(action).priority > 0) {
 427       /* total area needed to finish the action. Used in trace integration */
 428       total_area =
 429           (GENERIC_ACTION(action).remains) * sum_priority *
 430           GENERIC_ACTION(action).priority;
 431
 432       total_area /= cpu->power_peak;
 433
 434       GENERIC_ACTION(action).finish =
 435           surf_cpu_ti_solve_trace(cpu->avail_trace, now, total_area);
 436       /* verify which event will happen before (max_duration or finish time) */
 437       if ((GENERIC_ACTION(action).max_duration != NO_MAX_DURATION) &&
 438           (GENERIC_ACTION(action).start +
 439            GENERIC_ACTION(action).max_duration <
 440            GENERIC_ACTION(action).finish))
 441         min_finish = GENERIC_ACTION(action).start +
 442             GENERIC_ACTION(action).max_duration;
 443       else
 444         min_finish = GENERIC_ACTION(action).finish;
 445     } else {
 446       /* put the max duration time on heap */
 447       if (GENERIC_ACTION(action).max_duration != NO_MAX_DURATION)
 448         min_finish =
 449             (GENERIC_ACTION(action).start +
 450              GENERIC_ACTION(action).max_duration);
 451     }
 452     /* add in action heap */
 453     XBT_DEBUG("action(%p) index %d", action, action->index_heap);
 454     if (action->index_heap >= 0) {
 455       surf_action_cpu_ti_t heap_act =
 456           xbt_heap_remove(cpu_ti_action_heap, action->index_heap);
 457       if (heap_act != action)
 458         DIE_IMPOSSIBLE;
 459     }
 460     if (min_finish != NO_MAX_DURATION)
 461       xbt_heap_push(cpu_ti_action_heap, action, min_finish);
 462
 463     XBT_DEBUG
 464         ("Update finish time: Cpu(%s) Action: %p, Start Time: %lf Finish Time: %lf Max duration %lf",
 465          cpu->generic_resource.name, action, GENERIC_ACTION(action).start,
 466          GENERIC_ACTION(action).finish,
 467          GENERIC_ACTION(action).max_duration);
 468   }
 469 /* remove from modified cpu */
 470   xbt_swag_remove(cpu, cpu_ti_modified_cpu);
 471 #undef GENERIC_ACTION
 472 }
 473
 474 static double cpu_ti_share_resources(double now)
 475 {
 476   cpu_ti_t cpu, cpu_next;
 477   double min_action_duration = -1;
 478
 479 /* iterates over modified cpus to update share resources */
 480   xbt_swag_foreach_safe(cpu, cpu_next, cpu_ti_modified_cpu) {
 481     cpu_ti_update_action_finish_date(cpu, now);
 482   }
 483 /* get the min next event if heap not empty */
 484   if (xbt_heap_size(cpu_ti_action_heap) > 0)
 485     min_action_duration = xbt_heap_maxkey(cpu_ti_action_heap) - now;
 486
 487   XBT_DEBUG("Share resources, min next event date: %lf", min_action_duration);
 488
 489   return min_action_duration;
 490 }
 491
 492 static void cpu_ti_update_actions_state(double now, double delta)
 493 {
 494 #define GENERIC_ACTION(action) action->generic_action
 495   surf_action_cpu_ti_t action;
 496   while ((xbt_heap_size(cpu_ti_action_heap) > 0)
 497          && (xbt_heap_maxkey(cpu_ti_action_heap) <= now)) {
 498     action = xbt_heap_pop(cpu_ti_action_heap);
 499     XBT_DEBUG("Action %p: finish", action);
 500     GENERIC_ACTION(action).finish = surf_get_clock();
 501     /* set the remains to 0 due to precision problems when updating the remaining amount */
 502     GENERIC_ACTION(action).remains = 0;
 503     cpu_ti_action_state_set((surf_action_t) action, SURF_ACTION_DONE);
 504     /* update remaining amout of all actions */
 505     cpu_ti_update_remaining_amount(action->cpu, surf_get_clock());
 506   }
 507 #undef GENERIC_ACTION
 508 }
 509
 510 static void cpu_ti_update_resource_state(void *id,
 511                                          tmgr_trace_event_t event_type,
 512                                          double value, double date)
 513 {
 514   cpu_ti_t cpu = id;
 515   surf_action_cpu_ti_t action;
 516
 517   if (event_type == cpu->power_event) {
 518     tmgr_trace_t power_trace;
 519     surf_cpu_ti_tgmr_t trace;
 520     s_tmgr_event_t val;
 521
 522     XBT_DEBUG("Finish trace date: %lf value %lf date %lf", surf_get_clock(),
 523            value, date);
 524     /* update remaining of actions and put in modified cpu swag */
 525     cpu_ti_update_remaining_amount(cpu, date);
 526     xbt_swag_insert(cpu, cpu_ti_modified_cpu);
 527
 528     power_trace = cpu->avail_trace->power_trace;
 529     xbt_dynar_get_cpy(power_trace->event_list,
 530                       xbt_dynar_length(power_trace->event_list) - 1, &val);
 531     /* free old trace */
 532     surf_cpu_ti_free_tmgr(cpu->avail_trace);
 533     cpu->power_scale = val.value;
 534
 535     trace = xbt_new0(s_surf_cpu_ti_tgmr_t, 1);
 536     trace->type = TRACE_FIXED;
 537     trace->value = val.value;
 538     XBT_DEBUG("value %lf", val.value);
 539
 540     cpu->avail_trace = trace;
 541
 542     if (tmgr_trace_event_free(event_type))
 543       cpu->power_event = NULL;
 544
 545   } else if (event_type == cpu->state_event) {
 546     if (value > 0)
 547       cpu->state_current = SURF_RESOURCE_ON;
 548     else {
 549       cpu->state_current = SURF_RESOURCE_OFF;
 550
 551       /* put all action running on cpu to failed */
 552       xbt_swag_foreach(action, cpu->action_set) {
 553         if (surf_action_state_get((surf_action_t) action) ==
 554             SURF_ACTION_RUNNING
 555             || surf_action_state_get((surf_action_t) action) ==
 556             SURF_ACTION_READY
 557             || surf_action_state_get((surf_action_t) action) ==
 558             SURF_ACTION_NOT_IN_THE_SYSTEM) {
 559           action->generic_action.finish = date;
 560           cpu_ti_action_state_set((surf_action_t) action,
 561                                   SURF_ACTION_FAILED);
 562           if (action->index_heap >= 0) {
 563             surf_action_cpu_ti_t heap_act =
 564                 xbt_heap_remove(cpu_ti_action_heap, action->index_heap);
 565             if (heap_act != action)
 566               DIE_IMPOSSIBLE;
 567           }
 568         }
 569       }
 570     }
 571     if (tmgr_trace_event_free(event_type))
 572       cpu->state_event = NULL;
 573   } else {
 574     XBT_CRITICAL("Unknown event ! \n");
 575     xbt_abort();
 576   }
 577
 578   return;
 579 }
 580
 581 static surf_action_t cpu_ti_execute(void *cpu, double size)
 582 {
 583   surf_action_cpu_ti_t action = NULL;
 584   cpu_ti_t CPU = cpu;
 585
 586   XBT_IN("(%s,%g)", surf_resource_name(CPU), size);
 587   action =
 588       surf_action_new(sizeof(s_surf_action_cpu_ti_t), size, surf_cpu_model,
 589                       CPU->state_current != SURF_RESOURCE_ON);
 590   action->cpu = cpu;
 591   action->index_heap = -1;
 592
 593   xbt_swag_insert(CPU, cpu_ti_modified_cpu);
 594
 595   xbt_swag_insert(action, CPU->action_set);
 596
 597   action->suspended = 0;        /* Should be useless because of the
 598                                    calloc but it seems to help valgrind... */
 599
 600   XBT_OUT();
 601   return (surf_action_t) action;
 602 }
 603
 604 static void cpu_ti_action_update_index_heap(void *action, int i)
 605 {
 606   ((surf_action_cpu_ti_t) action)->index_heap = i;
 607 }
 608
 609 static surf_action_t cpu_ti_action_sleep(void *cpu, double duration)
 610 {
 611   surf_action_cpu_ti_t action = NULL;
 612
 613   if (duration > 0)
 614     duration = MAX(duration, MAXMIN_PRECISION);
 615
 616   XBT_IN("(%s,%g)", surf_resource_name(cpu), duration);
 617   action = (surf_action_cpu_ti_t) cpu_ti_execute(cpu, 1.0);
 618   action->generic_action.max_duration = duration;
 619   action->suspended = 2;
 620   if (duration == NO_MAX_DURATION) {
 621     /* Move to the *end* of the corresponding action set. This convention
 622        is used to speed up update_resource_state  */
 623     xbt_swag_remove(action, ((surf_action_t) action)->state_set);
 624     ((surf_action_t) action)->state_set =
 625         cpu_ti_running_action_set_that_does_not_need_being_checked;
 626     xbt_swag_insert(action, ((surf_action_t) action)->state_set);
 627   }
 628   XBT_OUT();
 629   return (surf_action_t) action;
 630 }
 631
 632 static void cpu_ti_action_suspend(surf_action_t action)
 633 {
 634   XBT_IN("(%p)", action);
 635   if (((surf_action_cpu_ti_t) action)->suspended != 2) {
 636     ((surf_action_cpu_ti_t) action)->suspended = 1;
 637     xbt_heap_remove(cpu_ti_action_heap,
 638                     ((surf_action_cpu_ti_t) action)->index_heap);
 639     xbt_swag_insert(ACTION_GET_CPU(action), cpu_ti_modified_cpu);
 640   }
 641   XBT_OUT();
 642 }
 643
 644 static void cpu_ti_action_resume(surf_action_t action)
 645 {
 646   XBT_IN("(%p)", action);
 647   if (((surf_action_cpu_ti_t) action)->suspended != 2) {
 648     ((surf_action_cpu_ti_t) action)->suspended = 0;
 649     xbt_swag_insert(ACTION_GET_CPU(action), cpu_ti_modified_cpu);
 650   }
 651   XBT_OUT();
 652 }
 653
 654 static int cpu_ti_action_is_suspended(surf_action_t action)
 655 {
 656   return (((surf_action_cpu_ti_t) action)->suspended == 1);
 657 }
 658
 659 static void cpu_ti_action_set_max_duration(surf_action_t action,
 660                                            double duration)
 661 {
 662   surf_action_cpu_ti_t ACT = (surf_action_cpu_ti_t) action;
 663   double min_finish;
 664
 665   XBT_IN("(%p,%g)", action, duration);
 666
 667   action->max_duration = duration;
 668
 669   if (duration >= 0)
 670     min_finish =
 671         (action->start + action->max_duration) <
 672         action->finish ? (action->start +
 673                           action->max_duration) : action->finish;
 674   else
 675     min_finish = action->finish;
 676
 677 /* add in action heap */
 678   if (ACT->index_heap >= 0) {
 679     surf_action_cpu_ti_t heap_act =
 680         xbt_heap_remove(cpu_ti_action_heap, ACT->index_heap);
 681     if (heap_act != ACT)
 682       DIE_IMPOSSIBLE;
 683   }
 684   xbt_heap_push(cpu_ti_action_heap, ACT, min_finish);
 685
 686   XBT_OUT();
 687 }
 688
 689 static void cpu_ti_action_set_priority(surf_action_t action,
 690                                        double priority)
 691 {
 692   XBT_IN("(%p,%g)", action, priority);
 693   action->priority = priority;
 694   xbt_swag_insert(ACTION_GET_CPU(action), cpu_ti_modified_cpu);
 695   XBT_OUT();
 696 }
 697
 698 static double cpu_ti_action_get_remains(surf_action_t action)
 699 {
 700   XBT_IN("(%p)", action);
 701   cpu_ti_update_remaining_amount((cpu_ti_t)
 702                                  ((surf_action_cpu_ti_t) action)->cpu,
 703                                  surf_get_clock());
 704   return action->remains;
 705   XBT_OUT();
 706 }
 707
 708 static e_surf_resource_state_t cpu_ti_get_state(void *cpu)
 709 {
 710   return ((cpu_ti_t) cpu)->state_current;
 711 }
 712
 713 static double cpu_ti_get_speed(void *cpu, double load)
 714 {
 715   return load * (((cpu_ti_t) cpu)->power_peak);
 716 }
 717
 718 /**
 719 * \brief Auxiliary function to update the CPU power scale.
 720 *
 721 *       This function uses the trace structure to return the power scale at the determined time a.
 722 * \param trace          Trace structure to search the updated power scale
 723 * \param a                              Time
 724 * \return CPU power scale
 725 */
 726 static double surf_cpu_ti_get_power_scale(surf_cpu_ti_tgmr_t trace,
 727                                           double a)
 728 {
 729   double reduced_a;
 730   int point;
 731   s_tmgr_event_t val;
 732
 733   reduced_a = a - floor(a / trace->last_time) * trace->last_time;
 734   point =
 735       surf_cpu_ti_binary_search(trace->trace->time_points, reduced_a, 0,
 736                                 trace->trace->nb_points - 1);
 737   xbt_dynar_get_cpy(trace->power_trace->event_list, point, &val);
 738   return val.value;
 739 }
 740
 741 static double cpu_ti_get_available_speed(void *cpu)
 742 {
 743   cpu_ti_t CPU = cpu;
 744   CPU->power_scale =
 745       surf_cpu_ti_get_power_scale(CPU->avail_trace, surf_get_clock());
 746 /* number between 0 and 1 */
 747   return CPU->power_scale;
 748 }
 749
 750 static void cpu_ti_finalize(void)
 751 {
 752   void **cpu;
 753   xbt_lib_cursor_t cursor;
 754   char *key;
 755
 756   xbt_lib_foreach(host_lib, cursor, key, cpu){
 757           if(cpu[SURF_CPU_LEVEL])
 758           {
 759                     cpu_ti_t CPU = cpu[SURF_CPU_LEVEL];
 760                     xbt_swag_free(CPU->action_set);
 761                     surf_cpu_ti_free_tmgr(CPU->avail_trace);
 762           }
 763   }
 764
 765   surf_model_exit(surf_cpu_model);
 766   surf_cpu_model = NULL;
 767
 768   xbt_swag_free
 769       (cpu_ti_running_action_set_that_does_not_need_being_checked);
 770   xbt_swag_free(cpu_ti_modified_cpu);
 771   cpu_ti_running_action_set_that_does_not_need_being_checked = NULL;
 772   xbt_heap_free(cpu_ti_action_heap);
 773 }
 774
 775 static void surf_cpu_ti_model_init_internal(void)
 776 {
 777   s_surf_action_t action;
 778   s_cpu_ti_t cpu;
 779
 780   surf_cpu_model = surf_model_init();
 781
 782   cpu_ti_running_action_set_that_does_not_need_being_checked =
 783       xbt_swag_new(xbt_swag_offset(action, state_hookup));
 784
 785   cpu_ti_modified_cpu =
 786       xbt_swag_new(xbt_swag_offset(cpu, modified_cpu_hookup));
 787
 788   surf_cpu_model->name = "CPU_TI";
 789
 790   surf_cpu_model->action_unref = cpu_ti_action_unref;
 791   surf_cpu_model->action_cancel = cpu_ti_action_cancel;
 792   surf_cpu_model->action_state_set = cpu_ti_action_state_set;
 793
 794   surf_cpu_model->model_private->resource_used = cpu_ti_resource_used;
 795   surf_cpu_model->model_private->share_resources = cpu_ti_share_resources;
 796   surf_cpu_model->model_private->update_actions_state =
 797       cpu_ti_update_actions_state;
 798   surf_cpu_model->model_private->update_resource_state =
 799       cpu_ti_update_resource_state;
 800   surf_cpu_model->model_private->finalize = cpu_ti_finalize;
 801
 802   surf_cpu_model->suspend = cpu_ti_action_suspend;
 803   surf_cpu_model->resume = cpu_ti_action_resume;
 804   surf_cpu_model->is_suspended = cpu_ti_action_is_suspended;
 805   surf_cpu_model->set_max_duration = cpu_ti_action_set_max_duration;
 806   surf_cpu_model->set_priority = cpu_ti_action_set_priority;
 807   surf_cpu_model->get_remains = cpu_ti_action_get_remains;
 808
 809   surf_cpu_model->extension.cpu.execute = cpu_ti_execute;
 810   surf_cpu_model->extension.cpu.sleep = cpu_ti_action_sleep;
 811
 812   surf_cpu_model->extension.cpu.get_state = cpu_ti_get_state;
 813   surf_cpu_model->extension.cpu.get_speed = cpu_ti_get_speed;
 814   surf_cpu_model->extension.cpu.get_available_speed =
 815       cpu_ti_get_available_speed;
 816   surf_cpu_model->extension.cpu.create_resource = cpu_ti_create_resource;
 817   surf_cpu_model->extension.cpu.add_traces = add_traces_cpu_ti;
 818
 819   cpu_ti_action_heap = xbt_heap_new(8, NULL);
 820   xbt_heap_set_update_callback(cpu_ti_action_heap,
 821                                cpu_ti_action_update_index_heap);
 822
 823 }
 824
 825 void surf_cpu_model_init_ti(const char *filename)
 826 {
 827   if (surf_cpu_model)
 828     return;
 829   surf_cpu_ti_model_init_internal();
 830   cpu_ti_define_callbacks(filename);
 831   xbt_dynar_push(model_list, &surf_cpu_model);
 832 }
 833
 834
 835 /**
 836 * \brief Integrate trace
 837 *
 838 * Wrapper around surf_cpu_integrate_trace_simple() to get
 839 * the cyclic effect.
 840 *
 841 * \param trace Trace structure.
 842 * \param a                      Begin of interval
 843 * \param b                      End of interval
 844 * \return the integrate value. -1 if an error occurs.
 845 */
 846 static double surf_cpu_ti_integrate_trace(surf_cpu_ti_tgmr_t trace,
 847                                           double a, double b)
 848 {
 849   double first_chunk;
 850   double middle_chunk;
 851   double last_chunk;
 852   int a_index, b_index;
 853
 854   if ((a < 0.0) || (a > b)) {
 855     XBT_CRITICAL
 856         ("Error, invalid integration interval [%.2f,%.2f]. You probably have a task executing with negative computation amount. Check your code.",
 857          a, b);
 858     xbt_abort();
 859   }
 860   if (a == b)
 861     return 0.0;
 862
 863   if (trace->type == TRACE_FIXED) {
 864     return ((b - a) * trace->value);
 865   }
 866
 867   if (ceil(a / trace->last_time) == a / trace->last_time)
 868     a_index = 1 + (int) (ceil(a / trace->last_time));
 869   else
 870     a_index = (int) (ceil(a / trace->last_time));
 871
 872   b_index = (int) (floor(b / trace->last_time));
 873
 874   if (a_index > b_index) {      /* Same chunk */
 875     return surf_cpu_ti_integrate_trace_simple(trace->trace,
 876                                               a - (a_index -
 877                                                    1) * trace->last_time,
 878                                               b -
 879                                               (b_index) *
 880                                               trace->last_time);
 881   }
 882
 883   first_chunk = surf_cpu_ti_integrate_trace_simple(trace->trace,
 884                                                    a - (a_index -
 885                                                         1) *
 886                                                    trace->last_time,
 887                                                    trace->last_time);
 888   middle_chunk = (b_index - a_index) * trace->total;
 889   last_chunk = surf_cpu_ti_integrate_trace_simple(trace->trace,
 890                                                   0.0,
 891                                                   b -
 892                                                   (b_index) *
 893                                                   trace->last_time);
 894
 895   XBT_DEBUG("first_chunk=%.2f  middle_chunk=%.2f  last_chunk=%.2f\n",
 896          first_chunk, middle_chunk, last_chunk);
 897
 898   return (first_chunk + middle_chunk + last_chunk);
 899 }
 900
 901 /**
 902  * \brief Auxiliary function to calculate the integral between a and b.
 903  *              It simply calculates the integral at point a and b and returns the difference
 904  *      between them.
 905  * \param trace         Trace structure
 906  * \param a                             Initial point
 907  * \param b     Final point
 908  * \return      Integral
 909 */
 910 static double surf_cpu_ti_integrate_trace_simple(surf_cpu_ti_trace_t trace,
 911                                                  double a, double b)
 912 {
 913   return surf_cpu_ti_integrate_trace_simple_point(trace,
 914                                                   b) -
 915       surf_cpu_ti_integrate_trace_simple_point(trace, a);
 916 }
 917
 918 /**
 919  * \brief Auxiliary function to calculate the integral at point a.
 920  * \param trace         Trace structure
 921  * \param a                             point
 922  * \return      Integral
 923 */
 924 static double surf_cpu_ti_integrate_trace_simple_point(surf_cpu_ti_trace_t
 925                                                        trace, double a)
 926 {
 927   double integral = 0;
 928   int ind;
 929   double a_aux = a;
 930   ind =
 931       surf_cpu_ti_binary_search(trace->time_points, a, 0,
 932                                 trace->nb_points - 1);
 933   integral += trace->integral[ind];
 934   XBT_DEBUG
 935       ("a %lf ind %d integral %lf ind + 1 %lf ind %lf time +1 %lf time %lf",
 936        a, ind, integral, trace->integral[ind + 1], trace->integral[ind],
 937        trace->time_points[ind + 1], trace->time_points[ind]);
 938   double_update(&a_aux, trace->time_points[ind]);
 939   if (a_aux > 0)
 940     integral +=
 941         ((trace->integral[ind + 1] -
 942           trace->integral[ind]) / (trace->time_points[ind + 1] -
 943                                    trace->time_points[ind])) * (a -
 944                                                                 trace->
 945                                                                 time_points
 946                                                                 [ind]);
 947   XBT_DEBUG("Integral a %lf = %lf", a, integral);
 948
 949   return integral;
 950 }
 951
 952 /**
 953 * \brief Calculate the time needed to execute "amount" on cpu.
 954 *
 955 * Here, amount can span multiple trace periods
 956 *
 957 * \param trace  CPU trace structure
 958 * \param a                              Initial time
 959 * \param amount Amount to be executed
 960 * \return       End time
 961 */
 962 static double surf_cpu_ti_solve_trace(surf_cpu_ti_tgmr_t trace, double a,
 963                                       double amount)
 964 {
 965   int quotient;
 966   double reduced_b;
 967   double reduced_amount;
 968   double reduced_a;
 969   double b;
 970
 971 /* Fix very small negative numbers */
 972   if ((a < 0.0) && (a > -EPSILON)) {
 973     a = 0.0;
 974   }
 975   if ((amount < 0.0) && (amount > -EPSILON)) {
 976     amount = 0.0;
 977   }
 978
 979 /* Sanity checks */
 980   if ((a < 0.0) || (amount < 0.0)) {
 981     XBT_CRITICAL
 982         ("Error, invalid parameters [a = %.2f, amount = %.2f]. You probably have a task executing with negative computation amount. Check your code.",
 983          a, amount);
 984     xbt_abort();
 985   }
 986
 987 /* At this point, a and amount are positive */
 988
 989   if (amount < EPSILON)
 990     return a;
 991
 992 /* Is the trace fixed ? */
 993   if (trace->type == TRACE_FIXED) {
 994     return (a + (amount / trace->value));
 995   }
 996
 997   XBT_DEBUG("amount %lf total %lf", amount, trace->total);
 998 /* Reduce the problem to one where amount <= trace_total */
 999   quotient = (int) (floor(amount / trace->total));
1000   reduced_amount = (trace->total) * ((amount / trace->total) -
1001                                      floor(amount / trace->total));
1002   reduced_a = a - (trace->last_time) * (int) (floor(a / trace->last_time));
1003
1004   XBT_DEBUG("Quotient: %d reduced_amount: %lf reduced_a: %lf", quotient,
1005          reduced_amount, reduced_a);
1006
1007 /* Now solve for new_amount which is <= trace_total */
1008 /*
1009          fprintf(stderr,"reduced_a = %.2f\n",reduced_a);
1010          fprintf(stderr,"reduced_amount = %.2f\n",reduced_amount);
1011  */
1012   reduced_b =
1013       surf_cpu_ti_solve_trace_somewhat_simple(trace, reduced_a,
1014                                               reduced_amount);
1015
1016 /* Re-map to the original b and amount */
1017   b = (trace->last_time) * (int) (floor(a / trace->last_time)) +
1018       (quotient * trace->last_time) + reduced_b;
1019   return b;
1020 }
1021
1022 /**
1023 * \brief Auxiliary function to solve integral
1024 *
1025 * Here, amount is <= trace->total
1026 * and a <=trace->last_time
1027 *
1028 */
1029 static double surf_cpu_ti_solve_trace_somewhat_simple(surf_cpu_ti_tgmr_t
1030                                                       trace, double a,
1031                                                       double amount)
1032 {
1033   double amount_till_end;
1034   double b;
1035
1036   XBT_DEBUG("Solve integral: [%.2f, amount=%.2f]", a, amount);
1037   amount_till_end =
1038       surf_cpu_ti_integrate_trace(trace, a, trace->last_time);
1039 /*
1040          fprintf(stderr,"amount_till_end=%.2f\n",amount_till_end);
1041  */
1042
1043   if (amount_till_end > amount) {
1044     b = surf_cpu_ti_solve_trace_simple(trace->trace, a, amount);
1045   } else {
1046     b = trace->last_time +
1047         surf_cpu_ti_solve_trace_simple(trace->trace, 0.0,
1048                                        amount - amount_till_end);
1049   }
1050   return b;
1051 }
1052
1053 /**
1054  * \brief Auxiliary function to solve integral.
1055  *      It returns the date when the requested amount of flops is available
1056  * \param trace         Trace structure
1057  * \param a                             Initial point
1058  * \param amount        Amount of flops
1059  * \return The date when amount is available.
1060 */
1061 static double surf_cpu_ti_solve_trace_simple(surf_cpu_ti_trace_t trace,
1062                                              double a, double amount)
1063 {
1064   double integral_a;
1065   int ind;
1066   double time;
1067   integral_a = surf_cpu_ti_integrate_trace_simple_point(trace, a);
1068   ind =
1069       surf_cpu_ti_binary_search(trace->integral, integral_a + amount, 0,
1070                                 trace->nb_points - 1);
1071   time = trace->time_points[ind];
1072   time +=
1073       (integral_a + amount -
1074        trace->integral[ind]) / ((trace->integral[ind + 1] -
1075                                  trace->integral[ind]) /
1076                                 (trace->time_points[ind + 1] -
1077                                  trace->time_points[ind]));
1078
1079   return time;
1080 }
1081
1082 /**
1083  * \brief Binary search in array.
1084  *      It returns the first point of the interval in which "a" is.
1085  * \param array         Array
1086  * \param a                             Value to search
1087  * \param low           Low bound to search in array
1088  * \param high          Upper bound to search in array
1089  * \return Index of point
1090 */
1091 static int surf_cpu_ti_binary_search(double *array, double a, int low,
1092                                      int high)
1093 {
1094   xbt_assert(low < high, "Wrong parameters: low (%d) should be smaller than"
1095       " high (%d)", low, high);
1096
1097   int mid;
1098   do {
1099     mid = low + (high - low) / 2;
1100     XBT_DEBUG("a %lf low %d high %d mid %d value %lf", a, low, high, mid,
1101         array[mid]);
1102
1103     if (array[mid] > a)
1104       high = mid;
1105     else
1106       low = mid;
1107   }
1108   while (low < high - 1);
1109
1110   return low;
1111 }