From baeda307c121812cf7501caefeb5c287c327f242 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Vladim=C3=ADr=20Vondru=C5=A1?= <mosra@centrum.cz>
Date: Wed, 27 Feb 2013 13:13:23 +0100
Subject: [PATCH] Math: do Matrix::trace() using diagonal() and Vector::sum().

Use (future) SIMD optimizations in Vector::sum() instead of doing it all
over again.
---
 src/Math/Matrix.h | 7 +------
 1 file changed, 1 insertion(+), 6 deletions(-)

diff --git a/src/Math/Matrix.h b/src/Math/Matrix.h
index b4f482391..6334038cc 100644
--- a/src/Math/Matrix.h
+++ b/src/Math/Matrix.h
@@ -100,12 +100,7 @@ template<std::size_t size, class T> class Matrix: public RectangularMatrix<size,
          * @f]
          */
         T trace() const {
-            T out(0);
-
-            for(std::size_t i = 0; i != size; ++i)
-                out += (*this)[i][i];
-
-            return out;
+            return this->diagonal().sum();
         }
 
         /** @brief %Matrix without given column and row */