Changeset 7315 for branches/FitnessLandscapeAnalysis/VRPProblemAnalyzer/Program.cs

Timestamp:

01/11/12 13:05:52 (12 years ago)

Author:

epitzer

Message:

#1696: Add VRP depot excentricity metrics

File:

• branches/FitnessLandscapeAnalysis/VRPProblemAnalyzer/Program.cs (modified) (6 diffs)

branches/FitnessLandscapeAnalysis/VRPProblemAnalyzer/Program.cs

@@ -1 +1 @@
 using System;
 using System.Collections.Generic;
+using System.Diagnostics;
 using System.Linq;
 using System.Text;
@@ -6 +7 @@
 
 namespace VRPProblemAnalyzer {
-  public static class Maths {
-    public static double StandardDeviation(this IEnumerable<double> data) {
-      double average = data.Average();
-      var individualDeviations = data.Select(num => Math.Pow(num - average, 2.0));
-      return Math.Sqrt(individualDeviations.Average());
-    }
-  }
   
-  class Program {
+  public class Program {
 
     private static double GetDistance(double[,] vertices, int source, int target) {
@@ -187 +181 @@
     }
 
+
     private static void Normalize(TSPLIBParser instance) {
       //normalize demands
@@ -195 +190 @@
       //normalize coordinates
       //-find bounds
-      double maxX = double.MinValue;
-      double maxY = double.MinValue;
-
-      int cities = instance.Vertices.Length / 2;
-      for (int i = 0; i < cities; i++) {
-        if (instance.Vertices[i, 0] > maxX)
-          maxX = instance.Vertices[i, 0];
-
-        if (instance.Vertices[i, 1] > maxY)
-          maxY = instance.Vertices[i, 1];
-      }
-
-      //-normalize
-      double factorX = 1;
-      double factorY = 1;
-      if (maxX > maxY) {
-        factorX = 1 / maxX;
-        factorY = (maxY / maxX) / maxY;
-      } else {
-        factorX = (maxX / maxY) / maxX;
-        factorY = 1 / maxY;
-      }
-
-      for (int i = 0; i < cities; i++) {
-        instance.Vertices[i, 0] *= factorX;
-        instance.Vertices[i, 1] *= factorY;
+      var cities = Utils.MatrixToPointList(instance.Vertices);
+      var minX = cities.Min(c => c.X);
+      var minY = cities.Min(c => c.Y);
+      var maxX = cities.Max(c => c.X);
+      var maxY = cities.Max(c => c.Y);
+      var rangeX = maxX - minX;
+      var rangeY = maxY - minY;
+      var factor = Math.Sqrt(0.5)/Math.Max(rangeX, rangeY);
+
+
+      for (int i = 0; i < instance.Vertices.GetLength(0); i++) {
+        instance.Vertices[i, 0] = (cities[i].X-minX)*factor;
+        instance.Vertices[i, 1] = (cities[i].Y-minY)*factor;
       }
     }
@@ -231 +213 @@
       var path = args[0];
       using (StreamWriter sw = new StreamWriter(Path.Combine(path, "analysis.csv"))) {
-        sw.WriteLine("Instance;Customers;Clustering;DistanceAvg;DistanceStdev;DemandAvg;DemandStdev");
+        sw.WriteLine("Instance;Customers;Clustering;DistanceAvg;DistanceStdev;DemandAvg;DemandStdev;GeographicExcentricity;DistanceExcentricity;DistanceDemandExcenetricity");
         
         string[] instances = Directory.GetFiles(path, "*.vrp", SearchOption.AllDirectories);
@@ -253 +235 @@
           sw.Write(GetDemandHeterogenity(parser));
           sw.Write(';');
+          sw.Write(DepotExcentricityCalculator.Geographic(parser.Vertices));
+          sw.Write(';');
+          sw.Write(DepotExcentricityCalculator.DistanceCentroid(parser.Vertices));
+          sw.Write(';');
+          sw.Write(DepotExcentricityCalculator.DemandDistanceCentroid(parser.Vertices, parser.Demands));
+          sw.Write(';');
           sw.WriteLine();
         }
       }
+      Console.WriteLine("Done. Press Enter...");
+      Console.ReadLine();
     }
   }