nagcpp/examples/ex_e02bb.cpp at master · numericalalgorithmsgroup/nagcpp · GitHub

History

44 lines (36 loc) · 1.27 KB

Raw

// example of using:

// nagcpp::fit::dim1_spline_eval (e02bb)

#include <iomanip>

#include <iostream>

#include "e02/nagcpp_e02bb.hpp"

int main(void) {

std::cout << "nagcpp::fit::dim1_spline_eval Example" << std::endl;

std::vector<double> lamda = {1.0, 1.0, 1.0, 1.0, 3.0, 6.0,

8.0, 9.0, 9.0, 9.0, 9.0};

std::vector<double> c = {1.0, 2.0, 4.0, 7.0, 6.0, 4.0, 3.0};

c.resize(lamda.size());

std::cout << std::fixed << std::setprecision(4);

// title for the results

std::cout << std::endl;

std::cout << " Value of" << std::endl;

std::cout << " R x cubic spline" << std::endl;

// evaluate the spline at m equally spaced points in the interval

// lamda[3] <= x <= lamda[ncap+3]

size_t m = 9;

size_t ncap = lamda.size() - 7;

double a = lamda[3];

double b = lamda[ncap + 3];

for (size_t r = 0; r < m; ++r) {

double x = ((m - r - 1) * a + r * b) / static_cast<double>(m - 1);

double s;

try {

nagcpp::fit::dim1_spline_eval(lamda, c, x, s);

} catch (nagcpp::error_handler::Exception &e) {

std::cout << e.msg << std::endl;

return 1;

}

std::cout << std::setw(3) << r << std::setw(14) << x << std::setw(14) << s

<< std::endl;

}

return 0;

}

nagcpp/examples/ex_e02bb.cpp at master · numericalalgorithmsgroup/nagcpp · GitHub

Uh oh!

Expand file tree

Breadcrumbs

ex_e02bb.cpp

Latest commit

History

Breadcrumbs

ex_e02bb.cpp

File metadata and controls