8 Commits

Author SHA1 Message Date
Jan Potocki b2517f51f2 Simplified if-else statement (TS) 2020-01-20 19:40:22 +01:00
Jan Potocki 20160689f4 Experimental: constraint for diversification (max 150% of current optimum length) 2020-01-20 19:33:40 +01:00
Jan Potocki 6cab940785 Fixed newline in output (TS) 2020-01-20 19:24:40 +01:00
Jan Potocki c10d296d82 Fixed TS crash when entire neighborhood is on tabu list 2020-01-20 19:23:47 +01:00
Jan Potocki 5d78916930 Added displaying timestamp for new optimal router 2020-01-15 19:43:00 +01:00
Jan Potocki b577fd6c13 New approach for global optimum sharing 2020-01-15 19:25:38 +01:00
Jan Potocki 74e1d0be8e Global shared optimum - almost ready 2020-01-15 19:04:22 +01:00
Jan Potocki 420ee4f275 Experimental: initial shared global optimum in TS 2020-01-15 18:22:38 +01:00
2 changed files with 27 additions and 7 deletions
+25 -5
View File
@@ -5,7 +5,6 @@
#include <queue> #include <queue>
#include <random> #include <random>
#include <thread> #include <thread>
#include <iostream> #include <iostream>
Graph::Graph() Graph::Graph()
@@ -426,7 +425,7 @@ std::vector<unsigned> Graph::travellingSalesmanTabuSearch(Graph &graph, unsigned
{ {
// ALGORYTM wielawotkowy oparty na metaheurystyce tabu search // ALGORYTM wielawotkowy oparty na metaheurystyce tabu search
// Pomocniczy kod uruchamiajacy watki wlasciwego algorytmu w najbardziej optymalny sposob // Pomocniczy kod uruchamiajacy watki wlasciwego algorytmu w najbardziej optymalny sposob
// Implementacja: Jan Potocki 2019-2020 // Implementacja: Jan Potocki 2019
std::vector<unsigned> startVertexVector; std::vector<unsigned> startVertexVector;
std::vector<std::thread> threadsVector; std::vector<std::thread> threadsVector;
@@ -492,13 +491,15 @@ void Graph::travellingSalesmanTabuSearchEngine(Graph &graph, unsigned tabuSteps,
// ALGORYTM oparty na metaheurystyce tabu search z dywersyfikacja i sasiedztwem typu swap // ALGORYTM oparty na metaheurystyce tabu search z dywersyfikacja i sasiedztwem typu swap
// Rdzen przeznaczony do uruchamiania jako jeden watek // Rdzen przeznaczony do uruchamiania jako jeden watek
// Projekt i implementacja: Jan Potocki 2017 // Projekt i implementacja: Jan Potocki 2017
// (refactoring 2019-2020) // (refactoring 2019)
Stopwatch onboardClock; Stopwatch onboardClock;
std::vector<unsigned> optimalRoute; // Tu bedziemy zapisywac optymalne (w danej chwili) rozwiazanie std::vector<unsigned> optimalRoute; // Tu bedziemy zapisywac optymalne (w danej chwili) rozwiazanie
int optimalRouteLength = -1; // -1 - bedziemy odtad uznawac, ze to jest nieskonczonosc ;-) int optimalRouteLength = -1; // -1 - bedziemy odtad uznawac, ze to jest nieskonczonosc ;-)
std::vector<unsigned> currentRoute; // Rozpatrywane rozwiazanie std::vector<unsigned> currentRoute; // Rozpatrywane rozwiazanie
// Wyznaczenie poczatkowego rozwiazania algorytmem zachlannym
//currentRoute = Graph::travellingSalesmanGreedy(graph);
currentRoute = startRoute; currentRoute = startRoute;
// Inicjalizacja glownej petli... // Inicjalizacja glownej petli...
@@ -565,6 +566,7 @@ void Graph::travellingSalesmanTabuSearchEngine(Graph &graph, unsigned tabuSteps,
// ...jezeli niespelnione - pomijamy ruch // ...jezeli niespelnione - pomijamy ruch
continue; continue;
if(nextRouteLength == -1 || nextRouteLength > neighbourRouteLength) if(nextRouteLength == -1 || nextRouteLength > neighbourRouteLength)
{ {
nextRouteLength = neighbourRouteLength; nextRouteLength = neighbourRouteLength;
@@ -599,7 +601,7 @@ void Graph::travellingSalesmanTabuSearchEngine(Graph &graph, unsigned tabuSteps,
stopCounter = 0; stopCounter = 0;
} }
// Synchronizacja globalnie najlepszej trasy // Synchronizacja globalnie najlepszej trasy (1)
globalOptimumMutex.lock(); globalOptimumMutex.lock();
if(globalOptimumLength == -1 || globalOptimumLength > nextRouteLength) if(globalOptimumLength == -1 || globalOptimumLength > nextRouteLength)
{ {
@@ -660,16 +662,34 @@ void Graph::travellingSalesmanTabuSearchEngine(Graph &graph, unsigned tabuSteps,
// Intensyfikacja przeszukiwania przez skrócenie kadencji // Intensyfikacja przeszukiwania przez skrócenie kadencji
// (jezeli w ostatnim przebiegu znaleziono nowe minimum) // (jezeli w ostatnim przebiegu znaleziono nowe minimum)
currentRoute = optimalRoute; currentRoute = optimalRoute;
currentTabuSteps = tabuSteps; /// 4; currentTabuSteps = tabuSteps / 4;
intensification = false; intensification = false;
// PEA 2 Plus // PEA 2 Plus
// Jan Potocki 2019 // Jan Potocki 2019
} }
else else
{ {
// Synchronizacja globalnie najlepszej trasy (2)
globalOptimumMutex.lock();
optimalRouteLength = globalOptimumLength;
optimalRoute = globalOptimum;
globalOptimumMutex.unlock();
// W innym przypadku wlasciwa dywersyfikacja przez wygenerowanie nowego // W innym przypadku wlasciwa dywersyfikacja przez wygenerowanie nowego
// rozwiazania startowego algorytmem hybrydowym losowo-zachlannym // rozwiazania startowego algorytmem hybrydowym losowo-zachlannym
// (ale nie znacznie gorszego niz juz znalezione)
int criticalLength = (15 * optimalRouteLength) / 10;
int currentRouteLength;
do
{
currentRoute = Graph::travellingSalesmanHybrid(graph); currentRoute = Graph::travellingSalesmanHybrid(graph);
currentRouteLength = 0;
for(int i = 1; i < currentRoute.size(); i++)
currentRouteLength += graph.getWeight(currentRoute.at(i - 1), currentRoute.at(i));
} while(currentRouteLength > criticalLength);
currentTabuSteps = tabuSteps; currentTabuSteps = tabuSteps;
intensification = false; intensification = false;
} }
+1 -1
View File
@@ -439,7 +439,6 @@ int main(int argc, char *argv[])
// Wyswietlenie trasy // Wyswietlenie trasy
unsigned distFromStart = 0; unsigned distFromStart = 0;
unsigned length = 0; unsigned length = 0;
cout << endl;
cout << route.at(0) << '\t' << length << '\t' << distFromStart << endl; cout << route.at(0) << '\t' << length << '\t' << distFromStart << endl;
for(int i = 1; i < route.size(); i++) for(int i = 1; i < route.size(); i++)
{ {
@@ -449,6 +448,7 @@ int main(int argc, char *argv[])
cout << route.at(i) << '\t' << length << '\t' << distFromStart << endl; cout << route.at(i) << '\t' << length << '\t' << distFromStart << endl;
} }
cout << endl;
cout << "Dlugosc trasy: " << distFromStart << endl; cout << "Dlugosc trasy: " << distFromStart << endl;
cout << endl; cout << endl;
cout << "Czas wykonania algorytmu [s]: " << clock.read() << endl; cout << "Czas wykonania algorytmu [s]: " << clock.read() << endl;