Method eval

src/Region/CubicSplineCurve.cpp:286–291

}

double CubicSplineCurve::eval(double a) const noexcept {
    const std::size_t i = locate(a);
    const double h = a_[i + 1] - a_[i];
    const double t = (a - a_[i]) / h;
    return segment_eval(t, h, b_[i], b_[i + 1], M_[i], M_[i + 1]);
}

double CubicSplineCurve::eval_da(double a) const noexcept {
    const std::size_t i = locate(a);