aboutsummaryrefslogtreecommitdiffhomepage
path: root/bench/btl/libs/eigen3/eigen3_interface.hh
diff options
context:
space:
mode:
Diffstat (limited to 'bench/btl/libs/eigen3/eigen3_interface.hh')
-rw-r--r--bench/btl/libs/eigen3/eigen3_interface.hh6
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);