From e6265da98033c6ebf4e939737b8de2cd4d38cbd6 Mon Sep 17 00:00:00 2001
From: Mara Bos <m-ou.se@m-ou.se>
Date: Tue, 16 Apr 2019 10:11:27 +0200
Subject: [PATCH] Implement Default for MatrixMN.

---
 src/base/array_storage.rs | 13 +++++++++++++
 src/base/matrix.rs        | 15 +++++++++++++++
 2 files changed, 28 insertions(+)

diff --git a/src/base/array_storage.rs b/src/base/array_storage.rs
index bebb8740..450127ba 100644
--- a/src/base/array_storage.rs
+++ b/src/base/array_storage.rs
@@ -48,6 +48,19 @@ where
 /// Renamed to [ArrayStorage].
 pub type MatrixArray<N, R, C> = ArrayStorage<N, R, C>;
 
+impl<N, R, C> Default for ArrayStorage<N, R, C>
+where
+    R: DimName,
+    C: DimName,
+    R::Value: Mul<C::Value>,
+    Prod<R::Value, C::Value>: ArrayLength<N>,
+    N: Default,
+{
+    fn default() -> Self {
+        ArrayStorage { data: Default::default() }
+    }
+}
+
 impl<N, R, C> Hash for ArrayStorage<N, R, C>
 where
     N: Hash,
diff --git a/src/base/matrix.rs b/src/base/matrix.rs
index 7119a8ad..1aa171e4 100644
--- a/src/base/matrix.rs
+++ b/src/base/matrix.rs
@@ -90,6 +90,21 @@ impl<N: Scalar, R: Dim, C: Dim, S: fmt::Debug> fmt::Debug for Matrix<N, R, C, S>
     }
 }
 
+impl<N, R, C, S> Default for Matrix<N, R, C, S>
+where
+    N: Scalar,
+    R: Dim,
+    C: Dim,
+    S: Default,
+{
+    fn default() -> Self {
+        Matrix {
+            data: Default::default(),
+            _phantoms: PhantomData,
+        }
+    }
+}
+
 #[cfg(feature = "serde-serialize")]
 impl<N, R, C, S> Serialize for Matrix<N, R, C, S>
 where