diff --git a/include/cglm-mat3.h b/include/cglm-mat3.h
index 88ce37a..ad0f544 100644
--- a/include/cglm-mat3.h
+++ b/include/cglm-mat3.h
@@ -148,34 +148,4 @@ glm_mat3_mulv(mat3 m, vec3 v, vec3 dest) {
   dest[2] = m[0][2] * v[0] + m[1][2] * v[1] + m[2][2] * v[2];
 }
 
-CGLM_INLINE
-void
-glm_mat3_print(mat3 matrix,
-               FILE * __restrict ostream) {
-  int i;
-  int j;
-
-#define m 3
-#define n 3
-
-  fprintf(ostream, "Matrix (float%dx%d):\n", m, n);
-
-  for (i = 0; i < m; i++) {
-    fprintf(ostream, "\t|");
-    for (j = 0; j < n; j++) {
-      fprintf(ostream, "%0.4f", matrix[j][i]);;
-
-      if (j != n - 1)
-        fprintf(ostream, "\t");
-    }
-
-    fprintf(ostream, "|\n");
-  }
-
-  fprintf(ostream, "\n");
-
-#undef m
-#undef n
-}
-
 #endif /* cglm_mat3_h */
diff --git a/include/cglm-util.h b/include/cglm-util.h
index 59c3a94..999946b 100644
--- a/include/cglm-util.h
+++ b/include/cglm-util.h
@@ -42,6 +42,36 @@ glm_mat4_print(mat4   matrix,
 #undef n
 }
 
+CGLM_INLINE
+void
+glm_mat3_print(mat3 matrix,
+               FILE * __restrict ostream) {
+  int i;
+  int j;
+
+#define m 3
+#define n 3
+
+  fprintf(ostream, "Matrix (float%dx%d):\n", m, n);
+
+  for (i = 0; i < m; i++) {
+    fprintf(ostream, "\t|");
+    for (j = 0; j < n; j++) {
+      fprintf(ostream, "%0.4f", matrix[j][i]);;
+
+      if (j != n - 1)
+        fprintf(ostream, "\t");
+    }
+
+    fprintf(ostream, "|\n");
+  }
+
+  fprintf(ostream, "\n");
+
+#undef m
+#undef n
+}
+
 CGLM_INLINE
 void
 glm_vec4_print(vec4 vec,
diff --git a/include/cglm.h b/include/cglm.h
index e5ed6c0..36e75b5 100644
--- a/include/cglm.h
+++ b/include/cglm.h
@@ -20,6 +20,7 @@ typedef vec4 versor;
 
 #include "cglm-vec.h"
 #include "cglm-mat.h"
+#include "cglm-mat3.h"
 #include "cglm-util.h"
 #include "cglm-affine.h"
 #include "cglm-cam.h"