Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Efficient linear max-minimization library : first version.
[simgrid.git] / src / surf / maxmin.c
1 /* Authors: Arnaud Legrand                                                  */
2
3 /* This program is free software; you can redistribute it and/or modify it
4    under the terms of the license (GNU LGPL) which comes with this package. */
5
6 #include "xbt/sysdep.h"
7 #include "xbt/error.h"
8 #include "maxmin_private.h"
9 #include <stdlib.h>
10
11 lmm_system_t lmm_system_new(void)
12 {
13   lmm_system_t l = NULL;
14   s_lmm_variable_t var;
15   s_lmm_constraint_t cnst;
16
17   l = xbt_new0(s_lmm_system_t, 1);
18
19   xbt_swag_init(&(l->variable_set),
20                 xbt_swag_offset(var, variable_set_hookup));
21   xbt_swag_init(&(l->constraint_set),
22                 xbt_swag_offset(cnst, constraint_set_hookup));
23
24   return l;
25 }
26
27 void lmm_system_free(lmm_system_t sys)
28 {
29   lmm_variable_t var = NULL;
30   lmm_constraint_t cnst = NULL;
31
32   while ((var = extract_variable(sys)))
33     lmm_var_free(var);
34
35   while ((cnst = extract_constraint(sys)))
36     lmm_cnst_free(cnst);
37
38   xbt_free(sys);
39 }
40
41 static void lmm_var_free(lmm_variable_t var)
42 {
43   int i;
44   lmm_element_t elem = NULL;
45
46   for (i = 0; i < var->cnsts_number; i++) {
47     elem = &var->cnsts[i];
48     xbt_swag_extract(elem, &(elem->constraint->element_set));
49   }
50   xbt_free(var->cnsts);
51   xbt_free(var);
52 }
53
54 static void lmm_cnst_free(lmm_constraint_t cnst)
55 {
56 /*   xbt_assert0(xbt_fifo_size(&(cnst->row)), */
57 /*               "This list should be empty!"); */
58   xbt_free(cnst);
59 }
60
61 lmm_constraint_t lmm_constraint_new(lmm_system_t sys,   /* void *id, */
62                                     FLOAT bound_value)
63 {
64   lmm_constraint_t cnst = NULL;
65   s_lmm_element_t elem;
66
67   cnst = xbt_new0(s_lmm_constraint_t, 1);
68 /*   cnst->id = id; */
69   xbt_swag_init(&(cnst->element_set),
70                 xbt_swag_offset(elem, element_set_hookup));
71
72   cnst->bound = bound_value;
73   cnst->usage = UNUSED_CONSTRAINT;
74   insert_constraint(sys, cnst);
75
76   return cnst;
77 }
78
79 void lmm_constraint_free(lmm_system_t sys, lmm_constraint_t cnst)
80 {
81   remove_constraint(sys, cnst);
82   lmm_cnst_free(cnst);
83 }
84
85 lmm_variable_t lmm_variable_new(lmm_system_t sys,       /* void *id, */
86                                 FLOAT weight_value, FLOAT bound,
87                                 int number_of_constraints)
88 {
89   lmm_variable_t var = NULL;
90
91   var = xbt_new0(s_lmm_variable_t, 1);
92 /*   var->id = id; */
93   var->cnsts = xbt_new0(s_lmm_element_t, number_of_constraints);
94   var->cnsts_size = number_of_constraints;
95   /* var->cnsts_number = 0; *//* Useless because of the calloc  */
96   var->weight_value = weight_value;
97   var->bound = bound;
98   var->value = -1;
99   insert_variable(sys, var);
100
101   return var;
102 }
103
104 void lmm_variable_free(lmm_system_t sys, lmm_variable_t var)
105 {
106   remove_variable(sys, var);
107   lmm_var_free(var);
108 }
109
110 void lmm_add_constraint(lmm_system_t sys, lmm_constraint_t cnst,
111                         lmm_variable_t var, FLOAT value)
112 {
113   lmm_element_t elem = NULL;
114
115 /*   xbt_assert0(var->cnsts_number<var->cnsts_size,"Too much constraint"); */
116
117   elem = &(var->cnsts[var->cnsts_number++]);
118
119   elem->value = value;
120   elem->constraint = cnst;
121
122   insert_elem_in_constraint(elem);
123
124 }
125
126 void lmm_solve(lmm_system_t sys)
127 {
128   void **id_table = NULL;
129
130   /* Compute Usage and store the variables that reach the maximum */
131
132   do {
133     /* Fix the variables that have to be */
134
135     /* Update Usage and store the variables that reach the maximum */
136
137   } while (1);
138
139 }