src/work/athos/lp/lp_sample.cc
changeset 1317 83f80464f111
child 1318 88edb143a87a
equal deleted inserted replaced
-1:000000000000 0:d29e439ead40
       
     1 //#include <stdio.h>
       
     2 //#include <stdlib.h>
       
     3 #include <iostream>
       
     4 #include <lemon/lp_glpk.h>
       
     5 using namespace lemon;
       
     6 
       
     7 int main()
       
     8 {     
       
     9   LpGlpk lp;
       
    10   typedef LpGlpk::Row Row;
       
    11   typedef LpGlpk::Col Col;
       
    12   //std::cout << "Hello" << std::endl;
       
    13 
       
    14 //   LPX *lp;
       
    15 //       int ia[1+1000], ja[1+1000];
       
    16 //       double ar[1+1000], Z, x1, x2, x3;
       
    17 // s1:   lp = lpx_create_prob();
       
    18 // s2:   lpx_set_prob_name(lp, "sample");
       
    19 // s3:   lpx_set_obj_dir(lp, LPX_MAX);
       
    20   lp.max();
       
    21 // s4:   lpx_add_rows(lp, 3);
       
    22 // s5:   lpx_set_row_name(lp, 1, "p");
       
    23 // s6:   lpx_set_row_bnds(lp, 1, LPX_UP, 0.0, 100.0);
       
    24 
       
    25 //  Row p = lp.addRow();
       
    26 
       
    27 // s7:   lpx_set_row_name(lp, 2, "q");
       
    28 // s8:   lpx_set_row_bnds(lp, 2, LPX_UP, 0.0, 600.0);
       
    29 
       
    30 //  Row q = lp.addRow();
       
    31 
       
    32 // s9:   lpx_set_row_name(lp, 3, "r");
       
    33 // s10:  lpx_set_row_bnds(lp, 3, LPX_UP, 0.0, 300.0);
       
    34 
       
    35 //  Row r = lp.addRow();
       
    36 
       
    37 // s11:  lpx_add_cols(lp, 3);
       
    38 // s12:  lpx_set_col_name(lp, 1, "x1");
       
    39 // s13:  lpx_set_col_bnds(lp, 1, LPX_LO, 0.0, 0.0);
       
    40 
       
    41 //  Col x1 = lp.addCol();
       
    42 
       
    43 // s14:  lpx_set_obj_coef(lp, 1, 10.0);
       
    44 // s15:  lpx_set_col_name(lp, 2, "x2");
       
    45 // s16:  lpx_set_col_bnds(lp, 2, LPX_LO, 0.0, 0.0);
       
    46 // s17:  lpx_set_obj_coef(lp, 2, 6.0);
       
    47 // s18:  lpx_set_col_name(lp, 3, "x3");
       
    48 // s19:  lpx_set_col_bnds(lp, 3, LPX_LO, 0.0, 0.0);
       
    49 // s20:  lpx_set_obj_coef(lp, 3, 4.0);
       
    50 // s21:  ia[1] = 1, ja[1] = 1, ar[1] =  1.0; /* a[1,1] =  1 */
       
    51 // s22:  ia[2] = 1, ja[2] = 2, ar[2] =  1.0; /* a[1,2] =  1 */
       
    52 // s23:  ia[3] = 1, ja[3] = 3, ar[3] =  1.0; /* a[1,3] =  1 */
       
    53 // s24:  ia[4] = 2, ja[4] = 1, ar[4] = 10.0; /* a[2,1] = 10 */
       
    54 // s25:  ia[5] = 3, ja[5] = 1, ar[5] =  2.0; /* a[3,1] =  2 */
       
    55 // s26:  ia[6] = 2, ja[6] = 2, ar[6] =  4.0; /* a[2,2] =  4 */
       
    56 // s27:  ia[7] = 3, ja[7] = 2, ar[7] =  2.0; /* a[3,2] =  2 */
       
    57 // s28:  ia[8] = 2, ja[8] = 3, ar[8] =  5.0; /* a[2,3] =  5 */
       
    58 // s29:  ia[9] = 3, ja[9] = 3, ar[9] =  6.0; /* a[3,3] =  6 */
       
    59 // s30:  lpx_load_matrix(lp, 9, ia, ja, ar);
       
    60 // s31:  lpx_simplex(lp);
       
    61 // s32:  Z = lpx_get_obj_val(lp);
       
    62 // s33:  x1 = lpx_get_col_prim(lp, 1);
       
    63 // s34:  x2 = lpx_get_col_prim(lp, 2);
       
    64 // s35:  x3 = lpx_get_col_prim(lp, 3);
       
    65 // s36:  printf("\nZ = %g; x1 = %g; x2 = %g; x3 = %g\n", Z, x1, x2, x3);
       
    66 // s37:  lpx_delete_prob(lp);
       
    67 //       return 0;
       
    68   return 0;
       
    69 }