+ assert(!lp);
+ lp= lpx_create_prob();
+ lpx_set_obj_dir(lp, LPX_MAX);
+ lpx_set_int_parm(lp, LPX_K_MSGLEV, DEBUGP(lp) ? 3 : 1);
+
+ if (DEBUGP(value)) {
+ lpx_set_prob_name(lp,(char*)"value_route");
+ lpx_set_obj_name(lp,(char*)"profit");
+ }
+
+ int legs= nislands-1;
+ int mass_constraints= setup_leg_constraints(max_mass, legs, "mass");
+ int volu_constraints= setup_leg_constraints(max_volu, legs, "volu");
+ int capi_constraints= setup_leg_constraints(max_capi, legs, "capi");
+
+ double delay_slot_loss_factor= 1.0;
+ for (s=0;
+ s<nislands;
+ s++, delay_slot_loss_factor *= LOSS_FACTOR_PER_DELAY_SLOT) {
+ int si= islands[s];
+
+ for (d=s; d<nislands; d++) {
+ int di= islands[d];
+ int already_d;
+ for (already_d=s+1; already_d<d; already_d++)
+ if (islands[already_d] == di)
+ /* visited this island already since we left s, uninteresting */
+ goto next_d;
+
+ if (d>s && di==si)
+ /* route has returned to si, no need to think more about s */
+ goto next_s;
+
+ /*----- actually add these trades to the LP problem -----*/
+
+ IslandPair *ip= ipair_get(islands[s], islands[d]);
+ TradesBlock *block= ip->trades;
+ int tradestodo= ip->ntrades;
+ if (!tradestodo)
+ goto next_d;
+
+ int inblock= 0;
+ int col= lpx_add_cols(lp,ip->ntrades);
+
+ double loss_factor= delay_slot_loss_factor * ip->distance_loss_factor;
+
+ while (tradestodo-- >0) {
+ if (inblock >= TRADES_PER_BLOCK) {
+ block= block->next;
+ inblock= 0;
+ }
+ Trade *t= &block->t[inblock++];
+
+ debugf(" TRADE %d#%d..%d#%d %d %d-%d\n",
+ si,s, di,d, t->commodid, t->src_price, t->dst_price);
+
+ nconstraint_rows=0;
+
+ avail_c(t, &iteps[s]->collect, t->src_price,"src", si, ss_ite_sell);
+ avail_c(t, &iteps[d]->deliver, t->dst_price,"dst", di, ss_ite_buy);
+
+ int leg;
+ for (leg=s; leg<d; leg++) {
+ add_leg_c(mass_constraints, leg, commodstable[t->commodid].mass);
+ add_leg_c(volu_constraints, leg, commodstable[t->commodid].volu);
+ add_leg_c(capi_constraints, leg, t->src_price);
+ }
+
+ double unit_profit= (t->dst_price - t->src_price) * loss_factor;
+ debugf(" unit profit %f\n", unit_profit);
+assert(unit_profit < 1e6);
+
+ lpx_set_col_bnds(lp, col, LPX_LO, 0, 0);
+ lpx_set_obj_coef(lp, col, unit_profit);
+ lpx_set_mat_col(lp, col, nconstraint_rows,
+ constraint_rows, constraint_coeffs);
+
+ if (DEBUGP(value)) {
+ char *name= masprintf("c%d_p%d_%d_p%d_%d",
+ t->commodid, s, t->src_price, d, t->dst_price);
+ lpx_set_col_name(lp, col, name);
+ free(name);
+ }
+
+ col++;
+ } /* while (tradestodo-- >0) */
+
+ /*----- that's done adding these trades to the LP problem -----*/
+
+ next_d:;
+ } /* for (d) */
+ next_s:;
+ } /* for (s) */
+
+ if (DEBUGP(lp))
+ lpx_write_cpxlp(lp, (char*)"/dev/stdout");
+
+ int ipr= lpx_interior(lp);
+ assert(ipr==LPX_E_OK);
+
+ if (DEBUGP(lp))
+ lpx_print_ips(lp, (char*)"/dev/stdout");
+
+ assert(lpx_ipt_status(lp) == LPX_T_OPT);
+ double profit= lpx_ipt_obj_val(lp);
+
+ lpx_delete_prob(lp);
+ lp= 0;
+
+ return profit;