Changeset 5648 for branches/QAP/HeuristicLab.Problems.QuadraticAssignment.Views/3.3/QAPView.cs

Timestamp:

03/10/11 01:49:10 (13 years ago)

Author:

abeham

Message:

#1330

• Unified QAP visualization in solution and problem view
• Fixed bug in gradient-descent gradient calculation when performing multidimensional scaling
• Extended QAPLIB parsers to cover some file format variations
• Added unit tests to check if all QAPLIB instances import without error
• Changed BestKnownSolution to be an OptionalValueParameter

Location:

branches/QAP

Files:

• . (modified) (1 prop)
• HeuristicLab.Problems.QuadraticAssignment.Views/3.3/QAPView.cs (modified) (14 diffs)

branches/QAP

@@ -1 +1 @@
 *.suo
+TestResults

@@ -1 +1 @@
 *.suo
+TestResults

branches/QAP/HeuristicLab.Problems.QuadraticAssignment.Views/3.3/QAPView.cs

@@ -35 +35 @@
     private Bitmap defaultBitmap;
 
+    #region Properties
     private DoubleMatrix distances;
     public DoubleMatrix Distances {
@@ -67 +68 @@
       }
     }
-
+    #endregion
+
+    #region Event Handling
     private void DeregisterDistancesEvents() {
       if (Distances != null) {
@@ -128 +131 @@
     }
 
+    private void pictureBox_SizeChanged(object sender, EventArgs e) {
+      SetupDefaultBitmap();
+      OnRedraw();
+    }
+    #endregion
+
     public QAPView() {
       InitializeComponent();
       SetupDefaultBitmap();
-      pictureBox.Resize += new EventHandler(pictureBox_Resize);
-    }
-
-    private void pictureBox_Resize(object sender, EventArgs e) {
-      SetupDefaultBitmap();
     }
 
@@ -141 +145 @@
       if (defaultBitmap != null) defaultBitmap.Dispose();
       defaultBitmap = new Bitmap(pictureBox.Width, pictureBox.Height);
-      using (Graphics g = Graphics.FromImage(defaultBitmap)) {
+      WriteCenteredTextToBitmap(ref defaultBitmap, "No visualization available");
+    }
+
+    private void WriteCenteredTextToBitmap(ref Bitmap bitmap, string text) {
+      using (Graphics g = Graphics.FromImage(bitmap)) {
         g.TextRenderingHint = TextRenderingHint.AntiAlias;
         g.SmoothingMode = SmoothingMode.AntiAlias;
 
-        string msg = "No visualization available";
-        SizeF strSize = g.MeasureString(msg, Font);
-        g.DrawString(msg, Font, Brushes.Black, (float)(pictureBox.Width - strSize.Width) / 2.0f, (float)(pictureBox.Height - strSize.Height) / 2.0f);
+        SizeF strSize = g.MeasureString(text, Font);
+        g.DrawString(text, Font, Brushes.Black, (float)(pictureBox.Width - strSize.Width) / 2.0f, (float)(pictureBox.Height - strSize.Height) / 2.0f);
       }
     }
@@ -161 +168 @@
     private void GenerateImage() {
       Bitmap newBitmap = null;
-      errorLabel.Text = "Error: -";
+      stressLabel.Text = "-";
       if (distancesRadioButton.Checked && Distances != null && Distances.Rows > 0
         && Distances.Rows == Distances.Columns) {
-        newBitmap = GenerateDistanceImage();
+        if (Distances.Rows > 30) {
+          newBitmap = new Bitmap(pictureBox.Width, pictureBox.Height);
+          WriteCenteredTextToBitmap(ref newBitmap, "Problem dimension too large for visualization.");
+        } else newBitmap = GenerateDistanceImage();
       } else if (weightsRadioButton.Checked && Weights != null && Weights.Rows > 0
         && Weights.Rows == Weights.Columns) {
-        newBitmap = GenerateWeightsImage();
+        if (Weights.Rows > 30) {
+          newBitmap = new Bitmap(pictureBox.Width, pictureBox.Height);
+          WriteCenteredTextToBitmap(ref newBitmap, "Problem dimension too large for visualization.");
+        } else newBitmap = GenerateWeightsImage();
       } else if (assignmentRadioButton.Checked
         && Assignment != null && Assignment.Length > 0
@@ -175 +188 @@
         && Distances.Rows == Distances.Columns
         && Assignment.Length == Weights.Rows
-        && Assignment.Length == Distances.Rows) {
+        && Assignment.Length == Distances.Rows
+        && Assignment.Validate()) {
         newBitmap = GenerateAssignmentImage();
       }
@@ -181 +195 @@
       pictureBox.Image = newBitmap != null ? newBitmap : defaultBitmap;
       if (bitmap != null) bitmap.Dispose();
-      bitmap = newBitmap;
-    }
-
+      if (newBitmap != null) bitmap = newBitmap;
+      else bitmap = null;
+    }
+
+    #region Draw distances
     private Bitmap GenerateDistanceImage() {
       if ((pictureBox.Width > 0) && (pictureBox.Height > 0)) {
         Bitmap newBitmap = new Bitmap(pictureBox.Width, pictureBox.Height);
 
-        double error;
-        DoubleMatrix coordinates = MultidimensionalScaling.Metric(distances, out error);
-        errorLabel.Text = "Error: " + error.ToString(CultureInfo.CurrentCulture.NumberFormat);
+        for (int i = 0; i < distances.Rows; i++) {
+          for (int j = i + 1; j < distances.Rows; j++)
+            if (distances[i, j] != distances[j, i]) {
+              WriteCenteredTextToBitmap(ref newBitmap, "Distance matrix is not symmetric");
+              return newBitmap;
+            }
+        }
+
+        double stress;
+        DoubleMatrix coordinates = MultidimensionalScaling.MetricByDistance(distances, out stress);
+        stressLabel.Text = stress.ToString("0.00", CultureInfo.CurrentCulture.NumberFormat);
         double xMin = double.MaxValue, yMin = double.MaxValue, xMax = double.MinValue, yMax = double.MinValue;
         double maxDistance = double.MinValue;
@@ -214 +238 @@
                                 newBitmap.Height - (border + ((int)((coordinates[i, 1] - yMin) * yStep))));
 
+        Random rand = new Random();
         using (Graphics graphics = Graphics.FromImage(newBitmap)) {
           graphics.SmoothingMode = SmoothingMode.AntiAlias;
           graphics.TextRenderingHint = TextRenderingHint.AntiAlias;
-          graphics.DrawString("Showing Locations spaced out according to their distances", Font, Brushes.Black, 5, 2);
+          graphics.DrawString("Showing locations spaced out according to their distances", Font, Brushes.Black, 5, 2);
 
           for (int i = 0; i < coordinates.Rows - 1; i++) {
@@ -233 +258 @@
                   caption = distances[i, j].ToString(CultureInfo.InvariantCulture.NumberFormat);
               }
+              if (!String.IsNullOrEmpty(caption)) {
+                double r = rand.NextDouble();
+                while (r < 0.2 || r > 0.8) r = rand.NextDouble();
+                float x = (float)(start.X + (end.X - start.X) * r + 5);
+                float y = (float)(start.Y + (end.Y - start.Y) * r + 5);
+                graphics.DrawString(caption, Font, Brushes.Black, x, y);
+              }
             }
           }
@@ -246 +278 @@
       return null;
     }
-
+    #endregion
+
+    #region Draw weights
     private Bitmap GenerateWeightsImage() {
       if ((pictureBox.Width > 0) && (pictureBox.Height > 0)) {
         Bitmap newBitmap = new Bitmap(pictureBox.Width, pictureBox.Height);
 
-        double error;
-        DoubleMatrix coordinates = MultidimensionalScaling.Classic(weights, out error);
-        errorLabel.Text = "Error: " + error.ToString(CultureInfo.CurrentCulture.NumberFormat);
+        double maxWeight = double.MinValue;
+        for (int i = 0; i < weights.Rows; i++)
+          for (int j = i + 1; j < weights.Rows; j++) {
+            if (weights[i, j] > maxWeight)
+              maxWeight = weights[i, j] + weights[j, i];
+
+            if (weights[i, j] != weights[j, i]) {
+              WriteCenteredTextToBitmap(ref newBitmap, "Weights matrix is not symmetric");
+              return newBitmap;
+            }
+          }
+
+        DoubleMatrix distances = new DoubleMatrix(weights.Rows, weights.Columns);
+        for (int i = 0; i < distances.Rows; i++)
+          for (int j = 0; j < distances.Columns; j++) {
+            distances[i, j] = maxWeight + 1 - weights[i, j];
+          }
+
+        double stress;
+        DoubleMatrix coordinates = MultidimensionalScaling.MetricByDistance(distances, out stress);
+        stressLabel.Text = stress.ToString("0.00", CultureInfo.CurrentCulture.NumberFormat);
+        double xMin = double.MaxValue, yMin = double.MaxValue, xMax = double.MinValue, yMax = double.MinValue;
+        for (int i = 0; i < coordinates.Rows; i++) {
+          if (xMin > coordinates[i, 0]) xMin = coordinates[i, 0];
+          if (yMin > coordinates[i, 1]) yMin = coordinates[i, 1];
+          if (xMax < coordinates[i, 0]) xMax = coordinates[i, 0];
+          if (yMax < coordinates[i, 1]) yMax = coordinates[i, 1];
+        }
+
+        int border = 20;
+        double xStep = xMax != xMin ? (pictureBox.Width - 2 * border) / (xMax - xMin) : 1;
+        double yStep = yMax != yMin ? (pictureBox.Height - 2 * border) / (yMax - yMin) : 1;
+
+        Point[] points = new Point[coordinates.Rows];
+        for (int i = 0; i < coordinates.Rows; i++)
+          points[i] = new Point(border + ((int)((coordinates[i, 0] - xMin) * xStep)),
+                                newBitmap.Height - (border + ((int)((coordinates[i, 1] - yMin) * yStep))));
+
+        Random rand = new Random();
+        using (Graphics graphics = Graphics.FromImage(newBitmap)) {
+          graphics.SmoothingMode = SmoothingMode.AntiAlias;
+          graphics.TextRenderingHint = TextRenderingHint.AntiAlias;
+          graphics.DrawString("Showing facilities spaced out according to their weights", Font, Brushes.Black, 5, 2);
+
+          for (int i = 0; i < coordinates.Rows - 1; i++) {
+            for (int j = i + 1; j < coordinates.Rows; j++) {
+              Point start = points[i], end = points[j];
+              string caption = String.Empty;
+              double d = Math.Max(distances[i, j], distances[j, i]);
+              double w = weights[i, j];
+              if (w > 0) {
+                float width = (float)Math.Ceiling(3.0 * w / maxWeight);
+                graphics.DrawLine(new Pen(Color.MediumBlue, width), start, end);
+                caption = w.ToString(CultureInfo.InvariantCulture.NumberFormat);
+              }
+              if (!String.IsNullOrEmpty(caption)) {
+                double r = rand.NextDouble();
+                while (r < 0.2 || r > 0.8) r = rand.NextDouble();
+                float x = (float)(start.X + (end.X - start.X) * r + 5);
+                float y = (float)(start.Y + (end.Y - start.Y) * r + 5);
+                graphics.DrawString(caption, Font, Brushes.Black, x, y);
+              }
+            }
+          }
+          for (int i = 0; i < points.Length; i++) {
+            Point p = new Point(points[i].X - 3, points[i].Y - 3);
+            graphics.FillRectangle(Brushes.Black, p.X, p.Y, 8, 8);
+            graphics.DrawString(i.ToString(), Font, Brushes.Black, p.X, p.Y + 10);
+          }
+        }
+        return newBitmap;
+      }
+      return null;
+    }
+    #endregion
+
+    #region Draw assignment
+    private Bitmap GenerateAssignmentImage() {
+      if ((pictureBox.Width > 0) && (pictureBox.Height > 0)) {
+        Bitmap newBitmap = new Bitmap(pictureBox.Width, pictureBox.Height);
+
+        for (int i = 0; i < distances.Rows; i++) {
+          for (int j = i + 1; j < distances.Rows; j++) {
+            if (distances[i, j] != distances[j, i]) {
+              WriteCenteredTextToBitmap(ref newBitmap, "Distance matrix is not symmetric");
+              return newBitmap;
+            }
+            if (weights[i, j] != weights[j, i]) {
+              WriteCenteredTextToBitmap(ref newBitmap, "Weights matrix is not symmetric");
+            }
+          }
+        }
+
+        double stress;
+        DoubleMatrix coordinates = MultidimensionalScaling.MetricByDistance(distances, out stress);
+        stressLabel.Text = stress.ToString("0.00", CultureInfo.CurrentCulture.NumberFormat);
         double xMin = double.MaxValue, yMin = double.MaxValue, xMax = double.MinValue, yMax = double.MinValue;
         double maxWeight = double.MinValue;
@@ -263 +390 @@
 
           for (int j = i + 1; j < coordinates.Rows; j++) {
-            if (weights[i, j] + weights[j, i] > maxWeight)
-              maxWeight = weights[i, j] + weights[j, i];
+            if (weights[i, j] > maxWeight) maxWeight = weights[i, j];
           }
         }
@@ -279 +405 @@
         Random rand = new Random();
         using (Graphics graphics = Graphics.FromImage(newBitmap)) {
-          graphics.SmoothingMode = SmoothingMode.AntiAlias;
-          graphics.TextRenderingHint = TextRenderingHint.AntiAlias;
-          graphics.DrawString("Showing Facilities spaced out according to their weights", Font, Brushes.Black, 5, 2);
-
-          for (int i = 0; i < coordinates.Rows - 1; i++) {
-            for (int j = i + 1; j < coordinates.Rows; j++) {
-              Point start = points[i], end = points[j];
+          graphics.SmoothingMode = System.Drawing.Drawing2D.SmoothingMode.AntiAlias;
+          graphics.TextRenderingHint = System.Drawing.Text.TextRenderingHint.AntiAlias;
+
+          for (int i = 0; i < assignment.Length - 1; i++) {
+            for (int j = i + 1; j < assignment.Length; j++) {
+              Point start, end;
               string caption = String.Empty;
-              double d = Math.Max(distances[i, j], distances[j, i]);
-              double w = weights[i, j] + weights[j, i];
+              double d = distances[i, j];
+              start = points[assignment[i]];
+              end = points[assignment[j]];
+              double w = weights[i, j];
               if (w > 0) {
-                float width = (float)Math.Ceiling(3.0 * w / maxWeight);
+                float width = (float)Math.Ceiling(4.0 * w / maxWeight);
                 graphics.DrawLine(new Pen(Color.MediumBlue, width), start, end);
                 caption = w.ToString(CultureInfo.InvariantCulture.NumberFormat);
@@ -303 +430 @@
             }
           }
+
           for (int i = 0; i < points.Length; i++) {
             Point p = new Point(points[i].X - 3, points[i].Y - 3);
@@ -313 +441 @@
       return null;
     }
-
-    private Bitmap GenerateAssignmentImage() {
-      return null;
-    }
+    #endregion
 
     private void CustomDispose(bool disposing) {

@@ -1 +1 @@
 *.suo
+TestResults