diff options
author | Gael Guennebaud <g.gael@free.fr> | 2014-04-17 21:01:45 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2014-04-17 21:01:45 +0200 |
commit | 9777a5ca60f0a82bb789f55912fd046ab7f3d15d (patch) | |
tree | 294378a09230f2abbf9f6da52234b9ce77c3005c /bench/btl/libs/blaze/blaze_interface.hh | |
parent | 9746396d1b8d039d3d0d6537ad477135e5e9d3f5 (diff) |
Various minor fixes in BTL
Diffstat (limited to 'bench/btl/libs/blaze/blaze_interface.hh')
-rw-r--r-- | bench/btl/libs/blaze/blaze_interface.hh | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/bench/btl/libs/blaze/blaze_interface.hh b/bench/btl/libs/blaze/blaze_interface.hh index 8020fef27..ed43ecdd4 100644 --- a/bench/btl/libs/blaze/blaze_interface.hh +++ b/bench/btl/libs/blaze/blaze_interface.hh @@ -84,25 +84,25 @@ public : X = (A*B); } -// static inline void transposed_matrix_matrix_product(const gene_matrix & A, const gene_matrix & B, gene_matrix & X, int N){ -// X = (trans(A)*trans(B)); -// } + static inline void transposed_matrix_matrix_product(const gene_matrix & A, const gene_matrix & B, gene_matrix & X, int N){ + X = A.transpose()*B.transpose(); + } -// static inline void ata_product(const gene_matrix & A, gene_matrix & X, int N){ -// X = (trans(A)*A); -// } + static inline void ata_product(const gene_matrix & A, gene_matrix & X, int N){ + X = (A.transpose()*A); + } -// static inline void aat_product(const gene_matrix & A, gene_matrix & X, int N){ -// X = (A*trans(A)); -// } + static inline void aat_product(const gene_matrix & A, gene_matrix & X, int N){ + X = (A*A.transpose()); + } static inline void matrix_vector_product(gene_matrix & A, gene_vector & B, gene_vector & X, int N){ X = (A*B); } -// static inline void atv_product(gene_matrix & A, gene_vector & B, gene_vector & X, int N){ -// X = (trans(A)*B); -// } + static inline void atv_product(gene_matrix & A, gene_vector & B, gene_vector & X, int N){ + X = (A.transpose()*B); + } static inline void axpy(const real coef, const gene_vector & X, gene_vector & Y, int N){ Y += coef * X; |