TSP_driver.cpp File Reference

#include "drivers/tsp/TSP_driver.h"
#include <string>
#include <sstream>
#include <vector>
#include <algorithm>
#include "cpp_common/Dmatrix.h"
#include "tsp/pgr_tsp.hpp"
#include "cpp_common/pgr_alloc.hpp"
#include "cpp_common/pgr_assert.h"

Include dependency graph for TSP_driver.cpp:

Go to the source code of this file.

Functions
void	do_pgr_tsp (Matrix_cell_t *distances, size_t total_distances, int64_t start_vid, int64_t end_vid, double initial_temperature, double final_temperature, double cooling_factor, int64_t tries_per_temperature, int64_t max_changes_per_temperature, int64_t max_consecutive_non_changes, bool randomize, double time_limit, General_path_element_t **return_tuples, size_t *return_count, char **log_msg, char **notice_msg, char **err_msg)

Function Documentation

do_pgr_tsp()

void do_pgr_tsp	(	Matrix_cell_t *	distances,
		size_t	total_distances,
		int64_t	start_vid,
		int64_t	end_vid,
		double	initial_temperature,
		double	final_temperature,
		double	cooling_factor,
		int64_t	tries_per_temperature,
		int64_t	max_changes_per_temperature,
		int64_t	max_consecutive_non_changes,
		bool	randomize,
		double	time_limit,
		General_path_element_t **	return_tuples,
		size_t *	return_count,
		char **	log_msg,
		char **	notice_msg,
		char **	err_msg
	)

Definition at line 45 of file TSP_driver.cpp.

                        {
    std::ostringstream log;
    std::ostringstream notice;
    std::ostringstream err;
 
    try {
        std::vector <Matrix_cell_t> data_costs(
                distances,
                distances + total_distances);
 
        pgrouting::tsp::Dmatrix costs(data_costs);
 
        if (!costs.has_no_infinity()) {
            err << "An Infinity value was found on the Matrix";
            *err_msg = pgr_msg(err.str().c_str());
            return;
        }
 
        if (!costs.is_symmetric()) {
            err << "A Non symmetric Matrix was given as input";
            *err_msg = pgr_msg(err.str().c_str());
            return;
        }
 
        double real_cost = -1;
 
        size_t idx_start = costs.has_id(start_vid) ?
            costs.get_index(start_vid) : 0;
 
        size_t idx_end = costs.has_id(end_vid) ?
            costs.get_index(end_vid) : 0;
 
        if (costs.has_id(start_vid)
                && costs.has_id(end_vid)
                && start_vid != end_vid) {
            /* An ending vertex needs to be by the starting vertex */
            real_cost = costs.distance(idx_start, idx_end);
            costs.set(idx_start, idx_end, 0);
        }
 
 
        log << "Processing Information \n"
            << "Initializing tsp class --->";
        pgrouting::tsp::TSP<pgrouting::tsp::Dmatrix> tsp(costs);
 
 
        log << " tsp.greedyInitial --->";
        tsp.greedyInitial(idx_start);
 
 
 
        log << " tsp.annealing --->";
        tsp.annealing(
                initial_temperature,
                final_temperature,
                cooling_factor,
                tries_per_temperature,
                max_changes_per_temperature,
                max_consecutive_non_changes,
                randomize,
                time_limit);
        log << " OK\n";
        log << tsp.get_log();
        log << tsp.get_stats();
 
        auto bestTour(tsp.get_tour());
 
        if (costs.has_id(start_vid)
                && costs.has_id(end_vid)
                && start_vid != end_vid) {
            costs.set(idx_start, idx_end, real_cost);
        }
 
 
        log << "\nBest cost reached = " << costs.tourCost(bestTour);
 
        auto start_ptr = std::find(
                bestTour.cities.begin(),
                bestTour.cities.end(),
                idx_start);
 
 
        std::rotate(
                bestTour.cities.begin(),
                start_ptr,
                bestTour.cities.end());
 
        if (costs.has_id(start_vid)
                && costs.has_id(end_vid)
                && start_vid != end_vid) {
            if (*(bestTour.cities.begin() + 1) == idx_end) {
                std::reverse(
                        bestTour.cities.begin() + 1,
                        bestTour.cities.end());
            }
        }
 
 
        std::vector< General_path_element_t > result;
        result.reserve(bestTour.cities.size() + 1);
        pgassert(bestTour.cities.size() == costs.size());
 
        bestTour.cities.push_back(bestTour.cities.front());
 
        auto prev_id = bestTour.cities.front();
        double agg_cost = 0;
        for (const auto &id : bestTour.cities) {
            if (id == prev_id) continue;
            General_path_element_t data;
            data.node = costs.get_id(prev_id);
            data.edge = static_cast<int64_t>(prev_id);
            data.cost = costs.distance(prev_id, id);
            data.agg_cost = agg_cost;
            result.push_back(data);
            agg_cost += data.cost;
            prev_id = id;
        }
 
        /* inserting the returning to starting point */
        {
            General_path_element_t data;
            data.node = costs.get_id(bestTour.cities.front());
            data.edge = static_cast<int64_t>(bestTour.cities.front());
            data.cost = costs.distance(prev_id, bestTour.cities.front());
            agg_cost += data.cost;
            data.agg_cost = agg_cost;
            result.push_back(data);
        }
 
        pgassert(result.size() == bestTour.cities.size());
        *return_count = bestTour.size();
        (*return_tuples) = pgr_alloc(result.size(), (*return_tuples));
 
        /* store the results */
        int seq = 0;
        for (const auto &row : result) {
            (*return_tuples)[seq] = row;
            ++seq;
        }
 
        *log_msg = log.str().empty()?
            *log_msg :
            pgr_msg(log.str().c_str());
        *notice_msg = notice.str().empty()?
            *notice_msg :
            pgr_msg(notice.str().c_str());
    } catch (AssertFailedException &except) {
        (*return_tuples) = pgr_free(*return_tuples);
        (*return_count) = 0;
        err << except.what();
        *err_msg = pgr_msg(err.str().c_str());
        *log_msg = pgr_msg(log.str().c_str());
    } catch (std::exception &except) {
        (*return_tuples) = pgr_free(*return_tuples);
        (*return_count) = 0;
        err << except.what();
        *err_msg = pgr_msg(err.str().c_str());
        *log_msg = pgr_msg(log.str().c_str());
    } catch(...) {
        (*return_tuples) = pgr_free(*return_tuples);
        (*return_count) = 0;
        err << "Caught unknown exception!";
        *err_msg = pgr_msg(err.str().c_str());
        *log_msg = pgr_msg(log.str().c_str());
    }
}

References General_path_element_t::agg_cost, pgrouting::tsp::TSP< MATRIX >::annealing(), General_path_element_t::cost, pgrouting::tsp::Dmatrix::distance(), General_path_element_t::edge, pgrouting::tsp::Dmatrix::get_id(), pgrouting::tsp::Dmatrix::get_index(), pgrouting::tsp::TSP< MATRIX >::get_log(), pgrouting::tsp::TSP< MATRIX >::get_stats(), pgrouting::tsp::TSP< MATRIX >::get_tour(), pgrouting::tsp::TSP< MATRIX >::greedyInitial(), pgrouting::tsp::Dmatrix::has_id(), pgrouting::tsp::Dmatrix::has_no_infinity(), pgrouting::tsp::Dmatrix::is_symmetric(), General_path_element_t::node, pgassert, pgr_alloc(), pgr_free(), pgr_msg(), pgrouting::tsp::Dmatrix::set(), pgrouting::tsp::Dmatrix::size(), pgrouting::tsp::Dmatrix::tourCost(), and AssertFailedException::what().

Referenced by process().