Try out std::simd backend of VecCore

Author: guitargeek

math/mathcore/inc/Math/Types.h

@@ -24,17 +24,15 @@
 #pragma GCC diagnostic pop
 #endif
 
+#define VECCORE_ENABLE_STD_SIMD
 #include <VecCore/VecCore>
+#undef VECCORE_ENABLE_STD_SIMD
 
 namespace ROOT {
 
 namespace Internal {
    using ScalarBackend = vecCore::backend::Scalar;
-#ifdef VECCORE_ENABLE_VC
-   using VectorBackend = vecCore::backend::VcVector;
-#else
-   using VectorBackend = vecCore::backend::Scalar;
-#endif
+   using VectorBackend = vecCore::backend::SIMDNative;
 }
    using Float_v  = typename Internal::VectorBackend::Float_v;
    using Double_v = typename Internal::VectorBackend::Double_v;

math/mathcore/inc/Math/Util.h

@@ -79,7 +79,7 @@ class TimingScope {
    inline T EvalLog(T x) {
       static const T epsilon = T(2.0 * std::numeric_limits<double>::min());
 #ifdef R__HAS_VECCORE
-      T logval = vecCore::Blend<T>(x <= epsilon, x / epsilon + std::log(epsilon) - T(1.0), std::log(x));
+      T logval = vecCore::Blend<T>(x <= epsilon, x / epsilon + vecCore::math::Log(epsilon) - T(1.0), vecCore::math::Log(x));
 #else
       T logval = x <= epsilon ? x / epsilon + std::log(epsilon) - T(1.0) : std::log(x);
 #endif