diff --git a/src/base/array_storage.rs b/src/base/array_storage.rs
index dc6b22ce..68da4c51 100644
--- a/src/base/array_storage.rs
+++ b/src/base/array_storage.rs
@@ -25,7 +25,6 @@ use crate::base::storage::{
     ContiguousStorage, ContiguousStorageMut, Owned, ReshapableStorage, Storage, StorageMut,
 };
 use crate::base::Scalar;
-use crate::SMatrix;
 
 /*
  *
@@ -300,14 +299,3 @@ where
         self.as_slice().iter().fold(0, |acc, e| acc + e.extent())
     }
 }
-
-// TODO: Where to put this impl block?
-impl<T, const R: usize, const C: usize> SMatrix<T, R, C> {
-    /// Creates a new statically-allocated matrix from the given [ArrayStorage].
-    #[inline(always)]
-    pub const fn from_array_storage(storage: ArrayStorage<T, R, C>) -> Self {
-        // This is sound because the row and column types are exactly the same as that of the
-        // storage, so there can be no mismatch
-        unsafe { Self::from_data_statically_unchecked(storage) }
-    }
-}
diff --git a/src/base/matrix.rs b/src/base/matrix.rs
index 3b4bc50c..f9b391ef 100644
--- a/src/base/matrix.rs
+++ b/src/base/matrix.rs
@@ -29,7 +29,7 @@ use crate::base::storage::{
     ContiguousStorage, ContiguousStorageMut, Owned, SameShapeStorage, Storage, StorageMut,
 };
 use crate::base::{Const, DefaultAllocator, OMatrix, OVector, Scalar, Unit};
-use crate::SimdComplexField;
+use crate::{ArrayStorage, DMatrix, DVector, Dynamic, SMatrix, SimdComplexField, VecStorage};
 
 /// A square matrix.
 pub type SquareMatrix<T, D, S> = Matrix<T, D, D, S>;
@@ -317,6 +317,47 @@ impl<T, R, C, S> Matrix<T, R, C, S> {
     }
 }
 
+impl<T, const R: usize, const C: usize> SMatrix<T, R, C> {
+    /// Creates a new statically-allocated matrix from the given [ArrayStorage].
+    ///
+    /// This method exists primarily as a workaround for the fact that `from_data` can not
+    /// work in `const fn` contexts.
+    #[inline(always)]
+    pub const fn from_array_storage(storage: ArrayStorage<T, R, C>) -> Self {
+        // This is sound because the row and column types are exactly the same as that of the
+        // storage, so there can be no mismatch
+        unsafe { Self::from_data_statically_unchecked(storage) }
+    }
+}
+
+// TODO: Consider removing/deprecating `from_vec_storage` once we are able to make
+// `from_data` const fn compatible
+impl<T> DMatrix<T> {
+    /// Creates a new heap-allocated matrix from the given [VecStorage].
+    ///
+    /// This method exists primarily as a workaround for the fact that `from_data` can not
+    /// work in `const fn` contexts.
+    pub const fn from_vec_storage(storage: VecStorage<T, Dynamic, Dynamic>) -> Self {
+        // This is sound because the dimensions of the matrix and the storage are guaranteed
+        // to be the same
+        unsafe { Self::from_data_statically_unchecked(storage) }
+    }
+}
+
+// TODO: Consider removing/deprecating `from_vec_storage` once we are able to make
+// `from_data` const fn compatible
+impl<T> DVector<T> {
+    /// Creates a new heap-allocated matrix from the given [VecStorage].
+    ///
+    /// This method exists primarily as a workaround for the fact that `from_data` can not
+    /// work in `const fn` contexts.
+    pub const fn from_vec_storage(storage: VecStorage<T, Dynamic, U1>) -> Self {
+        // This is sound because the dimensions of the matrix and the storage are guaranteed
+        // to be the same
+        unsafe { Self::from_data_statically_unchecked(storage) }
+    }
+}
+
 impl<T: Scalar, R: Dim, C: Dim, S: Storage<T, R, C>> Matrix<T, R, C, S> {
     /// Creates a new matrix with the given data.
     #[inline(always)]
diff --git a/src/base/vec_storage.rs b/src/base/vec_storage.rs
index f8d28769..29da127d 100644
--- a/src/base/vec_storage.rs
+++ b/src/base/vec_storage.rs
@@ -13,7 +13,6 @@ use crate::base::storage::{
 };
 use crate::base::{Scalar, Vector};
 
-use crate::{DMatrix, DVector};
 #[cfg(feature = "abomonation-serialize")]
 use abomonation::Abomonation;
 
@@ -411,21 +410,3 @@ impl<T> Extend<T> for VecStorage<T, Dynamic, U1> {
         self.nrows = Dynamic::new(self.data.len());
     }
 }
-
-impl<T> DMatrix<T> {
-    /// Creates a new heap-allocated matrix from the given [VecStorage].
-    pub const fn from_vec_storage(storage: VecStorage<T, Dynamic, Dynamic>) -> Self {
-        // This is sound because the dimensions of the matrix and the storage are guaranteed
-        // to be the same
-        unsafe { Self::from_data_statically_unchecked(storage) }
-    }
-}
-
-impl<T> DVector<T> {
-    /// Creates a new heap-allocated matrix from the given [VecStorage].
-    pub const fn from_vec_storage(storage: VecStorage<T, Dynamic, U1>) -> Self {
-        // This is sound because the dimensions of the matrix and the storage are guaranteed
-        // to be the same
-        unsafe { Self::from_data_statically_unchecked(storage) }
-    }
-}