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