#include <ThirdOrderDelayODE.h>

Public Member Functions
void	calculate (double tau, double inputCurrent, double dIdt, double d2Id2t, const double x, unsigned lenX, unsigned firstNode, double derivative) const

unsigned	getPadeColumn () const

unsigned	getPadeRow () const

unsigned	nParameters () const

void	setParameters (const double *pars, unsigned nPars)

	ThirdOrderDelayODE (unsigned, unsigned)

Private Attributes
double	a_

double	b_

double	c_

Detailed Description

Definition at line 9 of file ThirdOrderDelayODE.h.

Constructor & Destructor Documentation

◆ ThirdOrderDelayODE()

ThirdOrderDelayODE::ThirdOrderDelayODE	(	unsigned	,
		unsigned
	)

inline

Definition at line 11 of file ThirdOrderDelayODE.h.

11 : a_(1.0) {}

Member Function Documentation

◆ calculate()

void ThirdOrderDelayODE::calculate	(	double	tau,
		double	inputCurrent,
		double	dIdt,
		double	d2Id2t,
		const double *	x,
		unsigned	lenX,
		unsigned	firstNode,
		double *	derivative
	)		const

Definition at line 7 of file ThirdOrderDelayODE.cc.

                                                              {
   // Check input sanity
   if (lenX < firstNode + 3U)
     throw cms::Exception("In ThirdOrderDelayODE::calculate: insufficient number of variables");
   if (tau <= 0.0)
     throw cms::Exception("In ThirdOrderDelayODE::calculate: delay time is not positive");
   assert(x);
   assert(derivative);
  
   derivative[firstNode] = x[firstNode + 1];
   derivative[firstNode + 1] = x[firstNode + 2];
   derivative[firstNode + 2] =
       6.0 / a_ * (currentIn - x[firstNode] - c_ * tau * x[firstNode + 1] - b_ / 2.0 * tau * tau * x[firstNode + 2]) /
       tau / tau / tau;
 }

References a_, cms::cuda::assert(), b_, c_, funct::derivative(), Exception, metsig::tau, mitigatedMETSequence_cff::U, and x.

◆ getPadeColumn()

unsigned ThirdOrderDelayODE::getPadeColumn ( ) const

inline

Definition at line 23 of file ThirdOrderDelayODE.h.

23 { return 3U; }

References mitigatedMETSequence_cff::U.

◆ getPadeRow()

unsigned ThirdOrderDelayODE::getPadeRow ( ) const

inline

Definition at line 22 of file ThirdOrderDelayODE.h.

22 { return 0U; }

References mitigatedMETSequence_cff::U.

◆ nParameters()

unsigned ThirdOrderDelayODE::nParameters ( ) const

inline

Definition at line 24 of file ThirdOrderDelayODE.h.

24 { return 3U; }

References mitigatedMETSequence_cff::U.

◆ setParameters()

void ThirdOrderDelayODE::setParameters	(	const double *	pars,
		unsigned	nPars
	)

Definition at line 30 of file ThirdOrderDelayODE.cc.

                                                                                {
   assert(nPars == 3U);
   assert(pars);
   a_ = exp(pars[0]);
   b_ = exp(pars[1]);
   c_ = exp(pars[2]);
 }