+ 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;
+ 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_constraint(t, &iteps[s]->collect, t->src_price,t->src_qty,"src");
+ avail_constraint(t, &iteps[d]->deliver, t->dst_price,t->dst_qty,"dst");
+
+ 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);
+ }
+
+ lpx_set_col_bnds(lp, col, LPX_LO, 0, 0);
+ lpx_set_obj_coef(lp, col,
+ (t->dst_price - t->src_price) * loss_factor);
+ lpx_set_mat_col(lp, col, nconstraint_rows,
+ constraint_rows, constraint_coeffs);
+ if (DEBUGP(value)) {
+ char *name= masprintf("trade_commod%d_port%d_at%d_port%d_at%d",
+ t->commodid, s, t->src_price, d, t->dst_price);
+ lpx_set_col_name(lp, col, name);
+ free(name);
+ }
+ } /* while (tradestodo-- >0) */
+
+ /*----- that's done adding these trades to the LP problem -----*/
+
+ next_d:;
+ } /* for (d) */
+ next_s:;
+ } /* for (s) */
+
+ lpx_delete_prob(lp);
+ lp= 0;
+
+ obstack_free(&ites_obstack, free_to_reset_obstack);