#include <ScalarFunctionsImpl.h>

Inheritance diagram for CubicSplineScalarFunction:

Collaboration diagram for CubicSplineScalarFunction:

Public Member Functions
	CubicSplineScalarFunction (const std::vector< doublereal > &y_i, const std::vector< doublereal > &x_i, bool doNotExtrapolate=false)

virtual	~CubicSplineScalarFunction (void)

virtual doublereal	operator() (const doublereal x) const

virtual doublereal	ComputeDiff (const doublereal t, const integer order=1) const

Public Member Functions inherited from DifferentiableScalarFunction
virtual	~DifferentiableScalarFunction ()

Public Member Functions inherited from BasicScalarFunction
virtual	~BasicScalarFunction ()

Private Attributes
std::vector< doublereal >	Y_i

std::vector< doublereal >	X_i

std::vector< doublereal >	b

std::vector< doublereal >	c

std::vector< doublereal >	d

bool	doNotExtrapolate

Detailed Description

Definition at line 120 of file ScalarFunctionsImpl.h.

Constructor & Destructor Documentation

CubicSplineScalarFunction::CubicSplineScalarFunction	(	const std::vector< doublereal > &	y_i,
		const std::vector< doublereal > &	x_i,
		bool	doNotExtrapolate = `false`
	)

Definition at line 458 of file ScalarFunctionsImpl.cc.

References ASSERTMSGBREAK, b, c, d, MBDYN_EXCEPT_ARGS, spline(), X_i, and Y_i.

 : Y_i(y_i), X_i(x_i),
 doNotExtrapolate(doNotExtrapolate)
 {
         ASSERTMSGBREAK(Y_i.size() == X_i.size(),
                 "CubicSplineScalarFunction error, Y_i.size() != X_i.size()");
         std::vector<doublereal>::iterator xi, xe;
         xi = X_i.begin();
         xe = X_i.end() - 1;
         for (unsigned i = 0; xi != xe; ++xi, ++i) {
                 if (*xi >= *(xi + 1)) {
                         silent_cerr("CubicSplineScalarFunction error, "
                                 "X is not ordered: "
                                 "X[" << i << "]=" << *xi
                                 << " is not less than "
                                 "X[" << i + 1 << "]=" << *(xi + 1)
                                 << std::endl);
                         throw DataManager::ErrGeneric(MBDYN_EXCEPT_ARGS);
                 }
         }
         spline(X_i, Y_i, b, c, d);
 }

Here is the call graph for this function:

CubicSplineScalarFunction::~CubicSplineScalarFunction ( void )

virtual

Definition at line 484 of file ScalarFunctionsImpl.cc.

References NO_OP.

 {
         NO_OP;
 }

Member Function Documentation

doublereal CubicSplineScalarFunction::ComputeDiff	(	const doublereal	t,
		const integer	order = `1`
	)		const

virtual

Implements DifferentiableScalarFunction.

Definition at line 507 of file ScalarFunctionsImpl.cc.

References ASSERTMSGBREAK, b, c, d, operator()(), seval(), X_i, and Y_i.

 {
         ASSERTMSGBREAK(order >=0, "Error in CubicSplineScalarFunction::ComputeDiff, order<0");
         switch (order) {
         case 0: 
                 return this->operator()(x);
 
         case 1: 
                 return seval(x, X_i, Y_i, b, c, d, 1);
 
         case 2: 
                 return seval(x, X_i, Y_i, b, c, d, 2);
 
         case 3: 
                 return seval(x, X_i, Y_i, b, c, d, 3);
 
         default:
                 return 0.;
         }
 }

Here is the call graph for this function:

doublereal CubicSplineScalarFunction::operator() ( const doublereal x ) const

virtual

Implements DifferentiableScalarFunction.

Definition at line 490 of file ScalarFunctionsImpl.cc.

References b, c, d, doNotExtrapolate, seval(), X_i, and Y_i.

Referenced by ComputeDiff().

 {
         if (doNotExtrapolate) {
                 if (x <= X_i[0]) {
                         return Y_i[0];
                 }
 
                 int s = X_i.size() - 1;
                 if (x >= X_i[s]) {
                         return Y_i[s];
                 }
         }
 
         return seval(x, X_i, Y_i, b, c, d);
 }