1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/examples/sample.c Mon Dec 06 13:09:21 2010 +0100
1.3 @@ -0,0 +1,52 @@
1.4 +/* sample.c */
1.5 +
1.6 +#include <stdio.h>
1.7 +#include <stdlib.h>
1.8 +#include <glpk.h>
1.9 +
1.10 +int main(void)
1.11 +{ glp_prob *lp;
1.12 + int ia[1+1000], ja[1+1000];
1.13 + double ar[1+1000], z, x1, x2, x3;
1.14 +s1: lp = glp_create_prob();
1.15 +s2: glp_set_prob_name(lp, "sample");
1.16 +s3: glp_set_obj_dir(lp, GLP_MAX);
1.17 +s4: glp_add_rows(lp, 3);
1.18 +s5: glp_set_row_name(lp, 1, "p");
1.19 +s6: glp_set_row_bnds(lp, 1, GLP_UP, 0.0, 100.0);
1.20 +s7: glp_set_row_name(lp, 2, "q");
1.21 +s8: glp_set_row_bnds(lp, 2, GLP_UP, 0.0, 600.0);
1.22 +s9: glp_set_row_name(lp, 3, "r");
1.23 +s10: glp_set_row_bnds(lp, 3, GLP_UP, 0.0, 300.0);
1.24 +s11: glp_add_cols(lp, 3);
1.25 +s12: glp_set_col_name(lp, 1, "x1");
1.26 +s13: glp_set_col_bnds(lp, 1, GLP_LO, 0.0, 0.0);
1.27 +s14: glp_set_obj_coef(lp, 1, 10.0);
1.28 +s15: glp_set_col_name(lp, 2, "x2");
1.29 +s16: glp_set_col_bnds(lp, 2, GLP_LO, 0.0, 0.0);
1.30 +s17: glp_set_obj_coef(lp, 2, 6.0);
1.31 +s18: glp_set_col_name(lp, 3, "x3");
1.32 +s19: glp_set_col_bnds(lp, 3, GLP_LO, 0.0, 0.0);
1.33 +s20: glp_set_obj_coef(lp, 3, 4.0);
1.34 +s21: ia[1] = 1, ja[1] = 1, ar[1] = 1.0; /* a[1,1] = 1 */
1.35 +s22: ia[2] = 1, ja[2] = 2, ar[2] = 1.0; /* a[1,2] = 1 */
1.36 +s23: ia[3] = 1, ja[3] = 3, ar[3] = 1.0; /* a[1,3] = 1 */
1.37 +s24: ia[4] = 2, ja[4] = 1, ar[4] = 10.0; /* a[2,1] = 10 */
1.38 +s25: ia[5] = 3, ja[5] = 1, ar[5] = 2.0; /* a[3,1] = 2 */
1.39 +s26: ia[6] = 2, ja[6] = 2, ar[6] = 4.0; /* a[2,2] = 4 */
1.40 +s27: ia[7] = 3, ja[7] = 2, ar[7] = 2.0; /* a[3,2] = 2 */
1.41 +s28: ia[8] = 2, ja[8] = 3, ar[8] = 5.0; /* a[2,3] = 5 */
1.42 +s29: ia[9] = 3, ja[9] = 3, ar[9] = 6.0; /* a[3,3] = 6 */
1.43 +s30: glp_load_matrix(lp, 9, ia, ja, ar);
1.44 +s31: glp_simplex(lp, NULL);
1.45 +s32: z = glp_get_obj_val(lp);
1.46 +s33: x1 = glp_get_col_prim(lp, 1);
1.47 +s34: x2 = glp_get_col_prim(lp, 2);
1.48 +s35: x3 = glp_get_col_prim(lp, 3);
1.49 +s36: printf("\nz = %g; x1 = %g; x2 = %g; x3 = %g\n",
1.50 + z, x1, x2, x3);
1.51 +s37: glp_delete_prob(lp);
1.52 + return 0;
1.53 +}
1.54 +
1.55 +/* eof */