diff --git a/src/base/matrix.rs b/src/base/matrix.rs
index 91e47624..95893e22 100644
--- a/src/base/matrix.rs
+++ b/src/base/matrix.rs
@@ -274,6 +274,15 @@ impl<N: Scalar, R: Dim, C: Dim, S: Storage<N, R, C>> Matrix<N, R, C, S> {
         self.data.get_unchecked(irow, icol)
     }
 
+    /// Returns a pointer to the start of the matrix.
+    ///
+    /// If the matrix is not empty, this pointer is guaranteed to be aligned
+    /// and non-null.
+    #[inline]
+    pub fn as_ptr(&self) -> *const N {
+        self.data.ptr()
+    }
+
     /// Tests whether `self` and `rhs` are equal up to a given epsilon.
     ///
     /// See `relative_eq` from the `RelativeEq` trait for more details.
@@ -550,6 +559,15 @@ impl<N: Scalar, R: Dim, C: Dim, S: StorageMut<N, R, C>> Matrix<N, R, C, S> {
         MatrixIterMut::new(&mut self.data)
     }
 
+    /// Returns a mutable pointer to the start of the matrix.
+    ///
+    /// If the matrix is not empty, this pointer is guaranteed to be aligned
+    /// and non-null.
+    #[inline]
+    pub fn as_mut_ptr(&mut self) -> *mut N {
+        self.data.ptr_mut()
+    }
+
     /// Gets a mutable reference to the i-th element of this matrix.
     #[inline]
     pub unsafe fn get_unchecked_mut(&mut self, irow: usize, icol: usize) -> &mut N {