BTl: GMM++ LU is not a full pivoting LU

This commit is contained in:
Gael Guennebaud 2011-03-22 15:39:23 +01:00
parent 22c7609d72
commit da3f3586e0
2 changed files with 3 additions and 3 deletions

View File

@ -123,7 +123,7 @@ public :
gmm::lower_tri_solve(L, X, false);
}
static inline void lu_decomp(const gene_matrix & X, gene_matrix & R, int N){
static inline void partial_lu_decomp(const gene_matrix & X, gene_matrix & R, int N){
gmm::copy(X,R);
std::vector<int> ipvt(N);
gmm::lu_factor(R, ipvt);

View File

@ -20,7 +20,7 @@
#include "bench.hh"
#include "basic_actions.hh"
#include "action_hessenberg.hh"
#include "action_lu_decomp.hh"
#include "action_partial_lu.hh"
BTL_MAIN;
@ -40,7 +40,7 @@ int main()
bench<Action_trisolve<gmm_interface<REAL_TYPE> > >(MIN_MM,MAX_MM,NB_POINT);
//bench<Action_lu_solve<blitz_LU_solve_interface<REAL_TYPE> > >(MIN_LU,MAX_LU,NB_POINT);
bench<Action_lu_decomp<gmm_interface<REAL_TYPE> > >(MIN_MM,MAX_MM,NB_POINT);
bench<Action_partial_lu<gmm_interface<REAL_TYPE> > >(MIN_MM,MAX_MM,NB_POINT);
bench<Action_hessenberg<gmm_interface<REAL_TYPE> > >(MIN_MM,MAX_MM,NB_POINT);
bench<Action_tridiagonalization<gmm_interface<REAL_TYPE> > >(MIN_MM,MAX_MM,NB_POINT);