Genfun::RKIntegrator::RKFunction Class Reference

#include <RKIntegrator.hh>

Inheritance diagram for Genfun::RKIntegrator::RKFunction:

Public Member Functions
	RKFunction (RKData *data, unsigned int index)
virtual	~RKFunction ()
	RKFunction (const RKFunction &right)
virtual double	operator() (double argument) const override
virtual double	operator() (const Argument &a) const override
Public Member Functions inherited from Genfun::AbsFunction
	AbsFunction ()
	AbsFunction (const AbsFunction &right)
virtual	~AbsFunction ()
virtual unsigned int	dimensionality () const
virtual double	operator() (double argument) const =0
virtual double	operator() (const Argument &argument) const =0
virtual AbsFunction *	clone () const =0
virtual FunctionComposition	operator() (const AbsFunction &f) const
virtual ParameterComposition	operator() (const AbsParameter &p) const
Derivative	derivative (const Variable &v) const
Derivative	prime () const
virtual bool	hasAnalyticDerivative () const
virtual Derivative	partial (unsigned int) const

Detailed Description

Definition at line 140 of file RKIntegrator.hh.

Constructor & Destructor Documentation

RKFunction() [1/2]

Genfun::RKIntegrator::RKFunction::RKFunction	(	RKData *	data,
		unsigned int	index
	)

Definition at line 14 of file RKIntegrator.cc.

  :_data(data),
   _index(index)
{
  _data->ref();
}

~RKFunction()

Genfun::RKIntegrator::RKFunction::~RKFunction ( )

virtual

Definition at line 21 of file RKIntegrator.cc.

{
  _data->unref();
}

RKFunction() [2/2]

Genfun::RKIntegrator::RKFunction::RKFunction

(

const RKFunction &

right

)

Definition at line 26 of file RKIntegrator.cc.

  :AbsFunction(right),
   _data(right._data),
   _index(right._index)
{
  _data->ref();
}

Member Function Documentation

operator()() [1/2]

virtual double Genfun::RKIntegrator::RKFunction::operator() ( const Argument &  a ) const

inlineoverridevirtual

Implements Genfun::AbsFunction.

Definition at line 157 of file RKIntegrator.hh.

{return operator() (a[0]);}

operator()() [2/2]

double Genfun::RKIntegrator::RKFunction::operator() ( double  argument ) const

overridevirtual

Implements Genfun::AbsFunction.

Definition at line 35 of file RKIntegrator.cc.

                                                         {
  if (t<0) return 0;
  if (!_data->_locked) _data->lock();
 
  // Do this first, thereafter, just read the cache
  _data->recache();
  
 
  // If the cache is empty, make an entry for t=0;
  size_t nvar = _data->_startingValParameter.size();
  if (_data->_fx.empty()) {
    RKData::Data d((int)nvar);
    d.time=0;
    Argument arg((int)nvar);
    for (size_t f=0;f<nvar;f++) {
      d.variable[f]=_data->_startingValParameterCache[f];
      arg[(int)f]=d.variable[f];
    }
    _data->_fx.insert(d);
  }
 
  if (t==0) return (*_data->_fx.begin()).variable[_index];
 
  RKData::Data dt((int)nvar);
  dt.time=t;
  std::set<RKData::Data>::iterator l =_data->_fx.lower_bound(dt);
 
  // We did find an exact value (measure 0), just return it. 
  if (l!=_data->_fx.end() && (*l).time==t) {
    return (*l).variable[_index];
  }
 
  else {
    std::set<RKData::Data>::iterator u =_data->_fx.upper_bound(dt);
    
    while (u==_data->_fx.end()) {
      u--;
      RKData::Data newData((int)nvar);;
      _data->_stepper->step(_data,*u, newData, 0);
      _data->_fx.insert(l,newData);
      if (newData.time==t) return newData.variable[_index];
      u = _data->_fx.upper_bound(dt);
    }
    u--;
    _data->_stepper->step(_data,*u, dt, t);
    return dt.variable[_index];
  }
}

Referenced by operator()().

---

The documentation for this class was generated from the following files:

• RKIntegrator.hh
• RKIntegrator.cc