src/include/surf/maxmin.h

   1 /* Copyright (c) 2004-2013. 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 #ifndef _SURF_MAXMIN_H
   8 #define _SURF_MAXMIN_H
   9
  10 #include "portable.h"
  11 #include "xbt/misc.h"
  12 #include "surf/datatypes.h"
  13 #include <math.h>
  14
  15 extern double sg_maxmin_precision;
  16 #define MAXMIN_PRECISION sg_maxmin_precision
  17 static XBT_INLINE void double_update(double *variable, double value)
  18 {
  19   *variable -= value;
  20   if (*variable < MAXMIN_PRECISION)
  21     *variable = 0.0;
  22 }
  23
  24 static XBT_INLINE int double_positive(double value)
  25 {
  26   return (value > MAXMIN_PRECISION);
  27 }
  28
  29 static XBT_INLINE int double_equals(double value1, double value2)
  30 {
  31   return (fabs(value1 - value2) < MAXMIN_PRECISION);
  32 }
  33
  34 #ifdef __cplusplus
  35 extern "C" {
  36 #endif
  37
  38 XBT_PUBLIC(lmm_system_t) lmm_system_new(int selective_update);
  39 XBT_PUBLIC(void) lmm_system_free(lmm_system_t sys);
  40 void lmm_variable_disable(lmm_system_t sys, lmm_variable_t var);
  41
  42 XBT_PUBLIC(lmm_constraint_t) lmm_constraint_new(lmm_system_t sys, void *id,
  43                                                 double bound_value);
  44 void lmm_constraint_shared(lmm_constraint_t cnst);
  45 int lmm_constraint_is_shared(lmm_constraint_t cnst);
  46
  47 void lmm_constraint_free(lmm_system_t sys, lmm_constraint_t cnst);
  48
  49 double lmm_constraint_get_usage(lmm_constraint_t cnst);
  50
  51 XBT_PUBLIC(lmm_variable_t) lmm_variable_new(lmm_system_t sys, void *id,
  52                                             double weight_value,
  53                                             double bound,
  54                                             int number_of_constraints);
  55 XBT_PUBLIC(void) lmm_variable_free(lmm_system_t sys, lmm_variable_t var);
  56 XBT_PUBLIC(double) lmm_variable_getvalue(lmm_variable_t var);
  57 XBT_PUBLIC(double) lmm_variable_getbound(lmm_variable_t var);
  58
  59 XBT_PUBLIC(void) lmm_expand(lmm_system_t sys, lmm_constraint_t cnst,
  60                             lmm_variable_t var, double value);
  61 void lmm_expand_add(lmm_system_t sys, lmm_constraint_t cnst,
  62                     lmm_variable_t var, double value);
  63 void lmm_elem_set_value(lmm_system_t sys, lmm_constraint_t cnst,
  64                         lmm_variable_t var, double value);
  65
  66 lmm_constraint_t lmm_get_cnst_from_var(lmm_system_t sys,
  67                                        lmm_variable_t var, int num);
  68 double lmm_get_cnst_weight_from_var(lmm_system_t sys, lmm_variable_t var,
  69                                     int num);
  70 int lmm_get_number_of_cnst_from_var(lmm_system_t sys, lmm_variable_t var);
  71 lmm_variable_t lmm_get_var_from_cnst(lmm_system_t sys,
  72                                      lmm_constraint_t cnst,
  73                                      lmm_element_t * elem);
  74
  75 lmm_constraint_t lmm_get_first_active_constraint(lmm_system_t sys);
  76 lmm_constraint_t lmm_get_next_active_constraint(lmm_system_t sys,
  77                                                 lmm_constraint_t cnst);
  78 #ifdef HAVE_LATENCY_BOUND_TRACKING
  79 XBT_PUBLIC(int) lmm_is_variable_limited_by_latency(lmm_variable_t var);
  80 #endif
  81
  82 void *lmm_constraint_id(lmm_constraint_t cnst);
  83 void *lmm_variable_id(lmm_variable_t var);
  84
  85 void lmm_update(lmm_system_t sys, lmm_constraint_t cnst,
  86                 lmm_variable_t var, double value);
  87 void lmm_update_variable_bound(lmm_system_t sys, lmm_variable_t var,
  88                                double bound);
  89
  90
  91 XBT_PUBLIC(void) lmm_update_variable_weight(lmm_system_t sys,
  92                                             lmm_variable_t var,
  93                                             double weight);
  94 double lmm_get_variable_weight(lmm_variable_t var);
  95
  96 XBT_PUBLIC(void) lmm_update_constraint_bound(lmm_system_t sys,
  97                                              lmm_constraint_t cnst,
  98                                              double bound);
  99
 100 int lmm_constraint_used(lmm_system_t sys, lmm_constraint_t cnst);
 101
 102
 103 XBT_PUBLIC(void) lmm_solve(lmm_system_t sys);
 104
 105 XBT_PUBLIC(void) lagrange_solve(lmm_system_t sys);
 106 XBT_PUBLIC(void) bottleneck_solve(lmm_system_t sys);
 107
 108 /**
 109  * Default functions associated to the chosen protocol. When
 110  * using the lagrangian approach.
 111  */
 112
 113 XBT_PUBLIC(void) lmm_set_default_protocol_function(double (*func_f)
 114                                                     (lmm_variable_t var,
 115                                                      double x),
 116                                                    double (*func_fp)
 117                                                     (lmm_variable_t var,
 118                                                      double x),
 119                                                    double (*func_fpi)
 120                                                     (lmm_variable_t var,
 121                                                      double x));
 122
 123 XBT_PUBLIC(double func_reno_f) (lmm_variable_t var, double x);
 124 XBT_PUBLIC(double func_reno_fp) (lmm_variable_t var, double x);
 125 XBT_PUBLIC(double func_reno_fpi) (lmm_variable_t var, double x);
 126
 127 XBT_PUBLIC(double func_reno2_f) (lmm_variable_t var, double x);
 128 XBT_PUBLIC(double func_reno2_fp) (lmm_variable_t var, double x);
 129 XBT_PUBLIC(double func_reno2_fpi) (lmm_variable_t var, double x);
 130
 131 XBT_PUBLIC(double func_vegas_f) (lmm_variable_t var, double x);
 132 XBT_PUBLIC(double func_vegas_fp) (lmm_variable_t var, double x);
 133 XBT_PUBLIC(double func_vegas_fpi) (lmm_variable_t var, double x);
 134
 135 #ifdef __cplusplus
 136 }
 137 #endif
 138
 139 #endif                          /* _SURF_MAXMIN_H */