aboutsummaryrefslogtreecommitdiffhomepage
path: root/bench/btl
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2011-03-23 16:28:09 +0100
committerGravatar Gael Guennebaud <g.gael@free.fr>2011-03-23 16:28:09 +0100
commit3ef0da6efb537ab0647bd853a43411e251c56367 (patch)
treebed0be4f1f7f49e0235de610842ef1684039dcbe /bench/btl
parent816541d82c3b853c0c73b42f9e978f069afc43e1 (diff)
fix tridiagonalization action
Diffstat (limited to 'bench/btl')
-rw-r--r--bench/btl/actions/action_hessenberg.hh13
1 files changed, 8 insertions, 5 deletions
diff --git a/bench/btl/actions/action_hessenberg.hh b/bench/btl/actions/action_hessenberg.hh
index 5dd12a9bc..2100ebd89 100644
--- a/bench/btl/actions/action_hessenberg.hh
+++ b/bench/btl/actions/action_hessenberg.hh
@@ -139,11 +139,14 @@ public :
MESSAGE("Action_tridiagonalization Ctor");
// STL vector initialization
- typename Interface::stl_matrix tmp;
- init_matrix<pseudo_random>(tmp,_size);
- init_matrix<null_function>(X_stl,_size);
- STL_interface<typename Interface::real_type>::aat_product(tmp,X_stl,_size);
-
+ init_matrix<pseudo_random>(X_stl,_size);
+
+ for(int i=0; i<_size; ++i)
+ {
+ for(int j=0; j<i; ++j)
+ X_stl[i][j] = X_stl[j][i];
+ }
+
init_matrix<null_function>(C_stl,_size);
init_matrix<null_function>(resu_stl,_size);