Function root_sum_square

src/CPnumerics.cpp:6–12

#include <unsupported/Eigen/Polynomials>

double root_sum_square(const std::vector<double>& x) {
    double sum = 0;
    for (double i : x) {
        sum += pow(i, 2);
    }
    return sqrt(sum);
}
double interp1d(const std::vector<double>* x, const std::vector<double>* y, double x0) {
    std::size_t i = 0, L = 0, R = 0, M = 0;
    L = 0;