dvoisin/aro
Archived
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

TA GUEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEUUUUUUULE

Browse Source
This commit is contained in:
papush! 2019-10-22 20:05:01 +02:00
parent 3f980f1d09
commit 9cb04f58bd
2 changed files with 24 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
src/master.c
View File

@ -50,6 +50,9 @@ int perfect_rolls() {
int nb_coeffs = 4; int nb_coeffs = 4;
double y1, y2, y3, y4; double y1, y2, y3, y4;
int a[4] = {0}; int a[4] = {0};
glp_smcp params = {0};
glp_init_smcp(&params);
params.msg_lev = GLP_MSG_OFF;
while (1) { while (1) {
/* printf("nb_coeffs=%d\n", nb_coeffs); */ /* printf("nb_coeffs=%d\n", nb_coeffs); */
/* for (int i = 0; i < nb_coeffs + 4; i++) { */ /* for (int i = 0; i < nb_coeffs + 4; i++) { */
@ -61,19 +64,21 @@ int perfect_rolls() {
printf("col_num=%d\n", col_num); printf("col_num=%d\n", col_num);
/* solve problem */ /* solve problem */
glp_prob *tmp_prob = glp_create_prob(); glp_simplex(lp, &params);
glp_copy_prob(tmp_prob, lp, GLP_ON); y1 = glp_get_row_dual(lp, 1);
glp_simplex(tmp_prob, NULL); y2 = glp_get_row_dual(lp, 2);
y1 = glp_get_row_dual(tmp_prob, 1); y3 = glp_get_row_dual(lp, 3);
y2 = glp_get_row_dual(tmp_prob, 2); y4 = glp_get_row_dual(lp, 4);
y3 = glp_get_row_dual(tmp_prob, 3);
y4 = glp_get_row_dual(tmp_prob, 4);
column_generator(y1, y2, y3, y4, a, num_col_sp); column_generator(y1, y2, y3, y4, a, num_col_sp);
if (a[0] > 1) { if (a[0] > 1) {
/* for (int i = 1; i <= 4; i++) { */ fputs("Après résolution du sous-problème, on ajoute : ", stdout);
/* printf("a[%d]=%d\n", i, a[i]); */ int i = 1;
/* } */ for (i = 1; i < 4; i++) {
printf("%d, ", a[i]);
}
printf("%d\n", a[i]);
add_column(lp, 1.0); add_column(lp, 1.0);
for (int i = 1; i < 5; i++) { for (int i = 1; i < 5; i++) {
/* a[col_num, i] = a[i] */ /* a[col_num, i] = a[i] */

11
src/sub.c
View File

@ -45,8 +45,15 @@ void column_generator(int y1, int y2, int y3, int y4, int a[static 4], int num_c
glp_load_matrix(sous_prob, 4, ia, ja, ar); glp_load_matrix(sous_prob, 4, ia, ja, ar);
/* solve problem */ /* solve problem */
glp_simplex(sous_prob, NULL); glp_smcp params = {0};
glp_intopt(sous_prob, NULL); glp_init_smcp(&params);
params.msg_lev = GLP_MSG_OFF;
glp_iocp iocp = {0};
glp_init_iocp(&iocp);
iocp.msg_lev = GLP_MSG_OFF;
glp_simplex(sous_prob, &params);
glp_intopt(sous_prob, &iocp);
a[0] = glp_get_obj_val(sous_prob); a[0] = glp_get_obj_val(sous_prob);
a[1] = glp_mip_col_val(sous_prob, 1); a[1] = glp_mip_col_val(sous_prob, 1);