svn commit: r820200 - /incubator/hama/trunk/src/test/org/apache/hama/matrix/TestSingularValueDecomposition.java

[ed...@apache.org]

Author: edwardyoon
Date: Wed Sep 30 07:31:23 2009
New Revision: 820200

URL: http://svn.apache.org/viewvc?rev=820200&view=rev
Log:
I've added a few comments to TestSingularValueDecomposition.class

Modified:
    incubator/hama/trunk/src/test/org/apache/hama/matrix/TestSingularValueDecomposition.java

Modified: incubator/hama/trunk/src/test/org/apache/hama/matrix/TestSingularValueDecomposition.java
URL: http://svn.apache.org/viewvc/incubator/hama/trunk/src/test/org/apache/hama/matrix/TestSingularValueDecomposition.java?rev=820200&r1=820199&r2=820200&view=diff
==============================================================================
--- incubator/hama/trunk/src/test/org/apache/hama/matrix/TestSingularValueDecomposition.java (original)
+++ incubator/hama/trunk/src/test/org/apache/hama/matrix/TestSingularValueDecomposition.java Wed Sep 30 07:31:23 2009
@@ -39,6 +39,7 @@
   private static HamaConfiguration conf;
 
   // Let's assume the A = [4 0; 3-5]
+  private static double[][] matrixA = {{4, 0}, {3, -5}};
   // A'A = [25 -15; -15 25]
   private static double[][] values = { { 25, -15 }, { -15, 25 } };
   // Then, eigenvalues of A'A are 10, 40 
@@ -56,7 +57,7 @@
         m1 = new DenseMatrix(conf, 2, 2);
         for (int i = 0; i < 2; i++)
           for (int j = 0; j < 2; j++)
-            m1.set(i, j, values[i][j]);
+            m1.set(i, j, matrixA[i][j]);
       }
 
       protected void tearDown() {
@@ -66,11 +67,20 @@
     return setup;
   }
 
-  public void testLog() throws IOException {
-    // Find the eigen/singular values and vectors of A'A    
-    m1.jacobiEigenValue(1);
-    HTable table = m1.getHTable();
+  public void testEigenSingularValues() throws IOException {
+    Matrix aT = m1.transpose();
+    DenseMatrix aTa = (DenseMatrix) aT.mult(m1);
     
+    for(int i = 0; i < m1.getRows(); i++) {
+      for(int j = 0; j < m1.getRows(); j++) {
+        assertEquals(aTa.get(i, j), values[i][j]);
+      }
+    }
+    
+    // Find the eigen/singular values and vectors of A'A
+    aTa.jacobiEigenValue(1);
+    HTable table = aTa.getHTable();
+
     for(int x=0; x<2; x++) {
       double eigenvalue = BytesUtil.bytesToDouble(table.get(BytesUtil.getRowIndex(x), 
           Bytes.toBytes(JacobiEigenValue.EIVAL)).getValue());
@@ -78,6 +88,9 @@
       assertTrue(Math.abs(Math.pow(eigenvalue, 0.5) - singularvalues[x]) < .0000001);
     }
     
+    //TODO: need to compute the inverse of S, S(-1)
+    //TODO: need to find out the V(T)
+    
     // Therefore, U= AVS'1=[-0.8944 -0.4472; 0.4472 -0.8944]
     // A = USV'=[4 0; 3 -5]
   }