cubic_spline.cpp
Go to the documentation of this file.
1
12 /*****************************************************************************
13 ** Includes
14 *****************************************************************************/
15
17 #include "../../include/ecl/geometry/cubic_spline.hpp"
18
19 /*****************************************************************************
20 ** Namespaces
21 *****************************************************************************/
22
23 namespace ecl {
24
25 /*****************************************************************************
26 ** Implementation
27 *****************************************************************************/
28
29 double CubicSpline::operator()(const double &x) const {
31  int index = 0;
32  while ( x > discretised_domain[index+1] ) {
33  ++index;
34  }
35  return cubic_polynomials[index](x);
36 }
37
38 double CubicSpline::derivative(double x) const {
40  int index = 0;
41  while ( x > discretised_domain[index+1] ) {
42  ++index;
43  }
44  return cubic_polynomials[index].derivative()(x);
45 }
46
47 double CubicSpline::dderivative(double x) const {
49  int index = 0;
50  while ( x > discretised_domain[index+1] ) {
51  ++index;
52  }
53  return cubic_polynomials[index].derivative().derivative()(x);
54 }
55
56 } // namespace ecl
Array< CubicPolynomial > cubic_polynomials
Embedded control libraries.
#define LOC
double dderivative(double x) const
Spline second derivative.
double derivative(double x) const
Spline derivative.
#define ecl_assert_throw(expression, exception)
OutOfRangeError
Array< double > discretised_domain
double operator()(const double &x) const
Spline function.
reference back()
reference front()

ecl_geometry
Author(s): Daniel Stonier
autogenerated on Mon Feb 28 2022 22:18:49