diff --git a/unsupported/Eigen/src/FFT/ei_kissfft_impl.h b/unsupported/Eigen/src/FFT/ei_kissfft_impl.h
index 430953aeec7ead7bb5a306325da81df1d2ea166a..e2df9641ca38bcccac68fb85812d454e892e4503 100644
--- a/unsupported/Eigen/src/FFT/ei_kissfft_impl.h
+++ b/unsupported/Eigen/src/FFT/ei_kissfft_impl.h
@@ -29,6 +29,12 @@ struct kiss_cpx_fft
   {
     using numext::sin;
     using numext::cos;
+    typedef typename packet_traits<Scalar>::type PacketType;
+    int psize = unpacket_traits<PacketType>::size;
+    Scalar arr[unpacket_traits<PacketType>::size];
+    Scalar arr_s[unpacket_traits<PacketType>::size];
+    Scalar arr_c[unpacket_traits<PacketType>::size];
+
     m_inverse = inverse;
     m_twiddles.resize(nfft);
     double phinc =  0.25 * double(EIGEN_PI) / nfft;
@@ -37,6 +43,22 @@ struct kiss_cpx_fft
     if ((nfft&1)==0)
       m_twiddles[nfft/2] = Complex(Scalar(-1), Scalar(0));
     int i=1;
+
+    for (;(i+psize-1)*8<nfft;i+=psize)
+    {
+      for(int j=0;j<psize;j++){
+        arr[j] = (i+j)*8*phinc;
+      }
+      PacketType data = ploadu<PacketType>(arr);
+      PacketType c = pcos<PacketType>(data);
+      PacketType s = psin<PacketType>(data);
+      pstoreu<Scalar,PacketType>(arr_s, s);
+      pstoreu<Scalar,PacketType>(arr_c, c);
+      for(int j=0;j<psize;j++){
+        m_twiddles[i+j] = Complex(arr_c[j], arr_s[j]*flip);
+        m_twiddles[nfft-i-j] = Complex(arr_c[j], -1*arr_s[j]*flip);
+      }
+    }
     for (;i*8<nfft;++i)
     {
       Scalar c = Scalar(cos(i*8*phinc));
@@ -44,6 +66,22 @@ struct kiss_cpx_fft
       m_twiddles[i] = Complex(c, s*flip);
       m_twiddles[nfft-i] = Complex(c, -s*flip);
     }
+
+    for (;(i+psize-1)*4<nfft;i+=psize)
+    {
+      for(int j=0;j<psize;j++){
+        arr[j] = (2*nfft-8*(i+j))*phinc;
+      }
+      PacketType data = ploadu<PacketType>(arr);
+      PacketType c = pcos<PacketType>(data);
+      PacketType s = psin<PacketType>(data);
+      pstoreu<Scalar,PacketType>(arr_s, s);
+      pstoreu<Scalar,PacketType>(arr_c, c);
+      for(int j=0;j<psize;j++){
+        m_twiddles[i+j] = Complex(arr_s[j], arr_c[j]*flip);
+        m_twiddles[nfft-i-j] = Complex(arr_s[j], -1*arr_c[j]*flip);
+      }
+    }
     for (;i*4<nfft;++i)
     {
       Scalar c = Scalar(cos((2*nfft-8*i)*phinc));
@@ -51,6 +89,22 @@ struct kiss_cpx_fft
       m_twiddles[i] = Complex(s, c*flip);
       m_twiddles[nfft-i] = Complex(s, -c*flip);
     }
+
+    for (;(i+psize-1)*8<3*nfft;i+=psize)
+    {
+      for(int j=0;j<psize;j++){
+        arr[j] = (8*(i+j)-2*nfft)*phinc;
+      }
+      PacketType data = ploadu<PacketType>(arr);
+      PacketType c = pcos<PacketType>(data);
+      PacketType s = psin<PacketType>(data);
+      pstoreu<Scalar,PacketType>(arr_s, s);
+      pstoreu<Scalar,PacketType>(arr_c, c);
+      for(int j=0;j<psize;j++){
+        m_twiddles[i+j] = Complex(-1*arr_s[j], arr_c[j]*flip);
+        m_twiddles[nfft-i-j] = Complex(-1*arr_s[j], -1*arr_c[j]*flip);
+      }
+    }
     for (;i*8<3*nfft;++i)
     {
       Scalar c = Scalar(cos((8*i-2*nfft)*phinc));
@@ -58,6 +112,22 @@ struct kiss_cpx_fft
       m_twiddles[i] = Complex(-s, c*flip);
       m_twiddles[nfft-i] = Complex(-s, -c*flip);
     }
+
+    for (;(i+psize-1)*2<nfft;i+=psize)
+    {
+      for(int j=0;j<psize;j++){
+        arr[j] = (4*nfft-8*(i+j))*phinc;
+      }
+      PacketType data = ploadu<PacketType>(arr);
+      PacketType c = pcos<PacketType>(data);
+      PacketType s = psin<PacketType>(data);
+      pstoreu<Scalar,PacketType>(arr_s, s);
+      pstoreu<Scalar,PacketType>(arr_c, c);
+      for(int j=0;j<psize;j++){
+        m_twiddles[i+j] = Complex(-1*arr_c[j], arr_s[j]*flip);
+        m_twiddles[nfft-i-j] = Complex(-1*arr_c[j], -1*arr_s[j]*flip);
+      }
+    }   
     for (;i*2<nfft;++i)
     {
       Scalar c = Scalar(cos((4*nfft-8*i)*phinc));
@@ -65,6 +135,7 @@ struct kiss_cpx_fft
       m_twiddles[i] = Complex(-c, s*flip);
       m_twiddles[nfft-i] = Complex(-c, -s*flip);
     }
+
   }
 
   void factorize(int nfft)