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 SG_BEGIN_DECL()
  35
  36 XBT_PUBLIC(lmm_system_t) lmm_system_new(int selective_update);
  37 XBT_PUBLIC(void) lmm_system_free(lmm_system_t sys);
  38
  39 XBT_PUBLIC(lmm_constraint_t) lmm_constraint_new(lmm_system_t sys, void *id,
  40                                                 double bound_value);
  41 void lmm_constraint_shared(lmm_constraint_t cnst);
  42 int lmm_constraint_is_shared(lmm_constraint_t cnst);
  43
  44 void lmm_constraint_free(lmm_system_t sys, lmm_constraint_t cnst);
  45
  46 double lmm_constraint_get_usage(lmm_constraint_t cnst);
  47
  48 XBT_PUBLIC(lmm_variable_t) lmm_variable_new(lmm_system_t sys, void *id,
  49                                             double weight_value,
  50                                             double bound,
  51                                             int number_of_constraints);
  52 XBT_PUBLIC(void) lmm_variable_free(lmm_system_t sys, lmm_variable_t var);
  53 XBT_PUBLIC(double) lmm_variable_getvalue(lmm_variable_t var);
  54 XBT_PUBLIC(double) lmm_variable_getbound(lmm_variable_t var);
  55
  56 XBT_PUBLIC(void) lmm_shrink(lmm_system_t sys, lmm_constraint_t cnst,
  57                             lmm_variable_t var);
  58 XBT_PUBLIC(void) lmm_expand(lmm_system_t sys, lmm_constraint_t cnst,
  59                             lmm_variable_t var, double value);
  60 void lmm_expand_add(lmm_system_t sys, lmm_constraint_t cnst,
  61                     lmm_variable_t var, double value);
  62 void lmm_elem_set_value(lmm_system_t sys, lmm_constraint_t cnst,
  63                         lmm_variable_t var, double value);
  64
  65 lmm_constraint_t lmm_get_cnst_from_var(lmm_system_t sys,
  66                                        lmm_variable_t var, int num);
  67 double lmm_get_cnst_weight_from_var(lmm_system_t sys, lmm_variable_t var,
  68                                     int num);
  69 int lmm_get_number_of_cnst_from_var(lmm_system_t sys, lmm_variable_t var);
  70 lmm_variable_t lmm_get_var_from_cnst(lmm_system_t sys,
  71                                      lmm_constraint_t cnst,
  72                                      lmm_element_t * elem);
  73
  74 lmm_constraint_t lmm_get_first_active_constraint(lmm_system_t sys);
  75 lmm_constraint_t lmm_get_next_active_constraint(lmm_system_t sys,
  76                                                 lmm_constraint_t cnst);
  77 #ifdef HAVE_LATENCY_BOUND_TRACKING
  78 XBT_PUBLIC(int) lmm_is_variable_limited_by_latency(lmm_variable_t var);
  79 #endif
  80
  81 void *lmm_constraint_id(lmm_constraint_t cnst);
  82 void *lmm_variable_id(lmm_variable_t var);
  83
  84 void lmm_update(lmm_system_t sys, lmm_constraint_t cnst,
  85                 lmm_variable_t var, double value);
  86 void lmm_update_variable_bound(lmm_system_t sys, lmm_variable_t var,
  87                                double bound);
  88
  89
  90 XBT_PUBLIC(void) lmm_update_variable_weight(lmm_system_t sys,
  91                                             lmm_variable_t var,
  92                                             double weight);
  93 double lmm_get_variable_weight(lmm_variable_t var);
  94
  95 XBT_PUBLIC(void) lmm_update_constraint_bound(lmm_system_t sys,
  96                                              lmm_constraint_t cnst,
  97                                              double bound);
  98
  99 int lmm_constraint_used(lmm_system_t sys, lmm_constraint_t cnst);
 100
 101
 102 XBT_PUBLIC(void) lmm_solve(lmm_system_t sys);
 103
 104 XBT_PUBLIC(void) lagrange_solve(lmm_system_t sys);
 105 XBT_PUBLIC(void) bottleneck_solve(lmm_system_t sys);
 106
 107 /**
 108  * Default functions associated to the chosen protocol. When
 109  * using the lagrangian approach.
 110  */
 111
 112 XBT_PUBLIC(void) lmm_set_default_protocol_function(double (*func_f)
 113                                                     (lmm_variable_t var,
 114                                                      double x),
 115                                                    double (*func_fp)
 116                                                     (lmm_variable_t var,
 117                                                      double x),
 118                                                    double (*func_fpi)
 119                                                     (lmm_variable_t var,
 120                                                      double x));
 121
 122 XBT_PUBLIC(double func_reno_f) (lmm_variable_t var, double x);
 123 XBT_PUBLIC(double func_reno_fp) (lmm_variable_t var, double x);
 124 XBT_PUBLIC(double func_reno_fpi) (lmm_variable_t var, double x);
 125
 126 XBT_PUBLIC(double func_reno2_f) (lmm_variable_t var, double x);
 127 XBT_PUBLIC(double func_reno2_fp) (lmm_variable_t var, double x);
 128 XBT_PUBLIC(double func_reno2_fpi) (lmm_variable_t var, double x);
 129
 130 XBT_PUBLIC(double func_vegas_f) (lmm_variable_t var, double x);
 131 XBT_PUBLIC(double func_vegas_fp) (lmm_variable_t var, double x);
 132 XBT_PUBLIC(double func_vegas_fpi) (lmm_variable_t var, double x);
 133
 134 SG_END_DECL()
 135
 136 #endif                          /* _SURF_MAXMIN_H */