diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2016-12-02 11:32:03 +0100 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2016-12-02 11:32:03 +0100 |
| commit | d2718d662c9759c29ed98456ab4abb128fc1a614 (patch) | |
| tree | 3e9a590acc61e1c60e1740ed9124bc35d4edf913 /bench/btl/libs/STL/STL_interface.hh | |
| parent | 8c24723a09cd6b3852d2a5a0748f6a250d1d40f5 (diff) | |
Re-enable A^T*A action in BTL
Diffstat (limited to 'bench/btl/libs/STL/STL_interface.hh')
| -rw-r--r-- | bench/btl/libs/STL/STL_interface.hh | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/bench/btl/libs/STL/STL_interface.hh b/bench/btl/libs/STL/STL_interface.hh index ef4cc9233..16658c4ba 100644 --- a/bench/btl/libs/STL/STL_interface.hh +++ b/bench/btl/libs/STL/STL_interface.hh @@ -78,18 +78,18 @@ public : cible[i][j]=source[i][j]; } -// static inline void ata_product(const gene_matrix & A, gene_matrix & X, int N) -// { -// real somme; -// for (int j=0;j<N;j++){ -// for (int i=0;i<N;i++){ -// somme=0.0; -// for (int k=0;k<N;k++) -// somme += A[i][k]*A[j][k]; -// X[j][i]=somme; -// } -// } -// } + static inline void ata_product(const gene_matrix & A, gene_matrix & X, int N) + { + real somme; + for (int j=0;j<N;j++){ + for (int i=0;i<N;i++){ + somme=0.0; + for (int k=0;k<N;k++) + somme += A[i][k]*A[j][k]; + X[j][i]=somme; + } + } + } static inline void aat_product(const gene_matrix & A, gene_matrix & X, int N) { |