From 2946992ad4bc929d231282b477a23175cc1c0090 Mon Sep 17 00:00:00 2001
From: Georg Drenkhahn <Georg.Drenkhahn@gmail.com>
Date: Mon, 22 Sep 2014 17:59:02 +0200
Subject: [PATCH] Using StorageIndexType for loop assigning initial
 permutation. Adding assert for index overflow.

---
 Eigen/src/Core/PermutationMatrix.h | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/Eigen/src/Core/PermutationMatrix.h b/Eigen/src/Core/PermutationMatrix.h
index b67046144..d5113ce19 100644
--- a/Eigen/src/Core/PermutationMatrix.h
+++ b/Eigen/src/Core/PermutationMatrix.h
@@ -374,7 +374,9 @@ class PermutationMatrix : public PermutationBase<PermutationMatrix<SizeAtCompile
     PermutationMatrix(const Transpose<PermutationBase<Other> >& other)
       : m_indices(other.nestedPermutation().size())
     {
-      for (typename IndicesType::Index i=0; i<m_indices.size();++i) m_indices.coeffRef(other.nestedPermutation().indices().coeff(i)) = i;
+      eigen_internal_assert(m_indices.size() <= NumTraits<StorageIndexType>::highest());
+      for (StorageIndexType i=0; i<m_indices.size();++i)
+        m_indices.coeffRef(other.nestedPermutation().indices().coeff(i)) = i;
     }
     template<typename Lhs,typename Rhs>
     PermutationMatrix(internal::PermPermProduct_t, const Lhs& lhs, const Rhs& rhs)