Inheritance diagram for CubicBspline< T, GRIDTYPE, BC >:

Collaboration diagram for CubicBspline< T, GRIDTYPE, BC >:

Public Types
using	point_type = typename CubicBsplineGrid< T, GRIDTYPE, BC >::point_type

using	value_type = typename CubicBsplineGrid< T, GRIDTYPE, BC >::value_type

using	container_type = typename CubicBsplineGrid< T, GRIDTYPE, BC >::container_type

Public Member Functions
	CubicBspline ()
	default constructor More...

void	Init (point_type start, point_type end, const container_type &datain, bool closed)

void	Init (point_type start, point_type end, const container_type &datain, bool closed, T yp1, T ypn)

value_type	getValue (point_type x)

value_type	getDeriv (point_type x)

value_type	getDeriv2 (point_type x)

value_type	getDeriv3 (point_type x)

value_type	operator() (T x)

value_type	splint (T x)

value_type	splint (point_type x, value_type &dy, value_type &d2y)

value_type	interpolate (value_type p0, value_type p1, value_type p2, value_type p3, value_type &dy, value_type &d2y)

value_type	interpolate0 (value_type p0, value_type p1, value_type p2, value_type p3)

value_type	interpolate1 (value_type p0, value_type p1, value_type p2, value_type p3)

value_type	interpolate2 (value_type p0, value_type p1, value_type p2, value_type p3)

Public Attributes
int	i0
	index of current grid point More...

int	i1

int	i2

int	i3

value_type	OffSet
	constant shift More...

point_type	A [16]
	coefficients More...

point_type	dA [12]

point_type	d2A [8]

point_type	d3A [4]

container_type	P
	The control points. More...

Detailed Description

template<class T, unsigned GRIDTYPE, unsigned BC>
struct CubicBspline< T, GRIDTYPE, BC >

Definition at line 23 of file CubicBspline.h.

Member Typedef Documentation

◆ container_type

using container_type = typename CubicBsplineGrid<T, GRIDTYPE, BC>::container_type

Definition at line 27 of file CubicBspline.h.

◆ point_type

using point_type = typename CubicBsplineGrid<T, GRIDTYPE, BC>::point_type

Definition at line 25 of file CubicBspline.h.

◆ value_type

using value_type = typename CubicBsplineGrid<T, GRIDTYPE, BC>::value_type

Definition at line 26 of file CubicBspline.h.

Constructor & Destructor Documentation

◆ CubicBspline()

CubicBspline ( )

inline

default constructor

Initialize linear coefficients

Definition at line 46 of file CubicBspline.h.

References CubicBspline< T, GRIDTYPE, BC >::A, CubicBspline< T, GRIDTYPE, BC >::d2A, CubicBspline< T, GRIDTYPE, BC >::d3A, and CubicBspline< T, GRIDTYPE, BC >::dA.

                         : OffSet(0.0)
   {
     A[0]   = -1.0 / 6.0;
     A[1]   = 3.0 / 6.0;
     A[2]   = -3.0 / 6.0;
     A[3]   = 1.0 / 6.0;
     A[4]   = 3.0 / 6.0;
     A[5]   = -6.0 / 6.0;
     A[6]   = 3.0 / 6.0;
     A[7]   = 0.0 / 6.0;
     A[8]   = -3.0 / 6.0;
     A[9]   = 0.0 / 6.0;
     A[10]  = 3.0 / 6.0;
     A[11]  = 0.0 / 6.0;
     A[12]  = 1.0 / 6.0;
     A[13]  = 4.0 / 6.0;
     A[14]  = 1.0 / 6.0;
     A[15]  = 0.0 / 6.0;
     dA[0]  = -0.5;
     dA[1]  = 1.5;
     dA[2]  = -1.5;
     dA[3]  = 0.5;
     dA[4]  = 1.0;
     dA[5]  = -2.0;
     dA[6]  = 1.0;
     dA[7]  = 0.0;
     dA[8]  = -0.5;
     dA[9]  = 0.0;
     dA[10] = 0.5;
     dA[11] = 0.0;
     d2A[0] = -1.0;
     d2A[1] = 3.0;
     d2A[2] = -3.0;
     d2A[3] = 1.0;
     d2A[4] = 1.0;
     d2A[5] = -2.0;
     d2A[6] = 1.0;
     d2A[7] = 0.0;
     d3A[0] = -1.0;
     d3A[1] = 3.0;
     d3A[2] = -3.0;
     d3A[3] = 1.0;
   }

Member Function Documentation

◆ getDeriv()

value_type getDeriv ( point_type x )

inline

Definition at line 108 of file CubicBspline.h.

References CubicBspline< T, GRIDTYPE, BC >::i0, CubicBspline< T, GRIDTYPE, BC >::interpolate1(), CubicBspline< T, GRIDTYPE, BC >::OffSet, and CubicBspline< T, GRIDTYPE, BC >::P.

   {
     if (this->getGridPoint(x, i0))
       return interpolate1(P[i0], P[i0 + 1], P[i0 + 2], P[i0 + 3]);
     else
       return OffSet;
   }

◆ getDeriv2()

value_type getDeriv2 ( point_type x )

inline

Definition at line 116 of file CubicBspline.h.

References CubicBspline< T, GRIDTYPE, BC >::i0, CubicBspline< T, GRIDTYPE, BC >::interpolate2(), CubicBspline< T, GRIDTYPE, BC >::OffSet, and CubicBspline< T, GRIDTYPE, BC >::P.

   {
     if (this->getGridPoint(x, i0))
       return interpolate2(P[i0], P[i0 + 1], P[i0 + 2], P[i0 + 3]);
     else
       return OffSet;
   }

◆ getDeriv3()

value_type getDeriv3 ( point_type x )

inline

Definition at line 124 of file CubicBspline.h.

References CubicBspline< T, GRIDTYPE, BC >::d2A, CubicBspline< T, GRIDTYPE, BC >::i0, CubicBspline< T, GRIDTYPE, BC >::i1, CubicBspline< T, GRIDTYPE, BC >::i2, CubicBspline< T, GRIDTYPE, BC >::i3, CubicBspline< T, GRIDTYPE, BC >::OffSet, and CubicBspline< T, GRIDTYPE, BC >::P.

   {
     if (this->getGridPoint(x, i0))
       return GridDeltaInv * GridDeltaInv * GridDeltaInv *
           (tp[3] * (d2A[0] * P[i0] + d2A[1] * P[i1] + d2A[2] * P[i2] + d2A[3] * P[i3]));
     else
       return OffSet;
   }

◆ getValue()

value_type getValue ( point_type x )

inline

Definition at line 100 of file CubicBspline.h.

References CubicBspline< T, GRIDTYPE, BC >::i0, CubicBspline< T, GRIDTYPE, BC >::interpolate0(), CubicBspline< T, GRIDTYPE, BC >::OffSet, and CubicBspline< T, GRIDTYPE, BC >::P.

Referenced by CubicBspline< T, GRIDTYPE, BC >::operator()().

   {
     if (this->getGridPoint(x, i0))
       return interpolate0(P[i0], P[i0 + 1], P[i0 + 2], P[i0 + 3]);
     else
       return OffSet;
   }

◆ Init() [1/2]

void Init	(	point_type	start,
		point_type	end,
		const container_type &	datain,
		bool	closed
	)

inline

Definition at line 90 of file CubicBspline.h.

References CubicBspline< T, GRIDTYPE, BC >::P.

Referenced by CubicSplineBasisSet< RT >::reset().

   {
     this->spline(start, end, datain, P, closed);
   }

◆ Init() [2/2]

void Init	(	point_type	start,
		point_type	end,
		const container_type &	datain,
		bool	closed,
		T	yp1,
		T	ypn
	)

inline

Definition at line 95 of file CubicBspline.h.

References CubicBspline< T, GRIDTYPE, BC >::P.

   {
     this->spline(start, end, yp1, ypn, datain, P);
   }

◆ interpolate()

value_type interpolate	(	value_type	p0,
		value_type	p1,
		value_type	p2,
		value_type	p3,
		value_type &	dy,
		value_type &	d2y
	)

inline

Definition at line 169 of file CubicBspline.h.

Referenced by CubicBspline< T, GRIDTYPE, BC >::splint().

   {
     dy = GridDeltaInv *
         (tp[1] * (-0.5 * p0 + 1.5 * p1 - 1.5 * p2 + 0.5 * p3) + tp[2] * (p0 - 2.0 * p1 + p2) +
          tp[3] * (-0.5 * p0 + 0.5 * p2));
     d2y = GridDeltaInv2 * (tp[2] * (-p0 + 3.0 * p1 - 3.0 * p2 + p3) + tp[3] * (p0 - 2.0 * p1 + p2));
     const point_type onesixth = 1.0 / 6.0;
     return onesixth *
         (tp[0] * (-p0 + 3.0 * p1 - 3.0 * p2 + p3) + tp[1] * (3.0 * p0 - 6.0 * p1 + 3.0 * p2) +
          tp[2] * (-3.0 * p0 + 3.0 * p2) + tp[3] * (p0 + 4.0 * p1 + p2));
   }

◆ interpolate0()

value_type interpolate0	(	value_type	p0,
		value_type	p1,
		value_type	p2,
		value_type	p3
	)

inline

Definition at line 186 of file CubicBspline.h.

Referenced by CubicBspline< T, GRIDTYPE, BC >::getValue(), and CubicBspline< T, GRIDTYPE, BC >::splint().

   {
     const point_type onesixth = 1.0 / 6.0;
     return onesixth *
         (tp[0] * (-p0 + 3.0 * p1 - 3.0 * p2 + p3) + tp[1] * (3.0 * p0 - 6.0 * p1 + 3.0 * p2) +
          tp[2] * (-3.0 * p0 + 3.0 * p2) + tp[3] * (p0 + 4.0 * p1 + p2));
   }

◆ interpolate1()

value_type interpolate1	(	value_type	p0,
		value_type	p1,
		value_type	p2,
		value_type	p3
	)

inline

Definition at line 194 of file CubicBspline.h.

Referenced by CubicBspline< T, GRIDTYPE, BC >::getDeriv().

   {
     return GridDeltaInv *
         (tp[1] * (-0.5 * p0 + 1.5 * p1 - 1.5 * p2 + 0.5 * p3) + tp[2] * (p0 - 2.0 * p1 + p2) +
          tp[3] * (-0.5 * p0 + 0.5 * p2));
   }

◆ interpolate2()

value_type interpolate2	(	value_type	p0,
		value_type	p1,
		value_type	p2,
		value_type	p3
	)

inline

Definition at line 201 of file CubicBspline.h.

Referenced by CubicBspline< T, GRIDTYPE, BC >::getDeriv2().

   {
     return GridDeltaInv2 * (tp[2] * (-p0 + 3.0 * p1 - 3.0 * p2 + p3) + tp[3] * (p0 - 2.0 * p1 + p2));
   }

◆ operator()()

value_type operator() ( T x )

inline

Definition at line 133 of file CubicBspline.h.

References CubicBspline< T, GRIDTYPE, BC >::getValue().

133 { return getValue(x); }

CubicBspline::getValue

value_type getValue(point_type x)

Definition: CubicBspline.h:100

◆ splint() [1/2]

value_type splint ( T x )

inline

Definition at line 135 of file CubicBspline.h.

References CubicBspline< T, GRIDTYPE, BC >::i0, CubicBspline< T, GRIDTYPE, BC >::interpolate0(), CubicBspline< T, GRIDTYPE, BC >::OffSet, and CubicBspline< T, GRIDTYPE, BC >::P.

Referenced by CubicSplineBasisSet< RT >::df(), CubicSplineBasisSet< RT >::evaluate(), CubicSplineBasisSet< RT >::f(), and CubicSplineBasisSet< RT >::print().

   {
     if (this->getGridPoint(x, i0))
       return interpolate0(P[i0], P[i0 + 1], P[i0 + 2], P[i0 + 3]);
     else
       return OffSet;
   }

◆ splint() [2/2]

value_type splint	(	point_type	x,
		value_type &	dy,
		value_type &	d2y
	)

inline

Definition at line 143 of file CubicBspline.h.

References CubicBspline< T, GRIDTYPE, BC >::i0, CubicBspline< T, GRIDTYPE, BC >::interpolate(), CubicBspline< T, GRIDTYPE, BC >::OffSet, and CubicBspline< T, GRIDTYPE, BC >::P.

   {
     if (this->getGridPoint(x, i0))
     {
       return interpolate(P[i0], P[i0 + 1], P[i0 + 2], P[i0 + 3], dy, d2y);
     }
     else
     {
       dy  = 0.0;
       d2y = 0.0;
       return OffSet;
     }
     //Too slow
     //dy= GridDeltaInv *
     //  (tp[1]*(dA[0]*P[i0]+dA[1]*P[i1]+dA[ 2]*P[i2]+dA[ 3]*P[i3])+
     //   tp[2]*(dA[4]*P[i0]+dA[5]*P[i1]+dA[ 6]*P[i2]+dA[ 7]*P[i3])+
     //   tp[3]*(dA[8]*P[i0]+dA[9]*P[i1]+dA[10]*P[i2]+dA[11]*P[i3]));
     //d2y=GridDeltaInv2 *
     //  (tp[2]*(d2A[0]*P[i0]+d2A[1]*P[i1]+d2A[2]*P[i2]+d2A[3]*P[i3])+
     //   tp[3]*(d2A[4]*P[i0]+d2A[5]*P[i1]+d2A[6]*P[i2]+d2A[7]*P[i3]));
     //return
     //  tp[0]*(A[ 0]*P[i0]+A[ 1]*P[i1]+A[ 2]*P[i2]+A[ 3]*P[i3])+
     //  tp[1]*(A[ 4]*P[i0]+A[ 5]*P[i1]+A[ 6]*P[i2]+A[ 7]*P[i3])+
     //  tp[2]*(A[ 8]*P[i0]+A[ 9]*P[i1]+A[10]*P[i2]+A[11]*P[i3])+
     //  tp[3]*(A[12]*P[i0]+A[13]*P[i1]+A[14]*P[i2]+A[15]*P[i3]);
   }