Changeset 4566 in orxonox.OLD for orxonox/trunk/src/util/newmat/eigen.cpp

Timestamp:

Jun 9, 2005, 5:40:47 PM (19 years ago)

Author:

patrick

Message:

orxonox/trunk: made some tests with the newmat lib, now it should be no problem to solve all math problems in cd

File:

• orxonox/trunk/src/util/newmat/eigen.cpp (modified) (2 diffs)

orxonox/trunk/src/util/newmat/eigen.cpp

@@ -35 +35 @@
 int main(int argc, char **argv) {
 
+  cout << "start test: composing matrix" << endl;
+  
+  SymmetricMatrix C(3);
+
+  C(1,1) = 1;
+  C(1,2) = 4;
+  C(1,3) = 4;
+  C(2,1) = 4;  
+  C(2,2) = 2;
+  C(2,3) = 4;
+  C(3,1) = 4;
+  C(3,2) = 4;
+  C(3,3) = 3;
+  
+  cout << "matrix composed" << endl;
+
+
+  //SymmetricMatrix X;
+  //X << C;
+
+  cout << "The symmetrix matrix C" << endl;
+  cout << setw(5) << setprecision(0) << C << endl;
+
+  Matrix                        V(3,3); // for eigenvectors
+  DiagonalMatrix        D(3);   // for eigenvalues
+  
+  // the decomposition
+  Jacobi(C, D, V);
+
+
+  // Print the result
+  cout << "The eigenvalues matrix:" << endl;
+  cout << setw(10) << setprecision(5) << D << endl;
+  cout << "The eigenvectors matrix:" << endl;
+  cout << setw(10) << setprecision(5) << V << endl;
+
+  return 0;
+  /*
   int M = 3, N = 5;
   Matrix X(M,N); // Define an M x N general matrix
@@ -103 +141 @@
   cout << "solution for right hand side r2" << endl;
   cout << setw(10) << setprecision(5) << s2 << endl;
-        
+  */
+
   return 0;
 }