diff -r b3564d0e9c60 -r c5984e925384 src/work/athos/min_cost_flow.cc --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/src/work/athos/min_cost_flow.cc Tue May 25 12:31:18 2004 +0000 @@ -0,0 +1,110 @@ +#include +#include "test_tools.h" +#include +#include +//#include +//#include + +using namespace std; +using namespace hugo; + + + +bool passed = true; +/* +void check(bool rc, char *msg="") { + passed = passed && rc; + if(!rc) { + std::cerr << "Test failed! ("<< msg << ")" << std::endl; \ + + + } +} +*/ + + +int main() +{ + + typedef ListGraph::Node Node; + typedef ListGraph::Edge Edge; + + ListGraph graph; + + //Ahuja könyv példája + + Node s=graph.addNode(); + Node v1=graph.addNode(); + Node v2=graph.addNode(); + Node v3=graph.addNode(); + Node v4=graph.addNode(); + Node v5=graph.addNode(); + Node t=graph.addNode(); + + Edge s_v1=graph.addEdge(s, v1); + Edge v1_v2=graph.addEdge(v1, v2); + Edge s_v3=graph.addEdge(s, v3); + Edge v2_v4=graph.addEdge(v2, v4); + Edge v2_v5=graph.addEdge(v2, v5); + Edge v3_v5=graph.addEdge(v3, v5); + Edge v4_t=graph.addEdge(v4, t); + Edge v5_t=graph.addEdge(v5, t); + + + ListGraph::EdgeMap length(graph); + + length.set(s_v1, 6); + length.set(v1_v2, 4); + length.set(s_v3, 10); + length.set(v2_v4, 5); + length.set(v2_v5, 1); + length.set(v3_v5, 4); + length.set(v4_t, 8); + length.set(v5_t, 8); + + /* + ListGraph::EdgeMap capacity(graph); + + capacity.set(s_v1, 2); + capacity.set(v1_v2, 2); + capacity.set(s_v3, 1); + capacity.set(v2_v4, 1); + capacity.set(v2_v5, 1); + capacity.set(v3_v5, 1); + capacity.set(v4_t, 1); + capacity.set(v5_t, 2); + */ + + // ConstMap const1map(1); + std::cout << "Enhanced capacity scaling algorithm test (for the mincostflow problem)..." << std::endl; + + MinCostFlow< ListGraph, ListGraph::EdgeMap, ListGraph::NodeMap > + min_cost_flow_test(graph, length, supply_demand); + + int k=1; + + /* + check( min_cost_flow_test.run(s,t,k) == 1 && min_cost_flow_test.totalLength() == 19,"One path, total length should be 19"); + + check(min_cost_flow_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?"); + + k=2; + + check( min_cost_flow_test.run(s,t,k) == 2 && min_cost_flow_test.totalLength() == 41,"Two paths, total length should be 41"); + + check(min_cost_flow_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?"); + + + k=4; + + check( min_cost_flow_test.run(s,t,k) == 3 && min_cost_flow_test.totalLength() == 64,"Three paths, total length should be 64"); + + check(min_cost_flow_test.checkComplementarySlackness(), "Is the primal-dual solution pair really optimal?"); + + + cout << (passed ? "All tests passed." : "Some of the tests failed!!!") + << endl; + + return passed ? 0 : 1; + */ +}