test_DualMinCostFlow.cpp

Go to the documentation of this file.

 
#include <iostream>
#include <limbo/solvers/DualMinCostFlow.h>

void test(std::string const& filename, int alg)
{
    typedef limbo::solvers::LinearModel<int, int> model_type; 
    model_type optModel; 
    optModel.read(filename); 
 
    // print problem 
    optModel.print(std::cout); 
 
    // solve 
    limbo::solvers::MinCostFlowSolver<int, int>* minCostFlowSolver = NULL; 
    switch (alg)
    {
        case 0:
            minCostFlowSolver = new limbo::solvers::CostScaling<int, int>(); 
            break;
        case 1:
            minCostFlowSolver = new limbo::solvers::CapacityScaling<int, int>(); 
            break;
        case 2:
            minCostFlowSolver = new limbo::solvers::NetworkSimplex<int, int>(); 
            break;
        case 3:
        default:
            minCostFlowSolver = new limbo::solvers::CycleCanceling<int, int>(); 
            break; 
    }
    limbo::solvers::DualMinCostFlow<int, int> solver (&optModel); 
    limbo::solvers::SolverProperty status = solver(minCostFlowSolver);
    //limbo::solvers::SolverProperty status = solver();
    std::cout << "Problem solved " << limbo::solvers::toString(status) << "\n";
    delete minCostFlowSolver;
 
    // print solutions 
    optModel.printSolution(std::cout);
    // print problem 
    optModel.print(std::cout); 
 
    // print graph with solution information 
    solver.printGraph(true);
}

int main(int argc, char** argv)
{
    if (argc > 1)
    {
        int alg = 0; 
        if (argc > 2)
            alg = atoi(argv[2]);
        // test file API 
        test(argv[1], alg);
    }
    else 
        std::cout << "at least 1 argument required\n";
 
    return 0;
}