Changeset 9291f85 for molecuilder/src/unittests/gslmatrixsymmetricunittest.cpp

Timestamp:

Jan 10, 2010, 7:25:50 PM (16 years ago)

Author:

Frederik Heber <heber@…>

Children:

9cdbee

Parents:

7b991d

Message:

Added Determinant() function to class GSLMatrix.

• new function GSLMatrix::Determinant() which calculates the determinant for square matrices by LU decomposition.

Signed-off-by: Frederik Heber <heber@tabletINS.(none)>

File:

• molecuilder/src/unittests/gslmatrixsymmetricunittest.cpp (modified) (1 diff)

molecuilder/src/unittests/gslmatrixsymmetricunittest.cpp

@@ -247 +247 @@
   m->SetFromDoubleArray(array);
   CPPUNIT_ASSERT_EQUAL( true, m->IsPositiveDefinite() );
+
+  //determinant
+  m->SetIdentity();
+  CPPUNIT_ASSERT_EQUAL( 1., m->Determinant() );
+
+  m->SetZero();
+  CPPUNIT_ASSERT_EQUAL( 0., m->Determinant() );
+
+  m->Set( 0, 0, 1.);
+  m->Set( 1, 1, 1.);
+  m->Set( 2, 1, 1.);
+  CPPUNIT_ASSERT_EQUAL( 0., m->Determinant() );
+  
+  double array2[] = { 2., 0., 1.,
+                     -3., 1., 1.,
+                     1., 5.5, 1. };
+  m->SetFromDoubleArray(array2);
+  CPPUNIT_ASSERT_EQUAL( -26.5, m->Determinant() );
 };