| 1 | /* sample.c */ |
|---|
| 2 | |
|---|
| 3 | #include <stdio.h> |
|---|
| 4 | #include <stdlib.h> |
|---|
| 5 | #include <glpk.h> |
|---|
| 6 | |
|---|
| 7 | int main(void) |
|---|
| 8 | { glp_prob *lp; |
|---|
| 9 | int ia[1+1000], ja[1+1000]; |
|---|
| 10 | double ar[1+1000], z, x1, x2, x3; |
|---|
| 11 | s1: lp = glp_create_prob(); |
|---|
| 12 | s2: glp_set_prob_name(lp, "sample"); |
|---|
| 13 | s3: glp_set_obj_dir(lp, GLP_MAX); |
|---|
| 14 | s4: glp_add_rows(lp, 3); |
|---|
| 15 | s5: glp_set_row_name(lp, 1, "p"); |
|---|
| 16 | s6: glp_set_row_bnds(lp, 1, GLP_UP, 0.0, 100.0); |
|---|
| 17 | s7: glp_set_row_name(lp, 2, "q"); |
|---|
| 18 | s8: glp_set_row_bnds(lp, 2, GLP_UP, 0.0, 600.0); |
|---|
| 19 | s9: glp_set_row_name(lp, 3, "r"); |
|---|
| 20 | s10: glp_set_row_bnds(lp, 3, GLP_UP, 0.0, 300.0); |
|---|
| 21 | s11: glp_add_cols(lp, 3); |
|---|
| 22 | s12: glp_set_col_name(lp, 1, "x1"); |
|---|
| 23 | s13: glp_set_col_bnds(lp, 1, GLP_LO, 0.0, 0.0); |
|---|
| 24 | s14: glp_set_obj_coef(lp, 1, 10.0); |
|---|
| 25 | s15: glp_set_col_name(lp, 2, "x2"); |
|---|
| 26 | s16: glp_set_col_bnds(lp, 2, GLP_LO, 0.0, 0.0); |
|---|
| 27 | s17: glp_set_obj_coef(lp, 2, 6.0); |
|---|
| 28 | s18: glp_set_col_name(lp, 3, "x3"); |
|---|
| 29 | s19: glp_set_col_bnds(lp, 3, GLP_LO, 0.0, 0.0); |
|---|
| 30 | s20: glp_set_obj_coef(lp, 3, 4.0); |
|---|
| 31 | s21: ia[1] = 1, ja[1] = 1, ar[1] = 1.0; /* a[1,1] = 1 */ |
|---|
| 32 | s22: ia[2] = 1, ja[2] = 2, ar[2] = 1.0; /* a[1,2] = 1 */ |
|---|
| 33 | s23: ia[3] = 1, ja[3] = 3, ar[3] = 1.0; /* a[1,3] = 1 */ |
|---|
| 34 | s24: ia[4] = 2, ja[4] = 1, ar[4] = 10.0; /* a[2,1] = 10 */ |
|---|
| 35 | s25: ia[5] = 3, ja[5] = 1, ar[5] = 2.0; /* a[3,1] = 2 */ |
|---|
| 36 | s26: ia[6] = 2, ja[6] = 2, ar[6] = 4.0; /* a[2,2] = 4 */ |
|---|
| 37 | s27: ia[7] = 3, ja[7] = 2, ar[7] = 2.0; /* a[3,2] = 2 */ |
|---|
| 38 | s28: ia[8] = 2, ja[8] = 3, ar[8] = 5.0; /* a[2,3] = 5 */ |
|---|
| 39 | s29: ia[9] = 3, ja[9] = 3, ar[9] = 6.0; /* a[3,3] = 6 */ |
|---|
| 40 | s30: glp_load_matrix(lp, 9, ia, ja, ar); |
|---|
| 41 | s31: glp_simplex(lp, NULL); |
|---|
| 42 | s32: z = glp_get_obj_val(lp); |
|---|
| 43 | s33: x1 = glp_get_col_prim(lp, 1); |
|---|
| 44 | s34: x2 = glp_get_col_prim(lp, 2); |
|---|
| 45 | s35: x3 = glp_get_col_prim(lp, 3); |
|---|
| 46 | s36: printf("\nz = %g; x1 = %g; x2 = %g; x3 = %g\n", |
|---|
| 47 | z, x1, x2, x3); |
|---|
| 48 | s37: glp_delete_prob(lp); |
|---|
| 49 | return 0; |
|---|
| 50 | } |
|---|
| 51 | |
|---|
| 52 | /* eof */ |
|---|
