private AssignmentTwoResultList <GraphGenome> MultiStartLocalSearch(INeighborhood <GraphGenome> neighborhood)
{
//we moeten nog iets met ExperimentAmount doen...
//for(int i = 0; i < ExperimentAmount; ++i)
AssignmentTwoResultList <GraphGenome> results = new AssignmentTwoResultList <GraphGenome>("$MLS_{" + neighborhood.Name + "}$");
for (int j = 0; j < threads.Length; j++)
{
elapsedMilisecondsThreads[j] = new List <AssignmentTwoResults <GraphGenome> >();
int k = j;
int seed = main_random_source.Next();
int amount = OptimaAmount / threads.Length;
Thread T = new Thread(() => searchThread(k, amount, seed, neighborhood));
threads[j] = T;
}
Stopwatch sw = new Stopwatch();
sw.Start();
foreach (Thread T in threads)
{
T.Start();
}
while (this.threadsAreRunning())
{
;
}
sw.Stop();
GraphGenome optimum = null;
for (int j = 0; j < threads.Length; j++)
{
if (optimum == null || optimum.Fitness > solutionsThreads[j].Fitness)
{
optimum = solutionsThreads[j];
}
}
//List<AssignmentTwoResults> results = new List<AssignmentTwoResults>();
foreach (List <AssignmentTwoResults <GraphGenome> > L in elapsedMilisecondsThreads)
{
results.AddRange(L);
}
//Console.WriteLine("Best solution found: {0} in an average of {1} ms", optimum.Fitness, cpu_ticks.Average());
Console.WriteLine("Best solution found: {0} in an average of ... ms", optimum.Fitness); //, results.Sum() / OptimaAmount);
Console.WriteLine("Total time: {0} sec", sw.ElapsedMilliseconds / 1000f);
optimum.ToImage(String.Format("results/MLS[{0}]-{1}.bmp", neighborhood.Name, optimum.Fitness), 3000, 3000);
return(results);
}
private AssignmentTwoResultList <GraphGenome> IteratedLocalSearch(INeighborhood <GraphGenome> neighborhood, bool silent = false)
{
LocalSearch <GraphGenome> local_search =
new LocalSearch <GraphGenome>(500,
neighborhood,
new GraphComparer <GraphGenome>(), main_random_source);
IMutation <GraphGenome> pertubation = new ILSPertubation <GraphGenome>(main_random_source);
GraphGenome graph = new GraphGenome(500);
graph.CreateGraph("Graph500.txt");
GraphGenome optimum = null;
List <long> elapsedMilisecondsList = new List <long>();
AssignmentTwoResultList <GraphGenome> results = new AssignmentTwoResultList <GraphGenome>("$ILS_{" + neighborhood.Name + "}$");
for (int i = 0; i < ExperimentAmount; ++i)
{
AssignmentTwoResults <GraphGenome> res = new AssignmentTwoResults <GraphGenome>();
for (int j = 0; j < OptimaAmount; j++)
{
if (optimum == null)
{
graph = new GraphGenome(500);
graph.Generate(ref main_random_source);
}
else
{
graph = pertubation.Mutate(optimum);
}
Stopwatch sw = Stopwatch.StartNew();
GraphGenome inner_optimum = local_search.Search(graph);
sw.Stop();
elapsedMilisecondsList.Add(sw.ElapsedMilliseconds);
res.Add(inner_optimum, sw.ElapsedTicks);
if (!silent)
{
Console.WriteLine("Found {0} in {1} ticks or {2}ms. ({3}/{4})", inner_optimum.Fitness, sw.ElapsedTicks,
sw.ElapsedMilliseconds, j, OptimaAmount);
}
if (optimum == null || optimum.Fitness >= inner_optimum.Fitness)
{
optimum = inner_optimum;
}
optimum.FunctionEvaluations = 0;
}
results.Add(res);
}
//Console.WriteLine("Best solution found: {0} in an average of {1} ms", optimum.Fitness, cpu_ticks.Average());
Console.WriteLine("Best solution found: {0} in an average of {1} ms", optimum.Fitness, elapsedMilisecondsList.Sum() / OptimaAmount);
//Console.WriteLine("Total time: {0} sec", sw.ElapsedMilliseconds / 1000f);
optimum.ToImage(String.Format("results/ILS[{0}]-{1}.bmp", neighborhood.Name, optimum.Fitness), 3000, 3000);
//Console.ReadLine();
return(results);
}
