aboutsummaryrefslogtreecommitdiff
path: root/src/lapack.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'src/lapack.cpp')
-rw-r--r--src/lapack.cpp21
1 files changed, 0 insertions, 21 deletions
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) {