main.cpp

#include <iostream>
#include <vector>
#include <cmath>

#include "BSpline.h"

int main()
{
    const int nx = 20;

    std::vector<double> x(nx);
    std::vector<double> y(nx);

    const double pi = std::acos(-1.0);

    for (int i = 0; i < nx; ++i)
    {
        x[i] = static_cast<double>(i) / (nx - 1);
        y[i] = std::sin(2.0 * pi * x[i]);
    }

    double wl = 0.2;
    int bc_type = 2;
    int num_nodes = 20;

    BSpline<double> spline(
        x.data(),
        nx,
        y.data(),
        wl,
        bc_type,
        num_nodes
    );

    for (int i = 0; i <= 10; ++i)
    {
        double u = static_cast<double>(i) / 10.0;

        std::cout
            << "u = " << u
            << ", y = " << spline.evaluate(u)
            << ", dy = " << spline.slope(u)
            << std::endl;
    }

    return 0;
}