From 031f5a94196cd177467b6ca84ce308dabc8c24bf Mon Sep 17 00:00:00 2001 From: Pjotr Prins Date: Thu, 7 Dec 2017 10:17:25 +0000 Subject: More notes and removed unused LULndet --- src/lapack.cpp | 21 --------------------- 1 file changed, 21 deletions(-) (limited to 'src/lapack.cpp') diff --git a/src/lapack.cpp b/src/lapack.cpp index 26333b6..d15446b 100644 --- a/src/lapack.cpp +++ b/src/lapack.cpp @@ -306,27 +306,6 @@ double LULndet(const gsl_matrix *LU) { return gsl_linalg_LU_lndet((gsl_matrix *)LU); } -/* -double LULndet(gsl_matrix_float *LU) { - gsl_matrix *LU_double = gsl_matrix_alloc(LU->size1, LU->size2); - double d; - - // Copy float matrix to double. - for (size_t i = 0; i < LU->size1; i++) { - for (size_t j = 0; j < LU->size2; j++) { - gsl_matrix_set(LU_double, i, j, gsl_matrix_float_get(LU, i, j)); - } - } - - // LU decomposition. - d = gsl_linalg_LU_lndet(LU_double); - - // Free matrix - gsl_matrix_free(LU_double); - return d; -} -*/ - // LU solve. void LUSolve(const gsl_matrix *LU, const gsl_permutation *p, const gsl_vector *b, gsl_vector *x) { -- cgit v1.2.3