diff options
Diffstat (limited to 'bench/btl/libs')
-rw-r--r-- | bench/btl/libs/eigen3/eigen3_interface.hh | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/bench/btl/libs/eigen3/eigen3_interface.hh b/bench/btl/libs/eigen3/eigen3_interface.hh index 2fca393e9..31bcc1f93 100644 --- a/bench/btl/libs/eigen3/eigen3_interface.hh +++ b/bench/btl/libs/eigen3/eigen3_interface.hh @@ -195,16 +195,16 @@ public : } static inline void trisolve_lower_matrix(const gene_matrix & L, const gene_matrix& B, gene_matrix& X, int N){ - X = L.template triangularView<Lower>().solve(B); + X = L.template triangularView<Upper>().solve(B); } static inline void trmm(const gene_matrix & L, const gene_matrix& B, gene_matrix& X, int N){ - X = L.template triangularView<Lower>() * B; + X.noalias() = L.template triangularView<Lower>() * B; } static inline void cholesky(const gene_matrix & X, gene_matrix & C, int N){ C = X; - internal::llt_inplace<Lower>::blocked(C); + internal::llt_inplace<real,Lower>::blocked(C); //C = X.llt().matrixL(); // C = X; // Cholesky<gene_matrix>::computeInPlace(C); |