diff options
author | Gael Guennebaud <g.gael@free.fr> | 2011-03-23 16:28:09 +0100 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2011-03-23 16:28:09 +0100 |
commit | 3ef0da6efb537ab0647bd853a43411e251c56367 (patch) | |
tree | bed0be4f1f7f49e0235de610842ef1684039dcbe /bench/btl/actions | |
parent | 816541d82c3b853c0c73b42f9e978f069afc43e1 (diff) |
fix tridiagonalization action
Diffstat (limited to 'bench/btl/actions')
-rw-r--r-- | bench/btl/actions/action_hessenberg.hh | 13 |
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); |