Changeset 5871 for branches/QAP/HeuristicLab.Problems.QuadraticAssignment.Views

Timestamp:

03/29/11 16:54:52 (13 years ago)

Author:

abeham

Message:

#1330

• Fixed normalized stress value
• Changed method name and signature of MDS
• Adapted unit tests
• Adapted QAP View

File:

• branches/QAP/HeuristicLab.Problems.QuadraticAssignment.Views/3.3/QAPView.cs (modified) (4 diffs)

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

@@ -216 +216 @@
         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;
-            }
-        }
-
-        double stress;
-        DoubleMatrix coordinates = MultidimensionalScaling.MetricByDistance(distances, out stress);
-        stressLabel.Text = stress.ToString("0.00", CultureInfo.CurrentCulture.NumberFormat);
+        DoubleMatrix coordinates;
+        double stress = double.NaN;
+        try {
+          coordinates = MultidimensionalScaling.KruskalShepard(distances);
+          stress = MultidimensionalScaling.CalculateNormalizedStress(distances, coordinates);
+          stressLabel.Text = stress.ToString("0.00", CultureInfo.CurrentCulture.NumberFormat);
+        } catch {
+          WriteCenteredTextToBitmap(ref newBitmap, "Distance matrix is not symmetric");
+          stressLabel.Text = "-";
+          return newBitmap;
+        }
         double xMin = double.MaxValue, yMin = double.MaxValue, xMax = double.MinValue, yMax = double.MinValue;
         double maxDistance = double.MinValue;
@@ -302 +302 @@
             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;
-            }
           }
 
@@ -315 +310 @@
           }
 
-        double stress;
-        DoubleMatrix coordinates = MultidimensionalScaling.MetricByDistance(distances, out stress);
-        stressLabel.Text = stress.ToString("0.00", CultureInfo.CurrentCulture.NumberFormat);
+        DoubleMatrix coordinates;
+        double stress = double.NaN;
+        try {
+          coordinates = MultidimensionalScaling.KruskalShepard(distances);
+          stress = MultidimensionalScaling.CalculateNormalizedStress(distances, coordinates);
+          stressLabel.Text = stress.ToString("0.00", CultureInfo.CurrentCulture.NumberFormat);
+        } catch {
+          WriteCenteredTextToBitmap(ref newBitmap, "Weights matrix is not symmetric");
+          stressLabel.Text = "-";
+          return newBitmap;
+        }
         double xMin = double.MaxValue, yMin = double.MaxValue, xMax = double.MinValue, yMax = double.MinValue;
         for (int i = 0; i < coordinates.Rows; i++) {
@@ -382 +385 @@
             if (distances[i, j] != distances[j, i]) {
               WriteCenteredTextToBitmap(ref newBitmap, "Distance matrix is not symmetric");
+              stressLabel.Text = "-";
               return newBitmap;
             }
             if (weights[i, j] != weights[j, i]) {
               WriteCenteredTextToBitmap(ref newBitmap, "Weights matrix is not symmetric");
+              stressLabel.Text = "-";
+              return newBitmap;
             }
           }
         }
 
-        double stress;
-        DoubleMatrix coordinates = MultidimensionalScaling.MetricByDistance(distances, out stress);
-        stressLabel.Text = stress.ToString("0.00", CultureInfo.CurrentCulture.NumberFormat);
+        DoubleMatrix coordinates = null;
+        double stress = double.NaN;
+        try {
+          coordinates = MultidimensionalScaling.KruskalShepard(distances);
+          stress = MultidimensionalScaling.CalculateNormalizedStress(distances, coordinates);
+          stressLabel.Text = stress.ToString("0.00", CultureInfo.CurrentCulture.NumberFormat);
+        } catch {
+          WriteCenteredTextToBitmap(ref newBitmap, "Unknown error");
+          stressLabel.Text = "-";
+          return newBitmap;
+        }
+
         double xMin = double.MaxValue, yMin = double.MaxValue, xMax = double.MinValue, yMax = double.MinValue;
         double maxWeight = double.MinValue;