function.cpp

#include "function.h"
#include "sample_functions.h"

std::ostream& operator<<(std::ostream& out, const numerics::Polynomial& p) {
    uint32_t lead_exp = p.coeff.size() - 1;
    out << "p(x)=";
    for (uint32_t curr = 0; curr < lead_exp + 1; curr++) {
        out << p.coeff[curr] << "x^" << (lead_exp - curr);
        if (curr < lead_exp) {
            out << " + ";
        }
    }
    return out;
}

// Initialize sample functions
numerics::Parabola_2D numerics::functions::parabola_2d =
    numerics::Parabola_2D();
numerics::Periodic_Wave_2D numerics::functions::periodic_wave_2d =
    numerics::Periodic_Wave_2D();

numerics::PolyWrapper numerics::functions::parabola_1d =
    numerics::PolyWrapper({1., 0., -1.4});
numerics::PolyWrapper numerics::functions::higher_parabola_1d =
    numerics::PolyWrapper({1., 0.0, -2.3, 0., -1.4});
numerics::PolyWrapper numerics::functions::ref_polynomial_1d =
    numerics::PolyWrapper({1., -3., -1., 11., -17., 16., -4.});