examples/sample.c
changeset 1 c445c931472f
     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 */