diff --git a/matrix.hpp b/matrix.hpp
index 43668483..c77c25b7 100644
--- a/matrix.hpp
+++ b/matrix.hpp
@@ -225,6 +225,19 @@ namespace util {
     template <typename T> matrix<4,4,T> rotation    (T angle, vector<3,T> about);
 
 
+    template <size_t Rows, size_t Cols, typename ValueT>
+    bool
+    relatively_equal (const util::matrix<Rows,Cols,ValueT> &a,
+                      const util::matrix<Rows,Cols,ValueT> &b,
+                      const float percentage)
+    {
+        for (size_t r = 0; r < Rows; ++r)
+            if (!all (relatively_equal (a[r], b[r], percentage)))
+                return false;
+        return true;
+    }
+
+
     ///////////////////////////////////////////////////////////////////////////
     // Convert an affine rotation matrix to euler angles.
     //