From 67e582987792166ca15a8e0303968e78b5160626 Mon Sep 17 00:00:00 2001
From: Volodymyr Orlov <vorl@Echo.local>
Date: Wed, 25 Nov 2020 12:23:04 -0800
Subject: [PATCH] simplifies generic matrix.ab implementation

---
 src/linalg/high_order.rs | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/linalg/high_order.rs b/src/linalg/high_order.rs
index 359c4a1..493c737 100644
--- a/src/linalg/high_order.rs
+++ b/src/linalg/high_order.rs
@@ -19,7 +19,7 @@ pub trait HighOrderOperations<T: RealNumber>: BaseMatrix<T> {
     /// ```
     fn ab(&self, a_transpose: bool, b: &Self, b_transpose: bool) -> Self {
         match (a_transpose, b_transpose) {
-            (true, true) => self.transpose().matmul(&b.transpose()),
+            (true, true) => b.matmul(self).transpose(),
             (false, true) => self.matmul(&b.transpose()),
             (true, false) => self.transpose().matmul(b),
             (false, false) => self.matmul(b),