aboutsummaryrefslogtreecommitdiffhomepage
path: root/bench/btl/libs
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2010-02-26 12:51:20 +0100
committerGravatar Gael Guennebaud <g.gael@free.fr>2010-02-26 12:51:20 +0100
commitc05047d28e7d84d2a5312c3c958063d4415b6dcc (patch)
tree14de1065948bf4567cbece02966a6bde3216b08c /bench/btl/libs
parent3ac2b96a2f131e8162d39f0976cfb31b1a853237 (diff)
fix some BTL issues
Diffstat (limited to 'bench/btl/libs')
-rw-r--r--bench/btl/libs/eigen2/eigen2_interface.hh12
1 files changed, 6 insertions, 6 deletions
diff --git a/bench/btl/libs/eigen2/eigen2_interface.hh b/bench/btl/libs/eigen2/eigen2_interface.hh
index a8b5b884f..d8d1a607f 100644
--- a/bench/btl/libs/eigen2/eigen2_interface.hh
+++ b/bench/btl/libs/eigen2/eigen2_interface.hh
@@ -166,7 +166,7 @@ public :
}
static EIGEN_DONT_INLINE void rot(gene_vector & A, gene_vector & B, real c, real s, int N){
- ei_apply_rotation_in_the_plane(A, B, c, s);
+ ei_apply_rotation_in_the_plane(A, B, PlanarRotation<real>(c,s));
}
static inline void atv_product(gene_matrix & A, gene_vector & B, gene_vector & X, int N){
@@ -207,15 +207,15 @@ public :
}
static inline void lu_decomp(const gene_matrix & X, gene_matrix & C, int N){
+ C = X.fullPivLu().matrixLU();
+ }
+
+ static inline void partial_lu_decomp(const gene_matrix & X, gene_matrix & C, int N){
RowVectorXi piv(N);
int nb;
C = X;
ei_partial_lu_inplace(C,piv,nb);
- //C = X.lu().matrixLU();
- }
-
- static inline void partial_lu_decomp(const gene_matrix & X, gene_matrix & C, int N){
- C = X.partialPivLu().matrixLU();
+// C = X.partialPivLu().matrixLU();
}
static inline void tridiagonalization(const gene_matrix & X, gene_matrix & C, int N){