29 #ifndef RKIntegrator_h 30 #define RKIntegrator_h 1 63 const std::string & variableName=
"anon",
64 double defStartingValue=0.0,
65 double startingValueMin=0.0,
66 double startingValueMax=0.0);
72 double defStartingValue=0.0,
73 double startingValueMin=0.0,
74 double startingValueMax=0.0);
79 const RKFunction *
getFunction(
unsigned int i)
const;
96 std::vector<const RKFunction *> _fcn;
110 std::vector<double> variable;
111 mutable std::vector<double> firstDerivative;
114 Data(
int size): variable(size), firstDerivative(size), time(0) {}
123 std::vector<Parameter *> _startingValParameter;
124 std::vector<double> _startingValParameterCache;
126 std::vector <Parameter *> _controlParameter;
127 std::vector <double> _controlParameterCache;
129 std::vector<const AbsFunction *> _diffEqn;
136 friend class ImaginaryFriend;
156 virtual double operator ()(
double argument)
const;
157 virtual double operator ()(
const Argument &
a)
const {
return operator() (a[0]);}
166 const unsigned int _index;
179 double timeLimit=0)
const =0;
bool operator<(const HepRotation &r, const HepLorentzRotation <)
#define FUNCTION_OBJECT_DEF(classname)
const RKFunction * getFunction(unsigned int i) const
RKIntegrator(const RKStepper *stepper=NULL)
Parameter * createControlParameter(const std::string &variableName="anon", double defStartingValue=0.0, double startingValueMin=0.0, double startingValueMax=0.0)
Parameter * addDiffEquation(const AbsFunction *diffEquation, const std::string &variableName="anon", double defStartingValue=0.0, double startingValueMin=0.0, double startingValueMax=0.0)
bool operator==(const HepRotation &r, const HepLorentzRotation <)