diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2009-02-06 14:01:01 +0000 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2009-02-06 14:01:01 +0000 |
| commit | 19b035ee11093965387833ea7baece5e8f571602 (patch) | |
| tree | 18dc52723cc589ccce74eaaeee3000da86f40ffc /bench | |
| parent | cc90495e3091ef68927c8d3f1999fc618aec0c77 (diff) | |
s/cholesky/llt in precompiled lib and BTL
Diffstat (limited to 'bench')
| -rw-r--r-- | bench/benchCholesky.cpp | 7 | ||||
| -rw-r--r-- | bench/btl/libs/eigen2/eigen2_interface.hh | 2 |
2 files changed, 4 insertions, 5 deletions
diff --git a/bench/benchCholesky.cpp b/bench/benchCholesky.cpp index 6772b86cc..c722d6b98 100644 --- a/bench/benchCholesky.cpp +++ b/bench/benchCholesky.cpp @@ -117,15 +117,14 @@ __attribute__ ((noinline)) void benchLLT(const MatrixType& m) int main(int argc, char* argv[]) { -// const int dynsizes[] = {/*4,6,8,12,16,24,32,49,64,67,128,129,130,131,132,*/256,257,258,259,260,512,900,0}; + const int dynsizes[] = {4,6,8,16,24,32,49,64,128,256,512,900,0}; std::cout << "size no sqrt standard"; // #ifdef BENCH_GSL // std::cout << " GSL (standard + double + ATLAS) "; // #endif std::cout << "\n"; -// -// for (uint i=0; dynsizes[i]>0; ++i) -// benchLLT(Matrix<Scalar,Dynamic,Dynamic>(dynsizes[i],dynsizes[i])); + for (uint i=0; dynsizes[i]>0; ++i) + benchLLT(Matrix<Scalar,Dynamic,Dynamic>(dynsizes[i],dynsizes[i])); benchLLT(Matrix<Scalar,2,2>()); benchLLT(Matrix<Scalar,3,3>()); diff --git a/bench/btl/libs/eigen2/eigen2_interface.hh b/bench/btl/libs/eigen2/eigen2_interface.hh index 40d8f9c70..47fe58135 100644 --- a/bench/btl/libs/eigen2/eigen2_interface.hh +++ b/bench/btl/libs/eigen2/eigen2_interface.hh @@ -142,7 +142,7 @@ public : } static inline void cholesky(const gene_matrix & X, gene_matrix & C, int N){ - C = X.cholesky().matrixL(); + C = X.llt().matrixL(); // C = X; // Cholesky<gene_matrix>::computeInPlace(C); // Cholesky<gene_matrix>::computeInPlaceBlock(C); |