Changeset 9291f85 for molecuilder/src/gslmatrix.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/gslmatrix.cpp (modified) (1 diff)

molecuilder/src/gslmatrix.cpp

@@ -229 +229 @@
     return (gsl_linalg_cholesky_decomp (matrix) != GSL_EDOM);
 };
+
+
+/** Calculates the determinant of the matrix.
+ * if matrix is square, uses LU decomposition.
+ */
+double GSLMatrix::Determinant()
+{
+  int signum = 0;
+  assert (rows == columns && "Determinant can only be calculated for square matrices.");
+  gsl_permutation *p = gsl_permutation_alloc(rows);
+  gsl_linalg_LU_decomp(matrix, p, &signum);
+  gsl_permutation_free(p);
+  return gsl_linalg_LU_det(matrix, signum); 
+};
+