src/surf/surf_interface.cpp

   1 /* Copyright (c) 2004-2014. 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 "surf_private.h"
   8 #include "surf_interface.hpp"
   9 #include "network_interface.hpp"
  10 #include "cpu_interface.hpp"
  11 #include "workstation_interface.hpp"
  12 #include "vm_workstation_interface.hpp"
  13 #include "simix/smx_host_private.h"
  14 #include "surf_routing.hpp"
  15 #include "simgrid/sg_config.h"
  16 #include "mc/mc.h"
  17
  18 XBT_LOG_NEW_CATEGORY(surf, "All SURF categories");
  19 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(surf_kernel, surf,
  20                                 "Logging specific to SURF (kernel)");
  21
  22 /*********
  23  * Utils *
  24  *********/
  25
  26 /* This function is a pimple that we ought to fix. But it won't be easy.
  27  *
  28  * The surf_solve() function does properly return the set of actions that
  29  * changed. Instead, each model change a global data, and then the caller of
  30  * surf_solve must pick into these sets of action_failed and action_done.
  31  *
  32  * This was not clean but ok as long as we didn't had to restart the processes
  33  * when the resource comes back up.
  34  * We worked by putting sentinel actions on every resources we are interested
  35  * in, so that surf informs us if/when the corresponding resource fails.
  36  *
  37  * But this does not work to get Simix informed of when a resource comes back
  38  * up, and this is where this pimple comes. We have a set of resources that are
  39  * currently down and for which simix needs to know when it comes back up.
  40  * And the current function is called *at every simulation step* to sweep over
  41  * that set, searching for a resource that was turned back up in the meanwhile.
  42  * This is UGLY and slow.
  43  *
  44  * The proper solution would be to not rely on globals for the action_failed and
  45  * action_done swags. They must be passed as parameter by the caller (the
  46  * handling of these actions in simix may let you think that these two sets can
  47  * be merged, but their handling in SimDag induce the contrary unless this
  48  * simdag code can check by itself whether the action is done of failed -- seems
  49  * very doable, but yet more cleanup to do).
  50  *
  51  * Once surf_solve() is passed the set of actions that changed, you want to add
  52  * a new set of resources back up as parameter to this function. You also want
  53  * to add a boolean field "restart_watched" to each resource, and make sure that
  54  * whenever a resource with this field enabled comes back up, it's added to that
  55  * set so that Simix sees it and react accordingly. This would kill that need
  56  * for surf to call simix.
  57  *
  58  */
  59
  60 /*static void remove_watched_host(void *key)
  61 {
  62   xbt_dict_remove(watched_hosts_lib, *(char**)key);
  63 }*/
  64
  65 /*void surf_watched_hosts(void)
  66 {
  67   char *key;
  68   void *host;
  69   xbt_dict_cursor_t cursor;
  70   xbt_dynar_t hosts = xbt_dynar_new(sizeof(char*), NULL);
  71
  72   XBT_DEBUG("Check for host SURF_RESOURCE_ON on watched_hosts_lib");
  73   xbt_dict_foreach(watched_hosts_lib, cursor, key, host)
  74   {
  75     if(SIMIX_host_get_state((smx_host_t)host) == SURF_RESOURCE_ON){
  76       XBT_INFO("Restart processes on host: %s", SIMIX_host_get_name((smx_host_t)host));
  77       SIMIX_host_autorestart((smx_host_t)host);
  78       xbt_dynar_push_as(hosts, char*, key);
  79     }
  80     else
  81       XBT_DEBUG("See SURF_RESOURCE_OFF on host: %s",key);
  82   }
  83   xbt_dynar_map(hosts, remove_watched_host);
  84   xbt_dynar_free(&hosts);
  85 }*/
  86
  87 /* model_list_invoke contains only surf_workstation and surf_vm_workstation.
  88  * The callback functions of cpu_model and network_model will be called from
  89  * those of these workstation models. */
  90 xbt_dynar_t model_list = NULL; /* for destroying all models correctly */
  91 xbt_dynar_t model_list_invoke = NULL;  /* for invoking callbacks */
  92
  93 tmgr_history_t history = NULL;
  94 lmm_system_t maxmin_system = NULL;
  95 xbt_dynar_t surf_path = NULL;
  96 xbt_dynar_t host_that_restart = NULL;
  97 xbt_dict_t watched_hosts_lib;
  98
  99 surf_callback(void, void) surfExitCallbacks;
 100
 101 s_surf_model_description_t surf_plugin_description[] = {
 102     {"Energy", "Cpu energy consumption.", sg_energy_plugin_init},
 103      {NULL, NULL,  NULL}      /* this array must be NULL terminated */
 104 };
 105
 106 /* Don't forget to update the option description in smx_config when you change
 107    this */
 108 s_surf_model_description_t surf_network_model_description[] = {
 109   {"LV08",
 110    "Realistic network analytic model (slow-start modeled by multiplying latency by 10.4, bandwidth by .92; bottleneck sharing uses a payload of S=8775 for evaluating RTT). ",
 111    surf_network_model_init_LegrandVelho},
 112   {"Constant",
 113    "Simplistic network model where all communication take a constant time (one second). This model provides the lowest realism, but is (marginally) faster.",
 114    surf_network_model_init_Constant},
 115   {"SMPI",
 116    "Realistic network model specifically tailored for HPC settings (accurate modeling of slow start with correction factors on three intervals: < 1KiB, < 64 KiB, >= 64 KiB)",
 117    surf_network_model_init_SMPI},
 118   {"IB",
 119    "Realistic network model specifically tailored for HPC settings, with Infiniband contention model",
 120    surf_network_model_init_IB},
 121   {"CM02",
 122    "Legacy network analytic model (Very similar to LV08, but without corrective factors. The timings of small messages are thus poorly modeled).",
 123    surf_network_model_init_CM02},
 124 #ifdef HAVE_GTNETS
 125   {"GTNets",
 126    "Network pseudo-model using the GTNets simulator instead of an analytic model",
 127    surf_network_model_init_GTNETS},
 128 #endif
 129 #ifdef HAVE_NS3
 130   {"NS3",
 131    "Network pseudo-model using the NS3 tcp model instead of an analytic model",
 132   surf_network_model_init_NS3},
 133 #endif
 134   {"Reno",
 135    "Model from Steven H. Low using lagrange_solve instead of lmm_solve (experts only; check the code for more info).",
 136    surf_network_model_init_Reno},
 137   {"Reno2",
 138    "Model from Steven H. Low using lagrange_solve instead of lmm_solve (experts only; check the code for more info).",
 139    surf_network_model_init_Reno2},
 140   {"Vegas",
 141    "Model from Steven H. Low using lagrange_solve instead of lmm_solve (experts only; check the code for more info).",
 142    surf_network_model_init_Vegas},
 143   {NULL, NULL, NULL}      /* this array must be NULL terminated */
 144 };
 145
 146 s_surf_model_description_t surf_cpu_model_description[] = {
 147   {"Cas01",
 148    "Simplistic CPU model (time=size/power).",
 149    surf_cpu_model_init_Cas01},
 150   {NULL, NULL,  NULL}      /* this array must be NULL terminated */
 151 };
 152
 153 s_surf_model_description_t surf_workstation_model_description[] = {
 154   {"default",
 155    "Default workstation model. Currently, CPU:Cas01 and network:LV08 (with cross traffic enabled)",
 156    surf_workstation_model_init_current_default},
 157   {"compound",
 158    "Workstation model that is automatically chosen if you change the network and CPU models",
 159    surf_workstation_model_init_compound},
 160   {"ptask_L07", "Workstation model somehow similar to Cas01+CM02 but allowing parallel tasks",
 161    surf_workstation_model_init_ptask_L07},
 162   {NULL, NULL, NULL}      /* this array must be NULL terminated */
 163 };
 164
 165 s_surf_model_description_t surf_vm_workstation_model_description[] = {
 166   {"default",
 167    "Default vm workstation model.",
 168    surf_vm_workstation_model_init_HL13},
 169   {NULL, NULL, NULL}      /* this array must be NULL terminated */
 170 };
 171
 172 s_surf_model_description_t surf_optimization_mode_description[] = {
 173   {"Lazy",
 174    "Lazy action management (partial invalidation in lmm + heap in action remaining).",
 175    NULL},
 176   {"TI",
 177    "Trace integration. Highly optimized mode when using availability traces (only available for the Cas01 CPU model for now).",
 178     NULL},
 179   {"Full",
 180    "Full update of remaining and variables. Slow but may be useful when debugging.",
 181    NULL},
 182   {NULL, NULL, NULL}      /* this array must be NULL terminated */
 183 };
 184
 185 s_surf_model_description_t surf_storage_model_description[] = {
 186   {"default",
 187    "Simplistic storage model.",
 188    surf_storage_model_init_default},
 189   {NULL, NULL,  NULL}      /* this array must be NULL terminated */
 190 };
 191
 192 #ifdef CONTEXT_THREADS
 193 static xbt_parmap_t surf_parmap = NULL; /* parallel map on models */
 194 #endif
 195
 196 double NOW = 0;
 197 double *surf_mins = NULL; /* return value of share_resources for each model */
 198 int surf_min_index;       /* current index in surf_mins */
 199 double surf_min;               /* duration determined by surf_solve */
 200
 201 double surf_get_clock(void)
 202 {
 203   return NOW;
 204 }
 205
 206 #ifdef _XBT_WIN32
 207 # define FILE_DELIM "\\"
 208 #else
 209 # define FILE_DELIM "/"         /* FIXME: move to better location */
 210 #endif
 211
 212 FILE *surf_fopen(const char *name, const char *mode)
 213 {
 214   unsigned int cpt;
 215   char *path_elm = NULL;
 216   char *buff;
 217   FILE *file = NULL;
 218
 219   xbt_assert(name);
 220
 221   if (__surf_is_absolute_file_path(name))       /* don't mess with absolute file names */
 222     return fopen(name, mode);
 223
 224   /* search relative files in the path */
 225   xbt_dynar_foreach(surf_path, cpt, path_elm) {
 226     buff = bprintf("%s" FILE_DELIM "%s", path_elm, name);
 227     file = fopen(buff, mode);
 228     free(buff);
 229
 230     if (file)
 231       return file;
 232   }
 233   return NULL;
 234 }
 235
 236
 237 #ifndef MAX_DRIVE
 238 #define MAX_DRIVE 26
 239 #endif
 240
 241 #ifdef _XBT_WIN32
 242 #include <windows.h>
 243 static const char *disk_drives_letter_table[MAX_DRIVE] = {
 244   "A:\\",
 245   "B:\\",
 246   "C:\\",
 247   "D:\\",
 248   "E:\\",
 249   "F:\\",
 250   "G:\\",
 251   "H:\\",
 252   "I:\\",
 253   "J:\\",
 254   "K:\\",
 255   "L:\\",
 256   "M:\\",
 257   "N:\\",
 258   "O:\\",
 259   "P:\\",
 260   "Q:\\",
 261   "R:\\",
 262   "S:\\",
 263   "T:\\",
 264   "U:\\",
 265   "V:\\",
 266   "W:\\",
 267   "X:\\",
 268   "Y:\\",
 269   "Z:\\"
 270 };
 271 #endif
 272
 273 /*
 274  * Returns the initial path. On Windows the initial path is
 275  * the current directory for the current process in the other
 276  * case the function returns "./" that represents the current
 277  * directory on Unix/Linux platforms.
 278  */
 279
 280 const char *__surf_get_initial_path(void)
 281 {
 282
 283 #ifdef _XBT_WIN32
 284   unsigned i;
 285   char current_directory[MAX_PATH + 1] = { 0 };
 286   unsigned int len = GetCurrentDirectory(MAX_PATH + 1, current_directory);
 287   char root[4] = { 0 };
 288
 289   if (!len)
 290     return NULL;
 291
 292   strncpy(root, current_directory, 3);
 293
 294   for (i = 0; i < MAX_DRIVE; i++) {
 295     if (toupper(root[0]) == disk_drives_letter_table[i][0])
 296       return disk_drives_letter_table[i];
 297   }
 298
 299   return NULL;
 300 #else
 301   return "./";
 302 #endif
 303 }
 304
 305 /* The __surf_is_absolute_file_path() returns 1 if
 306  * file_path is a absolute file path, in the other
 307  * case the function returns 0.
 308  */
 309 int __surf_is_absolute_file_path(const char *file_path)
 310 {
 311 #ifdef _XBT_WIN32
 312   WIN32_FIND_DATA wfd = { 0 };
 313   HANDLE hFile = FindFirstFile(file_path, &wfd);
 314
 315   if (INVALID_HANDLE_VALUE == hFile)
 316     return 0;
 317
 318   FindClose(hFile);
 319   return 1;
 320 #else
 321   return (file_path[0] == '/');
 322 #endif
 323 }
 324
 325 /** Displays the long description of all registered models, and quit */
 326 void model_help(const char *category, s_surf_model_description_t * table)
 327 {
 328   int i;
 329   printf("Long description of the %s models accepted by this simulator:\n",
 330          category);
 331   for (i = 0; table[i].name; i++)
 332     printf("  %s: %s\n", table[i].name, table[i].description);
 333 }
 334
 335 int find_model_description(s_surf_model_description_t * table,
 336                            const char *name)
 337 {
 338   int i;
 339   char *name_list = NULL;
 340
 341   for (i = 0; table[i].name; i++)
 342     if (!strcmp(name, table[i].name)) {
 343       return i;
 344     }
 345   if (!table[0].name)
 346     xbt_die("No model is valid! This is a bug.");
 347   name_list = xbt_strdup(table[0].name);
 348   for (i = 1; table[i].name; i++) {
 349     name_list = (char *) xbt_realloc(name_list, strlen(name_list) + strlen(table[i].name) + 3);
 350     strcat(name_list, ", ");
 351     strcat(name_list, table[i].name);
 352   }
 353   xbt_die("Model '%s' is invalid! Valid models are: %s.", name, name_list);
 354   return -1;
 355 }
 356
 357 static XBT_INLINE void routing_asr_host_free(void *p)
 358 {
 359   delete ((RoutingEdgePtr) p);
 360 }
 361
 362 static XBT_INLINE void routing_asr_prop_free(void *p)
 363 {
 364   xbt_dict_t elm = (xbt_dict_t) p;
 365   xbt_dict_free(&elm);
 366 }
 367
 368 static XBT_INLINE void surf_cpu_free(void *r)
 369 {
 370   delete static_cast<CpuPtr>(r);
 371 }
 372
 373 static XBT_INLINE void surf_link_free(void *r)
 374 {
 375   delete static_cast<NetworkLinkPtr>(r);
 376 }
 377
 378 static XBT_INLINE void surf_workstation_free(void *r)
 379 {
 380   delete static_cast<WorkstationPtr>(r);
 381 }
 382
 383 static XBT_INLINE void surf_storage_free(void *r)
 384 {
 385   delete static_cast<StoragePtr>(r);
 386 }
 387
 388
 389 void sg_version(int *ver_major,int *ver_minor,int *ver_patch) {
 390   *ver_major = SIMGRID_VERSION_MAJOR;
 391   *ver_minor = SIMGRID_VERSION_MINOR;
 392   *ver_patch = SIMGRID_VERSION_PATCH;
 393 }
 394
 395 void surf_init(int *argc, char **argv)
 396 {
 397   XBT_DEBUG("Create all Libs");
 398   host_lib = xbt_lib_new();
 399   link_lib = xbt_lib_new();
 400   as_router_lib = xbt_lib_new();
 401   storage_lib = xbt_lib_new();
 402   storage_type_lib = xbt_lib_new();
 403   file_lib = xbt_lib_new();
 404   watched_hosts_lib = xbt_dict_new_homogeneous(NULL);
 405
 406   XBT_DEBUG("Add routing levels");
 407   ROUTING_HOST_LEVEL = xbt_lib_add_level(host_lib,routing_asr_host_free);
 408   ROUTING_ASR_LEVEL  = xbt_lib_add_level(as_router_lib,routing_asr_host_free);
 409   ROUTING_PROP_ASR_LEVEL = xbt_lib_add_level(as_router_lib,routing_asr_prop_free);
 410
 411   XBT_DEBUG("Add SURF levels");
 412   SURF_CPU_LEVEL = xbt_lib_add_level(host_lib,surf_cpu_free);
 413   SURF_WKS_LEVEL = xbt_lib_add_level(host_lib,surf_workstation_free);
 414   SURF_LINK_LEVEL = xbt_lib_add_level(link_lib,surf_link_free);
 415   SURF_STORAGE_LEVEL = xbt_lib_add_level(storage_lib,surf_storage_free);
 416
 417   xbt_init(argc, argv);
 418   if (!model_list)
 419     model_list = xbt_dynar_new(sizeof(ModelPtr), NULL);
 420   if (!model_list_invoke)
 421     model_list_invoke = xbt_dynar_new(sizeof(ModelPtr), NULL);
 422   if (!history)
 423     history = tmgr_history_new();
 424
 425   TRACE_add_start_function(TRACE_surf_alloc);
 426   TRACE_add_end_function(TRACE_surf_release);
 427
 428   sg_config_init(argc, argv);
 429
 430   if (MC_is_active())
 431     MC_memory_init();
 432 }
 433
 434 void surf_exit(void)
 435 {
 436   unsigned int iter;
 437   ModelPtr model = NULL;
 438
 439   TRACE_end();                  /* Just in case it was not called by the upper
 440                                  * layer (or there is no upper layer) */
 441
 442   sg_config_finalize();
 443
 444   xbt_dynar_free(&host_that_restart);
 445   xbt_dynar_free(&surf_path);
 446
 447   xbt_lib_free(&host_lib);
 448   xbt_lib_free(&link_lib);
 449   xbt_lib_free(&as_router_lib);
 450   xbt_lib_free(&storage_lib);
 451   xbt_lib_free(&storage_type_lib);
 452   xbt_lib_free(&file_lib);
 453   xbt_dict_free(&watched_hosts_lib);
 454
 455   xbt_dynar_foreach(model_list, iter, model)
 456     delete model;
 457   xbt_dynar_free(&model_list);
 458   xbt_dynar_free(&model_list_invoke);
 459   routing_exit();
 460
 461   surf_callback_emit(surfExitCallbacks);
 462
 463   if (maxmin_system) {
 464     lmm_system_free(maxmin_system);
 465     maxmin_system = NULL;
 466   }
 467   if (history) {
 468     tmgr_history_free(history);
 469     history = NULL;
 470   }
 471
 472 #ifdef CONTEXT_THREADS
 473   xbt_parmap_destroy(surf_parmap);
 474 #endif
 475
 476   xbt_free(surf_mins);
 477   surf_mins = NULL;
 478
 479   tmgr_finalize();
 480   surf_parse_lex_destroy();
 481   surf_parse_free_callbacks();
 482
 483   NOW = 0;                      /* Just in case the user plans to restart the simulation afterward */
 484 }
 485
 486 /*********
 487  * Model *
 488  *********/
 489
 490 Model::Model(const char *name)
 491   : p_maxminSystem(NULL)
 492   , p_name(name)
 493 {
 494   p_readyActionSet = new ActionList();
 495   p_runningActionSet = new ActionList();
 496   p_failedActionSet = new ActionList();
 497   p_doneActionSet = new ActionList();
 498
 499   p_modifiedSet = NULL;
 500   p_actionHeap = NULL;
 501   p_updateMechanism = UM_UNDEFINED;
 502   m_selectiveUpdate = 0;
 503 }
 504
 505 Model::~Model(){
 506   delete p_readyActionSet;
 507   delete p_runningActionSet;
 508   delete p_failedActionSet;
 509   delete p_doneActionSet;
 510 }
 511
 512 double Model::shareResources(double now)
 513 {
 514   //FIXME: set the good function once and for all
 515   if (p_updateMechanism == UM_LAZY)
 516     return shareResourcesLazy(now);
 517   else if (p_updateMechanism == UM_FULL)
 518     return shareResourcesFull(now);
 519   else
 520     xbt_die("Invalid cpu update mechanism!");
 521 }
 522
 523 double Model::shareResourcesLazy(double now)
 524 {
 525   ActionPtr action = NULL;
 526   double min = -1;
 527   double share;
 528
 529   XBT_DEBUG
 530       ("Before share resources, the size of modified actions set is %zd",
 531        p_modifiedSet->size());
 532
 533   lmm_solve(p_maxminSystem);
 534
 535   XBT_DEBUG
 536       ("After share resources, The size of modified actions set is %zd",
 537        p_modifiedSet->size());
 538
 539   while(!p_modifiedSet->empty()) {
 540     action = &(p_modifiedSet->front());
 541     p_modifiedSet->pop_front();
 542     int max_dur_flag = 0;
 543
 544     if (action->getStateSet() != p_runningActionSet)
 545       continue;
 546
 547     /* bogus priority, skip it */
 548     if (action->getPriority() <= 0 || action->getHat()==LATENCY)
 549       continue;
 550
 551     action->updateRemainingLazy(now);
 552
 553     min = -1;
 554     share = lmm_variable_getvalue(action->getVariable());
 555
 556     if (share > 0) {
 557       double time_to_completion;
 558       if (action->getRemains() > 0) {
 559         time_to_completion = action->getRemainsNoUpdate() / share;
 560       } else {
 561         time_to_completion = 0.0;
 562       }
 563       min = now + time_to_completion; // when the task will complete if nothing changes
 564     }
 565
 566     if ((action->getMaxDuration() != NO_MAX_DURATION)
 567         && (min == -1
 568             || action->getStartTime() +
 569             action->getMaxDuration() < min)) {
 570       min = action->getStartTime() +
 571           action->getMaxDuration();  // when the task will complete anyway because of the deadline if any
 572       max_dur_flag = 1;
 573     }
 574
 575
 576     XBT_DEBUG("Action(%p) corresponds to variable %d", action, action->getVariable()->id_int);
 577
 578     XBT_DEBUG("Action(%p) Start %f. May finish at %f (got a share of %f). Max_duration %f", action,
 579         action->getStartTime(), min, share,
 580         action->getMaxDuration());
 581
 582     if (min != -1) {
 583       action->heapUpdate(p_actionHeap, min, max_dur_flag ? MAX_DURATION : NORMAL);
 584       XBT_DEBUG("Insert at heap action(%p) min %f now %f", action, min,
 585                 now);
 586     } else DIE_IMPOSSIBLE;
 587   }
 588
 589   //hereafter must have already the min value for this resource model
 590   if (xbt_heap_size(p_actionHeap) > 0)
 591     min = xbt_heap_maxkey(p_actionHeap) - now;
 592   else
 593     min = -1;
 594
 595   XBT_DEBUG("The minimum with the HEAP %f", min);
 596
 597   return min;
 598 }
 599
 600 double Model::shareResourcesFull(double /*now*/) {
 601   THROW_UNIMPLEMENTED;
 602 }
 603
 604 double Model::shareResourcesMaxMin(ActionListPtr running_actions,
 605                           lmm_system_t sys,
 606                           void (*solve) (lmm_system_t))
 607 {
 608   ActionPtr action = NULL;
 609   double min = -1;
 610   double value = -1;
 611
 612   solve(sys);
 613
 614   ActionList::iterator it(running_actions->begin()), itend(running_actions->end());
 615   for(; it != itend ; ++it) {
 616     action = &*it;
 617     value = lmm_variable_getvalue(action->getVariable());
 618     if ((value > 0) || (action->getMaxDuration() >= 0))
 619       break;
 620   }
 621
 622   if (!action)
 623     return -1.0;
 624
 625   if (value > 0) {
 626     if (action->getRemains() > 0)
 627       min = action->getRemainsNoUpdate() / value;
 628     else
 629       min = 0.0;
 630     if ((action->getMaxDuration() >= 0) && (action->getMaxDuration() < min))
 631       min = action->getMaxDuration();
 632   } else
 633     min = action->getMaxDuration();
 634
 635
 636   for (++it; it != itend; ++it) {
 637         action = &*it;
 638     value = lmm_variable_getvalue(action->getVariable());
 639     if (value > 0) {
 640       if (action->getRemains() > 0)
 641         value = action->getRemainsNoUpdate() / value;
 642       else
 643         value = 0.0;
 644       if (value < min) {
 645         min = value;
 646         XBT_DEBUG("Updating min (value) with %p: %f", action, min);
 647       }
 648     }
 649     if ((action->getMaxDuration() >= 0) && (action->getMaxDuration() < min)) {
 650       min = action->getMaxDuration();
 651       XBT_DEBUG("Updating min (duration) with %p: %f", action, min);
 652     }
 653   }
 654   XBT_DEBUG("min value : %f", min);
 655
 656   return min;
 657 }
 658
 659 void Model::updateActionsState(double now, double delta)
 660 {
 661   if (p_updateMechanism == UM_FULL)
 662         updateActionsStateFull(now, delta);
 663   else if (p_updateMechanism == UM_LAZY)
 664         updateActionsStateLazy(now, delta);
 665   else
 666         xbt_die("Invalid cpu update mechanism!");
 667 }
 668
 669 void Model::updateActionsStateLazy(double /*now*/, double /*delta*/)
 670 {
 671  THROW_UNIMPLEMENTED;
 672 }
 673
 674 void Model::updateActionsStateFull(double /*now*/, double /*delta*/)
 675 {
 676   THROW_UNIMPLEMENTED;
 677 }
 678
 679 /************
 680  * Resource *
 681  ************/
 682
 683 Resource::Resource()
 684 : p_name(NULL), p_properties(NULL), p_model(NULL)
 685 {}
 686
 687 Resource::Resource(ModelPtr model, const char *name, xbt_dict_t props)
 688   : p_name(xbt_strdup(name)), p_properties(props), p_model(model)
 689   , m_running(true), m_stateCurrent(SURF_RESOURCE_ON)
 690 {}
 691
 692 Resource::Resource(ModelPtr model, const char *name, xbt_dict_t props, lmm_constraint_t constraint)
 693   : p_name(xbt_strdup(name)), p_properties(props), p_model(model)
 694   , m_running(true), m_stateCurrent(SURF_RESOURCE_ON), p_constraint(constraint)
 695 {}
 696
 697 Resource::Resource(ModelPtr model, const char *name, xbt_dict_t props, e_surf_resource_state_t stateInit)
 698   : p_name(xbt_strdup(name)), p_properties(props), p_model(model)
 699   , m_running(true), m_stateCurrent(stateInit)
 700 {}
 701
 702 Resource::~Resource() {
 703   xbt_free((void*)p_name);
 704   xbt_dict_free(&p_properties);
 705 }
 706
 707 e_surf_resource_state_t Resource::getState()
 708 {
 709   return m_stateCurrent;
 710 }
 711
 712 void Resource::setState(e_surf_resource_state_t state)
 713 {
 714   m_stateCurrent = state;
 715 }
 716
 717 bool Resource::isOn()
 718 {
 719   return m_running;
 720 }
 721
 722 void Resource::turnOn()
 723 {
 724   if (!m_running) {
 725     m_running = true;
 726   }
 727 }
 728
 729 void Resource::turnOff()
 730 {
 731   if (m_running) {
 732     m_running = false;
 733   }
 734 }
 735
 736 ModelPtr Resource::getModel() {
 737   return p_model;
 738 }
 739
 740 const char *Resource::getName() {
 741   return p_name;
 742 }
 743
 744 xbt_dict_t Resource::getProperties() {
 745   if (p_properties==NULL)
 746     p_properties = xbt_dict_new();
 747   return p_properties;
 748 }
 749
 750 lmm_constraint_t Resource::getConstraint() {
 751   return p_constraint;
 752 }
 753
 754 /**********
 755  * Action *
 756  **********/
 757
 758 const char *surf_action_state_names[6] = {
 759   "SURF_ACTION_READY",
 760   "SURF_ACTION_RUNNING",
 761   "SURF_ACTION_FAILED",
 762   "SURF_ACTION_DONE",
 763   "SURF_ACTION_TO_FREE",
 764   "SURF_ACTION_NOT_IN_THE_SYSTEM"
 765 };
 766
 767 void Action::initialize(ModelPtr model, double cost, bool failed,
 768                         lmm_variable_t var)
 769 {
 770   m_priority = 1.0;
 771   m_refcount = 1;
 772   m_remains = cost;
 773   m_maxDuration = NO_MAX_DURATION;
 774   m_finish = -1.0;
 775   m_failed = failed;
 776   m_start = surf_get_clock();
 777   m_cost = cost;
 778   p_model = model;
 779   p_data = NULL;
 780   p_variable = var;
 781   m_lastValue = 0;
 782   m_lastUpdate = 0;
 783   m_suspended = false;
 784   m_hat = NOTSET;
 785 }
 786
 787 Action::Action(ModelPtr model, double cost, bool failed)
 788 {
 789   initialize(model, cost, failed);
 790   p_category = NULL;
 791   p_stateHookup.prev = 0;
 792   p_stateHookup.next = 0;
 793   if (failed)
 794     p_stateSet = getModel()->getFailedActionSet();
 795   else
 796     p_stateSet = getModel()->getRunningActionSet();
 797
 798   p_stateSet->push_back(*this);
 799 }
 800
 801 Action::Action(ModelPtr model, double cost, bool failed, lmm_variable_t var)
 802 {
 803   initialize(model, cost, failed, var);
 804   p_category = NULL;
 805   p_stateHookup.prev = 0;
 806   p_stateHookup.next = 0;
 807   if (failed)
 808     p_stateSet = getModel()->getFailedActionSet();
 809   else
 810     p_stateSet = getModel()->getRunningActionSet();
 811
 812   p_stateSet->push_back(*this);
 813 }
 814
 815 Action::~Action() {
 816   xbt_free(p_category);
 817 }
 818
 819 void Action::finish() {
 820     m_finish = surf_get_clock();
 821 }
 822
 823 e_surf_action_state_t Action::getState()
 824 {
 825   if (p_stateSet ==  getModel()->getReadyActionSet())
 826     return SURF_ACTION_READY;
 827   if (p_stateSet ==  getModel()->getRunningActionSet())
 828     return SURF_ACTION_RUNNING;
 829   if (p_stateSet ==  getModel()->getFailedActionSet())
 830     return SURF_ACTION_FAILED;
 831   if (p_stateSet ==  getModel()->getDoneActionSet())
 832     return SURF_ACTION_DONE;
 833   return SURF_ACTION_NOT_IN_THE_SYSTEM;
 834 }
 835
 836 void Action::setState(e_surf_action_state_t state)
 837 {
 838   //surf_action_state_t action_state = &(action->model_type->states);
 839   XBT_IN("(%p,%s)", this, surf_action_state_names[state]);
 840   p_stateSet->erase(p_stateSet->iterator_to(*this));
 841   if (state == SURF_ACTION_READY)
 842     p_stateSet = getModel()->getReadyActionSet();
 843   else if (state == SURF_ACTION_RUNNING)
 844     p_stateSet = getModel()->getRunningActionSet();
 845   else if (state == SURF_ACTION_FAILED)
 846     p_stateSet = getModel()->getFailedActionSet();
 847   else if (state == SURF_ACTION_DONE)
 848     p_stateSet = getModel()->getDoneActionSet();
 849   else
 850     p_stateSet = NULL;
 851
 852   if (p_stateSet)
 853     p_stateSet->push_back(*this);
 854   XBT_OUT();
 855 }
 856
 857 double Action::getBound()
 858 {
 859   return (p_variable) ? lmm_variable_getbound(p_variable) : 0;
 860 }
 861
 862 void Action::setBound(double bound)
 863 {
 864   XBT_IN("(%p,%g)", this, bound);
 865   if (p_variable)
 866     lmm_update_variable_bound(getModel()->getMaxminSystem(), p_variable, bound);
 867
 868   if (getModel()->getUpdateMechanism() == UM_LAZY && getLastUpdate()!=surf_get_clock())
 869     heapRemove(getModel()->getActionHeap());
 870   XBT_OUT();
 871 }
 872
 873 double Action::getStartTime()
 874 {
 875   return m_start;
 876 }
 877
 878 double Action::getFinishTime()
 879 {
 880   /* keep the function behavior, some models (cpu_ti) change the finish time before the action end */
 881   return m_remains == 0 ? m_finish : -1;
 882 }
 883
 884 void Action::setData(void* data)
 885 {
 886   p_data = data;
 887 }
 888
 889 void Action::setCategory(const char *category)
 890 {
 891   XBT_IN("(%p,%s)", this, category);
 892   p_category = xbt_strdup(category);
 893   XBT_OUT();
 894 }
 895
 896 void Action::ref(){
 897   m_refcount++;
 898 }
 899
 900 void Action::setMaxDuration(double duration)
 901 {
 902   XBT_IN("(%p,%g)", this, duration);
 903   m_maxDuration = duration;
 904   if (getModel()->getUpdateMechanism() == UM_LAZY)      // remove action from the heap
 905     heapRemove(getModel()->getActionHeap());
 906   XBT_OUT();
 907 }
 908
 909 void Action::gapRemove() {}
 910
 911 void Action::setPriority(double priority)
 912 {
 913   XBT_IN("(%p,%g)", this, priority);
 914   m_priority = priority;
 915   lmm_update_variable_weight(getModel()->getMaxminSystem(), getVariable(), priority);
 916
 917   if (getModel()->getUpdateMechanism() == UM_LAZY)
 918         heapRemove(getModel()->getActionHeap());
 919   XBT_OUT();
 920 }
 921
 922 void Action::cancel(){
 923   setState(SURF_ACTION_FAILED);
 924   if (getModel()->getUpdateMechanism() == UM_LAZY) {
 925     if (actionLmmHook::is_linked())
 926       getModel()->getModifiedSet()->erase(getModel()->getModifiedSet()->iterator_to(*this));
 927     heapRemove(getModel()->getActionHeap());
 928   }
 929 }
 930
 931 int Action::unref(){
 932   m_refcount--;
 933   if (!m_refcount) {
 934     if (actionHook::is_linked())
 935       p_stateSet->erase(p_stateSet->iterator_to(*this));
 936     if (getVariable())
 937       lmm_variable_free(getModel()->getMaxminSystem(), getVariable());
 938     if (getModel()->getUpdateMechanism() == UM_LAZY) {
 939       /* remove from heap */
 940       heapRemove(getModel()->getActionHeap());
 941       if (actionLmmHook::is_linked())
 942         getModel()->getModifiedSet()->erase(getModel()->getModifiedSet()->iterator_to(*this));
 943     }
 944     delete this;
 945     return 1;
 946   }
 947   return 0;
 948 }
 949
 950 void Action::suspend()
 951 {
 952   XBT_IN("(%p)", this);
 953   if (m_suspended != 2) {
 954     lmm_update_variable_weight(getModel()->getMaxminSystem(), getVariable(), 0.0);
 955     m_suspended = 1;
 956     if (getModel()->getUpdateMechanism() == UM_LAZY)
 957       heapRemove(getModel()->getActionHeap());
 958   }
 959   XBT_OUT();
 960 }
 961
 962 void Action::resume()
 963 {
 964   XBT_IN("(%p)", this);
 965   if (m_suspended != 2) {
 966     lmm_update_variable_weight(getModel()->getMaxminSystem(), getVariable(), m_priority);
 967     m_suspended = 0;
 968     if (getModel()->getUpdateMechanism() == UM_LAZY)
 969       heapRemove(getModel()->getActionHeap());
 970   }
 971   XBT_OUT();
 972 }
 973
 974 bool Action::isSuspended()
 975 {
 976   return m_suspended == 1;
 977 }
 978 /* insert action on heap using a given key and a hat (heap_action_type)
 979  * a hat can be of three types for communications:
 980  *
 981  * NORMAL = this is a normal heap entry stating the date to finish transmitting
 982  * LATENCY = this is a heap entry to warn us when the latency is payed
 983  * MAX_DURATION =this is a heap entry to warn us when the max_duration limit is reached
 984  */
 985 void Action::heapInsert(xbt_heap_t heap, double key, enum heap_action_type hat)
 986 {
 987   m_hat = hat;
 988   xbt_heap_push(heap, this, key);
 989 }
 990
 991 void Action::heapRemove(xbt_heap_t heap)
 992 {
 993   m_hat = NOTSET;
 994   if (m_indexHeap >= 0) {
 995     xbt_heap_remove(heap, m_indexHeap);
 996   }
 997 }
 998
 999 void Action::heapUpdate(xbt_heap_t heap, double key, enum heap_action_type hat)
1000 {
1001   m_hat = hat;
1002   if (m_indexHeap >= 0) {
1003     xbt_heap_update(heap, m_indexHeap, key);
1004   }else{
1005     xbt_heap_push(heap, this, key);
1006   }
1007 }
1008
1009 /* added to manage the communication action's heap */
1010 void surf_action_lmm_update_index_heap(void *action, int i) {
1011   ((ActionPtr)action)->updateIndexHeap(i);
1012 }
1013
1014 void Action::updateIndexHeap(int i) {
1015   m_indexHeap = i;
1016 }
1017
1018 double Action::getRemains()
1019 {
1020   XBT_IN("(%p)", this);
1021   /* update remains before return it */
1022   if (getModel()->getUpdateMechanism() == UM_LAZY)      /* update remains before return it */
1023     updateRemainingLazy(surf_get_clock());
1024   XBT_OUT();
1025   return m_remains;
1026 }
1027
1028 double Action::getRemainsNoUpdate()
1029 {
1030   return m_remains;
1031 }
1032
1033 //FIXME split code in the right places
1034 void Action::updateRemainingLazy(double now)
1035 {
1036   double delta = 0.0;
1037
1038   if(getModel() == surf_network_model)
1039   {
1040     if (m_suspended != 0)
1041       return;
1042   }
1043   else
1044   {
1045     xbt_assert(p_stateSet == getModel()->getRunningActionSet(),
1046         "You're updating an action that is not running.");
1047
1048       /* bogus priority, skip it */
1049     xbt_assert(m_priority > 0,
1050         "You're updating an action that seems suspended.");
1051   }
1052
1053   delta = now - m_lastUpdate;
1054
1055   if (m_remains > 0) {
1056     XBT_DEBUG("Updating action(%p): remains was %f, last_update was: %f", this, m_remains, m_lastUpdate);
1057     double_update(&m_remains, m_lastValue * delta, sg_surf_precision*sg_maxmin_precision);
1058
1059     if (getModel() == surf_cpu_model_pm && TRACE_is_enabled()) {
1060       ResourcePtr cpu = static_cast<ResourcePtr>(lmm_constraint_id(lmm_get_cnst_from_var(getModel()->getMaxminSystem(), getVariable(), 0)));
1061       TRACE_surf_host_set_utilization(cpu->getName(), getCategory(), m_lastValue, m_lastUpdate, now - m_lastUpdate);
1062     }
1063     XBT_DEBUG("Updating action(%p): remains is now %f", this, m_remains);
1064   }
1065
1066   if(getModel() == surf_network_model)
1067   {
1068     if (m_maxDuration != NO_MAX_DURATION)
1069       double_update(&m_maxDuration, delta, sg_surf_precision);
1070
1071     //FIXME: duplicated code
1072     if ((m_remains <= 0) &&
1073         (lmm_get_variable_weight(getVariable()) > 0)) {
1074       finish();
1075       setState(SURF_ACTION_DONE);
1076       heapRemove(getModel()->getActionHeap());
1077     } else if (((m_maxDuration != NO_MAX_DURATION)
1078         && (m_maxDuration <= 0))) {
1079       finish();
1080       setState(SURF_ACTION_DONE);
1081       heapRemove(getModel()->getActionHeap());
1082     }
1083   }
1084
1085   m_lastUpdate = now;
1086   m_lastValue = lmm_variable_getvalue(getVariable());
1087 }