src/include/surf/maxmin.h

   1 /*      $Id$     */
   2
   3 /* Copyright (c) 2004 Arnaud Legrand. 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 #ifndef _SURF_MAXMIN_H
   9 #define _SURF_MAXMIN_H
  10
  11 #include <math.h>
  12 #include "portable.h"
  13 #include "xbt/misc.h"
  14 #include "surf/datatypes.h"
  15
  16 #define MAXMIN_PRECISION 0.00001
  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 XBT_PUBLIC(lmm_system_t) lmm_system_new(void);
  35 XBT_PUBLIC(void) lmm_system_free(lmm_system_t sys);
  36      void lmm_variable_disable(lmm_system_t sys, lmm_variable_t var);
  37
  38 XBT_PUBLIC(lmm_constraint_t) lmm_constraint_new(lmm_system_t sys, void *id,
  39                                                 double bound_value);
  40      void lmm_constraint_shared(lmm_constraint_t cnst);
  41      int lmm_constraint_is_shared(lmm_constraint_t cnst);
  42
  43      void lmm_constraint_free(lmm_system_t sys, lmm_constraint_t cnst);
  44
  45 XBT_PUBLIC(lmm_variable_t) lmm_variable_new(lmm_system_t sys, void *id,
  46                                             double weight_value,
  47                                             double bound,
  48                                             int number_of_constraints);
  49      void lmm_variable_free(lmm_system_t sys, lmm_variable_t var);
  50 XBT_PUBLIC(double) lmm_variable_getvalue(lmm_variable_t var);
  51 XBT_PUBLIC(double) lmm_variable_getbound(lmm_variable_t var);
  52
  53 XBT_PUBLIC(void) lmm_expand(lmm_system_t sys, lmm_constraint_t cnst,
  54                             lmm_variable_t var, double value);
  55      void lmm_expand_add(lmm_system_t sys, lmm_constraint_t cnst,
  56                          lmm_variable_t var, double value);
  57      void lmm_elem_set_value(lmm_system_t sys, lmm_constraint_t cnst,
  58                              lmm_variable_t var, double value);
  59
  60      lmm_constraint_t lmm_get_cnst_from_var(lmm_system_t sys,
  61                                             lmm_variable_t var, int num);
  62      int lmm_get_number_of_cnst_from_var(lmm_system_t sys,
  63                                          lmm_variable_t var);
  64      lmm_variable_t lmm_get_var_from_cnst(lmm_system_t sys,
  65                                           lmm_constraint_t cnst,
  66                                           lmm_element_t * elem);
  67
  68      lmm_constraint_t lmm_get_first_active_constraint(lmm_system_t sys);
  69      lmm_constraint_t lmm_get_next_active_constraint(lmm_system_t sys,
  70                                                      lmm_constraint_t cnst);
  71
  72      void *lmm_constraint_id(lmm_constraint_t cnst);
  73      void *lmm_variable_id(lmm_variable_t var);
  74
  75      void lmm_update(lmm_system_t sys, lmm_constraint_t cnst,
  76                      lmm_variable_t var, double value);
  77      void lmm_update_variable_bound(lmm_system_t sys, lmm_variable_t var,
  78                                     double bound);
  79
  80
  81 XBT_PUBLIC(void) lmm_update_variable_weight(lmm_system_t sys,
  82                                             lmm_variable_t var,
  83                                             double weight);
  84      double lmm_get_variable_weight(lmm_variable_t var);
  85
  86 XBT_PUBLIC(void) lmm_update_constraint_bound(lmm_system_t sys,
  87                                              lmm_constraint_t cnst,
  88                                              double bound);
  89
  90      int lmm_constraint_used(lmm_system_t sys, lmm_constraint_t cnst);
  91
  92
  93 XBT_PUBLIC(void) lmm_solve(lmm_system_t sys);
  94
  95 XBT_PUBLIC(void) lagrange_solve(lmm_system_t sys);
  96 XBT_PUBLIC(void) bottleneck_solve(lmm_system_t sys);
  97
  98 /**
  99  * Default functions associated to the chosen protocol. When
 100  * using the lagrangian approach.
 101  */
 102
 103 XBT_PUBLIC(void) lmm_set_default_protocol_function(double (*func_f)
 104                                                    (lmm_variable_t var, double x),
 105                                                    double (*func_fp)
 106                                                    (lmm_variable_t var, double x),
 107                                                    double (*func_fpi)
 108                                                    (lmm_variable_t var, double x));
 109
 110 XBT_PUBLIC(double func_reno_f) (lmm_variable_t var, double x);
 111 XBT_PUBLIC(double func_reno_fp) (lmm_variable_t var, double x);
 112 XBT_PUBLIC(double func_reno_fpi) (lmm_variable_t var, double x);
 113
 114 XBT_PUBLIC(double func_reno2_f) (lmm_variable_t var, double x);
 115 XBT_PUBLIC(double func_reno2_fp) (lmm_variable_t var, double x);
 116 XBT_PUBLIC(double func_reno2_fpi) (lmm_variable_t var, double x);
 117
 118 XBT_PUBLIC(double func_vegas_f) (lmm_variable_t var, double x);
 119 XBT_PUBLIC(double func_vegas_fp) (lmm_variable_t var, double x);
 120 XBT_PUBLIC(double func_vegas_fpi) (lmm_variable_t var, double x);
 121
 122
 123 #endif /* _SURF_MAXMIN_H */