tomoto 0.6.0 → 0.6.2

data/vendor/EigenRand/EigenRand/arch/NEON/PacketFilter.h

@@ -2,10 +2,10 @@
  * @file PacketFilter.h
  * @author bab2min (bab2min@gmail.com)
  * @brief
- * @version 0.4.1
- * @date 2022-08-13
+ * @version 0.5.1
+ * @date 2024-09-08
  *
- * @copyright Copyright (c) 2020-2021
+ * @copyright Copyright (c) 2020-2024
  *
  */
 

data/vendor/EigenRand/EigenRand/arch/NEON/RandUtils.h

@@ -2,10 +2,10 @@
  * @file RandUtils.h
  * @author bab2min (bab2min@gmail.com)
  * @brief
- * @version 0.4.1
- * @date 2022-08-13
+ * @version 0.5.1
+ * @date 2024-09-08
  *
- * @copyright Copyright (c) 2020-2021
+ * @copyright Copyright (c) 2020-2024
  *
  */
 
@@ -116,11 +116,29 @@ namespace Eigen
 			}
 		};
 
+#if EIGEN_ARCH_ARM64 && !EIGEN_APPLE_DOUBLE_NEON_BUG
+		template<typename Rng>
+		struct UniformRealUtils<Packet2d, Rng> : public RawbitsMaker<Packet4i, Rng>
+		{
+			EIGEN_STRONG_INLINE Packet2d zero_to_one(Rng& rng)
+			{
+				return pdiv((Packet2d)vcvtq_f64_s64(vreinterpretq_s64_s32(pand(this->rawbits(rng), vreinterpretq_s32_s64(vdupq_n_s64(0x7FFFFFFF))))),
+					pset1<Packet2d>(0x7FFFFFFF));
+			}
+
+			EIGEN_STRONG_INLINE Packet2d uniform_real(Rng& rng)
+			{
+				return bit_to_ur_double(this->rawbits(rng));
+			}
+		};
+
+#else
 		template<typename Gen, typename Urng, bool _mutable>
 		struct functor_traits<scalar_rng_adaptor<Gen, double, Urng, _mutable> >
 		{
 			enum { Cost = HugeCost, PacketAccess = 0, IsRepeatable = false };
 		};
+#endif
 	}
 }
 #endif

data/vendor/EigenRand/EigenRand/arch/SSE/MorePacketMath.h

@@ -2,10 +2,10 @@
  * @file MorePacketMath.h
  * @author bab2min (bab2min@gmail.com)
  * @brief
- * @version 0.4.1
- * @date 2022-08-13
+ * @version 0.5.1
+ * @date 2024-09-08
  *
- * @copyright Copyright (c) 2020-2021
+ * @copyright Copyright (c) 2020-2024
  *
  */
 
@@ -145,6 +145,12 @@ namespace Eigen
 			return _mm_cmpeq_epi32(a, b);
 		}
 
+		template<>
+		EIGEN_STRONG_INLINE Packet4f pcmpeq<Packet4f>(const Packet4f& a, const Packet4f& b)
+		{
+			return _mm_cmpeq_ps(a, b);
+		}
+
 		template<>
 		struct BitShifter<Packet4i>
 		{
@@ -414,7 +420,29 @@ namespace Eigen
 		{
 			return _psin(x);
 		}
+
+		template<> EIGEN_STRONG_INLINE bool predux_all(const Packet4f& x)
+		{
+			return _mm_movemask_ps(x) == 0x0F;
+		}
+
+		template<> EIGEN_STRONG_INLINE bool predux_all(const Packet4i& x)
+		{
+			return predux_all(_mm_castsi128_ps(x));
+		}
+
 #ifdef EIGENRAND_EIGEN_33_MODE
+
+		template<> EIGEN_STRONG_INLINE bool predux_any(const Packet4f& x)
+		{
+			return !!_mm_movemask_ps(x);
+		}
+
+		template<> EIGEN_STRONG_INLINE bool predux_any(const Packet4i& x)
+		{
+			return predux_any(_mm_castsi128_ps(x));
+		}
+
 		template<> EIGEN_DEFINE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS EIGEN_UNUSED
 			Packet2d plog<Packet2d>(const Packet2d& _x)
 		{
@@ -494,7 +522,7 @@ namespace Eigen
 			// negative arg will be NAN, 0 will be -INF
 			return pblendv(iszero_mask, minus_inf, _mm_or_pd(x, invalid_mask));
 		}
-	#endif
+#endif
 	}
 }
 

data/vendor/EigenRand/EigenRand/arch/SSE/PacketFilter.h

@@ -2,10 +2,10 @@
  * @file PacketFilter.h
  * @author bab2min (bab2min@gmail.com)
  * @brief
- * @version 0.4.1
- * @date 2022-08-13
+ * @version 0.5.1
+ * @date 2024-09-08
  *
- * @copyright Copyright (c) 2020-2021
+ * @copyright Copyright (c) 2020-2024
  *
  */
 

data/vendor/EigenRand/EigenRand/arch/SSE/RandUtils.h

@@ -2,10 +2,10 @@
  * @file RandUtils.h
  * @author bab2min (bab2min@gmail.com)
  * @brief
- * @version 0.4.1
- * @date 2022-08-13
+ * @version 0.5.1
+ * @date 2024-09-08
  *
- * @copyright Copyright (c) 2020-2021
+ * @copyright Copyright (c) 2020-2024
  *
  */
 

data/vendor/EigenRand/EigenRand/doc.h

@@ -12,7 +12,7 @@
  You can get 5~10 times speed by just replacing old Eigen's Random 
  or unvectorizable c++11 random number generators with EigenRand.
  
- EigenRand currently supports only x86-64 architecture (SSE, AVX, AVX2) and ARM64 NEON (experimental).
+ EigenRand currently supports only x86-64 architecture (SSE, AVX, AVX2) and ARM64 NEON.
 
  EigenRand is distributed under the MIT License.
 
@@ -85,7 +85,51 @@
 
  Every random distribution function has its corresponding -Like function.
 
- @section getting_started_4 Efficient Reusable Generator
+ @section getting_started_4 Vectorization over Parameters
+ EigenRand's random number generators typically accept scalar parameters. 
+ However, certain generators can generate random numbers efficiently for an array of parameters in an element-wise manner.
+ You can see the full list of distributions which support the vectorization over parameters at @link list_of_supported_distribution @endlink.
+ 
+ @code
+ #include <iostream>
+ #include <Eigen/Dense>
+ #include <EigenRand/EigenRand>
+
+ using namespace Eigen;
+
+ int main()
+ {
+   Rand::P8_mt19937_64 urng{ 42 };
+
+   ArrayXf a{ 10 }, b{ 10 }, c{ 10 };
+	 a << 1, 2, 3, 4, 5, 6, 7, 8, 9, 10;
+	 b << 10, 12, 14, 16, 18, 20, 22, 24, 26, 28;
+   
+   // You can use two array parameters.
+   // The shape of two parameters should be equal in this case.
+   c = Rand::uniformReal(urng, a, b);
+   std::cout << c << std::endl;
+   // c[0] is generated in the range [a[0], b[0]), 
+   // c[1] is generated in the range [a[1], b[1]) ...
+   
+   // Or you can provide one parameter as a scalar
+   // In this case, a scalar parameter is broadcast to the shape of the array parameter.
+   c = Rand::uniformReal(urng, -5, b);
+   std::cout << c << std::endl;
+   // c[0] is generated in the range [-5, b[0]), 
+   // c[1] is generated in the range [-5, b[1]) ...
+
+   c = Rand::uniformReal(urng, a, 11);
+   std::cout << c << std::endl;
+   // c[0] is generated in the range [a[0], 11), 
+   // c[1] is generated in the range [a[1], 11) ...
+   return 0;
+ }
+ @endcode
+
+
+
+ @section getting_started_5 Efficient Reusable Generator
  In the example above, functions, such as `Eigen::Rand::balancedLike`, `Eigen::Rand::normal` and so on, creates a generator internally each time to be called.
  If you want to generate random matrices from the same distribution, consider using Generator classes as following:
 
@@ -114,7 +158,7 @@
  }
  @endcode
 
- @section getting_started_5 Drawing samples from Multivariate Distribution
+ @section getting_started_6 Drawing samples from Multivariate Distribution
  EigenRand provides generators for some multivariate distributions.
 
  @code
@@ -169,33 +213,38 @@
  * 
  @section list_of_supported_distribution_1 Random Distributions for Real types
 
-| Function | Generator | Scalar Type | Description | Equivalent to |
-|:---:|:---:|:---:|:---:|:---:|
-| `Eigen::Rand::balanced` | `Eigen::Rand::BalancedGen` | float, double | generates real values in the [-1, 1] range | `Eigen::DenseBase<Ty>::Random` for floating point types |
-| `Eigen::Rand::beta` | `Eigen::Rand::BetaGen` | float, double | generates real values on a [beta distribution](https://en.wikipedia.org/wiki/Beta_distribution) |  |
-| `Eigen::Rand::cauchy` | `Eigen::Rand::CauchyGen` | float, double | generates real values on the [Cauchy distribution](https://en.wikipedia.org/wiki/Cauchy_distribution). | `std::cauchy_distribution` |
-| `Eigen::Rand::chiSquared` | `Eigen::Rand::ChiSquaredGen` | float, double | generates real values on a [chi-squared distribution](https://en.wikipedia.org/wiki/Chi-squared_distribution). | `std::chi_squared_distribution` |
-| `Eigen::Rand::exponential` | `Eigen::Rand::ExponentialGen` | float, double | generates real values on an [exponential distribution](https://en.wikipedia.org/wiki/Exponential_distribution). | `std::exponential_distribution` |
-| `Eigen::Rand::extremeValue` | `Eigen::Rand::ExtremeValueGen` | float, double | generates real values on an [extreme value distribution](https://en.wikipedia.org/wiki/Generalized_extreme_value_distribution). | `std::extreme_value_distribution` |
-| `Eigen::Rand::fisherF` | `Eigen::Rand::FisherFGen` | float, double | generates real values on the [Fisher's F distribution](https://en.wikipedia.org/wiki/F_distribution). | `std::fisher_f_distribution` |
-| `Eigen::Rand::gamma` | `Eigen::Rand::GammaGen` | float, double | generates real values on a [gamma distribution](https://en.wikipedia.org/wiki/Gamma_distribution). | `std::gamma_distribution` |
-| `Eigen::Rand::lognormal` | `Eigen::Rand::LognormalGen` | float, double | generates real values on a [lognormal distribution](https://en.wikipedia.org/wiki/Lognormal_distribution). | `std::lognormal_distribution` |
-| `Eigen::Rand::normal` | `Eigen::Rand::StdNormalGen`, `Eigen::Rand::NormalGen` | float, double | generates real values on a [normal distribution](https://en.wikipedia.org/wiki/Normal_distribution). | `std::normal_distribution` |
-| `Eigen::Rand::studentT` | `Eigen::Rand::StudentTGen` | float, double | generates real values on the [Student's t distribution](https://en.wikipedia.org/wiki/Student%27s_t-distribution). | `std::student_t_distribution` |
-| `Eigen::Rand::uniformReal` | `Eigen::Rand::StdUniformRealGen`, `Eigen::Rand::UniformRealGen` | float, double | generates real values in the `[-1, 0)` range. | `std::generate_canonical` |
-| `Eigen::Rand::weibull` | `Eigen::Rand::WeibullGen` | float, double | generates real values on the [Weibull distribution](https://en.wikipedia.org/wiki/Weibull_distribution). | `std::weibull_distribution` |
+| Function | Generator | Scalar Type | VoP | Description | Equivalent to |
+|:---:|:---:|:---:|:---:|:---:|:---:|
+| `Eigen::Rand::balanced` | `Eigen::Rand::BalancedGen` | float, double | yes | generates real values in the [-1, 1] range | `Eigen::DenseBase<Ty>::Random` for floating point types |
+| `Eigen::Rand::beta` | `Eigen::Rand::BetaGen` | float, double | | generates real values on a [beta distribution](https://en.wikipedia.org/wiki/Beta_distribution) |  |
+| `Eigen::Rand::cauchy` | `Eigen::Rand::CauchyGen` | float, double | yes | generates real values on the [Cauchy distribution](https://en.wikipedia.org/wiki/Cauchy_distribution). | `std::cauchy_distribution` |
+| `Eigen::Rand::chiSquared` | `Eigen::Rand::ChiSquaredGen` | float, double | | generates real values on a [chi-squared distribution](https://en.wikipedia.org/wiki/Chi-squared_distribution). | `std::chi_squared_distribution` |
+| `Eigen::Rand::exponential` | `Eigen::Rand::ExponentialGen` | float, double | yes | generates real values on an [exponential distribution](https://en.wikipedia.org/wiki/Exponential_distribution). | `std::exponential_distribution` |
+| `Eigen::Rand::extremeValue` | `Eigen::Rand::ExtremeValueGen` | float, double | yes | generates real values on an [extreme value distribution](https://en.wikipedia.org/wiki/Generalized_extreme_value_distribution). | `std::extreme_value_distribution` |
+| `Eigen::Rand::fisherF` | `Eigen::Rand::FisherFGen` | float, double | | generates real values on the [Fisher's F distribution](https://en.wikipedia.org/wiki/F_distribution). | `std::fisher_f_distribution` |
+| `Eigen::Rand::gamma` | `Eigen::Rand::GammaGen` | float, double | | generates real values on a [gamma distribution](https://en.wikipedia.org/wiki/Gamma_distribution). | `std::gamma_distribution` |
+| `Eigen::Rand::lognormal` | `Eigen::Rand::LognormalGen` | float, double | yes | generates real values on a [lognormal distribution](https://en.wikipedia.org/wiki/Lognormal_distribution). | `std::lognormal_distribution` |
+| `Eigen::Rand::normal` | `Eigen::Rand::StdNormalGen`, `Eigen::Rand::NormalGen` | float, double | yes | generates real values on a [normal distribution](https://en.wikipedia.org/wiki/Normal_distribution). | `std::normal_distribution` |
+| `Eigen::Rand::studentT` | `Eigen::Rand::StudentTGen` | float, double | yes | generates real values on the [Student's t distribution](https://en.wikipedia.org/wiki/Student%27s_t-distribution). | `std::student_t_distribution` |
+| `Eigen::Rand::uniformReal` | `Eigen::Rand::StdUniformRealGen`, `Eigen::Rand::UniformRealGen` | float, double | yes | generates real values in the `[0, 1)` range. | `std::generate_canonical` |
+| `Eigen::Rand::weibull` | `Eigen::Rand::WeibullGen` | float, double | yes | generates real values on the [Weibull distribution](https://en.wikipedia.org/wiki/Weibull_distribution). | `std::weibull_distribution` |
+
+* VoP indicates 'Vectorization over Parameters'. 
 
  @section list_of_supported_distribution_2 Random Distributions for Integer Types
 
-| Function | Generator | Scalar Type | Description | Equivalent to |
-|:---:|:---:|:---:|:---:|:---:|
-| `Eigen::Rand::binomial` | `Eigen::Rand::BinomialGen` | int | generates integers on a [binomial distribution](https://en.wikipedia.org/wiki/Binomial_distribution). | `std::binomial_distribution` |
-| `Eigen::Rand::discrete` | `Eigen::Rand::DiscreteGen` | int | generates random integers on a discrete distribution. | `std::discrete_distribution` |
-| `Eigen::Rand::geometric` | `Eigen::Rand::GeometricGen` | int | generates integers on a [geometric distribution](https://en.wikipedia.org/wiki/Geometric_distribution). | `std::geometric_distribution` |
-| `Eigen::Rand::negativeBinomial` | `Eigen::Rand::NegativeBinomialGen` | int | generates integers on a [negative binomial distribution](https://en.wikipedia.org/wiki/Negative_binomial_distribution). | `std::negative_binomial_distribution` |
-| `Eigen::Rand::poisson` | `Eigen::Rand::PoissonGen` | int | generates integers on the [Poisson distribution](https://en.wikipedia.org/wiki/Poisson_distribution). | `std::poisson_distribution` |
-| `Eigen::Rand::randBits` | `Eigen::Rand::RandbitsGen` | int | generates integers with random bits. | `Eigen::DenseBase<Ty>::Random` for integer types |
-| `Eigen::Rand::uniformInt` | `Eigen::Rand::UniformIntGen` | int | generates integers in the `[min, max]` range. | `std::uniform_int_distribution` |
+| Function | Generator | Scalar Type | VoP | Description | Equivalent to |
+|:---:|:---:|:---:|:---:|:---:|:---:|
+| `Eigen::Rand::bernoulli` | `Eigen::Rand::BernoulliGen` | int | yes | generates 0 or 1 on a [Bernoulli distribution](https://en.wikipedia.org/wiki/Bernoulli_distribution). | `std::bernoulli_distribution` |
+| `Eigen::Rand::binomial` | `Eigen::Rand::BinomialGen` | int | yes | generates integers on a [binomial distribution](https://en.wikipedia.org/wiki/Binomial_distribution). | `std::binomial_distribution` |
+| `Eigen::Rand::discrete` | `Eigen::Rand::DiscreteGen` | int | | generates random integers on a discrete distribution. | `std::discrete_distribution` |
+| `Eigen::Rand::geometric` | `Eigen::Rand::GeometricGen` | int | | generates integers on a [geometric distribution](https://en.wikipedia.org/wiki/Geometric_distribution). | `std::geometric_distribution` |
+| `Eigen::Rand::negativeBinomial` | `Eigen::Rand::NegativeBinomialGen` | int | | generates integers on a [negative binomial distribution](https://en.wikipedia.org/wiki/Negative_binomial_distribution). | `std::negative_binomial_distribution` |
+| `Eigen::Rand::poisson` | `Eigen::Rand::PoissonGen` | int | | generates integers on the [Poisson distribution](https://en.wikipedia.org/wiki/Poisson_distribution). | `std::poisson_distribution` |
+| `Eigen::Rand::randBits` | `Eigen::Rand::RandbitsGen` | int | | generates integers with random bits. | `Eigen::DenseBase<Ty>::Random` for integer types |
+| `Eigen::Rand::uniformInt` | `Eigen::Rand::UniformIntGen` | int | | generates integers in the `[min, max]` range. | `std::uniform_int_distribution` |
+
+* VoP indicates 'Vectorization over Parameters'. 
 
  @section list_of_distribution_3 Multivariate Random Distributions
 | Generator | Description | Equivalent to |
@@ -212,138 +261,40 @@
 |:---:|:---:|:---:|
 | `Eigen::Rand::Vmt19937_64` | a vectorized version of Mersenne Twister algorithm. It generates two 64bit random integers simultaneously with SSE2 and four integers with AVX2. | `std::mt19937_64` |
 | `Eigen::Rand::P8_mt19937_64` | a vectorized version of Mersenne Twister algorithm. Since it generates eight 64bit random integers simultaneously, the random values are the same regardless of architecture. | |
+
  * 
  * @page performance Performance
- * The following charts show the relative speed-up of EigenRand compared to Reference(C++ std or Eigen functions). Detailed results are below the charts.
-
- @section performance_1 Overview of Results at x86-64 Architecture
-
- \image html perf_no_vect.png
-
- \image html perf_sse2.png
-
- \image html perf_avx.png
-
- \image html perf_avx2.png
-
- \image html perf_mv_part1.png
-
- \image html perf_mv_part2.png
-
- * The following result is a measure of the time in seconds it takes to generate 1M random numbers. It shows the average of 20 times.
-
- @section performance_2 Overview of Results at ARM64 NEON (experimental)
-
- \image html perf_neon_v0.3.90.png
-
- \image html perf_mv_part1_neon_v0.3.90.png
-
- \image html perf_mv_part2_neon_v0.3.90.png
-
- * The following result is a measure of the time in seconds it takes to generate 1M random numbers. It shows the average of 20 times.
-
- @section performance_3 Intel(R) Xeon(R) Platinum 8171M CPU @ 2.60GHz (Ubuntu 16.04, gcc7.5)
-
-|  | C++ std (or Eigen) | EigenRand (No Vect.) | EigenRand (SSE2) | EigenRand (SSSE3) | EigenRand (AVX) | EigenRand (AVX2) |
-|---|---:|---:|---:|---:|---:|---:|
-| `balanced`* | 9.0 | 5.9 | 1.5 | 1.4 | 1.3 | 0.9 |
-| `balanced`(double)* | 8.7 | 6.4 | 3.3 | 2.9 | 1.7 | 1.7 |
-| `binomial(20, 0.5)` | 400.8 | 118.5 | 32.7 | 36.6 | 30.0 | 22.7 |
-| `binomial(50, 0.01)` | 71.7 | 22.5 | 7.7 | 8.3 | 7.9 | 6.6 |
-| `binomial(100, 0.75)` | 340.5 | 454.5 | 91.7 | 111.5 | 106.3 | 86.4 |
-| `cauchy` | 36.1 | 54.4 | 6.1 | 7.1 | 4.7 | 3.9 |
-| `chiSquared` | 80.5 | 249.5 | 64.6 | 58.0 | 29.4 | 28.8 |
-| `discrete`(int32) | - | 14.0 | 2.9 | 2.6 | 2.4 | 1.7 |
-| `discrete`(fp32) | - | 21.9 | 4.3 | 4.0 | 3.6 | 3.0 |
-| `discrete`(fp64) | 72.4 | 21.4 | 6.9 | 6.5 | 4.9 | 3.7 |
-| `exponential` | 31.0 | 25.3 | 5.5 | 5.3 | 3.3 | 2.9 |
-| `extremeValue` | 66.0 | 60.1 | 11.9 | 10.7 | 6.5 | 5.8 |
-| `fisherF(1, 1)` | 178.1 | 35.1 | 33.2 | 39.3 | 22.9 | 18.7 |
-| `fisherF(5, 5)` | 141.8 | 415.2 | 136.47 | 172.4 | 92.4 | 74.9 |
-| `gamma(0.2, 1)` | 207.8 | 211.4 | 54.6 | 51.2 | 26.9 | 27.0 |
-| `gamma(5, 3)` | 80.9 | 60.0 | 14.3 | 13.3 | 11.4 | 8.0 |
-| `gamma(10.5, 1)` | 81.1 | 248.6 | 63.3 | 58.5 | 29.2 | 28.4 |
-| `geometric` | 43.0 | 22.4 | 6.7 | 7.4 | 5.8 |  |
-| `lognormal` | 66.3 | 55.4 | 12.8 | 11.8 | 6.2 | 6.2 |
-| `negativeBinomial(10, 0.5)` | 312.0 | 301.4 | 82.9 | 100.6 | 95.3 | 77.9 |
-| `negativeBinomial(20, 0.25)` | 483.4 | 575.9 | 125.0 | 158.2 | 148.4 | 119.5 |
-| `normal(0, 1)` | 38.1 | 28.5 | 6.8 | 6.2 | 3.8 | 3.7 |
-| `normal(2, 3)` | 37.6 | 29.0 | 7.3 | 6.6 | 4.0 | 3.9 |
-| `poisson(1)` | 31.8 | 25.2 | 9.8 | 10.8 | 9.7 | 8.2 |
-| `poisson(16)` | 231.8 | 274.1 | 66.2 | 80.7 | 74.4 | 64.2 |
-| `randBits` | 5.2 | 5.4 | 1.4 | 1.3 | 1.1 | 1.0 |
-| `studentT(1)` | 122.7 | 120.1 | 15.3 | 19.2 | 12.6 | 9.4 |
-| `studentT(20)` | 102.2 | 111.1 | 15.4 | 19.2 | 12.2 | 9.4 |
-| `uniformInt(0~63)` | 22.4 | 4.7 | 1.7 | 1.6 | 1.4 | 1.1 |
-| `uniformInt(0~100k)` | 21.8 | 10.1 | 6.2 | 6.7 | 6.6 | 5.4 |
-| `uniformReal` | 12.9 | 5.7 | 1.4 | 1.2 | 1.4 | 0.7 |
-| `weibull` | 41.0 | 35.8 | 17.7 | 15.5 | 8.5 | 8.5 |
-
-* Since there is no equivalent class to `balanced` in C++11 std, we used Eigen::DenseBase::Random instead.
-
-|  | C++ std | EigenRand (No Vect.) | EigenRand (SSE2) | EigenRand (SSSE3) | EigenRand (AVX) | EigenRand (AVX2) |
-|---|---:|---:|---:|---:|---:|---:|
-| Mersenne Twister(int32) | 4.7 | 5.6 | 4.0 | 3.7 | 3.5 | 3.6 |
-| Mersenne Twister(int64) | 5.4 | 5.3 | 4.0 | 3.9 | 3.4 | 2.6 |
-
-|  | Python 3.6 + scipy 1.5.2 + numpy 1.19.2 | EigenRand (No Vect.) | EigenRand (SSE2) | EigenRand (SSSE3) | EigenRand (AVX) | EigenRand (AVX2) |
-|---|---:|---:|---:|---:|---:|---:|
-| `Dirichlet(4)` | 6.47 | 6.60 | 2.39 | 2.49 | 1.34 | 1.67 |
-| `Dirichlet(100)` | 75.95 | 189.97 | 66.60 | 72.11 | 38.86 | 34.98 |
-| `InvWishart(4)` | 140.18 | 7.62 | 4.21 | 4.54 | 3.58 | 3.39 |
-| `InvWishart(50)` | 1510.47 | 1737.4 | 697.39 | 733.69 | 604.59 | 554.006 |
-| `Multinomial(4, t=20)` | 3.32 | 4.12 | 0.95 | 1.06 | 1.00 | 1.03 |
-| `Multinomial(4, t=1000)` | 3.51 | 192.51 | 35.99 | 39.58 | 27.84 | 35.45 |
-| `Multinomial(100, t=20)` | 69.19 | 4.80 | 2.00 | 2.20 | 2.28 | 2.09 |
-| `Multinomial(100, t=1000)` | 139.74 | 179.43 | 49.48 | 56.19 | 40.78 | 43.18 |
-| `MvNormal(4)` | 2.32 | 0.96 | 0.36 | 0.37 | 0.25 | 0.30 |
-| `MvNormal(100)` | 49.09 | 57.18 | 17.17 | 18.51 | 10.82 | 11.03 |
-| `Wishart(4)` | 71.19 | 5.28 | 2.70 | 2.93 | 2.04 | 1.94 |
-| `Wishart(50)` | 1185.26 | 1360.49 | 492.91 | 517.44 | 359.03 | 324.60 |
-
- @section performance_4 AMD Ryzen 7 3700x CPU @ 3.60GHz (Windows 10, MSVC2017)
-
-|  | C++ std (or Eigen) | EigenRand (SSE2) | EigenRand (AVX) | EigenRand (AVX2) |
-|---|---:|---:|---:|---:|
-| `balanced`* | 20.8 | 1.9 | 2.0 | 1.4 |
-| `balanced`(double)* | 21.7 | 4.1 | 2.7 | 3.0 |
-| `binomial(20, 0.5)` | 416.0 | 27.7 | 28.9 | 29.1 |
-| `binomial(50, 0.01)` | 37.8 | 6.3 | 6.0 | 6.6 |
-| `binomial(100, 0.75)` | 309.1 | 72.4 | 66.0 | 67.0 |
-| `cauchy` | 42.2 | 4.8 | 5.1 | 2.7 |
-| `chiSquared` | 153.8 | 33.5 | 21.2 | 17.0 |
-| `discrete`(int32) | - | 2.4 | 2.3 | 2.5 |
-| `discrete`(fp32) | - | 2.6 | 2.3 | 3.5 |
-| `discrete`(fp64) | 55.8 | 5.1 | 4.7 | 4.3 |
-| `exponential` | 33.4 | 6.4 | 2.8 | 2.2 |
-| `extremeValue` | 39.4 | 7.8 | 4.6 | 4.0 |
-| `fisherF(1, 1)` | 103.9 | 25.3 | 14.9 | 11.7 |
-| `fisherF(5, 5)` | 295.7 | 85.5 | 58.3 | 44.8 |
-| `gamma(0.2, 1)` | 128.8 | 31.9 | 18.3 | 15.8 |
-| `gamma(5, 3)` | 156.1 | 9.7 | 8.0 | 5.0 |
-| `gamma(10.5, 1)` | 148.5 | 33.1 | 21.1 | 17.2 |
-| `geometric` | 27.1 | 6.6 | 4.3 | 4.1 |
-| `lognormal` | 104.0 | 6.6 | 4.7 | 3.5 |
-| `negativeBinomial(10, 0.5)` | 462.1 | 60.0 | 56.4 | 58.6 |
-| `negativeBinomial(20, 0.25)` | 357.6 | 84.5 | 80.6 | 78.4 |
-| `normal(0, 1)` | 48.8 | 4.2 | 3.7 | 2.3 |
-| `normal(2, 3)` | 48.8 | 4.5 | 3.8 | 2.4 |
-| `poisson(1)` | 46.4 | 7.9 | 7.4 | 8.2 |
-| `poisson(16)` | 192.4 | 43.2 | 40.4 | 40.9 |
-| `randBits` | 4.2 | 1.7 | 1.5 | 1.8 |
-| `studentT(1)` | 107.0 | 12.3 | 6.8 | 5.7 |
-| `studentT(20)` | 107.1 | 12.3 | 6.8 | 5.8 |
-| `uniformInt(0~63)` | 31.2 | 1.1 | 1.0 | 1.2 |
-| `uniformInt(0~100k)` | 27.7 | 5.6 | 5.6 | 5.4 |
-| `uniformReal` | 30.7 | 1.1 | 1.0 | 0.6 |
-| `weibull` | 46.5 | 10.6 | 6.4 | 5.2 |
-
-* Since there is no equivalent class to `balanced` in C++11 std, we used Eigen::DenseBase::Random instead.
-
-|  | C++ std | EigenRand (SSE2) | EigenRand (AVX) | EigenRand (AVX2) |
-|---|---:|---:|---:|---:|
-| Mersenne Twister(int32) | 5.0 | 3.4 | 3.4 | 3.3 |
-| Mersenne Twister(int64) | 5.1 | 3.9 | 3.9 | 3.3 |
+ * The following charts show the relative speed-up of EigenRand compared to references(equivalent functions of C++ std or Eigen for univariate distributions and Scipy for multivariate distributions).
+
+Since there is no equivalent class to `balanced` in C++11 std, we used Eigen::DenseBase::Random instead.
+
+Cases filled with orange are generators that are slower than reference functions.
+
+ @section performance_1 Windows 2019, MSVC 19.29.30147, Intel(R) Xeon(R) Platinum 8171M CPU, AVX2, Eigen 3.4.0
+ 
+ \image html perf_avx2_win.png width=80%
+ 
+ \image html perf_avx2_win_mv1.png width=80%
+ 
+ \image html perf_avx2_win_mv2.png width=80%
+
+ @section performance_2 Ubuntu 18.04, gcc 7.5.0, Intel(R) Xeon(R) Platinum 8370C CPU, AVX2, Eigen 3.4.0
+ 
+ \image html perf_avx2_ubu.png width=80%
+ 
+ \image html perf_avx2_ubu_mv1.png width=80%
+ 
+ \image html perf_avx2_ubu_mv2.png width=80%
+
+ @section performance_3 macOS Monterey 12.2.1, clang 13.1.6, Apple M1 Pro, NEON, Eigen 3.4.0
+ 
+ \image html perf_neon_mac.png width=80%
+ 
+ \image html perf_neon_mac_mv1.png width=80%
+ 
+ \image html perf_neon_mac_mv2.png width=80%
+
+ You can see the detailed numerical values used to plot the above charts on the <a href="https://github.com/bab2min/EigenRand/actions/workflows/release.yml" target="_blank">Action Results of GitHub repository</a>.
 
  * 
  */