[Math] Replace VecCore functions with std::simd functions

Author: guitargeek

hist/hist/inc/TF1.h

@@ -841,7 +841,7 @@ inline double TF1::EvalParVec(const Double_t *data, const Double_t *params)
       //    res = GetSave(x);
       return TMath::SignalingNaN();
    }
-   return vecCore::Get<ROOT::Double_v>(res, 0);
+   return res[0];
 }
 #endif
 

hist/hist/src/TFormula.cxx

@@ -3134,7 +3134,7 @@ Double_t TFormula::EvalPar(const Double_t *x,const Double_t *params) const
 
    if (fNdim == 0 || !x) {
       ROOT::Double_v ret =  DoEvalVec(nullptr, params);
-      return vecCore::Get( ret, 0 );
+      return ret[0];
    }
 
     // otherwise, regular Double_t inputs on a vectorized function
@@ -3150,15 +3150,15 @@ Double_t TFormula::EvalPar(const Double_t *x,const Double_t *params) const
          xvec[i] = x[i];
 
       ROOT::Double_v ans = DoEvalVec(xvec.data(), params);
-      return vecCore::Get(ans, 0);
+      return ans[0];
    }
    // allocating a vector is much slower (we do only for dim > 4)
    std::vector<ROOT::Double_v> xvec(fNdim);
    for (int i = 0; i < fNdim; i++)
       xvec[i] = x[i];
 
    ROOT::Double_v ans = DoEvalVec(xvec.data(), params);
-   return  vecCore::Get(ans, 0);
+   return ans[0];
 
 #else
    // this should never happen, because fVectorized can only be set true with
@@ -3393,16 +3393,16 @@ ROOT::Double_v TFormula::EvalParVec(const ROOT::Double_v *x, const Double_t *par
    if (gDebug)
       Info("EvalPar", "Function is not vectorized - converting ROOT::Double_v into Double_t and back");
 
-   const int vecSize = vecCore::VectorSize<ROOT::Double_v>();
+   const int vecSize = ROOT::Double_v::size();
    std::vector<Double_t>  xscalars(vecSize*fNdim);
 
    for (int i = 0; i < vecSize; i++)
       for (int j = 0; j < fNdim; j++)
-         xscalars[i*fNdim+j] = vecCore::Get(x[j],i);
+         xscalars[i * fNdim + j] = x[j][i];
 
    ROOT::Double_v answers(0.);
    for (int i = 0; i < vecSize; i++)
-      vecCore::Set(answers, i, DoEval(&xscalars[i*fNdim], params));
+      answers[i] = DoEval(&xscalars[i * fNdim], params);
 
    return answers;
 }

math/mathcore/inc/Fit/FitData.h

@@ -358,9 +358,9 @@ namespace ROOT {
          static unsigned VectorPadding(unsigned dataSize)
          {
             unsigned padding = 0;
-            unsigned modP = (dataSize) % vecCore::VectorSize<ROOT::Double_v>();
+            unsigned modP = (dataSize) % ROOT::Double_v::size();
             if (modP > 0)
-               padding = vecCore::VectorSize<ROOT::Double_v>() - modP;
+               padding = ROOT::Double_v::size() - modP;
             return padding;
          }
 #else