diff --git a/src/linalg/permutation_sequence.rs b/src/linalg/permutation_sequence.rs
index df7ad7dd5..cb738aaad 100644
--- a/src/linalg/permutation_sequence.rs
+++ b/src/linalg/permutation_sequence.rs
@@ -161,4 +161,10 @@ where
             -T::one()
         }
     }
+
+    /// Returns the permutation indices as a vector of tuples to the user.
+    #[inline]
+    pub fn permutation_indices(&self) -> Vec<(usize, usize)> {
+        self.ipiv.rows_range(..self.len).iter().cloned().collect()
+    }
 }