ThirdOrderDelayODE Class Reference

#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 ( unsigned  , unsigned    )

inline

Definition at line 11 of file ThirdOrderDelayODE.h.

: a_(1.0) {}

Member Function Documentation

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.

References a_, cms::cuda::assert(), b_, c_, and Exception.

                                                             {
  // 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;
}

unsigned ThirdOrderDelayODE::getPadeColumn ( ) const

inline

Definition at line 23 of file ThirdOrderDelayODE.h.

{ return 3U; }

unsigned ThirdOrderDelayODE::getPadeRow ( ) const

inline

Definition at line 22 of file ThirdOrderDelayODE.h.

{ return 0U; }

unsigned ThirdOrderDelayODE::nParameters ( ) const

inline

Definition at line 24 of file ThirdOrderDelayODE.h.

{ return 3U; }

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

Definition at line 30 of file ThirdOrderDelayODE.cc.

References a_, cms::cuda::assert(), b_, c_, and funct::exp().

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

Member Data Documentation

double ThirdOrderDelayODE::a_

private

Definition at line 30 of file ThirdOrderDelayODE.h.

Referenced by calculate(), and setParameters().

double ThirdOrderDelayODE::b_

private

Definition at line 31 of file ThirdOrderDelayODE.h.

Referenced by calculate(), and setParameters().

double ThirdOrderDelayODE::c_

private

Definition at line 32 of file ThirdOrderDelayODE.h.

Referenced by calculate(), and setParameters().