30 std::string StringFormat(const std::string& format, Args... args) {  in StringFormat()
31     auto size = std::snprintf(nullptr, 0, format.c_str(), args...);  in StringFormat()
33     std::vector<char> buffer(size + 1); // Add 1 for terminating null byte  in StringFormat()
34     std::snprintf(buffer.data(), buffer.size(), format.c_str(), args...);  in StringFormat()
35     return std::string(buffer.data(), size); // Exclude the terminating null byte  in StringFormat()
38 double polyEval(std::vector<double>& coeffs, double x) {  in polyEval()
48 std::vector<double> polyIntegrate(std::vector<double>& coeffs, double coi = 0.0) {  in polyIntegrate()
49     std::vector<double> integratedCoeffs(coeffs.size() + 1);  in polyIntegrate()
57 std::vector<double> polyFit(std::vector<double>& rates, std::vector<double>& qualities, int order) {  in polyFit()
71     std::vector<double> coeffs(order + 1);  in polyFit()
78 double getAvgImprovement(std::vector<double>& xA, std::vector<double>& yA, std::vector<double>& xB,  in getAvgImprovement()
79                          std::vector<double>& yB, int order) {  in getAvgImprovement()
80     std::vector<double> coeffsA = polyFit(xA, yA, order);  in getAvgImprovement()
81     std::vector<double> coeffsB = polyFit(xB, yB, order);  in getAvgImprovement()
82     std::vector<double> integratedCoeffsA = polyIntegrate(coeffsA);  in getAvgImprovement()
83     std::vector<double> integratedCoeffsB = polyIntegrate(coeffsB);  in getAvgImprovement()
84     double minX = std::max(*std::min_element(xA.begin(), xA.end()),  in getAvgImprovement()
85                            *std::min_element(xB.begin(), xB.end()));  in getAvgImprovement()
86     double maxX = std::min(*std::max_element(xA.begin(), xA.end()),  in getAvgImprovement()
87                            *std::max_element(xB.begin(), xB.end()));  in getAvgImprovement()
98     std::string msg;  in nativeGetBDRate()
106         std::vector<double> ratesA(len[0]);  in nativeGetBDRate()
108         std::vector<double> ratesB(len[0]);  in nativeGetBDRate()
110         std::vector<double> qualitiesA(len[0]);  in nativeGetBDRate()
112         std::vector<double> qualitiesB(len[0]);  in nativeGetBDRate()
116             ratesA[i] = std::log(ratesA[i]);  in nativeGetBDRate()
117             ratesB[i] = std::log(ratesB[i]);  in nativeGetBDRate()
128             bdRate = (std::exp(bdRate) - 1) * 100;  in nativeGetBDRate()