diff options
Diffstat (limited to 'bench/btl/libs/eigen2/eigen2_interface.hh')
-rw-r--r-- | bench/btl/libs/eigen2/eigen2_interface.hh | 6 |
1 files changed, 5 insertions, 1 deletions
diff --git a/bench/btl/libs/eigen2/eigen2_interface.hh b/bench/btl/libs/eigen2/eigen2_interface.hh index f3ba5e8af..39eaa3b65 100644 --- a/bench/btl/libs/eigen2/eigen2_interface.hh +++ b/bench/btl/libs/eigen2/eigen2_interface.hh @@ -17,7 +17,7 @@ // #ifndef EIGEN2_INTERFACE_HH #define EIGEN2_INTERFACE_HH -#include <cblas.h> +// #include <cblas.h> #include <Eigen/Core> #include <Eigen/Cholesky> #include <Eigen/LU> @@ -137,6 +137,10 @@ public : X = L.template marked<Lower>().solveTriangular(B); } + static inline void trisolve_lower_matrix(const gene_matrix & L, const gene_matrix& B, gene_matrix& X, int N){ + X = L.template marked<Lower>().solveTriangular(B); + } + static inline void cholesky(const gene_matrix & X, gene_matrix & C, int N){ C = X.cholesky().matrixL(); // C = X; |