tomoto 0.2.3 → 0.3.0

data/vendor/EigenRand/EigenRand/Dists/Discrete.h

@@ -2,10 +2,10 @@
  * @file Discrete.h
  * @author bab2min (bab2min@gmail.com)
  * @brief 
- * @version 0.3.0
- * @date 2020-10-07
+ * @version 0.3.3
+ * @date 2021-03-31
  *
- * @copyright Copyright (c) 2020
+ * @copyright Copyright (c) 2020-2021
  * 
  */
 
@@ -49,11 +49,11 @@ namespace Eigen
 				bitmask = o.bitmask;
 				if (msize)
 				{
-					arr = std::unique_ptr<_Precision[]>(new _Precision[1 << bitsize]);
-					alias = std::unique_ptr<_Size[]>(new _Size[1 << bitsize]);
+					arr = std::unique_ptr<_Precision[]>(new _Precision[(size_t)1 << bitsize]);
+					alias = std::unique_ptr<_Size[]>(new _Size[(size_t)1 << bitsize]);
 
-					std::copy(o.arr.get(), o.arr.get() + (1 << bitsize), arr.get());
-					std::copy(o.alias.get(), o.alias.get() + (1 << bitsize), alias.get());
+					std::copy(o.arr.get(), o.arr.get() + ((size_t)1 << bitsize), arr.get());
+					std::copy(o.alias.get(), o.alias.get() + ((size_t)1 << bitsize), alias.get());
 				}
 				return *this;
 			}
@@ -125,7 +125,7 @@ namespace Eigen
 					{
 						arr[under] = (_Precision)f[under];
 					}
-					alias[under] = over;
+					alias[under] = (_Size)over;
 					f[over] += f[under] - 1;
 					if (f[over] >= 1 || mm <= over)
 					{
@@ -152,7 +152,7 @@ namespace Eigen
 						{
 							arr[over] = 1;
 						}
-						alias[over] = over;
+						alias[over] = (_Size)over;
 					}
 				}
 
@@ -166,7 +166,7 @@ namespace Eigen
 					{
 						arr[under] = 1;
 					}
-					alias[under] = under;
+					alias[under] = (_Size)under;
 					for (under = mm; under < msize; ++under)
 					{
 						if (f[under] < 1)
@@ -179,7 +179,7 @@ namespace Eigen
 							{
 								arr[under] = 1;
 							}
-							alias[under] = under;
+							alias[under] = (_Size)under;
 						}
 					}
 				}
@@ -260,8 +260,9 @@ namespace Eigen
 				else
 				{
 					size_t bitcnt = bitsize;
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						_Scalar cands = (_Scalar)(rx & bitmask);
 						if (cands <= pdiff) return cands;
 						if (bitcnt + bitsize < 32)
@@ -293,8 +294,9 @@ namespace Eigen
 					auto& cm = Rand::detail::CompressMask<sizeof(Packet)>::get_inst();
 					auto plen = pset1<Packet>(pdiff + 1);
 					size_t bitcnt = bitsize;
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						// accept cands that only < plen
 						auto cands = pand(rx, pbitmask);
 						bool full = false;
@@ -304,7 +306,7 @@ namespace Eigen
 
 						if (bitcnt + bitsize < 32)
 						{
-							rx = psrl(rx, bitsize);
+							rx = psrl<-1>(rx, bitsize);
 							bitcnt += bitsize;
 						}
 						else
@@ -417,7 +419,7 @@ namespace Eigen
 					auto rx = randbits(std::forward<Rng>(rng));
 					auto albit = rx & alias_table.get_bitmask();
 					uint32_t alx = (uint32_t)(rx >> (sizeof(rx) * 8 - 31));
-					if (alx < alias_table.get_prob()[albit]) return albit;
+					if (alx < alias_table.get_prob()[albit]) return (_Scalar)albit;
 					return alias_table.get_alias()[albit];
 				}
 			}
@@ -451,7 +453,7 @@ namespace Eigen
 				{
 					auto rx = randbits.template packetOp<PacketType>(std::forward<Rng>(rng));
 					auto albit = pand(rx, pset1<PacketType>(alias_table.get_bitmask()));
-					auto c = pcmplt(psrl(rx, 1), pgather(alias_table.get_prob(), albit));
+					auto c = pcmplt(psrl<1>(rx), pgather(alias_table.get_prob(), albit));
 					ret = pblendv(c, albit, pgather(alias_table.get_alias(), albit));
 				}
 
@@ -673,6 +675,8 @@ namespace Eigen
 				}
 			}
 
+#ifdef EIGEN_VECTORIZE_NEON
+#else
 			template<typename Packet, typename Rng>
 			EIGEN_STRONG_INLINE const Packet packetOp(Rng&& rng)
 			{
@@ -681,7 +685,7 @@ namespace Eigen
 				if (!cdf.empty())
 				{
 					auto ret = pset1<Packet>(cdf.size());
-#ifdef EIGEN_VECTORIZE_AVX
+	#ifdef EIGEN_VECTORIZE_AVX
 					auto rx = ur.template packetOp<Packet4d>(std::forward<Rng>(rng));
 					for (auto& p : cdf)
 					{
@@ -689,7 +693,7 @@ namespace Eigen
 						auto r = combine_low32(c);
 						ret = padd(ret, r);
 					}
-#else
+	#else
 					auto rx1 = ur.template packetOp<DPacket>(rng),
 						rx2 = ur.template packetOp<DPacket>(rng);
 					for (auto& p : cdf)
@@ -697,25 +701,26 @@ namespace Eigen
 						auto pp = pset1<decltype(rx1)>(p);
 						ret = padd(ret, combine_low32(reinterpret_to_int(pcmplt(rx1, pp)), reinterpret_to_int(pcmplt(rx2, pp))));
 					}
-#endif
+	#endif
 					return ret;
 				}
 				else
 				{
-#ifdef EIGEN_VECTORIZE_AVX
+	#ifdef EIGEN_VECTORIZE_AVX
 					using RUtils = RawbitsMaker<Packet, Rng>;
 					auto albit = pand(RUtils{}.rawbits(rng), pset1<Packet>(alias_table.get_bitmask()));
 					auto c = reinterpret_to_int(pcmplt(ur.template packetOp<Packet4d>(rng), pgather(alias_table.get_prob(), _mm256_castsi128_si256(albit))));
 					return pblendv(combine_low32(c), albit, pgather(alias_table.get_alias(), albit));
-#else
+	#else
 					using RUtils = RawbitsMaker<Packet, Rng>;
 					auto albit = pand(RUtils{}.rawbits(rng), pset1<Packet>(alias_table.get_bitmask()));
 					auto c1 = reinterpret_to_int(pcmplt(ur.template packetOp<DPacket>(rng), pgather(alias_table.get_prob(), albit)));
 					auto c2 = reinterpret_to_int(pcmplt(ur.template packetOp<DPacket>(rng), pgather(alias_table.get_prob(), albit, true)));
 					return pblendv(combine_low32(c1, c2), albit, pgather(alias_table.get_alias(), albit));
-#endif
+	#endif
 				}
 			}
+#endif
 		};
 
 		template<typename> class BinomialGen;
@@ -777,8 +782,9 @@ namespace Eigen
 				{
 					_Scalar res;
 					double yx;
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						yx = std::tan(constant::pi * ur(rng));
 						res = (_Scalar)(sqrt_tmean * yx + mean);
 						if (res >= 0 && ur(rng) <= 0.9 * (1.0 + yx * yx)
@@ -800,8 +806,9 @@ namespace Eigen
 				{
 					Packet res = pset1<Packet>(0);
 					PacketType val = pset1<PacketType>(1), pne_mean = pset1<PacketType>(ne_mean);
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						val = pmul(val, ur.template packetOp<PacketType>(rng));
 						auto c = reinterpret_to_int(pcmplt(pne_mean, val));
 						if (pmovemask(c) == 0) break;
@@ -817,15 +824,16 @@ namespace Eigen
 						pmean = pset1<PacketType>(mean),
 						plog_mean = pset1<PacketType>(log_mean),
 						pg1 = pset1<PacketType>(g1);
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						PacketType fres, yx, psin, pcos;
 						psincos(pmul(ppi, ur.template packetOp<PacketType>(rng)), psin, pcos);
 						yx = pdiv(psin, pcos);
 						fres = ptruncate(padd(pmul(psqrt_tmean, yx), pmean));
 
 						auto p1 = pmul(padd(pmul(yx, yx), pset1<PacketType>(1)), pset1<PacketType>(0.9));
-						auto p2 = pexp(psub(psub(pmul(fres, plog_mean), plgamma(padd(fres, pset1<PacketType>(1)))), pg1));
+						auto p2 = pexp(psub(psub(pmul(fres, plog_mean), plgamma_approx(padd(fres, pset1<PacketType>(1)))), pg1));
 
 						auto c1 = pcmple(pset1<PacketType>(0), fres);
 						auto c2 = pcmple(ur.template packetOp<PacketType>(rng), pmul(p1, p2));
@@ -852,7 +860,7 @@ namespace Eigen
 
 			PoissonGen<_Scalar> poisson;
 			_Scalar trials;
-			double p, small_p, g1, sqrt_v, log_small_p, log_small_q;
+			double p = 0, small_p = 0, g1 = 0, sqrt_v = 0, log_small_p = 0, log_small_q = 0;
 
 		public:
 			using Scalar = _Scalar;
@@ -903,8 +911,9 @@ namespace Eigen
 				}
 				else
 				{
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						double ys;
 						ys = std::tan(constant::pi * poisson.ur(rng));
 						res = (_Scalar)(sqrt_v * ys + poisson.mean);
@@ -954,8 +963,9 @@ namespace Eigen
 						plog_small_p = pset1<PacketType>(log_small_p),
 						plog_small_q = pset1<PacketType>(log_small_q),
 						pg1 = pset1<PacketType>(g1);
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						PacketType fres, ys, psin, pcos;
 						psincos(pmul(ppi, poisson.ur.template packetOp<PacketType>(rng)), psin, pcos);
 						ys = pdiv(psin, pcos);
@@ -964,8 +974,8 @@ namespace Eigen
 						auto p1 = pmul(pmul(pset1<PacketType>(1.2), psqrt_v), padd(pset1<PacketType>(1), pmul(ys, ys)));
 						auto p2 = pexp(
 							padd(padd(psub(
-								psub(pg1, plgamma(padd(fres, pset1<PacketType>(1)))),
-								plgamma(psub(padd(ptrials, pset1<PacketType>(1)), fres))
+								psub(pg1, plgamma_approx(padd(fres, pset1<PacketType>(1)))),
+								plgamma_approx(psub(padd(ptrials, pset1<PacketType>(1)), fres))
 							), pmul(fres, plog_small_p)), pmul(psub(ptrials, fres), plog_small_q))
 						);
 
@@ -1503,6 +1513,17 @@ namespace Eigen
 			};
 		}
 	}
+
+#ifdef EIGEN_VECTORIZE_NEON
+	namespace internal
+	{
+		template<typename _Scalar, typename Urng, bool _mutable>
+		struct functor_traits<scalar_rng_adaptor<Rand::DiscreteGen<_Scalar, double>, _Scalar, Urng, _mutable> >
+		{
+			enum { Cost = HugeCost, PacketAccess = 0, IsRepeatable = false };
+		};
+	}
+#endif
 }
 
 #endif

data/vendor/EigenRand/EigenRand/Dists/GammaPoisson.h

@@ -2,10 +2,10 @@
  * @file GammaPoisson.h
  * @author bab2min (bab2min@gmail.com)
  * @brief
- * @version 0.3.0
- * @date 2020-10-07
+ * @version 0.3.3
+ * @date 2021-03-31
  *
- * @copyright Copyright (c) 2020
+ * @copyright Copyright (c) 2020-2021
  *
  */
 
@@ -58,7 +58,6 @@ namespace Eigen
 			EIGEN_STRONG_INLINE const Packet packetOp(Rng&& rng)
 			{
 				using namespace Eigen::internal;
-				using ur_base = UniformRealGen<float>;
 				using PacketType = decltype(reinterpret_to_float(std::declval<Packet>()));
 
 				auto mean = gamma.template packetOp<PacketType>(rng);
@@ -66,8 +65,9 @@ namespace Eigen
 				PacketType val = pset1<PacketType>(1), pne_mean = pexp(pnegate(mean));
 				if (pmovemask(pcmplt(pset1<PacketType>(12), mean)) == 0)
 				{
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						val = pmul(val, ur.template packetOp<PacketType>(rng));
 						auto c = reinterpret_to_int(pcmplt(pne_mean, val));
 						if (pmovemask(c) == 0) break;
@@ -81,16 +81,17 @@ namespace Eigen
 					const PacketType ppi = pset1<PacketType>(constant::pi),
 						psqrt_tmean = psqrt(pmul(pset1<PacketType>(2), mean)),
 						plog_mean = plog(mean),
-						pg1 = psub(pmul(mean, plog_mean), plgamma(padd(mean, pset1<PacketType>(1))));
-					while (1)
+						pg1 = psub(pmul(mean, plog_mean), plgamma_approx(padd(mean, pset1<PacketType>(1))));
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						PacketType fres, yx, psin, pcos;
 						psincos(pmul(ppi, ur.template packetOp<PacketType>(rng)), psin, pcos);
 						yx = pdiv(psin, pcos);
 						fres = ptruncate(padd(pmul(psqrt_tmean, yx), mean));
 
 						auto p1 = pmul(padd(pmul(yx, yx), pset1<PacketType>(1)), pset1<PacketType>(0.9));
-						auto p2 = pexp(psub(psub(pmul(fres, plog_mean), plgamma(padd(fres, pset1<PacketType>(1)))), pg1));
+						auto p2 = pexp(psub(psub(pmul(fres, plog_mean), plgamma_approx(padd(fres, pset1<PacketType>(1)))), pg1));
 
 						auto c1 = pcmple(pset1<PacketType>(0), fres);
 						auto c2 = pcmple(ur.template packetOp<PacketType>(rng), pmul(p1, p2));

data/vendor/EigenRand/EigenRand/Dists/NormalExp.h

@@ -2,10 +2,10 @@
  * @file NormalExp.h
  * @author bab2min (bab2min@gmail.com)
  * @brief 
- * @version 0.3.0
- * @date 2020-10-07
+ * @version 0.3.3
+ * @date 2021-03-31
  * 
- * @copyright Copyright (c) 2020
+ * @copyright Copyright (c) 2020-2021
  * 
  */
 
@@ -27,7 +27,7 @@ namespace Eigen
 		{
 			static_assert(std::is_floating_point<_Scalar>::value, "normalDist needs floating point types.");
 			bool valid = false;
-			UniformRealGen<_Scalar> ur;
+			StdUniformRealGen<_Scalar> ur;
 			
 		public:
 			using Scalar = _Scalar;
@@ -44,8 +44,9 @@ namespace Eigen
 				valid = true;
 
 				_Scalar v1, v2, sx;
-				while (1)
+				for (int _i = 0; ; ++_i)
 				{
+					EIGENRAND_CHECK_INFINITY_LOOP();
 					v1 = 2 * ur(rng) - 1;
 					v2 = 2 * ur(rng) - 1;
 					sx = v1 * v1 + v2 * v2;
@@ -186,7 +187,7 @@ namespace Eigen
 		class StudentTGen : public GenBase<StudentTGen<_Scalar>, _Scalar>
 		{
 			static_assert(std::is_floating_point<_Scalar>::value, "studentT needs floating point types.");
-			UniformRealGen<_Scalar> ur;
+			StdUniformRealGen<_Scalar> ur;
 			_Scalar n;
 
 		public:
@@ -213,8 +214,9 @@ namespace Eigen
 			{
 				using namespace Eigen::internal;
 				_Scalar v1, v2, sx;
-				while (1)
+				for (int _i = 0; ; ++_i)
 				{
+					EIGENRAND_CHECK_INFINITY_LOOP();
 					v1 = 2 * ur(rng) - 1;
 					v2 = 2 * ur(rng) - 1;
 					sx = v1 * v1 + v2 * v2;
@@ -238,8 +240,7 @@ namespace Eigen
 					psub(pexp(pmul(plog(u1), pset1<Packet>(-2 / n))), pset1<Packet>(1))
 				));
 				auto theta = pmul(pset1<Packet>(2 * constant::pi), u2);
-				Packet sintheta, costheta;
-
+				//Packet sintheta, costheta;
 				//psincos(theta, sintheta, costheta);
 				return pmul(radius, psin(theta));
 			}
@@ -257,7 +258,7 @@ namespace Eigen
 		{
 			friend GammaGen<_Scalar>;
 			static_assert(std::is_floating_point<_Scalar>::value, "expDist needs floating point types.");
-			UniformRealGen<_Scalar> ur;
+			StdUniformRealGen<_Scalar> ur;
 			_Scalar lambda = 1;
 
 		public:
@@ -325,7 +326,7 @@ namespace Eigen
 			GammaGen(_Scalar _alpha = 1, _Scalar _beta = 1)
 				: alpha{ _alpha }, beta{ _beta }
 			{
-				px = constant::e / (alpha + constant::e);
+				px = (_Scalar)(constant::e / (alpha + constant::e));
 				sqrt = std::sqrt(2 * alpha - 1);
 			}
 
@@ -342,8 +343,9 @@ namespace Eigen
 				if (alpha < 1)
 				{
 					_Scalar ux, vx, xx, qx;
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						ux = expon.ur(rng);
 						vx = expon.ur.nzur_scalar(rng);
 
@@ -380,8 +382,9 @@ namespace Eigen
 					return -beta * std::log(yx);
 				}
 
-				while (1)
+				for (int _i = 0; ; ++_i)
 				{
+					EIGENRAND_CHECK_INFINITY_LOOP();
 					_Scalar yx, xx;
 					yx = std::tan(constant::pi * expon.ur(rng));
 					xx = sqrt * yx + alpha - 1;
@@ -404,8 +407,9 @@ namespace Eigen
 				RUtils ru;
 				if (alpha < 1)
 				{
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						Packet ux = ru.uniform_real(rng);
 						Packet vx = ru.nonzero_uniform_real(rng);
 
@@ -447,8 +451,9 @@ namespace Eigen
 				}
 				else
 				{
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						Packet alpha_1 = pset1<Packet>(alpha - 1);
 						Packet ys, yc;
 						psincos(pmul(pset1<Packet>(constant::pi), ru.uniform_real(rng)), ys, yc);
@@ -482,7 +487,7 @@ namespace Eigen
 		class WeibullGen : public GenBase<WeibullGen<_Scalar>, _Scalar>
 		{
 			static_assert(std::is_floating_point<_Scalar>::value, "weilbullDist needs floating point types.");
-			UniformRealGen<_Scalar> ur;
+			StdUniformRealGen<_Scalar> ur;
 			_Scalar a = 1, b = 1;
 
 		public:
@@ -531,7 +536,7 @@ namespace Eigen
 		class ExtremeValueGen : public GenBase<ExtremeValueGen<_Scalar>, _Scalar>
 		{
 			static_assert(std::is_floating_point<_Scalar>::value, "extremeValueDist needs floating point types.");
-			UniformRealGen<_Scalar> ur;
+			StdUniformRealGen<_Scalar> ur;
 			_Scalar a = 0, b = 1;
 
 		public:
@@ -623,7 +628,7 @@ namespace Eigen
 		class CauchyGen : public GenBase<CauchyGen<_Scalar>, _Scalar>
 		{
 			static_assert(std::is_floating_point<_Scalar>::value, "cauchyDist needs floating point types.");
-			UniformRealGen<_Scalar> ur; 
+			StdUniformRealGen<_Scalar> ur;
 			_Scalar a = 0, b = 1;
 
 		public:
@@ -680,7 +685,7 @@ namespace Eigen
 			friend FisherFGen<_Scalar>;
 			static_assert(std::is_floating_point<_Scalar>::value, "betaDist needs floating point types.");
 			int cache_rest_cnt = 0;
-			UniformRealGen<_Scalar> ur;
+			StdUniformRealGen<_Scalar> ur;
 			_Scalar a, b;
 			GammaGen<_Scalar> gd1, gd2;
 
@@ -711,8 +716,9 @@ namespace Eigen
 				if (a < 1 && b < 1)
 				{
 					_Scalar x, p1, p2;
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						p1 = std::pow(ur(rng), 1 / a);
 						p2 = std::pow(ur(rng), 1 / b);
 						x = p1 + p2;
@@ -735,8 +741,9 @@ namespace Eigen
 				{
 					auto& cm = Rand::detail::CompressMask<sizeof(Packet)>::get_inst();
 					Packet x, p1, p2;
-					while (1)
+					for (int _i = 0; ; ++_i)
 					{
+						EIGENRAND_CHECK_INFINITY_LOOP();
 						p1 = pexp(pmul(plog(ur.template packetOp<Packet>(rng)), pset1<Packet>(1 / a)));
 						p2 = pexp(pmul(plog(ur.template packetOp<Packet>(rng)), pset1<Packet>(1 / b)));
 						x = padd(p1, p2);

data/vendor/EigenRand/EigenRand/EigenRand

@@ -2,10 +2,10 @@
  * @file EigenRand
  * @author bab2min (bab2min@gmail.com)
  * @brief 
- * @version 0.3.0
- * @date 2020-10-07
+ * @version 0.4.0
+ * @date 2021-09-17
  * 
- * @copyright Copyright (c) 2020
+ * @copyright Copyright (c) 2020-2021
  * 
  */
 
@@ -13,7 +13,12 @@
 #define EIGENRAND_EIGENRAND_H
 
 #include <Eigen/Dense>
-#include <EigenRand/Macro.h>
-#include <EigenRand/Core.h>
 
-#endif
+#include <Eigen/src/Core/util/DisableStupidWarnings.h>
+
+#include "Macro.h"
+#include "Core.h"
+
+#include <Eigen/src/Core/util/ReenableStupidWarnings.h>
+
+#endif

data/vendor/EigenRand/EigenRand/Macro.h

@@ -2,10 +2,10 @@
  * @file Macro.h
  * @author bab2min (bab2min@gmail.com)
  * @brief 
- * @version 0.3.0
- * @date 2020-10-07
+ * @version 0.4.0
+ * @date 2021-09-28
  * 
- * @copyright Copyright (c) 2020
+ * @copyright Copyright (c) 2020-2021
  * 
  */
 
@@ -13,12 +13,18 @@
 #define EIGENRAND_MACRO_H
 
 #define EIGENRAND_WORLD_VERSION 0
-#define EIGENRAND_MAJOR_VERSION 3
-#define EIGENRAND_MINOR_VERSION 2
+#define EIGENRAND_MAJOR_VERSION 4
+#define EIGENRAND_MINOR_VERSION 0
 
-#if EIGEN_VERSION_AT_LEAST(3,3,7)
+#if EIGEN_VERSION_AT_LEAST(3,3,10)
+#define EIGENRAND_EIGEN_34_MODE
+#elif EIGEN_VERSION_AT_LEAST(3,3,4)
+#define EIGENRAND_EIGEN_33_MODE
+#endif
+
+#if EIGEN_VERSION_AT_LEAST(3,3,4) && !EIGEN_VERSION_AT_LEAST(3,4,1)
 #else
-#error At least Eigen 3.3.7 is required.
+#error Eigen 3.3.4 ~ 3.4.0 is required.
 #endif
 
 #endif