/usr/include/CLHEP/GenericFunctions/ExtendedButcherTableau.hh is in libclhep-dev 2.1.4.1-1.1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 | #ifndef _ExtendedButcherTableau_h_
#define _ExtendedButcherTableau_h_
// This class defines a ExtendedButcher Tableau, which completely specifies
// an *embedded* Runge-Kutte integration scheme. ExtendedButcher Tableau
// are described in Numerical Methods for Ordinary Differential Equations,
// John Wiley & sons, West Sussex England.
//
// General form is :
//
// c|A
// ---
// |b^T
// |bHat^T
//
// where A is a matrix and b, bHat, and c are column vectors.
//
// The ExtendedButcher Tableau Class presents itself as an empty structure
// that the user has to fill up. One can blithely fill write into
// any element of A, b, bHat or c. Space is automatically allocated.
#include <vector>
#include <string>
namespace Genfun {
class ExtendedButcherTableau {
public:
// Constructor:
inline ExtendedButcherTableau(const std::string &name, unsigned int order, unsigned int orderHat);
// Returns the name:
inline const std::string & name() const;
// Returns the order of the main formula
inline unsigned int order() const;
// Returns the order of the controlling formula
inline unsigned int orderHat() const;
// Returns the number of steps:
inline unsigned int nSteps() const;
// Write access to elements:
inline double & A(unsigned int i, unsigned int j);
inline double & b(unsigned int i);
inline double & bHat(unsigned int i);
inline double & c(unsigned int i);
// Read access to elements (inline for speed)
inline const double & A(unsigned int i, unsigned int j) const;
inline const double & b(unsigned int i) const;
inline const double & bHat(unsigned int i) const;
inline const double & c(unsigned int i) const;
private:
std::vector< std::vector<double> > _A;
std::vector<double> _b;
std::vector<double> _bHat;
std::vector<double> _c;
std::string _name;
unsigned int _order;
unsigned int _orderHat;
};
class HeunEulerXtTableau: public ExtendedButcherTableau {
// Constructor:
public:
inline HeunEulerXtTableau();
};
class BogackiShampineXtTableau: public ExtendedButcherTableau {
// Constructor:
public:
inline BogackiShampineXtTableau();
};
class FehlbergRK45F2XtTableau: public ExtendedButcherTableau {
// Constructor:
public:
inline FehlbergRK45F2XtTableau();
};
class CashKarpXtTableau: public ExtendedButcherTableau {
// Constructor:
public:
inline CashKarpXtTableau();
};
}
inline std::ostream & operator << (std::ostream & o, const Genfun::ExtendedButcherTableau & b);
#include "CLHEP/GenericFunctions/ExtendedButcherTableau.icc"
#endif
|