diff --git a/src/cgmath/matrix.rs b/src/cgmath/matrix.rs
index 02ba5e1..e5620ed 100644
--- a/src/cgmath/matrix.rs
+++ b/src/cgmath/matrix.rs
@@ -417,13 +417,13 @@ pub trait Matrix
     fn is_symmetric(&self) -> bool;
 }
 
-impl<S: BaseFloat> Add<Matrix2<S>, Matrix2<S>> for Matrix2<S> { #[inline] fn add(&self, other: &Matrix2<S>) -> Matrix2<S> { build(|i| self.i(i).add_v(other.i(i))) } }
-impl<S: BaseFloat> Add<Matrix3<S>, Matrix3<S>> for Matrix3<S> { #[inline] fn add(&self, other: &Matrix3<S>) -> Matrix3<S> { build(|i| self.i(i).add_v(other.i(i))) } }
-impl<S: BaseFloat> Add<Matrix4<S>, Matrix4<S>> for Matrix4<S> { #[inline] fn add(&self, other: &Matrix4<S>) -> Matrix4<S> { build(|i| self.i(i).add_v(other.i(i))) } }
+impl<S: BaseFloat> Add<Matrix2<S>, Matrix2<S>> for Matrix2<S> { #[inline] fn add(&self, other: &Matrix2<S>) -> Matrix2<S> { self.add_m(other) } }
+impl<S: BaseFloat> Add<Matrix3<S>, Matrix3<S>> for Matrix3<S> { #[inline] fn add(&self, other: &Matrix3<S>) -> Matrix3<S> { self.add_m(other) } }
+impl<S: BaseFloat> Add<Matrix4<S>, Matrix4<S>> for Matrix4<S> { #[inline] fn add(&self, other: &Matrix4<S>) -> Matrix4<S> { self.add_m(other) } }
 
-impl<S: BaseFloat> Sub<Matrix2<S>, Matrix2<S>> for Matrix2<S> { #[inline] fn sub(&self, other: &Matrix2<S>) -> Matrix2<S> { build(|i| self.i(i).sub_v(other.i(i))) } }
-impl<S: BaseFloat> Sub<Matrix3<S>, Matrix3<S>> for Matrix3<S> { #[inline] fn sub(&self, other: &Matrix3<S>) -> Matrix3<S> { build(|i| self.i(i).sub_v(other.i(i))) } }
-impl<S: BaseFloat> Sub<Matrix4<S>, Matrix4<S>> for Matrix4<S> { #[inline] fn sub(&self, other: &Matrix4<S>) -> Matrix4<S> { build(|i| self.i(i).sub_v(other.i(i))) } }
+impl<S: BaseFloat> Sub<Matrix2<S>, Matrix2<S>> for Matrix2<S> { #[inline] fn sub(&self, other: &Matrix2<S>) -> Matrix2<S> { self.sub_m(other) } }
+impl<S: BaseFloat> Sub<Matrix3<S>, Matrix3<S>> for Matrix3<S> { #[inline] fn sub(&self, other: &Matrix3<S>) -> Matrix3<S> { self.sub_m(other) } }
+impl<S: BaseFloat> Sub<Matrix4<S>, Matrix4<S>> for Matrix4<S> { #[inline] fn sub(&self, other: &Matrix4<S>) -> Matrix4<S> { self.sub_m(other) } }
 
 impl<S: BaseFloat> Neg<Matrix2<S>> for Matrix2<S> { #[inline] fn neg(&self) -> Matrix2<S> { build(|i| self.i(i).neg()) } }
 impl<S: BaseFloat> Neg<Matrix3<S>> for Matrix3<S> { #[inline] fn neg(&self) -> Matrix3<S> { build(|i| self.i(i).neg()) } }
@@ -433,9 +433,9 @@ impl<S: BaseFloat> Zero for Matrix2<S> { #[inline] fn zero() -> Matrix2<S> { Mat
 impl<S: BaseFloat> Zero for Matrix3<S> { #[inline] fn zero() -> Matrix3<S> { Matrix3::zero() } #[inline] fn is_zero(&self) -> bool { *self == zero() } }
 impl<S: BaseFloat> Zero for Matrix4<S> { #[inline] fn zero() -> Matrix4<S> { Matrix4::zero() } #[inline] fn is_zero(&self) -> bool { *self == zero() } }
 
-impl<S: BaseFloat> Mul<Matrix2<S>, Matrix2<S>> for Matrix2<S> { #[inline] fn mul(&self, other: &Matrix2<S>) -> Matrix2<S> { build(|i| self.i(i).mul_v(other.i(i))) } }
-impl<S: BaseFloat> Mul<Matrix3<S>, Matrix3<S>> for Matrix3<S> { #[inline] fn mul(&self, other: &Matrix3<S>) -> Matrix3<S> { build(|i| self.i(i).mul_v(other.i(i))) } }
-impl<S: BaseFloat> Mul<Matrix4<S>, Matrix4<S>> for Matrix4<S> { #[inline] fn mul(&self, other: &Matrix4<S>) -> Matrix4<S> { build(|i| self.i(i).mul_v(other.i(i))) } }
+impl<S: BaseFloat> Mul<Matrix2<S>, Matrix2<S>> for Matrix2<S> { #[inline] fn mul(&self, other: &Matrix2<S>) -> Matrix2<S> { self.mul_m(other) } }
+impl<S: BaseFloat> Mul<Matrix3<S>, Matrix3<S>> for Matrix3<S> { #[inline] fn mul(&self, other: &Matrix3<S>) -> Matrix3<S> { self.mul_m(other) } }
+impl<S: BaseFloat> Mul<Matrix4<S>, Matrix4<S>> for Matrix4<S> { #[inline] fn mul(&self, other: &Matrix4<S>) -> Matrix4<S> { self.mul_m(other) } }
 
 impl<S: BaseFloat> One for Matrix2<S> { #[inline] fn one() -> Matrix2<S> { Matrix2::identity() } }
 impl<S: BaseFloat> One for Matrix3<S> { #[inline] fn one() -> Matrix3<S> { Matrix3::identity() } }
diff --git a/src/test/matrix.rs b/src/test/matrix.rs
index 8b036a4..d95a787 100644
--- a/src/test/matrix.rs
+++ b/src/test/matrix.rs
@@ -136,6 +136,7 @@ fn test_add_m() {
     let mut mut_a = matrix2::A;
     mut_a.add_self_m(&matrix2::B);
     assert_eq!(mut_a, matrix2::A.add_m(&matrix2::B));
+    assert_eq!(mut_a, matrix2::A + matrix2::B);
 
     // Matrix3
     assert_eq!(matrix3::A.add_m(&matrix3::B),
@@ -145,6 +146,7 @@ fn test_add_m() {
     let mut mut_a = matrix3::A;
     mut_a.add_self_m(&matrix3::B);
     assert_eq!(mut_a, matrix3::A.add_m(&matrix3::B));
+    assert_eq!(mut_a, matrix3::A + matrix3::B);
 
     // Matrix4
     assert_eq!(matrix4::A.add_m(&matrix4::B),
@@ -155,6 +157,7 @@ fn test_add_m() {
     let mut mut_a = matrix4::A;
     mut_a.add_self_m(&matrix4::B);
     assert_eq!(mut_a, matrix4::A.add_m(&matrix4::B));
+    assert_eq!(mut_a, matrix4::A + matrix4::B);
 }
 
 #[test]
@@ -166,6 +169,7 @@ fn test_sub_m() {
     let mut mut_a = matrix2::A;
     mut_a.sub_self_m(&matrix2::B);
     assert_eq!(mut_a, matrix2::A.sub_m(&matrix2::B));
+    assert_eq!(matrix2::A.sub_m(&matrix2::B), matrix2::A - matrix2::B);
 
     // Matrix3
     assert_eq!(matrix3::A.sub_m(&matrix3::B),
@@ -175,6 +179,7 @@ fn test_sub_m() {
     let mut mut_a = matrix3::A;
     mut_a.sub_self_m(&matrix3::B);
     assert_eq!(mut_a, matrix3::A.sub_m(&matrix3::B));
+    assert_eq!(matrix3::A.sub_m(&matrix3::B), matrix3::A - matrix3::B);
 
     // Matrix4
     assert_eq!(matrix4::A.sub_m(&matrix4::B),
@@ -185,6 +190,7 @@ fn test_sub_m() {
     let mut mut_a = matrix4::A;
     mut_a.sub_self_m(&matrix4::B);
     assert_eq!(mut_a, matrix4::A.sub_m(&matrix4::B));
+    assert_eq!(matrix4::A.sub_m(&matrix4::B), matrix4::A - matrix4::B);
 }
 
 #[test]
@@ -208,6 +214,10 @@ fn test_mul_m() {
                             110.0, 254.0, 398.0, 542.0,
                             120.0, 280.0, 440.0, 600.0,
                             130.0, 306.0, 482.0, 658.0));
+
+    assert_eq!(matrix2::A.mul_m(&matrix2::B), matrix2::A * matrix2::B);
+    assert_eq!(matrix3::A.mul_m(&matrix3::B), matrix3::A * matrix3::B);
+    assert_eq!(matrix4::A.mul_m(&matrix4::B), matrix4::A * matrix4::B);
 }
 
 #[test]