Use affinity_partitioner

Author: TLCFEM

patch_size.patch

@@ -1,5 +1,5 @@
 diff --git a/Eigen/src/SVD/BDCSVD.h b/Eigen/src/SVD/BDCSVD.h
-index 6b85d1d5a..81aa7fbf5 100644
+index 6b85d1d5a..f7d023c7a 100644
 --- a/Eigen/src/SVD/BDCSVD.h
 +++ b/Eigen/src/SVD/BDCSVD.h
 @@ -30,6 +30,10 @@
@@ -8,7 +8,7 @@ index 6b85d1d5a..81aa7fbf5 100644
 
 +#include <tbb/parallel_reduce.h>
 +
-+tbb::affinity_partitioner ap;
++static tbb::affinity_partitioner ap;
 +
  #ifdef EIGEN_BDCSVD_DEBUG_VERBOSE
  #include <iostream>

src/Gauss.hpp

@@ -28,8 +28,6 @@ inline mpreal MP_PI_HALF{MP_PI / 2};
 
 extern Config config;
 
-extern tbb::affinity_partitioner ap;
-
 class Evaluation {
     const mpreal ONE = mpreal(1, config.precision_bits);
     const mpreal TWO = mpreal(2, config.precision_bits);
@@ -68,6 +66,8 @@ class Evaluation {
 };
 
 class LegendrePolynomial {
+    static tbb::affinity_partitioner ap;
+
     const mpreal ONE = mpreal(1, config.precision_bits);
     const mpreal TWO = mpreal(2, config.precision_bits);
 
@@ -105,6 +105,8 @@ class LegendrePolynomial {
     [[nodiscard]] mpreal weight(const int i) const { return this->_w[i]; }
 };
 
+tbb::affinity_partitioner LegendrePolynomial::ap{};
+
 class Quadrature {
     LegendrePolynomial poly;