math/example/formatter_b_spline_output.cpp at c8c07a5325aad3963fcc8aa2162cb6d5d38eaba7 · boostorg/math · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
//  (C) Copyright John Maddock 2019.
//  Use, modification and distribution are subject to the
//  Boost Software License, Version 1.0. (See accompanying file
//  LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

#include <boost/math/tools/formatting.hpp>
#include <iostream>
#include <vector>
#include <boost/math/interpolators/cardinal_cubic_b_spline.hpp>
#include <boost/math/interpolators/cardinal_quadratic_b_spline.hpp>
#include <boost/math/interpolators/cardinal_quintic_b_spline.hpp>

int main() {
   std::vector<double> v(10);
   for (size_t i = 0; i < v.size(); ++i)
   {
      v[i] = static_cast<double>(5 + 6 * i);
   }

   double step = 0.1;
   double a = 5;
   {
      std::cout << "Quadratic B Spline:\n\n";
      boost::math::interpolators::cardinal_quadratic_b_spline<double> spline(v.data(), v.size(), a, step);

      boost::math::tools::text_printer printer(std::cout);
      printer << spline << std::endl << std::endl;

      boost::math::tools::latex_printer lprinter(std::cout);
      lprinter << spline << std::endl << std::endl;

      boost::math::tools::html_printer hprinter(std::cout);
      hprinter << spline << std::endl << std::endl;

      boost::math::tools::docbook_printer dprinter(std::cout);
      dprinter << spline << std::endl << std::endl;
   }
   {
      std::cout << "Cubic B Spline:\n\n";
      boost::math::interpolators::cardinal_cubic_b_spline<double> spline(v.data(), v.size(), a, step);

      boost::math::tools::text_printer printer(std::cout);
      printer << spline << std::endl << std::endl;

      boost::math::tools::latex_printer lprinter(std::cout);
      lprinter << spline << std::endl << std::endl;

      boost::math::tools::html_printer hprinter(std::cout);
      hprinter << spline << std::endl << std::endl;

      boost::math::tools::docbook_printer dprinter(std::cout);
      dprinter << spline << std::endl << std::endl;
   }
   {
      std::cout << "Quintic B Spline:\n\n";
      boost::math::interpolators::cardinal_quintic_b_spline<double> spline(v.data(), v.size(), a, step);

      boost::math::tools::text_printer printer(std::cout);
      printer << spline << std::endl << std::endl;

      boost::math::tools::latex_printer lprinter(std::cout);
      lprinter << spline << std::endl << std::endl;

      boost::math::tools::html_printer hprinter(std::cout);
      hprinter << spline << std::endl << std::endl;

      boost::math::tools::docbook_printer dprinter(std::cout);
      dprinter << spline << std::endl << std::endl;
   }
   return 0;
}