Changeset 438

Timestamp:

07/28/09 15:53:00 (15 years ago)

Author:

vbarta

Message:

testing matrix inversion

Location:

library/tests

Files:

• square_mat_stress.cpp (modified) (1 diff)
• square_mat_test.cpp (modified) (1 diff)

library/tests/square_mat_stress.cpp

@@ -72 +72 @@
         cout << "opupdt: " << elapsed << " s" << endl;
 
+        TMatrix invmat(sz);
+
+        tt.tic();
+        sqmat.inv(invmat);
+        elapsed = tt.toc();
+
+        if (!fast) {
+            mat invA = inv(A);
+            CHECK_CLOSE(invA, invmat.to_mat(), epsilon);
+        }
+
+        cout << "inv: " << elapsed << " s" << endl;
+
         sz *= 7;
     }

library/tests/square_mat_test.cpp

@@ -45 +45 @@
     res = A + w * outer_product(v, v);
     CHECK_CLOSE(res, sqmat2.to_mat(), epsilon);
+
+    TMatrix invmat(sz);
+    sqmat.inv(invmat);
+    mat invA = inv(A);
+    CHECK_CLOSE(invA, invmat.to_mat(), epsilon);
 
     double d = det(A);