chiark / gitweb /
prepare cplex input file for glpsol
[ypp-sc-tools.main.git] / yarrg / web / routetrade
index e6aff44eece5329da0bccc817e641ac7597f7d92..1924503b9ce706b888c2450de30403ba19f96a8b 100644 (file)
@@ -100,16 +100,17 @@ my $stmt= "
           OR   ", @flow_conds)."
        )
          AND   buy.price > sell.price
-       GROUP BY commodname, commods.commodid,
-               org_id, org_price, dst_id, dst_price
+       GROUP BY commods.commodid, org_id, org_price, dst_id, dst_price
        ORDER BY org_name, dst_name, profit DESC, commodname,
                 org_price, dst_price DESC
      ";
 
 my $sth= $dbh->prepare($stmt);
 $sth->execute(@query_params);
+my @flows;
 
 </%perl>
+
 % if ($qa->{'debug'}) {
 <pre>
 <% $stmt |h %>
@@ -117,7 +118,64 @@ $sth->execute(@query_params);
 </pre>
 % }
 
-<& dumpqueryresults, sth =>$sth &>
+% {
+<& dumpqueryresults:start, sth => $sth &>
+%   my $flow;
+%   while ($flow= $sth->fetchrow_hashref()) {
+%      $flow->{Ix}= @flows;
+%      $flow->{Var}= "f$flow->{Ix}";
+%      push @flows, $flow;
+<& dumpqueryresults:row, sth => $sth, row => $flow &>
+%   }
+<& dumpqueryresults:end &>
+% }
+
+<%perl>
+
+my $cplex= "
+Maximize
+
+  totalprofit:
+                  ".(join " +
+                  ", map { "$_->{profit} $_->{Var}" } @flows)."
+
+Subject To
+";
+
+my %avail_csts;
+foreach my $flow (@flows) {
+       foreach my $od (qw(org dst)) {
+               my $cstname= join '_',
+                       'avail',
+                       $flow->{'commodid'},
+                       $od,
+                       $flow->{"${od}_id"},
+                       $flow->{"${od}_price"};
+               push @{ $avail_csts{$cstname}{Flows} }, $flow->{Var};
+               $avail_csts{$cstname}{Qty}= $flow->{"${od}_qty"};
+       }
+}
+foreach my $cstname (sort keys %avail_csts) {
+       my $c= $avail_csts{$cstname};
+       $cplex .= "
+   ".  sprintf("%-30s","$cstname:")." ".
+       join("+", @{ $c->{Flows} }).
+       " <= ".$c->{Qty}."\n";
+}
+
+$cplex.= "
+Bounds
+        ".(join "
+        ", map { "$_->{Var} >= 0" } @flows)."
+
+End
+";
+# glpsol --cpxlp /dev/stdin <t.cplex -o /dev/stdout
+
+</%perl>
+<pre>
+<% $cplex |h %>
+</pre>
 
 <%init>
 use CommodsWeb;