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;